On the causes, consequences, and avoidance of PCR duplicates: Towards a theory of library complexity.

Library preparation protocols for most sequencing technologies involve PCR amplification of the template DNA, which open the possibility that a given template DNA molecule is sequenced multiple times. Reads arising from this phenomenon, known as PCR duplicates, inflate the cost of sequencing and can jeopardize the reliability of affected experiments. Despite the pervasiveness of this artefact, our understanding of its causes and of its impact on downstream statistical analyses remains essentially empirical. Here, we develop a general quantitative model of amplification distortions in sequencing data sets, which we leverage to investigate the factors controlling the occurrence of PCR duplicates. We show that the PCR duplicate rate is determined primarily by the ratio between library complexity and sequencing depth, and that amplification noise (including in its dependence on the number of PCR cycles) only plays a secondary role for this artefact. We confirm our predictions using new and published RAD-seq libraries and provide a method to estimate library complexity and amplification noise in any data set containing PCR duplicates. We discuss how amplification-related artefacts impact downstream analyses, and in particular genotyping accuracy. The proposed framework unites the numerous observations made on PCR duplicates and will be useful to experimenters of all sequencing technologies where DNA availability is a concern.

Genome-wide parallelism underlies contemporary adaptation in urban lizards.

Urbanization drastically transforms landscapes, resulting in fragmentation, degradation, and the loss of local biodiversity. Yet, urban environments also offer opportunities to observe rapid evolutionary change in wild populations that survive and even thrive in these novel habitats. In many ways, cities represent replicated "natural experiments" in which geographically separated populations adaptively respond to similar selection pressures over rapid evolutionary timescales. Little is known, however, about the genetic basis of adaptive phenotypic differentiation in urban populations nor the extent to which phenotypic parallelism is reflected at the genomic level with signatures of parallel selection. Here, we analyzed the genomic underpinnings of parallel urban-associated phenotypic change in Anolis cristatellus, a small-bodied neotropical lizard found abundantly in both urbanized and forested environments. We show that phenotypic parallelism in response to parallel urban environmental change is underlain by genomic parallelism and identify candidate loci across the Anolis genome associated with this adaptive morphological divergence. Our findings point to polygenic selection on standing genetic variation as a key process to effectuate rapid morphological adaptation. Identified candidate loci represent several functions associated with skeletomuscular development, morphology, and human disease. Taken together, these results shed light on the genomic basis of complex morphological adaptations, provide insight into the role of contingency and determinism in adaptation to novel environments, and underscore the value of urban environments to address fundamental evolutionary questions.

Humanization of wildlife gut microbiota in urban environments.

Urbanization is rapidly altering Earth's environments, demanding investigation of the impacts on resident wildlife. Here, we show that urban populations of coyotes (Canis latrans), crested anole lizards (Anolis cristatellus), and white-crowned sparrows (Zonotrichia leucophrys) acquire gut microbiota constituents found in humans, including gut bacterial lineages associated with urbanization in humans. Comparisons of urban and rural wildlife and human populations revealed significant convergence of gut microbiota among urban populations relative to rural populations. All bacterial lineages overrepresented in urban wildlife relative to rural wildlife and differentially abundant between urban and rural humans were also overrepresented in urban humans relative to rural humans. Remarkably, the bacterial lineage most overrepresented in urban anoles was a Bacteroides sequence variant that was also the most significantly overrepresented in urban human populations. These results indicate parallel effects of urbanization on human and wildlife gut microbiota and suggest spillover of bacteria from humans into wildlife in cities.

Vertebrate species, such as reptiles, birds or mammals, harbour distinct communities of microbes in their digestive systems. These miniature ecosystems ­ also known as microbiomes ­ are unique to each owner and species, reflecting their diverse lifestyles and evolutionary history. Urbanisation can disrupt these delicate intestinal communities. Humans and other animals living in cities have different gut microbes to their counterparts living in rural areas. And captive species in homes and zoos often acquire human gut bacteria in their digestive systems, which can lead to health problems in these animals. So far, it has been unclear whether such a humanization of gut bacteria also affects wild animals living in and around cities. To investigate this further, Dillard et al. compared the gut microbes of wild reptiles, birds, and mammals living in close contact with humans in North America, such as coyotes, crested anole lizards and white-crowned sparrows. DNA sequencing showed that in urban environments, the composition of gut bacteria living in all three wildlife species resembled the ones in humans. The types of bacteria overrepresented in the guts of urban humans were also overrepresented in urban wildlife. This suggests that urbanization can affect the composition of gut bacteria in wildlife species by disrupting or replacing portions of their microbiome. The reason for this pattern is unclear. It is possible that humans might be sharing their gut microbes directly with city animals, or that a human-like diet is causing the change. Given the role that gut microbes play in health and disease, it is important to find out whether these changes cause the animals any harm.

Genetic variation in haemoglobin is associated with evolved changes in breathing in high-altitude deer mice.

Physiological systems often have emergent properties but the effects of genetic variation on physiology are often unknown, which presents a major challenge to understanding the mechanisms of phenotypic evolution. We investigated whether genetic variants in haemoglobin (Hb) that contribute to high-altitude adaptation in deer mice (Peromyscus maniculatus) are associated with evolved changes in the control of breathing. We created F2 inter-population hybrids of highland and lowland deer mice to test for phenotypic associations of α- and ß-globin variants on a mixed genetic background. Hb genotype had expected effects on Hb-O2 affinity that were associated with differences in arterial O2 saturation in hypoxia. However, high-altitude genotypes were also associated with breathing phenotypes that should contribute to enhancing O2 uptake in hypoxia. Mice with highland α-globin exhibited a more effective breathing pattern, with highland homozygotes breathing deeper but less frequently across a range of inspired O2, and this difference was comparable to the evolved changes in breathing pattern in deer mouse populations native to high altitude. The ventilatory response to hypoxia was augmented in mice that were homozygous for highland ß-globin. The association of globin variants with variation in breathing phenotypes could not be recapitulated by acute manipulation of Hb-O2 affinity, because treatment with efaproxiral (a synthetic drug that acutely reduces Hb-O2 affinity) had no effect on breathing in normoxia or hypoxia. Therefore, adaptive variation in Hb may have unexpected effects on physiology in addition to the canonical function of this protein in circulatory O2 transport.

Physiological costs of undocumented human migration across the southern United States border.

Political, economic, and climatic upheaval can result in mass human migration across extreme terrain in search of more humane living conditions, exposing migrants to environments that challenge human tolerance. An empirical understanding of the biological stresses associated with these migrations will play a key role in the development of social, political, and medical strategies for alleviating adverse effects and risk of death. We model physiological stress associated with undocumented migration across a commonly traversed section of the southern border of the United States and find that locations of migrant death are disproportionately clustered within regions of greatest predicted physiological stress (evaporative water loss). Minimum values of estimated evaporative water loss were sufficient to cause severe dehydration and associated proximate causes of mortality. Integration of future climate predictions into models increased predicted physiological costs of migration by up to 34.1% over the next 30 years.

Ivory poaching and the rapid evolution of tusklessness in African elephants.

Understanding the evolutionary consequences of wildlife exploitation is increasingly important as harvesting becomes more efficient. We examined the impacts of ivory poaching during the Mozambican Civil War (1977 to 1992) on the evolution of African savanna elephants (Loxodonta africana) in Gorongosa National Park. Poaching resulted in strong selection that favored tusklessness amid a rapid population decline. Survey data revealed tusk-inheritance patterns consistent with an X chromosome­linked dominant, male-lethal trait. Whole-genome scans implicated two candidate genes with known roles in mammalian tooth development (AMELX and MEP1a), including the formation of enamel, dentin, cementum, and the periodontium. One of these loci (AMELX) is associated with an X-linked dominant, male-lethal syndrome in humans that diminishes the growth of maxillary lateral incisors (homologous to elephant tusks). This study provides evidence for rapid, poaching-mediated selection for the loss of a prominent anatomical trait in a keystone species.

Selection on adaptive and maladaptive gene expression plasticity during thermal adaptation to urban heat islands.

Phenotypic plasticity enables a single genotype to produce multiple phenotypes in response to environmental variation. Plasticity may play a critical role in the colonization of novel environments, but its role in adaptive evolution is controversial. Here we suggest that rapid parallel regulatory adaptation of Anolis lizards to urban heat islands is due primarily to selection for reduced and/or reversed heat-induced plasticity that is maladaptive in urban thermal conditions. We identify evidence for polygenic selection across genes of the skeletal muscle transcriptome associated with heat tolerance. Forest lizards raised in common garden conditions exhibit heat-induced changes in expression of these genes that largely correlate with decreased heat tolerance, consistent with maladaptive regulatory response to high-temperature environments. In contrast, urban lizards display reduced gene expression plasticity after heat challenge in common garden and a significant increase in gene expression change that is congruent with greater heat tolerance, a putatively adaptive state in warmer urban environments. Genes displaying maladaptive heat-induced plasticity repeatedly show greater genetic divergence between urban and forest habitats than those displaying adaptive plasticity. These results highlight the role of selection against maladaptive regulatory plasticity during rapid adaptive modification of complex systems in the wild.

The adaptive benefit of evolved increases in hemoglobin-O2 affinity is contingent on tissue O2 diffusing capacity in high-altitude deer mice.

BACKGROUND: Complex organismal traits are often the result of multiple interacting genes and sub-organismal phenotypes, but how these interactions shape the evolutionary trajectories of adaptive traits is poorly understood. We examined how functional interactions between cardiorespiratory traits contribute to adaptive increases in the capacity for aerobic thermogenesis (maximal O2 consumption, V̇O2max, during acute cold exposure) in high-altitude deer mice (Peromyscus maniculatus). We crossed highland and lowland deer mice to produce F2 inter-population hybrids, which expressed genetically based variation in hemoglobin (Hb) O2 affinity on a mixed genetic background. We then combined physiological experiments and mathematical modeling of the O2 transport pathway to examine the links between cardiorespiratory traits and V̇O2max. RESULTS: Physiological experiments revealed that increases in Hb-O2 affinity of red blood cells improved blood oxygenation in hypoxia but were not associated with an enhancement in V̇O2max. Sensitivity analyses performed using mathematical modeling showed that the influence of Hb-O2 affinity on V̇O2max in hypoxia was contingent on the capacity for O2 diffusion in active tissues. CONCLUSIONS: These results suggest that increases in Hb-O2 affinity would only have adaptive value in hypoxic conditions if concurrent with or preceded by increases in tissue O2 diffusing capacity. In high-altitude deer mice, the adaptive benefit of increasing Hb-O2 affinity is contingent on the capacity to extract O2 from the blood, which helps resolve controversies about the general role of hemoglobin function in hypoxia tolerance.

X chromosome-dependent disruption of placental regulatory networks in hybrid dwarf hamsters.

Embryonic development in mammals is highly sensitive to changes in gene expression within the placenta. The placenta is also highly enriched for genes showing parent-of-origin or imprinted expression, which is predicted to evolve rapidly in response to parental conflict. However, little is known about the evolution of placental gene expression, or if divergence of placental gene expression plays an important role in mammalian speciation. We used crosses between two species of dwarf hamsters (Phodopus sungorus and Phodopus campbelli) to examine the genetic and regulatory underpinnings of severe placental overgrowth in their hybrids. Using quantitative genetic mapping and mitochondrial substitution lines, we show that overgrowth of hybrid placentas was primarily caused by genetic differences on the maternally inherited P. sungorus X chromosome. Mitochondrial interactions did not contribute to abnormal hybrid placental development, and there was only weak correspondence between placental disruption and embryonic growth. Genome-wide analyses of placental transcriptomes from the parental species and first- and second-generation hybrids revealed a central group of co-expressed X-linked and autosomal genes that were highly enriched for maternally biased expression. Expression of this gene network was strongly correlated with placental size and showed widespread misexpression dependent on epistatic interactions with X-linked hybrid incompatibilities. Collectively, our results indicate that the X chromosome is likely to play a prominent role in the evolution of placental gene expression and the accumulation of hybrid developmental barriers between mammalian species.

Parallel selection on thermal physiology facilitates repeated adaptation of city lizards to urban heat islands.

Only recently have we begun to understand the ecological and evolutionary effects of urbanization on species, with studies revealing drastic impacts on community composition, gene flow, behaviour, morphology and physiology. However, our understanding of how adaptive evolution allows species to persist, and even thrive, in urban landscapes is still nascent. Here, we examine phenotypic, genomic and regulatory impacts of urbanization on a widespread lizard, the Puerto Rican crested anole (Anolis cristatellus). We find that urban lizards endure higher environmental temperatures and display greater heat tolerance than their forest counterparts. A single non-synonymous polymorphism within a protein synthesis gene (RARS) is associated with heat tolerance plasticity within urban heat islands and displays parallel signatures of selection in cities. Additionally, we identify groups of differentially expressed genes between habitats showing elevated genetic divergence in multiple urban-forest comparisons. These genes display evidence of adaptive regulatory evolution within cities and disproportionately cluster within regulatory modules associated with heat tolerance. This study provides evidence of temperature-mediated selection in urban heat islands with repeatable impacts on physiological evolution at multiple levels of biological hierarchy.

Physiological and regulatory underpinnings of geographic variation in reptilian cold tolerance across a latitudinal cline.

Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CTmin ), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co-expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity-limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.

Winter storms drive rapid phenotypic, regulatory, and genomic shifts in the green anole lizard.

Extreme environmental perturbations offer opportunities to observe the effects of natural selection in wild populations. During the winter of 2013-2014, the southeastern United States endured an extreme cold event. We used thermal performance, transcriptomics, and genome scans to measure responses of lizard populations to storm-induced selection. We found significant increases in cold tolerance at the species' southern limit. Gene expression in southern survivors shifted toward patterns characteristic of northern populations. Comparing samples before and after the extreme winter, 14 genomic regions were differentiated in the surviving southern population; four also exhibited signatures of local adaptation across the latitudinal gradient and implicate genes involved in nervous system function. Together, our results suggest that extreme winter events can rapidly produce strong selection on natural populations at multiple biological levels that recapitulate geographic patterns of local adaptation.

Out of Florida: mtDNA reveals patterns of migration and Pleistocene range expansion of the Green Anole lizard (Anolis carolinensis).

Anolis carolinensis is an emerging model species and the sole member of its genus native to the United States. Considerable morphological and physiological variation has been described in the species, and the recent sequencing of its genome makes it an attractive system for studies of genome variation. To inform future studies of molecular and phenotypic variation within A. carolinensis, a rigorous account of intraspecific population structure and relatedness is needed. Here, we present the most extensive phylogeographic study of this species to date. Phylogenetic analyses of mitochondrial DNA sequence data support the previous hypothesis of a western Cuban origin of the species. We found five well-supported, geographically distinct mitochondrial haplotype clades throughout the southeastern United States. Most Florida populations fall into one of three divergent clades, whereas the vast majority of populations outside Florida belong to a single, shallowly diverged clade. Genetic boundaries do not correspond to major rivers, but may reflect effects of Pleistocene glaciation events and the Appalachian Mountains on migration and expansion of the species. Phylogeographic signal should be examined using nuclear loci to complement these findings.