Hybrid speciation in sparrows II: a role for sex chromosomes?

Homoploid hybrid speciation in animals is poorly understood, mainly because of the scarcity of well-documented cases. Here, we present the results of a multilocus sequence analysis on the house sparrow (Passer domesticus), Spanish sparrow (P. hispaniolensis) and their proposed hybrid descendant, the Italian sparrow (P. italiae). The Italian sparrow is shown to be genetically intermediate between the house sparrow and Spanish sparrow, exhibiting genealogical discordance and a mosaic pattern of alleles derived from either of the putative parental species. The average variation on the Z chromosome was significantly reduced compared with autosomal variation in the putative parental species, the house sparrow and Spanish sparrow. Additionally, divergence between the two species was elevated on the Z chromosome relative to the autosomes. This pattern of variation and divergence is consistent with reduced introgression of Z-linked genes and/or a faster-Z effect (increased rate of adaptive divergence on the Z). F(ST) -outlier tests were consistent with the faster-Z hypothesis: two of five Z-linked loci (CHD1Z and PLAA) were identified as candidates for being subject to positive, divergent selection in the putative parental species. Interestingly, the two latter genes showed a mosaic pattern in the (hybrid) Italian sparrow; that is, the Italian sparrow was found to be fixed for Spanish sparrow alleles at CHD1Z and to mainly have house sparrow alleles at PLAA. Preliminary evidence presented in this study thus suggests that sex chromosomes may play a significant role in this case of homoploid hybrid speciation.

Nuclear and mitochondrial phylogeography of the European fire-bellied toads Bombina bombina and Bombina variegata supports their independent histories

Exact location and number of glacial refugia still remain unclear for many European cold-blooded terrestrial vertebrates. We performed a fine-scaled multilocus phylogeographic analysis of two Bombina species combining mitochondrial variation of 950 toads from 385 sites and nuclear genes (Rag-1, Ncx-1) from a subset of samples to reconstruct their colonization and contemporary variation patterns. We identified the lowlands northwest of the Black Sea and the Carpathians to be important refugial areas for B. bombina and B. variegata, respectively. This result emphasizes the importance of Central European refugia for ectothermic terrestrial species, far north of the Mediterranean areas regarded as exclusive glacial refugia for the animals. Additional refugia for B. variegata have been located in the southern Apennines and Balkans. In contrast, no evidence for the importance of other east European plains as refugial regions has been found. The distribution of mtDNA and Ncx-1 variation suggests the presence of local refugia near the Black Sea for B. bombina; however, coalescent simulations did not allow to distinguish whether one or two refugia were present in the region. Strong genetic drift apparently accompanied postglacial expansions reducing diversity in the colonization areas. Extended sampling, coupled with the multilocus isolation with migration analysis, revealed a limited and geographically restricted gene flow from the Balkan to Carpathian populations of B. variegata. However, despite proximity of inferred B. bombina and B. variegata refugia, gene exchange between them was not detected.

Hybridization is a recurrent evolutionary stimulus in wild yeast speciation

Hybridization can result in reproductively isolated and phenotypically distinct lineages that evolve as independent hybrid species. How frequently hybridization leads to speciation remains largely unknown. Here we examine the potential recurrence of hybrid speciation in the wild yeast Saccharomyces paradoxus in North America, which comprises two endemic lineages SpB and SpC, and an incipient hybrid species, SpC*. Using whole-genome sequences from more than 300 strains, we uncover the hybrid origin of another group, SpD, that emerged from hybridization between SpC* and one of its parental species, the widespread SpB. We show that SpD has the potential to evolve as a novel hybrid species, because it displays phenotypic novelties that include an intermediate transcriptome profile, and partial reproductive isolation with its most abundant sympatric parental species, SpB. Our findings show that repetitive cycles of divergence and hybridization quickly generate diversity and reproductive isolation, providing the raw material for speciation by hybridization.

Hybridization can contribute to diversity from the genomic to the species level. Here, Eberlein, Hénault et al. investigate genomic, transcriptomic and phenotypic variation among wild lineages of the yeast Saccharomyces paradoxus and suggest that an incipient species has formed by recurrent hybridization.

Constraint and adaptation in newt toll-like receptor genes

Acute die-offs of amphibian populations worldwide have been linked to the emergence of viral and fungal diseases. Inter and intraspecific immunogenetic differences may influence the outcome of infection. Toll-like receptors (TLRs) are an essential component of innate immunity and also prime acquired defenses. We report the first comprehensive assessment of TLR gene variation for urodele amphibians. The Lissotriton newt TLR repertoire includes representatives of 13 families and is compositionally most similar to that of the anuran Xenopus. Both ancient and recent gene duplications have occurred in urodeles, bringing the total number of TLR genes to at least 21. Purifying selection has predominated the evolution of newt TLRs in both long (∼70 Ma) and medium (∼18 Ma) timescales. However, we find evidence for both purifying and positive selection acting on TLRs in two recently diverged (2-5 Ma) allopatric evolutionary lineages (Lissotriton montandoni and L. vulgaris graecus). Overall, both forms of selection have been stronger in L. v. graecus, while constraint on most TLR genes in L. montandoni appears relaxed. The differences in selection regimes are unlikely to be biased by demographic effects because these were controlled by means of a historical demographic model derived from an independent data set of 62 loci. We infer that TLR genes undergo distinct trajectories of adaptive evolution in closely related amphibian lineages, highlight the potential of TLRs to capture the signatures of different assemblages of pathogenic microorganisms, and suggest differences between lineages in the relative roles of innate and acquired immunity.

Molecular Inversion Probes for targeted resequencing in non-model organisms

Applications that require resequencing of hundreds or thousands of predefined genomic regions in numerous samples are common in studies of non-model organisms. However few approaches at the scale intermediate between multiplex PCR and sequence capture methods are available. Here we explored the utility of Molecular Inversion Probes (MIPs) for the medium-scale targeted resequencing in a non-model system. Markers targeting 112 bp of exonic sequence were designed from transcriptome of Lissotriton newts. We assessed performance of 248 MIP markers in a sample of 85 individuals. Among the 234 (94.4%) successfully amplified markers 80% had median coverage within one order of magnitude, indicating relatively uniform performance; coverage uniformity across individuals was also high. In the analysis of polymorphism and segregation within family, 77% of 248 tested MIPs were confirmed as single copy Mendelian markers. Genotyping concordance assessed using replicate samples exceeded 99%. MIP markers for targeted resequencing have a number of advantages: high specificity, high multiplexing level, low sample requirement, straightforward laboratory protocol, no need for preparation of genomic libraries and no ascertainment bias. We conclude that MIP markers provide an effective solution for resequencing targets of tens or hundreds of kb in any organism and in a large number of samples.

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata): Inference from transcriptome-wide coalescence analyses

Ancient origins, profound ecological divergence, and extensive hybridization make the fire‐bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small‐scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes.

Balancing selection and introgression of newt immune-response genes

The importance of interspecific introgression as a source of adaptive variation is increasingly recognized. Theory predicts that beneficial genetic variants cross species boundaries easily even when interspecific hybridization is rare and gene flow is strongly constrained throughout the genome. However, it remains unclear whether certain classes of genes are particularly prone to adaptive introgression. Genes affected by balancing selection (BS) may constitute such a class, because forms of BS that favour novel, initially rare alleles, should facilitate introgression. We tested this hypothesis in hybridizing newts by comparing 13 genes with signatures of BS, in particular an excess of common non-synonymous polymorphisms, to the genomic background (154 genes). Parapatric hybridizing taxa were less differentiated in BS candidate genes than more closely related allopatric lineages, while the opposite was observed in the control genes. Coalescent and forward simulations that explored neutral and BS scenarios under isolation and migration showed that processes other than differential gene flow are unlikely to account for this pattern. We conclude that BS, probably involving a form of novel allele advantage, promotes introgression. This mechanism may be a source of adaptively relevant variation in hybridizing species over prolonged periods.

The neutral rate of whole-genome duplication varies among yeast species and their hybrids

Hybridization and polyploidization are powerful mechanisms of speciation. Hybrid speciation often coincides with whole-genome duplication (WGD) in eukaryotes. This suggests that WGD may allow hybrids to thrive by increasing fitness, restoring fertility and/or increasing access to adaptive mutations. Alternatively, it has been suggested that hybridization itself may trigger WGD. Testing these models requires quantifying the rate of WGD in hybrids without the confounding effect of natural selection. Here we show, by measuring the spontaneous rate of WGD of more than 1300 yeast crosses evolved under relaxed selection, that some genotypes or combinations of genotypes are more prone to WGD, including some hybrids between closely related species. We also find that higher WGD rate correlates with higher genomic instability and that WGD increases fertility and genetic variability. These results provide evidence that hybridization itself can promote WGD, which in turn facilitates the evolution of hybrids.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

The distributions of the six species constituting the smooth newt species complex (Lissotriton vulgaris sensu lato and L. montandoni) – an addition to the New Atlas of Amphibians and Reptiles of Europe

The ‘smooth newt’, the taxon traditionally referred to as Lissotriton vulgaris, consists of multiple morphologically distinct taxa. Given the uncertainty concerning the validity and rank of these taxa, L. vulgaris sensu lato has often been treated as a single, polytypic species. A recent study, driven by genetic data, proposed to recognize five species, L. graecus, L. kosswigi, L. lantzi, L. schmidtleri and a more restricted L. vulgaris. The Carpathian newt L. montandoni was confirmed to be a closely related sister species. We propose to refer to this collective of six Lissotriton species as the smooth newt or Lissotriton vulgaris species complex. Guided by comprehensive genomic data from throughout the range of the smooth newt species complex we 1) delineate the distribution ranges, 2) provide a distribution database, and 3) produce distribution maps according to the format of the New Atlas of Amphibians and Reptiles of Europe, for the six constituent species. This allows us to 4) highlight regions where more research is needed to determine the position of contact zones.

Heterogeneous Mutation Rates and Spectra in Yeast Hybrids

Mutation rates and spectra vary between species and among populations. Hybridization can contribute to this variation, but its role remains poorly understood. Estimating mutation rates requires controlled conditions where the effect of natural selection can be minimized. One way to achieve this is through mutation accumulation experiments coupled with genome sequencing. Here, we investigate 400 mutation accumulation lines initiated from 11 genotypes spanning intralineage, interlineage, and interspecific crosses of the yeasts Saccharomyces paradoxus and S. cerevisiae and propagated for 770 generations. We find significant differences in mutation rates and spectra among crosses, which are not related to the level of divergence of parental strains but are specific to some genotype combinations. Differences in number of generations and departures from neutrality play a minor role, whereas polyploidy and loss of heterozygosity impact mutation rates in some of the hybrid crosses in an opposite way.

Cross-feeding affects the target of resistance evolution to an antifungal drug

Pathogenic fungi are a cause of growing concern. Developing an efficient and safe antifungal is challenging because of the similar biological properties of fungal and host cells. Consequently, there is an urgent need to better understand the mechanisms underlying antifungal resistance to prolong the efficacy of current molecules. A major step in this direction would be to be able to predict or even prevent the acquisition of resistance. We leverage the power of experimental evolution to quantify the diversity of paths to resistance to the antifungal 5-fluorocytosine (5-FC), commercially known as flucytosine. We generated hundreds of independent 5-FC resistant mutants derived from two genetic backgrounds from wild isolates of Saccharomyces cerevisiae. Through automated pin-spotting, whole-genome and amplicon sequencing, we identified the most likely causes of resistance for most strains. Approximately a third of all resistant mutants evolved resistance through a pleiotropic drug response, a potentially novel mechanism in response to 5-FC, marked by cross-resistance to fluconazole. These cross-resistant mutants are characterized by a loss of respiration and a strong tradeoff in drug-free media. For the majority of the remaining two thirds, resistance was acquired through loss-of-function mutations in FUR1, which encodes an important enzyme in the metabolism of 5-FC. We describe conditions in which mutations affecting this particular step of the metabolic pathway are favored over known resistance mutations affecting a step upstream, such as the well-known target cytosine deaminase encoded by FCY1. This observation suggests that ecological interactions may dictate the identity of resistance hotspots.

Similarities in biological processes can be used to bridge ecology and molecular biology

Much of the research in biology aims to understand the origin of diversity. Naturally, ecological diversity was the first object of study, but we now have the necessary tools to probe diversity at molecular scales. The inherent differences in how we study diversity at different scales caused the disciplines of biology to be organized around these levels, from molecular biology to ecology. Here, we illustrate that there are key properties of each scale that emerge from the interactions of simpler components and that these properties are often shared across different levels of organization. This means that ideas from one level of organization can be an inspiration for novel hypotheses to study phenomena at another level. We illustrate this concept with examples of events at the molecular level that have analogs at the organismal or ecological level and vice versa. Through these examples, we illustrate that biological processes at different organization levels are governed by general rules. The study of the same phenomena at different scales could enrich our work through a multidisciplinary approach, which should be a staple in the training of future scientists.

Independent Evolution Has Led to Distinct Genomic Signatures in Dutch Elm Disease-Causing Fungi and Other Vascular Wilts-Causing Fungal Pathogens

Vascular wilts are important diseases caused by plant pathogenic fungi that result in the rapid death of their plant hosts. This is due to a systemic defense mechanism whereby the plant induces the compartmentalization of the infected vascular system in order to reduce the propagation of the fungus. The ascomycete class Sordariomycetes contains several species that cause vascular wilts in diverse plant hosts, and they can be classified into four taxonomic orders. The genetic mechanisms of pathogenesis have already been investigated in Fusarium and Verticillium species, but they have not yet been compared with other well-known wilt-causing species, especially fungi causing oak wilt or Dutch elm disease (DED). Here we analyzed 20 whole genome assemblies of wilt-causing fungi together with 56 other species using phylogenetic approaches to trace expansions and contractions of orthologous gene families and gene classes related to pathogenicity. We found that the wilt-causing pathogens evolved seven times, experiencing the largest fold changes in different classes of genes almost every time. However, some similarities exist across groups of wilt pathogens, particularly in Microascales and Ophiostomatales, and these include the common gains and losses of genes that make up secondary metabolite clusters (SMC). DED pathogens do not experience large-scale gene expansions, with most of the gene classes, except for some SMC families, reducing in number. We also found that gene family expansions in the most recent common ancestors of wilt pathogen groups are enriched for carbohydrate metabolic processes. Our study shows that wilt-causing species evolve primarily through distinct changes in their repertoires of pathogenicity-related genes and that there is the potential importance of carbohydrate metabolism genes for regulating osmosis in those pathogens that penetrate the plant vascular system.

Hybrid adaptation is hampered by Haldane's sieve

Hybrids between species exhibit plastic genomic architectures that could foster or slow down their adaptation. When challenged to evolve in an environment containing a UV mimetic drug, yeast hybrids have reduced adaptation rates compared to parents. We find that hybrids and their parents converge onto similar molecular mechanisms of adaptation by mutations in pleiotropic transcription factors, but at a different pace. After 100 generations, mutations in these genes tend to be homozygous in the parents but heterozygous in the hybrids. We hypothesize that a lower rate of loss of heterozygosity (LOH) in hybrids could limit fitness gain. Using genome editing, we first demonstrate that mutations display incomplete dominance, requiring homozygosity to show full impact and to entirely circumvent Haldane's sieve, which favors the fixation of dominant mutations. Second, tracking mutations in earlier generations confirmed a different rate of LOH in hybrids. Together, these findings show that Haldane's sieve slows down adaptation in hybrids, revealing an intrinsic constraint of hybrid genomic architecture that can limit the role of hybridization in adaptive evolution.

Most azole resistance mutations in the Candida albicans drug target confer cross-resistance without intrinsic fitness cost

Azole antifungals are the main drugs used to treat fungal infections. Amino acid substitutions in the drug target Erg11 (Cyp51) are a common resistance mechanism in pathogenic yeasts. How many and which mutations confer resistance is, however, largely unknown. Here we measure the impact of nearly 4,000 amino acid variants of Candida albicans Erg11 on the susceptibility to six clinical azoles. This was achieved by deep mutational scanning of CaErg11 expressed in Saccharomyces cerevisiae. We find that a large fraction of mutations lead to resistance (33%), most resistance mutations confer cross-resistance (88%) and only a handful of resistance mutations show a significant fitness cost (9%). Our results reveal that resistance to azoles can arise through a large set of mutations and this will probably lead to azole pan-resistance, with little evolutionary compromise. This resource will help inform treatment choices in clinical settings and guide the development of new drugs.

Genomic investigations of successful invasions: the picture emerging from recent studies

Invasion biology aims to identify traits and mechanisms that contribute to successful invasions, while also providing general insights into the mechanisms underlying population expansion and adaptation to rapid climate and habitat changes. Certain phenotypic attributes have been linked to successful invasions, and the role of genetics has been critical in understanding adaptation of invasive species. Nevertheless, a comprehensive summary evaluating the most common evolutionary mechanisms associated with successful invasions across species and environments is still lacking. Here we present a systematic review of studies since 2015 that have applied genomic tools to investigate mechanisms of successful invasions across different organisms. We examine demographic patterns such as changes in genomic diversity at the population level, the presence of genetic bottlenecks and gene flow in the invasive range. We review mechanisms of adaptation such as selection from standing genetic variation and de novo mutations, hybridisation and introgression, all of which can have an impact on invasion success. This comprehensive review of recent articles on the genomic diversity of invasive species led to the creation of a searchable database to provide researchers with an accessible resource. Analysis of this database allowed quantitative assessment of demographic and adaptive mechanisms acting in invasive species. A predominant role of admixture in increasing levels of genetic diversity enabling molecular adaptation in novel habitats is the most important finding of our study. The "genetic paradox" of invasive species was not validated in genomic data across species and ecosystems. Even though the presence of genetic drift and bottlenecks is commonly reported upon invasion, a large reduction in genomic diversity is rarely observed. Any decrease in genetic diversity is often relatively mild and almost always restored via gene flow between different invasive populations. The fact that loci under selection are frequently detected suggests that adaptation to novel habitats on a molecular level is not hindered. The above findings are confirmed herein for the first time in a semi-quantitative manner by molecular data. We also point to gaps and potential improvements in the design of studies of mechanisms driving rapid molecular adaptation in invasive populations. These include the scarcity of comprehensive studies that include sampling from multiple native and invasive populations, identification of invasion sources, longitudinal population sampling, and the integration of fitness measures into genomic analyses. We also note that the potential of whole genome studies is often not exploited fully in predicting invasive potential. Comparative genomic studies identifying genome features promoting invasions are underrepresented despite their potential for use as a tool in invasive species control.