Tracing histoplasmosis genomic epidemiology and species occurrence across the USA

ABSTRACTHistoplasmosis is an endemic mycosis in North America frequently reported along the Ohio and Mississippi River Valleys, although autochthonous cases occur in non-endemic areas. In the United States, the disease is provoked by two genetically distinct clades of Histoplasma capsulatum sensu lato, Histoplasma mississippiense (Nam1) and H. ohiense (Nam2). To bridge the molecular epidemiological gap, we genotyped 93 Histoplasma isolates (62 novel genomes) including clinical, environmental, and veterinarian samples from a broader geographical range by whole-genome sequencing, followed by evolutionary and species niche modelling analyses. We show that histoplasmosis is caused by two major lineages, H. ohiense and H. mississippiense; with sporadic cases caused by H. suramericanum in California and Texas. While H. ohiense is prevalent in eastern states, H. mississipiense was found to be prevalent in the central and western portions of the United States, but also geographically overlapping in some areas suggesting that these species might co-occur. Species Niche Modelling revealed that H. ohiense thrives in places with warmer and drier conditions, while H. mississippiense is endemic to areas with cooler temperatures and more precipitation. In addition, we predicted multiple areas of secondary contact zones where the two species co-occur, potentially facilitating gene exchange and hybridization. This study provides the most comprehensive understanding of the genomic epidemiology of histoplasmosis in the USA and lays a blueprint for the study of invasive fungal diseases.

Demographic history and the efficacy of selection in the globally invasive mosquito Aedes aegypti

Aedes aegypti is the main vector species of yellow fever, dengue, zika and chikungunya. The species is originally from Africa but has experienced a spectacular expansion in its geographic range to a large swath of the world, the demographic effects of which have remained largely understudied. In this report, we examine whole-genome sequences from 6 countries in Africa, North America, and South America to investigate the demographic history of the spread of Ae. aegypti into the Americas its impact on genomic diversity. In the Americas, we observe patterns of strong population structure consistent with relatively low (but probably non-zero) levels of gene flow but occasional long-range dispersal and/or recolonization events. We also find evidence that the colonization of the Americas has resulted in introduction bottlenecks. However, while each sampling location shows evidence of a past population contraction and subsequent recovery, our results suggest that the bottlenecks in America have led to a reduction in genetic diversity of only ~35% relative to African populations, and the American samples have retained high levels of genetic diversity (expected heterozygosity of ~0.02 at synonymous sites) and have experienced only a minor reduction in the efficacy of selection. These results evoke the image of an invasive species that has expanded its range with remarkable genetic resilience in the face of strong eradication pressure.

Aedes albopictus is present in the lowlands of southern Zambia

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and three rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on mtDNA barcoding. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to understanding the potential public health implications of this disease vector in southern Africa.

The spread of Aedes albopictus (Diptera: Culicidae) in the islands of São Tomé and Príncipe

The mosquito Aedes albopictus (Diptera: Culicidae) is a vector species of the causal agents of Dengue, yellow fever, and Zika among other diseases pathogens. The species originated in Southeast Asia and has spread widely and rapidly in the last century. The species has been reported in localities from the Gulf of Guinea since the early 2000s, but systematic sampling has been scant. We sampled Ae. albopictus twice, in 2013 and 2023 across the altitudinal gradient in São Tomé and found that the species was present in all sampled years at altitudes up to 680 m. We also found some evidence of increases in proportional representation compared to Ae. aegypti over time. We report the presence of the species in Príncipe for the first time, suggesting that the range of Ae. albopictus is larger than previously thought. Finally, we use bioclimatic niche modeling to infer the potential range of Ae. albopictus and infer that the species has the potential to spread across a large portion of São Tomé and Príncipe. Our results suggest that Ae. albopictus has established itself as a resident species of the islands of the Gulf of Guinea and should be incorporated into the list of potential vectors that need to be surveyed and controlled.

Phenotypic characterization of cryptic species in the fungal pathogen Histoplasma

Histoplasmosis is an endemic mycosis that often presents as a respiratory infection in immunocompromised patients. Hundreds of thousands of new infections are reported annually around the world. The etiological agent of the disease, Histoplasma, is a dimorphic fungus commonly found in the soil where it grows as mycelia. Humans can become infected by Histoplasma through inhalation of its spores (conidia) or mycelial particles. The fungi transitions into the yeast phase in the lungs at 37°C. Once in the lungs, yeast cells reside and proliferate inside alveolar macrophages. We have previously described that Histoplasma is composed of at least five cryptic species that differ genetically, and assigned new names to the lineages. Here we evaluated multiple phenotypic characteristics of 12 strains from five phylogenetic species of Histoplasma to identify phenotypic traits that differentiate between these species: H. capsulatum sensu stricto, H. ohiense, H. mississippiense, H. suramericanum, and an African lineage. We report diagnostic traits for two species. The other three species can be identified by a combination of traits. Our results suggest that 1) there are significant phenotypic differences among the cryptic species of Histoplasma, and 2) that those differences can be used to positively distinguish those species in a clinical setting and for further study of the evolution of this fungal pathogen.

Ancestral neural circuits potentiate the origin of a female sexual behavior

Courtship interactions are remarkably diverse in form and complexity among species. How neural circuits evolve to encode new behaviors that are functionally integrated into these dynamic social interactions is unknown. Here we report a recently originated female sexual behavior in the island endemic Drosophila species D. santomea, where females signal receptivity to male courtship songs by spreading their wings, which in turn promotes prolonged songs in courting males. Copulation success depends on this female signal and correlates with males' ability to adjust his singing in such a social feedback loop. Functional comparison of sexual circuitry across species suggests that a pair of descending neurons, which integrates male song stimuli and female internal state to control a conserved female abdominal behavior, drives wing spreading in D. santomea. This co-option occurred through the refinement of a pre-existing, plastic circuit that can be optogenetically activated in an outgroup species. Combined, our results show that the ancestral potential of a socially-tuned key circuit node to engage the wing motor program facilitates the expression of a new female behavior in appropriate sensory and motivational contexts. More broadly, our work provides insights into the evolution of social behaviors, particularly female behaviors, and the underlying neural mechanisms.

Aedes albopictus is present in the lowlands of southern Zambia

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and Zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and two rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two rural sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on CO1 sequence from mtDNA. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to the process of determining the potential public health implications of this disease vector in Central Africa.

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.

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.

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.

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.

A need for standardized reporting of introgression: Insights from studies across eukaryotes

With the rise of affordable next-generation sequencing technology, introgression-or the exchange of genetic materials between taxa-has become widely perceived to be a ubiquitous phenomenon in nature. Although this claim is supported by several keystone studies, no thorough assessment of the frequency of introgression across eukaryotes in nature has been performed to date. In this manuscript, we aim to address this knowledge gap by examining patterns of introgression across eukaryotes. We collated a single statistic, Patterson's D, which can be used as a test for introgression across 123 studies to further assess how taxonomic group, divergence time, and sequencing technology influence reports of introgression. Overall, introgression has mostly been measured in plants and vertebrates, with less attention given to the rest of the Eukaryotes. We find that the most frequently used metrics to detect introgression are difficult to compare across studies and even more so across biological systems due to differences in study effort, reporting standards, and methodology. Nonetheless, our analyses reveal several intriguing patterns, including the observation that differences in sequencing technologies may bias values of Patterson's D and that introgression may differ throughout the course of the speciation process. Together, these results suggest the need for a unified approach to quantifying introgression in natural communities and highlight important areas of future research that can be better assessed once this unified approach is met.

Phylogenetic climatic niche conservatism in sandflies (Diptera: Phlebotominae) and their relatives

Phylogenetic niche conservatism is a pattern in which closely related species are more similar than distant relatives in their niche-related traits. Species in the family Psychodidae show notable diversity in climatic niche, and present an opportunity to test for phylogenetic niche conservatism, which is as yet rarely studied in insects. Some species (in the subfamily Phlebotominae) transmit Leishmania parasites, responsible for the disease leishmaniasis, and their geographic range has been systematically characterized. Psychodid genus ranges can be solely tropical, confined to the temperate zones, or span both. We obtained observation site data, and associated climate data, for 234 psychodid species to understand which aspects of climate most closely predict distribution. Temperature and seasonality are strong determinants of species occurrence within the clade. Next, we built a phylogeny of Psychodidae, and found a positive relationship between pairwise genetic distance and climate niche differentiation, which indicates strong niche conservatism. This result is also supported by strong phylogenetic signals of metrics of climate differentiation. Finally, we used ancestral trait reconstruction to infer the tropicality (i.e., proportion of latitudinal range in the tropics minus the proportion of the latitudinal range in temperate areas) of ancestral species, and counted transitions to and from tropicality states. We find that tropical and temperate species produced almost entirely tropical and temperate descendant species, respectively. Taken together, our results imply that climate niches in psychodids are strongly predicted by phylogeny, and represent a formal test of a key prediction of phylogenetic niche conservatism in a clade with implications for human health.

Diversity of cave Phlebotomines (Diptera: Psychodidae) from a Colombian cave

Sandflies are vector species of Leishmania, among many other pathogens, with a global distribution and a variety of ecological niches. Previous samplings have found that karstic formations (i.e., caves, grottos, and folds formed by the erosion of limestone) serve as a natural habitat to sandfly species. The majority of samplings of cave sandfly diversity have occurred in Brazil and to date none have studied the species composition in a cave in the Northern Andes. We collected sandflies in the Cave "Los Guácharos", in the state of Antioquia, Colombia. The sampling was carried out during two consecutive nights in September 2019. CDC-type light traps were installed inside the cavern and in other surrounding karst systems (caves, rock-breaks, and folds). In total, we identified 17 species of sandfly from the cave and surrounding karst systems, including a new record for Colombia (Bichromomyia olmeca), and provide the first karstic reports for four other species (Lutzomyia gomezi, Lutzomyia hartmanni, Pintomyia ovallesi, and Psychodopygus panamensis). We then used the results of our survey and published literature to test two hypotheses. First, that sandfly diversity in Neotropical caves is richest nearer to the equator, and second that there is a phylogenetic signal of karstic habitat use in sandflies. Counter to our predictions, we found no evidence that diversity follows a latitudinal gradient. Further, we find no evidence of a phylogenetic signal of karstic habitat use, instead finding that the use of caves likely evolved multiple times across several genera. Our results highlight the importance of a wide sampling to understand the natural habitat of sandflies and other disease vectors.

The evolutionary trajectory of drosophilid walking

Neural circuits must both execute the behavioral repertoire of individuals and account for behavioral variation across species. Understanding how this variation emerges over evolutionary time requires large-scale phylogenetic comparisons of behavioral repertoires. Here, we describe the evolution of walking in fruit flies by capturing high-resolution, unconstrained movement from 13 species and 15 strains of drosophilids. We find that walking can be captured in a universal behavior space, the structure of which is evolutionarily conserved. However, the occurrence of and transitions between specific movements have evolved rapidly, resulting in repeated convergent evolution in the temporal structure of locomotion. Moreover, a meta-analysis demonstrates that many behaviors evolve more rapidly than other traits. Thus, the architecture and physiology of locomotor circuits can execute precise individual movements in one species and simultaneously support rapid evolutionary changes in the temporal ordering of these modular elements across clades.

Widespread introgression across a phylogeny of 155 Drosophila genomes

Genome-scale sequence data have invigorated the study of hybridization and introgression, particularly in animals. However, outside of a few notable cases, we lack systematic tests for introgression at a larger phylogenetic scale across entire clades. Here, we leverage 155 genome assemblies from 149 species to generate a fossil-calibrated phylogeny and conduct multilocus tests for introgression across 9 monophyletic radiations within the genus Drosophila. Using complementary phylogenomic approaches, we identify widespread introgression across the evolutionary history of Drosophila. Mapping gene-tree discordance onto the phylogeny revealed that both ancient and recent introgression has occurred across most of the 9 clades that we examined. Our results provide the first evidence of introgression occurring across the evolutionary history of Drosophila and highlight the need to continue to study the evolutionary consequences of hybridization and introgression in this genus and 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.

Genomic Diversity Analysis Reveals a Strong Population Structure in Histoplasma capsulatum LAmA (Histoplasma suramericanum)

Histoplasmosis is a severe mycotic disease affecting thousands of immunocompetent and immunocompromised individuals with high incidence in Latin America, where the disease agents are Histoplasma capsulatum and Histoplasma suramericanum. In this work, we used whole-genome sequencing to infer the species diversity and the population structure of H. suramericanum in South America. We find evidence for strong population structure and little admixture within the species. Genome-level phylogenetic trees indicate the existence of at least three different discrete populations. We recovered the existence of a previously identified population, LAmB, and confirm that it is highly differentiated along the whole genome. We also find that H. suramericanum is composed of two populations, one in Northern South America, and another in the southern portion of the continent. Moreover, one of the lineages from the southern population is endemic to Rio de Janeiro and there was no association with clinical data and species isolated from patients with histoplasmosis. Our results point out the need to characterize the symptomatology of histoplasmosis caused by different species and lineages of Histoplasma spp.

P-elements strengthen reproductive isolation within the Drosophila simulans species complex

Determining mechanisms that underlie reproductive isolation (RI) is key to understanding how species boundaries are maintained in nature. Transposable elements (TEs) are ubiquitous across eukaryotic genomes. However, the role of TEs in modulating the strength of RI between species is poorly understood. Several species of Drosophila have been found to harbor P-elements (PEs), yet only D. simulans is known to be currently polymorphic for their presence in wild populations. PEs can cause RI between PE-containing (P) and PE-lacking (M) lineages of the same species. However, it is unclear whether they also contribute to the magnitude of RI between species. Here, we use the simulans species complex to assess whether differences in PE status between D. simulans and its sister species, which do not harbor PEs, contribute to multiple barriers to gene flow between species. We show that crosses involving a P D. simulans father and an M mother from a sister species exhibit lower F1 female fecundity than crosses involving an M D. simulans father and an M sister-species mother. We also find that another TE, I-element, might play a minor role in determining the frequency of dysgenesis between species. Our results suggest that the presence of PEs in a species can strengthen isolation from its sister species, providing evidence that TEs can play a role in RI.

Highly contiguous assemblies of 101 drosophilid genomes

Over 100 years of studies in Drosophila melanogaster and related species in the genus Drosophila have facilitated key discoveries in genetics, genomics, and evolution. While high-quality genome assemblies exist for several species in this group, they only encompass a small fraction of the genus. Recent advances in long-read sequencing allow high-quality genome assemblies for tens or even hundreds of species to be efficiently generated. Here, we utilize Oxford Nanopore sequencing to build an open community resource of genome assemblies for 101 lines of 93 drosophilid species encompassing 14 species groups and 35 sub-groups. The genomes are highly contiguous and complete, with an average contig N50 of 10.5 Mb and greater than 97% BUSCO completeness in 97/101 assemblies. We show that Nanopore-based assemblies are highly accurate in coding regions, particularly with respect to coding insertions and deletions. These assemblies, along with a detailed laboratory protocol and assembly pipelines, are released as a public resource and will serve as a starting point for addressing broad questions of genetics, ecology, and evolution at the scale of hundreds of species.

Mitochondrial genomes of the human pathogens Coccidioides immitis and Coccidioides posadasii

Fungal mitochondrial genomes encode genes involved in crucial cellular processes, such as oxidative phosphorylation and mitochondrial translation, and the molecule has been used as a molecular marker for population genetics studies. Coccidioides immitis and C. posadasii are endemic fungal pathogens that cause coccidioidomycosis in arid regions across both American continents. To date, approximately 150 Coccidioides isolates have been sequenced to infer patterns of variation in nuclear genomes. However, less attention has been given to the mitochondrial genomes of Coccidioides. In this report, we describe the assembly and annotation of mitochondrial reference genomes for two representative strains of C. posadasii and C. immitis, as well as assess population variation among 77 selected genomes. The sizes of the circular-mapping molecules are 68.2 Kb in C. immitis and 75.1 Kb in C. posadasii. We identify 14 mitochondrial protein-coding genes common to most fungal mitochondria, which are largely syntenic across different populations and species of Coccidioides. Both Coccidioides species are characterized by a large number of group I and II introns, harboring twice the number of elements as compared to closely related Onygenales. The introns contain complete or truncated ORFs with high similarity to homing endonucleases of the LAGLIDADG and GIY-YIG families. Phylogenetic comparisons of mitochondrial and nuclear genomes show extensive phylogenetic discordance suggesting that the evolution of the two types of genetic material is not identical. This work represents the first assessment of mitochondrial genomes among isolates of both species of Coccidioides, and provides a foundation for future functional work.

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.

Pure species discriminate against hybrids in the Drosophila melanogaster species subgroup

Introgression, the exchange of alleles between species, is a common event in nature. This transfer of alleles between species must happen through fertile hybrids. Characterizing the traits that cause defects in hybrids illuminates how and when gene flow is expected to occur. Inviability and sterility are extreme examples of fitness reductions but are not the only type of defects in hybrids. Some traits specific to hybrids are more subtle but are important to determine their fitness. In this report, we study whether F1 hybrids between two species pairs of Drosophila are as attractive as the parental species. We find that in both species pairs, the sexual attractiveness of the F1 hybrids is reduced and that pure species discriminate strongly against them. We also find that the cuticular hydrocarbon (CHC) profile of the female hybrids is intermediate between the parental species. Perfuming experiments show that modifying the CHC profile of the female hybrids to resemble pure species improves their chances of mating. Our results show that behavioral discrimination against hybrids might be an important component of the persistence of species that can hybridize.

Comparative studies on speciation: 30 years since Coyne and Orr

Understanding the processes of population divergence and speciation remains a core question in evolutionary biology. For nearly a hundred years evolutionary geneticists have characterized reproductive isolation (RI) mechanisms and specific barriers to gene flow required for species formation. The seminal work of Coyne and Orr provided the first comprehensive comparative analysis of speciation. By combining phylogenetic hypotheses and species range data with estimates of genetic divergence and multiple mechanisms of RI across Drosophila, Coyne and Orr's influential meta-analyses answered fundamental questions and motivated new analyses that continue to push the field forward today. Now 30 years later, we revisit the five questions addressed by Coyne and Orr, identifying results that remain well supported and others that seem less robust with new data. We then consider the future of speciation research, with emphasis on areas where novel methods and data motivate potential progress. While the literature remains biased towards Drosophila and other model systems, we are enthusiastic about the future of the field.

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.

Who are we now? A demographic assessment of three evolution societies

Scientific societies have the potential to catalyze support for communities that have been historically excluded from science. Many of these societies have formed committees to propose and administer initiatives to promote the career and well-being of their members, with a special emphasis on racial and ethnic minorities. Yet, these societies are rarely armed with data to inform their proposals. Three of the evolution societies (American Society of Naturalists, "ASN"; Society of Systematic Biologists, "SSB"; Society for the Study of Evolution, "SSE") have also formed Diversity, Equity, and Inclusion committees in the last few years. As a first step in determining the needs of the societies, these committees collected data on the demographic characteristics of the societies' constituents by surveying the attendants of the Evolution 2019 meeting. Here, we report the proportions for different demographic groups in attendance at the meeting and compare these proportions to the demographics of recipients of Ph.D. degrees either in evolutionary biology or in the broader life sciences, as well as population demographics of the USA. Our results indicate that historically excluded groups are still underrepresented across US-based evolutionary biology professional societies. We explore whether demographic composition differs at different professional stages and find that representation for women and LGBTQ+ members decreases as the career stage progresses. We also find some evidence for heterogeneity across societies in terms of racial composition. Finally, we discuss the caveats and limitations of our procedures. Our results will serve to inform future efforts to collect demographic data at the society levels, which should in turn be used to design and implement evidence-based initiatives for inclusion and equity. This report should be a starting point for systematic efforts to characterize the ever-changing representation in evolutionary biology and to work toward the inclusion of all groups.

Paracoccidioides Genomes Reflect High Levels of Species Divergence and Little Interspecific Gene Flow

The fungus Paracoccidioides is a prevalent human pathogen endemic to South America. The genus is composed of five species. In this report, we use 37 whole-genome sequences to study the allocation of genetic variation in Paracoccidioides We tested three genome-wide predictions of advanced speciation, namely, that all species should be reciprocally monophyletic, that species pairs should be highly differentiated along the whole genome, and that there should be low rates of interspecific gene exchange. We find support for these three hypotheses. Species pairs with older divergences show no evidence of gene exchange, while more recently diverged species pairs show evidence of modest rates of introgression. Our results indicate that as divergence progresses, species boundaries become less porous among Paracoccidioides species. Our results suggest that species in Paracoccidioides are at different stages along the divergence continuum.IMPORTANCEParacoccidioides is the causal agent of a systemic mycosis in Latin America. Most of the inference of the evolutionary history of Paracoccidioides has used only a few molecular markers. In this report, we evaluate the extent of genome divergence among Paracoccidioides species and study the possibility of interspecific gene exchange. We find that all species are highly differentiated. We also find that the amount of gene flow between species is low and in some cases is even completely absent in spite of geographic overlap. Our study constitutes a systematic effort to identify species boundaries in fungal pathogens and to determine the extent of gene exchange among fungal species.

A novel Sporothrix brasiliensis genomic variant in Midwestern Brazil: evidence for an older and wider sporotrichosis epidemic

Sporotrichosis is a subcutaneous infection caused by fungi from the genus Sporothrix. It is transmitted by inoculation of infective particles found in plant-contaminated material or diseased animals, characterizing the classic sapronotic and emerging zoonotic transmission, respectively. Since 1998, southeastern Brazil has experienced a zoonotic sporotrichosis epidemic caused by S. brasiliensis, centred in the state of Rio de Janeiro. Our observation of feline sporotrichosis cases in Brasília (Midwestern Brazil), around 900 km away from Rio de Janeiro, led us to question whether the epidemic caused by S. brasiliensis has spread from the epicentre in Rio de Janeiro, emerged independently in the two locations, or if the disease has been present and unrecognized in Midwestern Brazil. A retrospective analysis of 91 human and 4 animal cases from Brasília, ranging from 1993 to 2018, suggests the occurrence of both sapronotic and zoonotic transmission. Molecular typing of the calmodulin locus identified S. schenckii as the agent in two animals and all seven human patients from which we were able to recover clinical isolates. In two other animals, the disease was caused by S. brasiliensis. Whole-genome sequence typing of seven Sporothrix spp. strains from Brasília and Rio de Janeiro suggests that S. brasiliensis isolates from Brasília are genetically distinct from those obtained at the epicentre of the outbreak in Rio de Janeiro, both in phylogenomic and population genomic analyses. The two S. brasiliensis populations seem to have separated between 2.2 and 3.1 million years ago, indicating independent outbreaks or that the zoonotic S. brasiliensis outbreak might have started earlier and be more widespread in South America than previously recognized.

Environmental and Genetic Contributions to Imperfect wMel-Like Wolbachia Transmission and Frequency Variation

Maternally transmitted Wolbachia bacteria infect about half of all insect species. They usually show imperfect maternal transmission and often produce cytoplasmic incompatibility (CI). Irrespective of CI, Wolbachia frequencies tend to increase when rare only if they benefit host fitness. Several Wolbachia, including wMel that infects Drosophila melanogaster, cause weak or no CI and persist at intermediate frequencies. On the island of São Tomé off West Africa, the frequencies of wMel-like Wolbachia infecting Drosophila yakuba (wYak) and Drosophila santomea (wSan) fluctuate, and the contributions of imperfect maternal transmission, fitness effects, and CI to these fluctuations are unknown. We demonstrate spatial variation in wYak frequency and transmission on São Tomé. Concurrent field estimates of imperfect maternal transmission do not predict spatial variation in wYak frequencies, which are highest at high altitudes where maternal transmission is the most imperfect. Genomic and genetic analyses provide little support for D. yakuba effects on wYak transmission. Instead, rearing at cool temperatures reduces wYak titer and increases imperfect transmission to levels observed on São Tomé. Using mathematical models of Wolbachia frequency dynamics and equilibria, we infer that temporally variable imperfect transmission or spatially variable effects on host fitness and reproduction are required to explain wYak frequencies. In contrast, spatially stable wSan frequencies are plausibly explained by imperfect transmission, modest fitness effects, and weak CI. Our results provide insight into causes of wMel-like frequency variation in divergent hosts. Understanding this variation is crucial to explain Wolbachia spread and to improve wMel biocontrol of human disease in transinfected mosquito systems.

The importance of intrinsic postzygotic barriers throughout the speciation process

Intrinsic postzygotic barriers can play an important and multifaceted role in speciation, but their contribution is often thought to be reserved to the final stages of the speciation process. Here, we review how intrinsic postzygotic barriers can contribute to speciation, and how this role may change through time. We outline three major contributions of intrinsic postzygotic barriers to speciation. (i) reduction of gene flow: intrinsic postzygotic barriers can effectively reduce gene exchange between sympatric species pairs. We discuss the factors that influence how effective incompatibilities are in limiting gene flow. (ii) early onset of species boundaries via rapid evolution: intrinsic postzygotic barriers can evolve between recently diverged populations or incipient species, thereby influencing speciation relatively early in the process. We discuss why the early origination of incompatibilities is expected under some biological models, and detail how other (and often less obvious) incompatibilities may also serve as important barriers early on in speciation. (iii) reinforcement: intrinsic postzygotic barriers can promote the evolution of subsequent reproductive isolation through processes such as reinforcement, even between relatively recently diverged species pairs. We incorporate classic and recent empirical and theoretical work to explore these three facets of intrinsic postzygotic barriers, and provide our thoughts on recent challenges and areas in the field in which progress can be made. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

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.

Paternally Inherited P-Element Copy Number Affects the Magnitude of Hybrid Dysgenesis in Drosophila simulans and D. melanogaster

Transposable elements (TEs) are repetitive regions of DNA that are able to self-replicate and reinsert themselves throughout host genomes. Since the discovery of TEs, a prevalent question has been whether increasing TE copy number has an effect on the fitness of their hosts. P-elements (PEs) in Drosophila are a well-studied TE that has strong phenotypic effects. When a female without PEs (M) is crossed to a male with them (P), the resulting females are often sterile, a phenomenon called hybrid dysgenesis (HD). Here, we used short- and long-read sequencing to infer the number of PEs in the genomes of dozens of isofemale lines from two Drosophila species and measured whether the magnitude of HD was correlated with the number of PEs in the paternal genome. Consistent with previous reports, we find evidence for a positive correlation between the paternal PE copy number and the magnitude of HD in progeny from ♀M × ♂ P crosses for both species. Other crosses are not affected by the number of PE copies. We also find that the correlation between the strength of HD and PE copy number differs between species, which suggests that there are genetic differences that might make some genomes more resilient to the potentially deleterious effects of TEs. Our results suggest that PE copy number interacts with other factors in the genome and the environment to cause HD and that the importance of these interactions is species specific.

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.

Rapid and Predictable Evolution of Admixed Populations Between Two Drosophila Species Pairs

The consequences of hybridization are varied, ranging from the origin of new lineages, introgression of some genes between species, to the extinction of one of the hybridizing species. We generated replicate admixed populations between two pairs of sister species of Drosophila: D. simulans and D. mauritiana; and D. yakuba and D. santomea Each pair consisted of a continental species and an island endemic. The admixed populations were maintained by random mating in discrete generations for over 20 generations. We assessed morphological, behavioral, and fitness-related traits from each replicate population periodically, and sequenced genomic DNA from the populations at generation 20. For both pairs of species, species-specific traits and their genomes regressed to those of the continental species. A few alleles from the island species persisted, but they tended to be proportionally rare among all sites in the genome and were rarely fixed within the populations. This paucity of alleles from the island species was particularly pronounced on the X-chromosome. These results indicate that nearly all foreign genes were quickly eliminated after hybridization and that selection against the minor species genome might be similar across experimental replicates.

Wolbachia Acquisition by Drosophila yakuba-Clade Hosts and Transfer of Incompatibility Loci Between Distantly Related Wolbachia

Maternally transmitted Wolbachia infect about half of insect species, yet the predominant mode(s) of Wolbachia acquisition remains uncertain. Species-specific associations could be old, with Wolbachia and hosts codiversifying (i.e., cladogenic acquisition), or relatively young and acquired by horizontal transfer or introgression. The three Drosophila yakuba-clade hosts [(D. santomea, D. yakuba) D. teissieri] diverged ∼3 MYA and currently hybridize on the West African islands Bioko and São Tomé. Each species is polymorphic for nearly identical Wolbachia that cause weak cytoplasmic incompatibility (CI)-reduced egg hatch when uninfected females mate with infected males. D. yakuba-clade Wolbachia are closely related to wMel, globally polymorphic in D. melanogaster We use draft Wolbachia and mitochondrial genomes to demonstrate that D. yakuba-clade phylogenies for Wolbachia and mitochondria tend to follow host nuclear phylogenies. However, roughly half of D. santomea individuals, sampled both inside and outside of the São Tomé hybrid zone, have introgressed D. yakuba mitochondria. Both mitochondria and Wolbachia possess far more recent common ancestors than the bulk of the host nuclear genomes, precluding cladogenic Wolbachia acquisition. General concordance of Wolbachia and mitochondrial phylogenies suggests that horizontal transmission is rare, but varying relative rates of molecular divergence complicate chronogram-based statistical tests. Loci that cause CI in wMel are disrupted in D. yakuba-clade Wolbachia; but a second set of loci predicted to cause CI are located in the same WO prophage region. These alternative CI loci seem to have been acquired horizontally from distantly related Wolbachia, with transfer mediated by flanking Wolbachia-specific ISWpi1 transposons.

Reinforcement

Reinforcement, or the strengthening of reproductive isolation in response to hybridization costs, is a case in which natural selection directly contributes to the origin of new species. From its conception, reinforcement has enjoyed alternate periods of enthusiasm and rejection, in which powerful verbal arguments identified obstacles that were subsequently solved by novel theory. Empirical data has provided some of the strongest arguments for the plausibility of reinforcement and for its frequency in nature, indicating that it is likely to be a common aspect of species divergence and the completion of speciation.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change

Diverse traits often covary between species [1-3]. The possibility that a single mutation could contribute to the evolution of several characters between species [3] is rarely investigated as relatively few cases are dissected at the nucleotide level. Drosophila santomea has evolved additional sex comb sensory teeth on its legs and has lost two sensory bristles on its genitalia. We present evidence that a single nucleotide substitution in an enhancer of the scute gene contributes to both changes. The mutation alters a binding site for the Hox protein Abdominal-B in the developing genitalia, leading to bristle loss, and for another factor in the developing leg, leading to bristle gain. Our study suggests that morphological evolution between species can occur through a single nucleotide change affecting several sexually dimorphic traits. VIDEO ABSTRACT.

A Maladaptive Combination of Traits Contributes to the Maintenance of a Drosophila Hybrid Zone

Drosophila teissieri and D. yakuba diverged approximately 3 mya and are thought to share a large, ancestral, African range [1-3]. These species now co-occur in parts of continental Africa and in west Africa on the island of Bioko [1, 4]. While D. yakuba is a human commensal, D. teissieri seems to be associated with Parinari fruits, restricting its range to forests [4-6]. Genome data indicate introgression, despite no evidence of contemporary hybridization. Here we report the discovery of D. yakuba-D. teissieri hybrids at the interface of secondary forests and disturbed, open habitats on Bioko. We demonstrate that hybrids are the F₁ progeny of D. yakuba females and D. teissieri males. At high temperatures like those found on Bioko, D. teissieri females are generally less receptive to mating, and in combination with temperature effects on egg lay and egg-to-adult viability, this decreases the potential for gene flow between female D. teissieri and male D. yakuba relative to the reciprocal cross. Field and laboratory experiments demonstrate that F₁ hybrids have a maladaptive combination of D. yakuba behavior and D. teissieri physiology, generating additional barriers to gene flow. Nevertheless, analysis of introgressed and non-introgressed regions of the genome indicate that, while rare, gene flow is relatively recent. Our observations identify precise intrinsic and extrinsic factors that, along with hybrid male sterility, limit gene flow and maintain these species. These data contribute to a growing body of literature that suggests the Gulf of Guinea may be a hotspot for hybridization.

Recent admixture between species of the fungal pathogen Histoplasma

Hybridization between species of pathogens has the potential to speed evolution of virulence by providing the raw material for adaptation through introgression or by assembling new combinations of virulence traits. Fungal diseases are a source high morbidity, and remain difficult to treat. Yet the frequency of hybridization between fungal species has rarely been explored, and the functional role of introgressed alleles remains largely unknown. Histoplasma mississippiense and H. ohiense are sympatric throughout their range in North America and have distinct virulence strategies, making them an ideal system to examine the role introgression may play in fungal pathogens. We identified introgressed tracts in the genomes of a sample of H. mississippiense and H. ohiense isolates. We found strong evidence in each species for recent admixture, but introgressed alleles were present at low frequencies, suggesting that they were deleterious. Consistent with this, coding and regulatory sequences were strongly depleted within introgressed regions, whereas intergenic regions were enriched, indicating that functional introgressed alleles were frequently deleterious in their new genomic context. Surprisingly, we found only two isolates with substantial admixture: the H. mississippiense and H. ohiense genomic reference strains, WU24 and G217B, respectively. Our results show that recent admixture has occurred, that it is frequently deleterious and that conclusions based on studies of the H. mississippiense and H. ohiense type strains should be revisited with more representative samples from the genus.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

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.

Understanding the extent to which species or diverged populations hybridize in nature is crucially important if we are to understand the speciation process. Accordingly numerous research groups have developed methodology for finding the genetic evidence of such introgression. In this report we develop a supervised machine learning approach for uncovering loci which have introgressed across species boundaries. We show that our method, FILET, has greater accuracy and power than competing methods in discovering introgression, and in addition can detect the directionality associated with the gene flow between species. Using whole genome sequences from Drosophila simulans and Drosophila sechellia we show that FILET discovers quite extensive introgression between these species that has occurred mostly from D. simulans to D. sechellia. Our work highlights the complex process of speciation even within a well-studied system and points to the growing importance of supervised machine learning in population genetics.