Interspecific introgressive origin of genomic diversity in the house mouse

Phylogenetic incongruence and the evolutionary origins of cardenolide-resistant forms of Na(+) ,K(+) -ATPase in Danaus butterflies

Many distantly related insect species are specialized feeders of cardenolide-containing host plants such as milkweed (Asclepias spp.). Previous studies have revealed frequent, parallel substitution of a functionally important amino acid substitution (N122H) in the alpha subunit of Na(+) ,K(+) -ATPase in a number of these species. This substitution facilitates the ability of these insects to feed on their toxic hosts and sequester cardenolides for their own use in defense. Among milkweed butterflies of the genus Danaus, the previously established phylogeny for this group suggests that N122H arose independently and fixed in two distinct lineages. We reevaluate this conclusion by examining Danaus phylogenetic relationships using >400 orthologous gene sequences assembled from transcriptome data. Our results indicate that the three Danaus species known to harbor the N122H substitution are more closely related than previously thought, consistent with a single, common origin for N122H. However, we also find evidence of both incomplete lineage sorting and post-speciation genetic exchange among these butterfly species, raising the possibility of collateral evolution of cardenolide-insensitivity in this species group.

A genomic perspective on hybridization and speciation

Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.

A beak size locus in Darwin's finches facilitated character displacement during a drought

Ecological character displacement is a process of morphological divergence that reduces competition for limited resources. We used genomic analysis to investigate the genetic basis of a documented character displacement event in Darwin's finches on Daphne Major in the Galápagos Islands: The medium ground finch diverged from its competitor, the large ground finch, during a severe drought. We discovered a genomic region containing the HMGA2 gene that varies systematically among Darwin's finch species with different beak sizes. Two haplotypes that diverged early in the radiation were involved in the character displacement event: Genotypes associated with large beak size were at a strong selective disadvantage in medium ground finches (selection coefficient s = 0.59). Thus, a major locus has apparently facilitated a rapid ecological diversification in the adaptive radiation of Darwin's finches.

Reticulate evolutionary history and extensive introgression in mosquito species revealed by phylogenetic network analysis

The role of hybridization and subsequent introgression has been demonstrated in an increasing number of species. Recently, Fontaine et al. (Science, 347, 2015, 1258524) conducted a phylogenomic analysis of six members of the Anopheles gambiae species complex. Their analysis revealed a reticulate evolutionary history and pointed to extensive introgression on all four autosomal arms. The study further highlighted the complex evolutionary signals that the co-occurrence of incomplete lineage sorting (ILS) and introgression can give rise to in phylogenomic analyses. While tree-based methodologies were used in the study, phylogenetic networks provide a more natural model to capture reticulate evolutionary histories. In this work, we reanalyse the Anopheles data using a recently devised framework that combines the multispecies coalescent with phylogenetic networks. This framework allows us to capture ILS and introgression simultaneously, and forms the basis for statistical methods for inferring reticulate evolutionary histories. The new analysis reveals a phylogenetic network with multiple hybridization events, some of which differ from those reported in the original study. To elucidate the extent and patterns of introgression across the genome, we devise a new method that quantifies the use of reticulation branches in the phylogenetic network by each genomic region. Applying the method to the mosquito data set reveals the evolutionary history of all the chromosomes. This study highlights the utility of 'network thinking' and the new insights it can uncover, in particular in phylogenomic analyses of large data sets with extensive gene tree incongruence.

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery-reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996-2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make-up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.

Inferring Orthologs: Open Questions and Perspectives

With the increasing number of sequenced genomes and their comparisons, the detection of orthologs is crucial for reliable functional annotation and evolutionary analyses of genes and species. Yet, the dynamic remodeling of genome content through gain, loss, transfer of genes, and segmental and whole-genome duplication hinders reliable orthology detection. Moreover, the lack of direct functional evidence and the questionable quality of some available genome sequences and annotations present additional difficulties to assess orthology. This article reviews the existing computational methods and their potential accuracy in the high-throughput era of genome sequencing and anticipates open questions in terms of methodology, reliability, and computation. Appropriate taxon sampling together with combination of methods based on similarity, phylogeny, synteny, and evolutionary knowledge that may help detecting speciation events appears to be the most accurate strategy. This review also raises perspectives on the potential determination of orthology throughout the whole species phylogeny.

Mapping the genomic architecture of adaptive traits with interspecific introgressive origin: a coalescent-based approach

Recent studies of eukaryotes including human and Neandertal, mice, and butterflies have highlighted the major role that interspecific introgression has played in adaptive trait evolution. A common question arises in each case: what is the genomic architecture of the introgressed traits? One common approach that can be used to address this question is association mapping, which looks for genotypic markers that have significant statistical association with a trait. It is well understood that sample relatedness can be a confounding factor in association mapping studies if not properly accounted for. Introgression and other evolutionary processes (e.g., incomplete lineage sorting) typically introduce variation among local genealogies, which can also differ from global sample structure measured across all genomic loci. In contrast, state-of-the-art association mapping methods assume fixed sample relatedness across the genome, which can lead to spurious inference. We therefore propose a new association mapping method called Coal-Map, which uses coalescent-based models to capture local genealogical variation alongside global sample structure. Using simulated and empirical data reflecting a range of evolutionary scenarios, we compare the performance of Coal-Map against EIGENSTRAT, a leading association mapping method in terms of its popularity, power, and type I error control. Our empirical data makes use of hundreds of mouse genomes for which adaptive interspecific introgression has recently been described. We found that Coal-Map's performance is comparable or better than EIGENSTRAT in terms of statistical power and false positive rate. Coal-Map's performance advantage was greatest on model conditions that most closely resembled empirically observed scenarios of adaptive introgression. These conditions had: (1) causal SNPs contained in one or a few introgressed genomic loci and (2) varying rates of gene flow - from high rates to very low rates where incomplete lineage sorting dominated as a primary cause of local genealogical variation.

Selective sweeps versus introgression - population genetic dynamics of the murine leukemia virus receptor Xpr1 in wild populations of the house mouse (Mus musculus)

Background

The interaction between viruses and their receptors in the host can be expected to lead to an evolutionary arms race resulting in cycles of rapid adaptations. We focus here on the receptor gene Xpr1 (xenotropic and polytropic retrovirus receptor 1) for murine leukemia viruses (MLVs). In a previous screen for selective sweeps in mouse populations we discovered that a population from Germany was almost monomorphic for Xpr1 haplotypes, while a population from France was polymorphic.

Results

Here we analyze Xpr1 sequences and haplotypes from a broad sample of wild mouse populations of two subspecies, M. m. domesticus and M. m. musculus, to trace the origins of this distinctive polymorphism pattern. We show that the high polymorphism in the population in France is caused by a relatively recent invasion of a haplotype from a population in Iran, rather than a selective sweep in Germany. The invading haplotype codes for a novel receptor variant, which has itself undergone a recent selective sweep in the Iranian population.

Conclusions

Our data support a scenario in which Xpr1 is frequently subject to positive selection, possibly as a response to resistance development against recurrently emerging infectious viruses. During such an infection cycle, receptor variants that may convey viral resistance can be captured from another population and quickly introgress into populations actively dealing with the infectious virus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0528-5) contains supplementary material, which is available to authorized users.

Divergence and gene flow among Darwin's finches: A genome-wide view of adaptive radiation driven by interspecies allele sharing

A recent analysis of the genomes of Darwin's finches revealed extensive interspecies allele sharing throughout the history of the radiation and identified a key locus responsible for morphological evolution in this group. The radiation of Darwin's finches on the Galápagos archipelago has long been regarded as an iconic study system for field ecology and evolutionary biology. Coupled with an extensive history of field work, these latest findings affirm the increasing acceptance of introgressive hybridization, or gene flow between species, as a significant contributor to adaptive evolution. Here, we review and discuss these findings in relation to both classical work on Darwin's finches and contemporary work showing similar evolutionary signatures in other biological systems. The continued unification of genomic data with field biology promises to further elucidate the molecular basis of adaptation in Darwin's finches and well beyond.

Gene flow between Drosophila yakuba and Drosophila santomea in subunit V of cytochrome c oxidase: A potential case of cytonuclear cointrogression

Introgression is the effective exchange of genetic information between species through natural hybridization. Previous genetic analyses of the Drosophila yakuba—D. santomea hybrid zone showed that the mitochondrial genome of D. yakuba had introgressed into D. santomea and completely replaced its native form. Since mitochondrial proteins work intimately with nuclear‐encoded proteins in the oxidative phosphorylation (OXPHOS) pathway, we hypothesized that some nuclear genes in OXPHOS cointrogressed along with the mitochondrial genome. We analyzed nucleotide variation in the 12 nuclear genes that form cytochrome c oxidase (COX) in 33 Drosophila lines. COX is an OXPHOS enzyme composed of both nuclear‐ and mitochondrial‐encoded proteins and shows evidence of cytonuclear coadaptation in some species. Using maximum‐likelihood methods, we detected significant gene flow from D. yakuba to D. santomea for the entire COX complex. Interestingly, the signal of introgression is concentrated in the three nuclear genes composing subunit V, which shows population migration rates significantly greater than the background level of introgression in these species. The detection of introgression in three proteins that work together, interact directly with the mitochondrial‐encoded core, and are critical for early COX assembly suggests this could be a case of cytonuclear cointrogression.

Insights into mammalian biology from the wild house mouse Mus musculus

The house mouse, Mus musculus, was established in the early 1900s as one of the first genetic model organisms owing to its short generation time, comparatively large litters, ease of husbandry, and visible phenotypic variants. For these reasons and because they are mammals, house mice are well suited to serve as models for human phenotypes and disease. House mice in the wild consist of at least three distinct subspecies and harbor extensive genetic and phenotypic variation both within and between these subspecies. Wild mice have been used to study a wide range of biological processes, including immunity, cancer, male sterility, adaptive evolution, and non-Mendelian inheritance. Despite the extensive variation that exists among wild mice, classical laboratory strains are derived from a limited set of founders and thus contain only a small subset of this variation. Continued efforts to study wild house mice and to create new inbred strains from wild populations have the potential to strengthen house mice as a model system.

DOI:http://dx.doi.org/10.7554/eLife.05959.001