Improving the estimation of genetic distances from Next-Generation Sequencing data

The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau

Tibetan sheep are the most common and widespread domesticated animals on the Qinghai-Tibetan Plateau (QTP) and have played an essential role in the permanent human occupation of this high-altitude region. However, the precise timing, route, and process of sheep pastoralism in the QTP region remain poorly established, and little is known about the underlying genomic changes that occurred during the process. Here, we investigate the genomic variation in Tibetan sheep using whole-genome sequences, single nucleotide polymorphism arrays, mitochondrial DNA, and Y-chromosomal variants in 986 samples throughout their distribution range. We detect strong signatures of selection in genes involved in the hypoxia and ultraviolet signaling pathways (e.g., HIF-1 pathway and HBB and MITF genes) and in genes associated with morphological traits such as horn size and shape (e.g., RXFP2). We identify clear signals of argali (Ovis ammon) introgression into sympatric Tibetan sheep, covering 5.23-5.79% of their genomes. The introgressed genomic regions are enriched in genes related to oxygen transportation system, sensory perception, and morphological phenotypes, in particular the genes HBB and RXFP2 with strong signs of adaptive introgression. The spatial distribution of genomic diversity and demographic reconstruction of the history of Tibetan sheep show a stepwise pattern of colonization with their initial spread onto the QTP from its northeastern part ∼3,100 years ago, followed by further southwest expansion to the central QTP ∼1,300 years ago. Together with archeological evidence, the date and route reveal the history of human expansions on the QTP by the Tang-Bo Ancient Road during the late Holocene. Our findings contribute to a depth understanding of early pastoralism and the local adaptation of Tibetan sheep as well as the late-Holocene human occupation of the QTP.

Host genetic determinants of the gut microbiota of wild mice

Identifying a common set of genes that mediate host-microbial interactions across populations and species of mammals has broad relevance for human health and animal biology. However, the genetic basis of the gut microbial composition in natural populations remains largely unknown outside of humans. Here, we used wild house mouse populations as a model system to ask three major questions: (a) Does host genetic relatedness explain interindividual variation in gut microbial composition? (b) Do population differences in the microbiota persist in a common environment? (c) What are the host genes associated with microbial richness and the relative abundance of bacterial genera? We found that host genetic distance is a strong predictor of the gut microbial composition as characterized by 16S amplicon sequencing. Using a common garden approach, we then identified differences in microbial composition between populations that persisted in a shared laboratory environment. Finally, we used exome sequencing to associate host genetic variants with microbial diversity and relative abundance of microbial taxa in wild mice. We identified 20 genes that were associated with microbial diversity or abundance including a macrophage-derived cytokine (IL12a) that contained three nonsynonymous mutations. Surprisingly, we found a significant overrepresentation of candidate genes that were previously associated with microbial measurements in humans. The homologous genes that overlapped between wild mice and humans included genes that have been associated with traits related to host immunity and obesity in humans. Gene-bacteria associations identified in both humans and wild mice suggest some commonality to the host genetic determinants of gut microbial composition across mammals.

Local adaptation does not lead to genome-wide differentiation in lava flow lizards

Adaptation can occur with or without genome-wide differentiation. If adaptive loci are linked to traits involved in reproductive isolation, genome-wide divergence is likely, and speciation is possible. However, adaptation can also lead to phenotypic differentiation without genome-wide divergence if levels of ongoing gene flow are high. Here, we use the replicated occurrence of melanism in lava flow lizards to assess the relationship between local adaptation and genome-wide differentiation. We compare patterns of phenotypic and genomic divergence among lava flow and nonlava populations for three lizard species and three lava flows in the Chihuahuan Desert. We find that local phenotypic adaptation (melanism) is not typically accompanied by genome-wide differentiation. Specifically, lava populations do not generally exhibit greater divergence from nonlava populations than expected by geography alone, regardless of whether the lava formation is 5,000 or 760,000 years old. We also infer that gene flow between lava and nonlava populations is ongoing in all lava populations surveyed. Recent work in the isolation by environment and ecological speciation literature suggests that environmentally driven genome-wide differentiation is common in nature. However, local adaptation may often simply be local adaptation rather than an early stage of ecological speciation.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Killer whale genomes reveal a complex history of recurrent admixture and vicariance

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

Signatures of Selection on Standing Genetic Variation Underlie Athletic and Navigational Performance in Racing Pigeons

Racing pigeons have been selectively bred to find their way home quickly over what are often extremely long distances. This breed is of substantial commercial value and is also an excellent avian model to gain empirical insights into the evolution of traits associated with flying performance and spatial orientation. Here, we investigate the molecular basis of the superior athletic and navigational capabilities of racing pigeons using whole-genome and RNA sequencing data. We inferred multiple signatures of positive selection distributed across the genome of racing pigeons. The strongest signature overlapped the CASK gene, a gene implicated in the formation of neuromuscular junctions. However, no diagnostic alleles were found between racing pigeons and other breeds, and only a small proportion of highly differentiated variants were exclusively detected in racing pigeons. We can thus conclude that very few individual genetic changes, if any, are either strictly necessary or sufficient for superior athletics and navigation. Gene expression analysis between racing and nonracing breeds revealed modest differences in muscle (213) and brain (29). These transcripts, however, showed only slightly elevated levels of genetic differentiation between the two groups, suggesting that most differential expression is not causative but likely a consequence of alterations in regulatory networks. Our results show that the unique suite of traits that enable fast flight, long endurance, and accurate navigation in racing pigeons, do not result from few loci acting as master switches but likely from a polygenic architecture that leveraged standing genetic variation available at the onset of the breed formation.

The complex geography of domestication of the African rice Oryza glaberrima

While the domestication history of Asian rice has been extensively studied, details of the evolution of African rice remain elusive. The inner Niger delta has been suggested as the center of origin but molecular data to support this hypothesis is lacking. Here, we present a comprehensive analysis of the evolutionary and domestication history of African rice. By analyzing whole genome re-sequencing data from 282 individuals of domesticated African rice Oryza glaberrima and its progenitor O. barthii, we hypothesize a non-centric (i.e. multiregional) domestication origin for African rice. Our analyses showed genetic structure within O. glaberrima that has a geographical association. Furthermore, we have evidence that the previously hypothesized O. barthii progenitor populations in West Africa have evolutionary signatures similar to domesticated rice and carried causal domestication mutations, suggesting those progenitors were either mislabeled or may actually represent feral wild-domesticated hybrids. Phylogeographic analysis of genes involved in the core domestication process suggests that the origins of causal domestication mutations could be traced to wild progenitors in multiple different locations in West and Central Africa. In addition, measurements of panicle threshability, a key early domestication trait for seed shattering, were consistent with the gene phylogeographic results. We suggest seed non-shattering was selected from multiple genotypes, possibly arising from different geographical regions. Based on our evidence, O. glaberrima was not domesticated from a single centric location but was a result of diffuse process where multiple regions contributed key alleles for different domestication traits.

For many crops it is not clear how they were domesticated from their wild progenitors. Transition from a wild to domesticated state required a series of genetic changes, and studying the evolutionary origin of these domestication-causing mutations are key to understanding the domestication origins of a crop. Moreover, population comparisons provide insight into the relationship between wild and cultivated populations and the evolutionary history of domestication. In this study, we investigated the domestication history of Oryza glaberrima, a rice species that was domesticated in West Africa independent from the Asian rice species O. sativa. Using genome-wide data from a large sample of domesticated and wild African rice samples we did not find evidence that supported the established domestication model for O. glaberrima—a single domestication origin. Rather, our evidence suggests the domestication process for African rice was initiated in multiple regions of West Africa, caused potentially by the local environments and cultivation preference of people. Hence domestication of African rice was a multi-regional process.

Codweb: Whole-genome sequencing uncovers extensive reticulations fueling adaptation among Atlantic, Arctic, and Pacific gadids

Introgressive hybridization creates networks of genetic relationships across species. Among marine fish of the Gadidae family, Pacific cod and walleye pollock are separate invasions of an Atlantic cod ancestor into the Pacific. Cods are ecological success stories, and their ecologies allow them to support the largest fisheries of the world. The enigmatic walleye pollock differs morphologically, behaviorally, and ecologically from its relatives, representing a niche shift. Here, we apply whole-genome sequencing to Pacific, Arctic, and Atlantic gadids and reveal extensive introgression among them with the ABBA-BABA test and pseudolikelihood phylogenetic network analysis. We propose that walleye pollock resulted from extensive adaptive introgression or homoploid hybrid speciation. The path of evolution of these taxa is more web than a tree. Their ability to invade and expand into new habitats and become ecologically successful may depend on genes acquired through adaptive introgression or hybrid speciation.

Phenotypic and genetic diversity in aposematic Malagasy poison frogs (genus Mantella)

Intraspecific color variation has long fascinated evolutionary biologists. In species with bright warning coloration, phenotypic diversity is particularly compelling because many factors, including natural and sexual selection, contribute to intraspecific variation. To better understand the causes of dramatic phenotypic variation in Malagasy poison frogs, we quantified genetic structure and color and pattern variation across three closely related species, Mantella aurantiaca, Mantella crocea, and Mantella milotympanum. Although our restriction site-associated DNA (RAD) sequencing approach identified clear genetic clusters, they do not align with current species designations, which has important conservation implications for these imperiled frogs. Moreover, our results suggest that levels of intraspecific color variation within this group have been overestimated, while species diversity has been underestimated. Within major genetic clusters, we observed distinct patterns of variation including: populations that are phenotypically similar yet genetically distinct, populations where phenotypic and genetic breaks coincide, and populations that are genetically similar but have high levels of within-population phenotypic variation. We also detected admixture between two of the major genetic clusters. Our study suggests that several mechanisms-including hybridization, selection, and drift-are contributing to phenotypic diversity. Ultimately, our work underscores the need for a reevaluation of how polymorphic and polytypic populations and species are classified, especially in aposematic organisms.

MobiSeq: De novo SNP discovery in model and non-model species through sequencing the flanking region of transposable elements

In recent years, the availability of reduced representation library (RRL) methods has catalysed an expansion of genome-scale studies to characterize both model and non-model organisms. Most of these methods rely on the use of restriction enzymes to obtain DNA sequences at a genome-wide level. These approaches have been widely used to sequence thousands of markers across individuals for many organisms at a reasonable cost, revolutionizing the field of population genomics. However, there are still some limitations associated with these methods, in particular the high molecular weight DNA required as starting material, the reduced number of common loci among investigated samples, and the short length of the sequenced site-associated DNA. Here, we present MobiSeq, a RRL protocol exploiting simple laboratory techniques, that generates genomic data based on PCR targeted enrichment of transposable elements and the sequencing of the associated flanking region. We validate its performance across 103 DNA extracts derived from three mammalian species: grey wolf (Canis lupus), red deer complex (Cervus sp.) and brown rat (Rattus norvegicus). MobiSeq enables the sequencing of hundreds of thousands loci across the genome and performs SNP discovery with relatively low rates of clonality. Given the ease and flexibility of MobiSeq protocol, the method has the potential to be implemented for marker discovery and population genomics across a wide range of organisms-enabling the exploration of diverse evolutionary and conservation questions.

Current geography masks dynamic history of gene flow during speciation in northern Australian birds

Genome divergence is greatly influenced by gene flow during early stages of speciation. As populations differentiate, geographic barriers can constrain gene flow and so affect the dynamics of divergence and speciation. Current geography, specifically disjunction and continuity of ranges, is often used to predict the historical gene flow during the divergence process. We test this prediction in eight meliphagoid bird species complexes codistributed in four regions. These regions are separated by known biogeographical barriers across northern Australia and Papua New Guinea. We find that bird populations currently separated by terrestrial habitat barriers within Australia and marine barriers between Australia and Papua New Guinea have a range of divergence levels and probability of gene flow not associated with current range connectivity. Instead, geographic distance and historical range connectivity better predict divergence and probability of gene flow. In this dynamic environmental context, we also find support for a nonlinear decrease of the probability of gene flow during the divergence process. The probability of gene flow initially decreases gradually after a certain level of divergence is reached. Its decrease then accelerates until the probability is close to zero. This implies that although geographic connectivity may have more of an effect early in speciation, other factors associated with higher divergence may play a more important role in influencing gene flow midway through and later in speciation. Current geographic connectivity may then mislead inferences regarding potential for gene flow during speciation under a complex and dynamic history of geographic and reproductive isolation.

Multiple Origin but Single Domestication Led to Oryza sativa

The domestication scenario that led to Asian rice (Oryza sativa) is a contentious topic. Here, we have reanalyzed a previously published large-scale wild and domesticated rice data set, which was also analyzed by two studies but resulted in two contrasting domestication models. We suggest that the analysis of false-positive selective sweep regions and phylogenetic analysis of concatenated genomic regions may have been the sources that contributed to the different results. In the end, our result indicates that Asian rice originated from multiple wild progenitor subpopulations; however, de novo domestication appears to have occurred only once and the domestication alleles were transferred between rice subpopulations through introgression.

Ancient genomes revisit the ancestry of domestic and Przewalski’s horses

The Eneolithic Botai culture of the Central Asian steppes provides the earliest archaeological evidence for horse husbandry, ~5500 years ago, but the exact nature of early horse domestication remains controversial. We generated 42 ancient-horse genomes, including 20 from Botai. Compared to 46 published ancient- and modern-horse genomes, our data indicate that Przewalski’s horses are the feral descendants of horses herded at Botai and not truly wild horses. All domestic horses dated from ~4000 years ago to present only show ~2.7% of Botai-related ancestry. This indicates that a massive genomic turnover underpins the expansion of the horse stock that gave rise to modern domesticates, which coincides with large-scale human population expansions during the Early Bronze Age.

Ancient genomic changes associated with domestication of the horse

The domestication of the horse was a seminal event in human cultural evolution. Libradoet al.obtained genome sequences from 14 horses from the Bronze and Iron Ages, about 2000 to 4000 years ago, soon after domestication. They identified variants determining coat color and genes selected during the domestication process. They could also see evidence of admixture with archaic horses and the demography of the domestication process, which included the accumulation of deleterious variants. The horse appears to have undergone a different type of domestication process than animals that were domesticated simply for food.

Science, this issue p.442