Humans sweat to cool their bodies and have by far the highest eccrine sweat gland density among primates. Humans’ high eccrine gland density has long been recognized as a hallmark human evolutionary adaptation, but its genetic basis has been unknown. In humans, expression of the Engrailed 1 (EN1) transcription factor correlates with the onset of eccrine gland formation. In mice, regulation of ectodermal En1 expression is a major determinant of natural variation in eccrine gland density between strains, and increased En1 expression promotes the specification of more eccrine glands. Here, we show that regulation of EN1 has evolved specifically on the human lineage to promote eccrine gland formation. Using comparative genomics and validation of ectodermal enhancer activity in mice, we identified a human EN1 skin enhancer, hECE18. We showed that multiple epistatically interacting derived substitutions in the human ECE18 enhancer increased its activity compared with nonhuman ape orthologs in cultured keratinocytes. Repression of hECE18 in human cultured keratinocytes specifically attenuated EN1 expression, indicating this element positively regulates EN1 in this context. In a humanized enhancer knock-in mouse, hECE18 increased developmental En1 expression in the skin to induce the formation of more eccrine glands. Our study uncovers a genetic basis contributing to the evolution of one of the most singular human adaptations and implicates multiple interacting mutations in a single enhancer as a mechanism for human evolutionary change.

Domestic dogs have experienced population bottlenecks, recent inbreeding, and strong artificial selection. These processes have simplified the genetic architecture of complex traits, allowed deleterious variation to persist, and increased both identity-by-descent (IBD) segments and runs of homozygosity (ROH). As such, dogs provide an excellent model for examining how these evolutionary processes influence disease. We assembled a dataset containing 4,414 breed dogs, 327 village dogs, and 380 wolves genotyped at 117,288 markers and data for clinical and morphological phenotypes. Breed dogs have an enrichment of IBD and ROH, relative to both village dogs and wolves, and we use these patterns to show that breed dogs have experienced differing severities of bottlenecks in their recent past. We then found that ROH burden is associated with phenotypes in breed dogs, such as lymphoma. We next test the prediction that breeds with greater ROH have more disease alleles reported in the Online Mendelian Inheritance in Animals (OMIA). Surprisingly, the number of causal variants identified correlates with the popularity of that breed rather than the ROH or IBD burden, suggesting an ascertainment bias in OMIA. Lastly, we use the distribution of ROH across the genome to identify genes with depletions of ROH as potential hotspots for inbreeding depression and find multiple exons where ROH are never observed. Our results suggest that inbreeding has played a large role in shaping genetic and phenotypic variation in dogs and that future work on understudied breeds may reveal new disease-causing variation.

Maintaining internal salt and water balance in response to fluctuating external conditions is essential for animal survival. This is particularly true for insects as their high surface-to-volume ratio makes them highly susceptible to osmotic stress. However, the cellular and hormonal mechanisms that mediate the systemic control of osmotic homeostasis in beetles (Coleoptera), the largest group of insects, remain largely unidentified. Here, we demonstrate that eight neurons in the brain of the red flour beetle Tribolium castaneum respond to internal changes in osmolality by releasing diuretic hormone (DH) 37 and DH47—homologs of vertebrate corticotropin-releasing factor (CRF) hormones—to control systemic water balance. Knockdown of the gene encoding the two hormones (Urinate, Urn8) reduces Malpighian tubule secretion and restricts organismal fluid loss, whereas injection of DH37 or DH47 reverses these phenotypes. We further identify a CRF-like receptor, Urinate receptor (Urn8R), which is exclusively expressed in a functionally unique secondary cell in the beetle tubules, as underlying this response. Activation of Urn8R increases K⁺ secretion, creating a lumen-positive transepithelial potential that drives fluid secretion. Together, these data show that beetle Malpighian tubules operate by a fundamentally different mechanism than those of other insects. Finally, we adopt a fluorescent labeling strategy to identify the evolutionary origin of this unusual tubule architecture, revealing that it evolved in the last common ancestor of the higher beetle families. Our work thus uncovers an important homeostatic program that is key to maintaining osmotic control in beetles, which evolved parallel to the radiation of the “advanced” beetle lineages.

Oral venom systems evolved multiple times in numerous vertebrates enabling the exploitation of unique predatory niches. Yet how and when they evolved remains poorly understood. Up to now, most research on venom evolution has focused strictly on the toxins. However, using toxins present in modern day animals to trace the origin of the venom system is difficult, since they tend to evolve rapidly, show complex patterns of expression, and were incorporated into the venom arsenal relatively recently. Here we focus on gene regulatory networks associated with the production of toxins in snakes, rather than the toxins themselves. We found that overall venom gland gene expression was surprisingly well conserved when compared to salivary glands of other amniotes. We characterized the “metavenom network,” a network of ∼3,000 nonsecreted housekeeping genes that are strongly coexpressed with the toxins, and are primarily involved in protein folding and modification. Conserved across amniotes, this network was coopted for venom evolution by exaptation of existing members and the recruitment of new toxin genes. For instance, starting from this common molecular foundation, Heloderma lizards, shrews, and solenodon, evolved venoms in parallel by overexpression of kallikreins, which were common in ancestral saliva and induce vasodilation when injected, causing circulatory shock. Derived venoms, such as those of snakes, incorporated novel toxins, though still rely on hypotension for prey immobilization. These similarities suggest repeated cooption of shared molecular machinery for the evolution of oral venom in mammals and reptiles, blurring the line between truly venomous animals and their ancestors.

Dive capacities of air-breathing vertebrates are dictated by onboard O₂ stores, suggesting that physiologic specialization of diving birds such as penguins may have involved adaptive changes in convective O₂ transport. It has been hypothesized that increased hemoglobin (Hb)-O₂ affinity improves pulmonary O₂ extraction and enhances the capacity for breath-hold diving. To investigate evolved changes in Hb function associated with the aquatic specialization of penguins, we integrated comparative measurements of whole-blood and purified native Hb with protein engineering experiments based on site-directed mutagenesis. We reconstructed and resurrected ancestral Hb representing the common ancestor of penguins and the more ancient ancestor shared by penguins and their closest nondiving relatives (order Procellariiformes, which includes albatrosses, shearwaters, petrels, and storm petrels). These two ancestors bracket the phylogenetic interval in which penguin-specific changes in Hb function would have evolved. The experiments revealed that penguins evolved a derived increase in Hb-O₂ affinity and a greatly augmented Bohr effect (i.e., reduced Hb-O₂ affinity at low pH). Although an increased Hb-O₂ affinity reduces the gradient for O₂ diffusion from systemic capillaries to metabolizing cells, this can be compensated by a concomitant enhancement of the Bohr effect, thereby promoting O₂ unloading in acidified tissues. We suggest that the evolved increase in Hb-O₂ affinity in combination with the augmented Bohr effect maximizes both O₂ extraction from the lungs and O₂ unloading from the blood, allowing penguins to fully utilize their onboard O₂ stores and maximize underwater foraging time.

Species in a shared environment tend to evolve similar adaptations under the influence of their phylogenetic context. Using snowfinches, a monophyletic group of passerine birds (Passeridae), we study the relative roles of ancestral and species-specific adaptations to an extreme high-elevation environment, the Qinghai–Tibet Plateau. Our ancestral trait reconstruction shows that the ancestral snowfinch occupied high elevations and had a larger body mass than most nonsnowfinches in Passeridae. Subsequently, this phenotypic adaptation diversified in the descendant species. By comparing high-quality genomes from representatives of the three phylogenetic lineages, we find that about 95% of genes under positive selection in the descendant species are different from those in the ancestor. Consistently, the biological functions enriched for these species differ from those of their ancestor to various degrees (semantic similarity values ranging from 0.27 to 0.5), suggesting that the three descendant species have evolved divergently from the initial adaptation in their common ancestor. Using a functional assay to a highly selective gene, DTL, we demonstrate that the nonsynonymous substitutions in the ancestor and descendant species have improved the repair capacity of ultraviolet-induced DNA damage. The repair kinetics of the DTL gene shows a twofold to fourfold variation across the ancestor and the descendants. Collectively, this study reveals an exceptional case of adaptive evolution to high-elevation environments, an evolutionary process with an initial adaptation in the common ancestor followed by adaptive diversification of the descendant species.

Despite their close genetic relatedness, apes and African and Asian monkeys (AAMs) differ in their susceptibility to severe bacterial and viral infections that are important causes of human disease. Such differences between humans and other primates are thought to be a result, at least in part, of interspecies differences in immune response to infection. However, because of the lack of comparative functional data across species, it remains unclear in what ways the immune systems of humans and other primates differ. Here, we report the whole-genome transcriptomic responses of ape species (human and chimpanzee) and AAMs (rhesus macaque and baboon) to bacterial and viral stimulation. We find stark differences in the responsiveness of these groups, with apes mounting a markedly stronger early transcriptional response to both viral and bacterial stimulation, altering the transcription of ∼40% more genes than AAMs. Additionally, we find that genes involved in the regulation of inflammatory and interferon responses show the most divergent early transcriptional responses across primates and that this divergence is attenuated over time. Finally, we find that relative to AAMs, apes engage a much less specific immune response to different classes of pathogens during the early hours of infection, up-regulating genes typical of anti-viral and anti-bacterial responses regardless of the nature of the stimulus. Overall, these findings suggest apes exhibit increased sensitivity to bacterial and viral immune stimulation, activating a broader array of defense molecules that may be beneficial for early pathogen killing at the potential cost of increased energy expenditure and tissue damage.

The rapid rise of antibiotic resistance, combined with the increasing cost and difficulties to develop new antibiotics, calls for treatment strategies that enable more sustainable antibiotic use. The development of such strategies, however, is impeded by the lack of suitable experimental approaches that allow testing their effects under realistic epidemiological conditions. Here, we present an approach to compare the effect of alternative multidrug treatment strategies in vitro using a robotic liquid-handling platform. We use this framework to study resistance evolution and spread implementing epidemiological population dynamics for treatment, transmission, and patient admission and discharge, as may be observed in hospitals. We perform massively parallel experimental evolution over up to 40 d and complement this with a computational model to infer the underlying population-dynamical parameters. We find that in our study, combination therapy outperforms monotherapies, as well as cycling and mixing, in minimizing resistance evolution and maximizing uninfecteds, as long as there is no influx of double resistance into the focal treated community.

Mitochondrial protein import requires outer membrane receptors that evolved independently in different lineages. Here we used quantitative proteomics and in vitro binding assays to investigate the substrate preferences of ATOM46 and ATOM69, the two mitochondrial import receptors of Trypanosoma brucei. The results show that ATOM46 prefers presequence-containing, hydrophilic proteins that lack transmembrane domains (TMDs), whereas ATOM69 prefers presequence-lacking, hydrophobic substrates that have TMDs. Thus, the ATOM46/yeast Tom20 and the ATOM69/yeast Tom70 pairs have similar substrate preferences. However, ATOM46 mainly uses electrostatic, and Tom20 hydrophobic, interactions for substrate binding. In vivo replacement of T. brucei ATOM46 by yeast Tom20 did not restore import. However, replacement of ATOM69 by the recently discovered Tom36 receptor of Trichomonas hydrogenosomes, while not allowing for growth, restored import of a large subset of trypanosomal proteins that lack TMDs. Thus, even though ATOM69 and Tom36 share the same domain structure and topology, they have different substrate preferences. The study establishes complementation experiments, combined with quantitative proteomics, as a highly versatile and sensitive method to compare in vivo preferences of protein import receptors. Moreover, it illustrates the role determinism and contingencies played in the evolution of mitochondrial protein import receptors.

Through the secretion of “public goods” molecules, microbes cooperatively exploit their habitat. This is known as a major driver of the functioning of microbial communities, including in human disease. Understanding why microbial species cooperate is therefore crucial to achieve successful microbial community management, such as microbiome manipulation. A leading explanation is that of Hamilton’s inclusive-fitness framework. A cooperator can indirectly transmit its genes by helping the reproduction of an individual carrying similar genes. Therefore, all else being equal, as relatedness among individuals increases, so should cooperation. However, the predictive power of relatedness, particularly in microbes, is surrounded by controversy. Using phylogenetic comparative analyses across the full diversity of the human gut microbiota and six forms of cooperation, we find that relatedness is predictive of the cooperative gene content evolution in gut-microbe genomes. Hence, relatedness is predictive of cooperation over broad microbial taxonomic levels that encompass variation in other life-history and ecology details. This supports the generality of Hamilton’s central insights and the relevance of relatedness as a key parameter of interest to advance microbial predictive and engineering science.

Gene retroposition is known to contribute to patterns of gene evolution and adaptations. However, possible negative effects of gene retroposition remain largely unexplored since most previous studies have focused on between-species comparisons where negatively selected copies are mostly not observed, as they are quickly lost from populations. Here, we show for natural house mouse populations that the primary rate of retroposition is orders of magnitude higher than the long-term rate. Comparisons with single-nucleotide polymorphism distribution patterns in the same populations show that most retroposition events are deleterious. Transcriptomic profiling analysis shows that new retroposed copies become easily subject to transcription and have an influence on the expression levels of their parental genes, especially when transcribed in the antisense direction. Our results imply that the impact of retroposition on the mutational load has been highly underestimated in natural populations. This has additional implications for strategies of disease allele detection in humans.

Many microorganisms with high prevalence in host populations are beneficial to the host and maintained by specialized transmission mechanisms. Although microbial promotion of host fitness and specificity of the associations undoubtedly enhance microbial prevalence, it is an open question whether these symbiotic traits are also a prerequisite for the evolutionary origin of prevalent microbial taxa. To address this issue, we investigate how processes without positive microbial effects on host fitness or host choice can influence the prevalence of certain microbes in a host population. Specifically, we develop a theoretical model to assess the conditions under which particular microbes can become enriched in animal hosts even when they are not providing a specific benefit to a particular host. We find increased prevalence of specific microbes in a host when both show some overlap in their lifecycles, and especially when both share dispersal routes across a patchy habitat distribution. Our results emphasize that host enrichment per se is not a reliable indicator of beneficial host–microbe interactions. The resulting increase in time spent associated with a host may nevertheless give rise to new selection conditions, which can favor microbial adaptations toward a host-associated lifestyle, and, thus, it could be the foundation for subsequent evolution of mutually beneficial coevolved symbioses.

Mapping landscape connectivity is important for controlling invasive species and disease vectors. Current landscape genetics methods are often constrained by the subjectivity of creating resistance surfaces and the difficulty of working with interacting and correlated environmental variables. To overcome these constraints, we combine the advantages of a machine-learning framework and an iterative optimization process to develop a method for integrating genetic and environmental (e.g., climate, land cover, human infrastructure) data. We validate and demonstrate this method for the Aedes aegypti mosquito, an invasive species and the primary vector of dengue, yellow fever, chikungunya, and Zika. We test two contrasting metrics to approximate genetic distance and find Cavalli-Sforza–Edwards distance (CSE) performs better than linearized F_ST. The correlation (R) between the model’s predicted genetic distance and actual distance is 0.83. We produce a map of genetic connectivity for Ae. aegypti’s range in North America and discuss which environmental and anthropogenic variables are most important for predicting gene flow, especially in the context of vector control.

Antagonistic interactions between the sexes are important drivers of evolutionary divergence. Interlocus sexual conflict is generally described as a conflict between alleles at two interacting loci whose identity and genomic location are arbitrary, but with opposite fitness effects in each sex. We build on previous theory by suggesting that when loci under interlocus sexual conflict are located on the sex chromosomes it can lead to cycles of antagonistic coevolution between them and therefore between the sexes. We tested this hypothesis by performing experimental crosses using Drosophila melanogaster where we reciprocally exchanged the sex chromosomes between five allopatric wild-type populations in a round-robin design. Disrupting putatively coevolved sex chromosome pairs resulted in increased male reproductive success in 16 of 20 experimental populations (10 of which were individually significant), but also resulted in lower offspring egg-to-adult viability that affected both male and female fitness. After 25 generations of experimental evolution these sexually antagonistic fitness effects appeared to be resolved. To formalize our hypothesis, we developed population genetic models of antagonistic coevolution using fitness expressions based on our empirical results. Our model predictions support the conclusion that antagonistic coevolution between the sex chromosomes is plausible under the fitness effects observed in our experiments. Together, our results lend both empirical and theoretical support to the idea that cycles of antagonistic coevolution can occur between sex chromosomes and illustrate how this process, in combination with autosomal coadaptation, may drive genetic and phenotypic divergence between populations.

Predator specialization has often been considered an evolutionary “dead end” due to the constraints associated with the evolution of morphological and functional optimizations throughout the organism. However, in some predators, these changes are localized in separate structures dedicated to prey capture. One of the most extreme cases of this modularity can be observed in siphonophores, a clade of pelagic colonial cnidarians that use tentilla (tentacle side branches armed with nematocysts) exclusively for prey capture. Here we study how siphonophore specialists and generalists evolve, and what morphological changes are associated with these transitions. To answer these questions, we: a) Measured 29 morphological characters of tentacles from 45 siphonophore species, b) mapped these data to a phylogenetic tree, and c) analyzed the evolutionary associations between morphological characters and prey-type data from the literature. Instead of a dead end, we found that siphonophore specialists can evolve into generalists, and that specialists on one prey type have directly evolved into specialists on other prey types. Our results show that siphonophore tentillum morphology has strong evolutionary associations with prey type, and suggest that shifts between prey types are linked to shifts in the morphology, mode of evolution, and evolutionary correlations of tentilla and their nematocysts. The evolutionary history of siphonophore specialization helps build a broader perspective on predatory niche diversification via morphological innovation and evolution. These findings contribute to understanding how specialization and morphological evolution have shaped present-day food webs.

The central question in the origin of life is to understand how structure can emerge from randomness. The Eigen theory of replication states, for sequences that are copied one base at a time, that the replication fidelity has to surpass an error threshold to avoid that replicated specific sequences become random because of the incorporated replication errors [M. Eigen, Naturwissenschaften 58 (10), 465–523 (1971)]. Here, we showed that linking short oligomers from a random sequence pool in a templated ligation reaction reduced the sequence space of product strands. We started from 12-mer oligonucleotides with two bases in all possible combinations and triggered enzymatic ligation under temperature cycles. Surprisingly, we found the robust creation of long, highly structured sequences with low entropy. At the ligation site, complementary and alternating sequence patterns developed. However, between the ligation sites, we found either an A-rich or a T-rich sequence within a single oligonucleotide. Our modeling suggests that avoidance of hairpins was the likely cause for these two complementary sequence pools. What emerged was a network of complementary sequences that acted both as templates and substrates of the reaction. This self-selecting ligation reaction could be restarted by only a few majority sequences. The findings showed that replication by random templated ligation from a random sequence input will lead to a highly structured, long, and nonrandom sequence pool. This is a favorable starting point for a subsequent Darwinian evolution searching for higher catalytic functions in an RNA world scenario.

Changes to feeding structures are a fundamental component of the vertebrate transition from water to land. Classically, this event has been characterized as a shift from an aquatic, suction-based mode of prey capture involving cranial kinesis to a biting-based feeding system utilizing a rigid skull capable of capturing prey on land. Here we show that a key intermediate, Tiktaalik roseae, was capable of cranial kinesis despite significant restructuring of the skull to facilitate biting and snapping. Lateral sliding joints between the cheek and dermal skull roof, as well as independent mobility between the hyomandibula and palatoquadrate, enable the suspensorium of T. roseae to expand laterally in a manner similar to modern alligator gars and polypterids. This movement can expand the spiracular and opercular cavities during feeding and respiration, which would direct fluid through the feeding apparatus. Detailed analysis of the sutural morphology of T. roseae suggests that the ability to laterally expand the cheek and palate was maintained during the fish-to-tetrapod transition, implying that limited cranial kinesis was plesiomorphic to the earliest limbed vertebrates. Furthermore, recent kinematic studies of feeding in gars demonstrate that prey capture with lateral snapping can synergistically combine both biting and suction, rather than trading off one for the other. A “gar-like” stage in early tetrapod evolution might have been an important intermediate step in the evolution of terrestrial feeding systems by maintaining suction-generation capabilities while simultaneously elaborating a mechanism for biting-based prey capture.