The fitness of an introgressing haplotype changes over the course of divergence and depends on its size and genomic location.

The genomic era has made clear that introgression, or the movement of genetic material between species, is a common feature of evolution. Examples of both adaptive and deleterious introgression exist in a variety of systems. What is unclear is how the fitness of an introgressing haplotype changes as species diverge or as the size of the introgressing haplotype changes. In a simple model, we show that introgression may more easily occur into parts of the genome which have not diverged heavily from a common ancestor. The key insight is that alleles from a shared genetic background are likely to have positive epistatic interactions, increasing the fitness of a larger introgressing block. In regions of the genome where few existing substitutions are disrupted, this positive epistasis can be larger than incompatibilities with the recipient genome. Further, we show that early in the process of divergence, introgression of large haplotypes can be favored more than introgression of individual alleles. This model is consistent with observations of a positive relationship between recombination rate and introgression frequency across the genome; however, it generates several novel predictions. First, the model suggests that the relationship between recombination rate and introgression may not exist, or may be negative, in recently diverged species pairs. Furthermore, the model suggests that introgression that replaces existing derived variation will be more deleterious than introgression at sites carrying ancestral variants. These predictions are tested in an example of introgression in Drosophila melanogaster, with some support for both. Finally, the model provides a potential alternative explanation to asymmetry in the direction of introgression, with expectations of higher introgression from rapidly diverged populations into slowly evolving ones.

Strong Positive Selection in Aedes aegypti and the Rapid Evolution of Insecticide Resistance.

Aedes aegypti vectors the pathogens that cause dengue, yellow fever, Zika virus, and chikungunya and is a serious threat to public health in tropical regions. Decades of work has illuminated many aspects of Ae. aegypti's biology and global population structure and has identified insecticide resistance genes; however, the size and repetitive nature of the Ae. aegypti genome have limited our ability to detect positive selection in this mosquito. Combining new whole genome sequences from Colombia with publicly available data from Africa and the Americas, we identify multiple strong candidate selective sweeps in Ae. aegypti, many of which overlap genes linked to or implicated in insecticide resistance. We examine the voltage-gated sodium channel gene in three American cohorts and find evidence for successive selective sweeps in Colombia. The most recent sweep encompasses an intermediate-frequency haplotype containing four candidate insecticide resistance mutations that are in near-perfect linkage disequilibrium with one another in the Colombian sample. We hypothesize that this haplotype may continue to rapidly increase in frequency and perhaps spread geographically in the coming years. These results extend our knowledge of how insecticide resistance has evolved in this species and add to a growing body of evidence suggesting that Ae. aegypti has an extensive genomic capacity to rapidly adapt to insecticide-based vector control.

Little effect of the tan locus on pigmentation in female hybrids between Drosophila santomea and D. melanogaster.

Previous work on Drosophila santomea suggested that its absence of abdominal pigmentation, compared to the other darkly pigmented species, is based on mutations in the cis-regulatory region of tan, inactivating the expression of that gene in the abdomen of D. santomea males and females. Our discovery that D. santomea males can produce viable hybrids when mated to D. melanogaster females enables us to use the armamentarium of genetic tools in the latter species to study the genetic basis of this interspecific difference in pigmentation. Hybridization tests using D. melanogaster deficiencies that include tan show no evidence that this locus is involved in the lighter pigmentation of D. santomea females; rather, the pigmentation difference appears to involve at least four other loci in the region. Earlier results implicating tan may have been based on a type of transgenic analysis that can give misleading results about the genes involved in an evolutionary change.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations.

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Isotopic niches do not follow the expectations of niche conservatism in the bird genus Cinclodes.

Phenotypic traits are expected to be more similar among closely related species than among species that diverged long ago (all else being equal). This pattern, known as phylogenetic niche conservatism, also applies to traits that are important to determine the niche of species. To test this hypothesis on ecological niches, we analysed isotopic data from 254 museum study skins from 12 of the 16 species of the bird genus Cinclodes and measured stable isotope ratios for four different elements: carbon, nitrogen, hydrogen and oxygen. We find that all traits, measured individually, or as a composite measurement, lack any phylogenetic signal, which in turn suggests a high level of lability in ecological niches. We compared these metrics to the measurements of morphological traits in the same genus and found that isotopic niches are uniquely evolutionarily labile compared to other traits. Our results suggest that, in Cinclodes, the realized niche evolves much faster than expected by the constraints of phylogenetic history and poses the question of whether this is a general pattern across the tree of life.

An Indian lineage of Histoplasma with strong signatures of differentiation and selection.

Histoplasma, a genus of dimorphic fungi, is the etiological agent of histoplasmosis, a pulmonary disease widespread across the globe. Whole genome sequencing has revealed that the genus harbors a previously unrecognized diversity of cryptic species. To date, studies have focused on Histoplasma isolates collected in the Americas with little knowledge of the genomic variation from other localities. In this report, we report the existence of a well-differentiated lineage of Histoplasma occurring in the Indian subcontinent. The group is differentiated enough to satisfy the requirements of a phylogenetic species, as it shows extensive genetic differentiation along the whole genome and has little evidence of gene exchange with other Histoplasma species. Next, we leverage this genetic differentiation to identify genetic changes that are unique to this group and that have putatively evolved through rapid positive selection. We found that none of the previously known virulence factors have evolved rapidly in the Indian lineage but find evidence of strong signatures of selection on other alleles potentially involved in clinically-important phenotypes. Our work serves as an example of the importance of correctly identifying species boundaries to understand the extent of selection in the evolution of pathogenic lineages. IMPORTANCE: Whole genome sequencing has revolutionized our understanding of microbial diversity, including human pathogens. In the case of fungal pathogens, a limiting factor in understanding the extent of their genetic diversity has been the lack of systematic sampling. In this piece, we show the results of a collection in the Indian subcontinent of the pathogenic fungus Histoplasma, the causal agent of a systemic mycosis. We find that Indian samples of Histoplasma form a distinct clade which is highly differentiated from other Histoplasma species. We also show that the genome of this lineage shows unique signals of natural selection. This work exemplifies how the combination of a robust sampling along with population genetics, and phylogenetics can reveal the precise genetic changes that differentiate lineages of fungal pathogens.

Single nucleotide polymorphisms and chromosomal copy number variation may impact the Sporothrix brasiliensis antifungal susceptibility and sporotrichosis clinical outcomes.

Feline-transmitted sporotrichosis has garnered attention due to the recent high incidence and the lack of efficient control in the epicenter of the epidemic, Rio de Janeiro, Brazil. Sporothrix brasiliensis is the major pathogen involved in feline-to-human sporotrichosis in Brazil and displays more virulent genotypes than the closely related species S. schenckii. Over the last two decades, several reports of antifungal-resistant strains have emerged. Sequencing and comparison analysis of the outbreak strains allowed us to observe that the azole non-wild-type S. brasiliensis strain CFP 1054 had significant chromosomal variations compared to wild-type strains. One of these variants includes a region of 231 Kb containing 75 duplicated genes, which were overrepresented for lipid and isoprenoid metabolism. We also identified an additional strain (CFP 1055) that was resistant to itraconazole and amphotericin B, which had a single nucleotide polymorphism in the tac1 gene. The patients infected with these two strains showed protracted clinical course and sequelae. Even though our sample size is modest, these results suggest the possibility of identifying specific point mutations and large chromosomal duplications potentially associated with antifungal resistance and clinical outcomes of sporotrichosis.

Genetic Diversity and Thermal Performance in Invasive and Native Populations of African Fig Flies.

During biological invasions, invasive populations can suffer losses of genetic diversity that are predicted to negatively impact their fitness/performance. Despite examples of invasive populations harboring lower diversity than conspecific populations in their native range, few studies have linked this lower diversity to a decrease in fitness. Using genome sequences, we show that invasive populations of the African fig fly, Zaprionus indianus, have less genetic diversity than conspecific populations in their native range and that diversity is proportionally lower in regions of the genome experiencing low recombination rates. This result suggests that selection may have played a role in lowering diversity in the invasive populations. We next use interspecific comparisons to show that genetic diversity remains relatively high in invasive populations of Z. indianus when compared with other closely related species. By comparing genetic diversity in orthologous gene regions, we also show that the genome-wide landscape of genetic diversity differs between invasive and native populations of Z. indianus indicating that invasion not only affects amounts of genetic diversity but also how that diversity is distributed across the genome. Finally, we use parameter estimates from thermal performance curves for 13 species of Zaprionus to show that Z. indianus has the broadest thermal niche of measured species, and that performance does not differ between invasive and native populations. These results illustrate how aspects of genetic diversity in invasive species can be decoupled from measures of fitness, and that a broad thermal niche may have helped facilitate Z. indianus's range expansion.

Recurrent Collection of Drosophila melanogaster from Wild African Environments and Genomic Insights into Species History.

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

Temperature-Dependent Competitive Outcomes between the Fruit Flies Drosophila santomea and Drosophila yakuba.

AbstractChanges in temperature associated with climate change can alter species' distributions, drive adaptive evolution, and in some cases cause extinction. Research has tended to focus on the direct effects of temperature, but changes in temperature can also have indirect effects on populations and species. Here, we test whether temperature can indirectly affect the fitness of Drosophila santomea and Drosophila yakuba by altering the nature of interspecific competition. We show that when raised in isolation, both D. santomea and D. yakuba display similar variation in relative fitness across temperatures of 18°, 22°, and 25°C. However, D. santomea has higher fitness than D. yakuba when experiencing interspecific competition at 18°C, while the inverse is true at 25°C. Patterns of fitness across thermal and competitive environments therefore indicate that the outcome of interspecific competition varies with temperature. We then use a coexistence experiment to show that D. santomea is rapidly (within eight generations) extirpated when maintained with D. yakuba at 25°C. By contrast, D. santomea remains as (or more) abundant than D. yakuba over the course of ∼10 generations when maintained at 18°C. Our results provide an example of how the thermal environment can affect interspecific competition and suggest that some species may become more prone to extinction under scenarios of climate change through indirect effects of the thermal environment on competitive advantages between species.

Genomic signatures of admixture and selection are shared among populations of Zaprionus indianus across the western hemisphere.

Introduced species have become an increasingly common component of biological communities around the world. A central goal in invasion biology is therefore to identify the demographic and evolutionary factors that underlie successful introductions. Here we use whole genome sequences, collected from populations in the native and introduced range of the African fig fly, Zaprionus indianus, to quantify genetic relationships among them, identify potential sources of the introductions, and test for selection at different spatial scales. We find that geographically widespread populations in the western hemisphere are genetically more similar to each other than to lineages sampled across Africa, and that these populations share a mixture of alleles derived from differentiated African lineages. Using patterns of allele-sharing and demographic modelling we show that Z. indinaus have undergone a single expansion across the western hemisphere with admixture between African lineages predating this expansion. We also find support for selection that is shared across populations in the western hemisphere, and in some cases, with a subset of African populations. This suggests either that parallel selection has acted across a large part of Z. indianus's introduced range; or, more parsimoniously, that Z. indianus has experienced selection early on during (or prior-to) its expansion into the western hemisphere. We suggest that the range expansion of Z. indianus has been facilitated by admixture and selection, and that management of this invasion could focus on minimizing future admixture by controlling the movement of individuals within this region rather than between the western and eastern hemisphere.

Genetic divergence and the number of hybridizing species affect the path to homoploid hybrid speciation.

Hybridization is often maladaptive and in some instances has led to the loss of biodiversity. However, hybridization can also promote speciation, such as during homoploid hybrid speciation, thereby generating biodiversity. Despite examples of homoploid hybrid species, the importance of hybridization as a speciation mechanism is still widely debated, and we lack a general understanding of the conditions most likely to generate homoploid hybrid species. Here we show that the level of genetic divergence between hybridizing species has a large effect on the probability that their hybrids evolve reproductive isolation. We find that populations of hybrids formed by parental species with intermediate levels of divergence were more likely to mate assortatively, and discriminate against their parental species, than those generated from weakly or strongly diverged parental species. Reproductive isolation was also found between hybrid populations, suggesting differential sorting of parental traits across populations. Finally, hybrid populations derived from three species were more likely to evolve reproductive isolation than those derived from two species, supporting arguments that hybridization-supplied genetic diversity can lead to the evolution of novel "adaptive systems" and promote speciation. Our results illustrate when we expect hybridization and admixture to promote hybrid speciation. Whether homoploid hybrid speciation is a common speciation mechanism in general remains an outstanding empirical question.

Supervised machine learning reveals introgressed loci in the genomes of Drosophila simulans and D. sechellia.

Hybridization and gene flow between species appears to be common. Even though it is clear that hybridization is widespread across all surveyed taxonomic groups, the magnitude and consequences of introgression are still largely unknown. Thus it is crucial to develop the statistical machinery required to uncover which genomic regions have recently acquired haplotypes via introgression from a sister population. We developed a novel machine learning framework, called FILET (Finding Introgressed Loci via Extra-Trees) capable of revealing genomic introgression with far greater power than competing methods. FILET works by combining information from a number of population genetic summary statistics, including several new statistics that we introduce, that capture patterns of variation across two populations. We show that FILET is able to identify loci that have experienced gene flow between related species with high accuracy, and in most situations can correctly infer which population was the donor and which was the recipient. Here we describe a data set of outbred diploid Drosophila sechellia genomes, and combine them with data from D. simulans to examine recent introgression between these species using FILET. Although we find that these populations may have split more recently than previously appreciated, FILET confirms that there has indeed been appreciable recent introgression (some of which might have been adaptive) between these species, and reveals that this gene flow is primarily in the direction of D. simulans to D. sechellia.

Genomic diversity of the human pathogen Paracoccidioides across the South American continent.

Paracoccidioidomycosis (PCM) is a life-threatening systemic mycosis widely reported in the Gran Chaco ecosystem. The disease is caused by different species from the genus Paracoccidioides, which are all endemic to South and Central America. Here, we sequenced and analyzed 31 isolates of Paracoccidioides across South America, with particular focus on isolates from Argentina and Paraguay. The de novo sequenced isolates were compared with publicly available genomes. Phylogenetics and population genomics revealed that PCM in Argentina and Paraguay is caused by three distinct Paracoccidioides genotypes, P. brasiliensis (S1a and S1b) and P. restrepiensis (PS3). P. brasiliensis S1a isolates from Argentina are frequently associated with chronic forms of the disease. Our results suggest the existence of extensive molecular polymorphism among Paracoccidioides species, and provide a framework to begin to dissect the connection between genotypic differences in the pathogen and the clinical outcomes of the disease.

Fine scale mapping of genomic introgressions within the Drosophila yakuba clade.

The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. However, the amount of interspecific genetic exchange (introgression) mediated by hybridization remains largely unknown, although recent progress in genome sequencing has made measuring introgression more tractable. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of São Tomé, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. In this report, we quantify the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). We sequenced the genomes of 86 individuals from all three species. We also developed and applied a new statistical framework, using a hidden Markov approach, to identify introgression. We found that introgression has occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, we find that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). We additionally show that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, our results indicate that genetic exchange between these species is currently rare.

The Rate of Evolution of Postmating-Prezygotic Reproductive Isolation in Drosophila.

Reproductive isolation is an intrinsic aspect of species formation. For that reason, the identification of the precise isolating traits, and the rates at which they evolve, is crucial to understanding how species originate and persist. Previous work has measured the rates of evolution of prezygotic and postzygotic barriers to gene flow, yet no systematic analysis has studied the rates of evolution of postmating-prezygotic (PMPZ) barriers. We measured the magnitude of two barriers to gene flow that act after mating occurs but before fertilization. We also measured the magnitude of a premating barrier (female mating rate in nonchoice experiments) and two postzygotic barriers (hybrid inviability and hybrid sterility) for all pairwise crosses of all nine known extant species within the melanogaster subgroup. Our results indicate that PMPZ isolation evolves faster than hybrid inviability but slower than premating isolation. Next, we partition postzygotic isolation into different components and find that, as expected, hybrid sterility evolves faster than hybrid inviability. These results lend support for the hypothesis that, in Drosophila, reproductive isolation mechanisms (RIMs) that act early in reproduction (or in development) tend to evolve faster than those that act later in the reproductive cycle. Finally, we tested whether there was evidence for reinforcing selection at any RIM. We found no evidence for generalized evolution of reproductive isolation via reinforcement which indicates that there is no pervasive evidence of this evolutionary process. Our results indicate that PMPZ RIMs might have important evolutionary consequences in initiating speciation and in the persistence of new species.

Fungal species boundaries in the genomics era.

Genomic data has opened new possibilities to understand how organisms change over time, and could enable the discovery of previously undescribed species. Although taxonomy used to be based on phenotypes, molecular data has frequently revealed that morphological traits are insufficient to describe biodiversity. Genomics holds the promise of revealing even more genetic discontinuities, but the parameters on how to describe species from genomic data remain unclear. Fungi have been a successful case in which the use of molecular markers has uncovered the existence of genetic boundaries where no crosses are possible. In this minireview, we highlight recent advances, propose a set of standards to use genomic sequences to uncover species boundaries, point out potential pitfalls, and present possible future research directions.

Species boundaries in the human pathogen Paracoccidioides.

The use of molecular taxonomy for identifying recently diverged species has transformed the study of speciation in fungi. The pathogenic fungus Paracoccidioides spp has been hypothesized to be composed of five phylogenetic species, four of which compose the brasiliensis species complex. Nuclear gene genealogies support this divergence scenario, but mitochondrial loci do not; while all species from the brasiliensis complex are differentiated at nuclear coding loci, they are not at mitochondrial loci. We addressed the source of this incongruity using 11 previously published gene fragments, 10 newly-sequenced nuclear non-coding loci, and 10 microsatellites. We hypothesized and further demonstrated that the mito-nuclear incongruence in the brasiliensis species complex results from interspecific hybridization and mitochondrial introgression, a common phenomenon in eukaryotes. Additional population genetic analyses revealed possible nuclear introgression but much less than that seen in the mitochondrion. Our results are consistent with a divergence scenario of secondary contact and subsequent mitochondrial introgression despite the continued persistence of species boundaries. We also suggest that yeast morphology slightly-but significantly-differs across all five Paracoccidioides species and propose to elevate four of these phylogenetic species to formally described taxonomic species.

Fine mapping of dominant X-linked incompatibility alleles in Drosophila hybrids.

Sex chromosomes have a large effect on reproductive isolation and play an important role in hybrid inviability. In Drosophila hybrids, X-linked genes have pronounced deleterious effects on fitness in male hybrids, which have only one X chromosome. Several studies have succeeded at locating and identifying recessive X-linked alleles involved in hybrid inviability. Nonetheless, the density of dominant X-linked alleles involved in interspecific hybrid viability remains largely unknown. In this report, we study the effects of a panel of small fragments of the D. melanogaster X-chromosome carried on the D. melanogaster Y-chromosome in three kinds of hybrid males: D. melanogaster/D. santomea, D. melanogaster/D. simulans and D. melanogaster/D. mauritiana. D. santomea and D. melanogaster diverged over 10 million years ago, while D. simulans (and D. mauritiana) diverged from D. melanogaster over 3 million years ago. We find that the X-chromosome from D. melanogaster carries dominant alleles that are lethal in mel/san, mel/sim, and mel/mau hybrids, and more of these alleles are revealed in the most divergent cross. We then compare these effects on hybrid viability with two D. melanogaster intraspecific crosses. Unlike the interspecific crosses, we found no X-linked alleles that cause lethality in intraspecific crosses. Our results reveal the existence of dominant alleles on the X-chromosome of D. melanogaster which cause lethality in three different interspecific hybrids. These alleles only cause inviability in hybrid males, yet have little effect in hybrid females. This suggests that X-linked elements that cause hybrid inviability in males might not do so in hybrid females due to differing sex chromosome interactions.