Genomic signatures of isolation, hybridization, and selection during speciation of island finches

Sister species occurring sympatrically on islands are rare and offer unique opportunities to understand how speciation can proceed in the face of gene flow. The São Tomé grosbeak is a massive-billed, 'giant' finch endemic to the island of São Tomé in the Gulf of Guinea, where it has diverged from its co-occurring sister species the Príncipe seedeater, an average-sized finch that also inhabits two neighbouring islands. Here, we show that the grosbeak carries a large number of unique alleles different from all three Príncipe seedeater populations, but also shares many alleles with the sympatric São Tomé population of the seedeater, a genomic signature signifying divergence in isolation as well as subsequent introgressive hybridization. Furthermore, genomic segments that remain unique to the grosbeak are situated close to genes, including genes that determine bill morphology, suggesting the preservation of adaptive variation through natural selection during divergence with gene flow. This study reveals a complex speciation process whereby genetic drift, introgression, and selection during periods of isolation and secondary contact all have shaped the diverging genomes of these sympatric island endemic finches.

Counteracting Forces of Introgressive Hybridization and Interspecific Competition Shape the Morphological Traits of Cryptic Lberian Eptesicus Bats

Cryptic species that coexist in sympatry are likely to simultaneously experience strong competition and hybridization. The first phenomenon would lead to character displacement, whereas the second can potentially promote morphological similarity through adaptive introgression. The main goal of this work was to investigate the effect of introgressive hybridization on the morphology of cryptic Iberian Eptesicus bats when facing counteracting evolutionary forces from interspecific competition. We found substantial overlap both in dentition and in wing morphology traits, though mainly in individuals in sympatry. The presence of hybrids contributes to a fifth of this overlap, with hybrids showing traits with intermediate morphometry. Thus, introgressive hybridization may contribute to species adaptation to trophic and ecological space responding directly to the macro-habitats characteristics of the sympatric zone and to local prey availability. On the other hand, fur shade tended to be browner and brighter in hybrids than parental species. Colour differences could result from partitioning of resources as an adaptation to environmental factors such as roost and microhabitats. We argue that a balance between adaptive introgression and niche partitioning shapes species interactions with the environment through affecting morphological traits under selection.

Molecular Structure and Hybridization Patterns of Abramis brama × Rutilus rutilus Hybrids from Modrac Reservoir, Bosnia and Herzegovina

Interspecific hybridization in the Cyprinidae family has been recorded worldwide, with Abramis brama (bream) and Rutilus rutilus (roach) as one of the often-reported hybridizing pairs. The only account of such an event in Bosnia and Herzegovina has been in Modrac Reservoir. Using morphological and molecular markers, the presence of hybrids was surveyed, the hybridization direction was determined and the hybrid group structure in this ecosystem was evaluated. Our findings confirmed unhindered natural hybridization between roach and bream in Modrac Reservoir. Over 50% of the hybrid specimens were classified as F2 hybrids by the NewHybrids software, while the rest were categorized as pure parental form, making it the first such finding in Europe. The analysis of mitochondrial cytochrome b showed that 90% of hybrid individuals were of bream maternal origin. The hybrid group expressed higher mean values of observed heterozygosity and gene diversity than both parental species. Signs of introgressive hybridization between parental species were detected. The hybrid zone of Modrac Reservoir appears to follow the intermediate or “flat” hybrid model based on the balanced distribution of parental and hybrid genotypes. Further investigation is needed to elucidate the factors that enable the survival and mating success of post-F1 individuals.

Distribution of tiger salamanders in northern Sonora, Mexico: comparison of sampling methods and possible implications for an endangered subspecies

Many aquatic species in the arid USA-Mexico borderlands region are imperiled, but limited information on distributions and threats often hinders management. To provide information on the distribution of the Western Tiger Salamander (Ambystoma mavortium), including the USA-federally endangered Sonoran Tiger Salamander (Ambystoma mavortium stebbinsi), we used traditional (seines, dip-nets) and modern (environmental DNA [eDNA]) methods to sample 91 waterbodies in northern Sonora, Mexico, during 2015-2018. The endemic Sonoran Tiger Salamander is threatened by introgressive hybridization and potential replacement by another sub-species of the Western Tiger Salamander, the non-native Barred Tiger Salamander (A. m. mavortium). Based on occupancy models that accounted for imperfect detection, eDNA sampling provided a similar detection probability (0.82 [95% CI: 0.56-0.94]) as seining (0.83 [0.46-0.96]) and much higher detection than dip-netting (0.09 [0.02-0.23]). Volume of water filtered had little effect on detection, possibly because turbid sites had greater densities of salamanders. Salamanders were estimated to occur at 51 sites in 3 river drainages in Sonora. These results indicate tiger salamanders are much more widespread in northern Sonora than previously documented, perhaps aided by changes in land and water management practices. However, because the two subspecies of salamanders cannot be reliably distinguished based on morphology or eDNA methods that are based on mitochondrial DNA, we are uncertain if we detected only native genotypes or if we documented recent invasion of the area by the non-native sub-species. Thus, there is an urgent need for methods to reliably distinguish the subspecies so managers can identify appropriate interventions.

Genetic diversity analysis of the invasive gall pest Leptocybe invasa (Hymenoptera: Apodemidae) from China

Leptocybe invasa Fisher et LaSalle is a global invasive pest that seriously damages Eucalyptus plants. Studying the genetic diversity, genetic structure and introgression hybridization of L. invasa in China is of great significance for clarifying the breeding strategy, future invasion and diffusion trends of L. invasa in China and developing scientific prevention and control measures. Genetic diversity and phylogenetic analyses of 320 L. invasa female adults from 14 geographic populations in China were conducted using 10 polymorphic microsatellite loci (SSRs) and mitochondrial DNA cytochrome oxidase I gene sequences (COIs). (1) The Bayesian phylogenetic tree and haplotype network diagram showed that only haplotype Hap3 existed in L. invasa lineage B in China, while haplotypes Hap1 and Hap2 existed in lineage A, among which haplotype Hap2 was found for the first time. The nucleotide and haplotype diversities of lineage A were higher than those of lineage B. (2) The SSR genetic diversity of the Wuzhou Guangxi, Ganzhou Jiangxi and Panzhihua Sichuan populations was higher than that of the other 11 populations, and the SSR genetic diversity of lineage A was higher than that of lineage B. (3) The AMOVA analysis of mitochondrial COI data showed that 75.55% of the variation was among populations, and 99.86% of the variation was between lineages, while the AMOVA analysis of nuclear SSR data showed that 35.26% of the variation was among populations, and 47.04% of the variation was between lineages. There were obvious differences in the sources of variation between the COI and SSR data. (4) The optimal K value of COI and SSR data in structure analysis was 2, and PCoA analysis also divided the dataset into two obvious categories. The UPMGA phylogenetic tree based on SSR data clustered 14 geographic species into two groups. The results of genetic structure analysis supported the existence of two lineages, A and B, in China. (5) Structural analysis showed that there was obvious introgressive hybridization in Wuzhou Guangxi, Ganzhou Jiangxi, Panzhihua Sichuan and other populations. These results suggest that lineage introgressive hybridization has occurred in the L. invasa population in China. The introgressive hybridization degree and genetic diversity of lineage A are obviously higher than those of lineage B. Lineage introgressive hybridization may be the driving force for further L. invasa invasion and diffusion in China in the future.

Genomics Facilitates Evaluation and Monitoring of McCloud River Redband Trout ( Oncorhynchus mykiss stonei )

The McCloud River Redband Trout (MRRT; Oncorhynchus mykiss stonei ) is a unique subspecies of rainbow trout that inhabits the isolated Upper McCloud River of Northern California. A major threat to MRRT is introgressive hybridization with non-native rainbow trout from historical stocking and contemporary unauthorized introductions . To help address this concern, we collected RAD-sequencing data on 308 total individuals from MRRT and other California O. mykiss populations and examined population structure using Principal Component and admixture analyses. Our results are consistent with previous studies; we found that populations of MRRT in Sheepheaven, Swamp, Edson, and Moosehead creeks are nonintrogressed. Additionally, we saw no evidence of introgression in Dry Creek, and suggest further investigation to determine if it can be considered a core MRRT conservation population. Sheepheaven Creek was previously thought to be the sole historical lineage of MRRT, but our analysis identified three: Sheepheaven, Edson, and Dry creeks, all of which should be preserved. Finally, we discovered diagnostic and polymorphic SNP markers for monitoring introgression and genetic diversity in MRRT. Collectively, our results provide a valuable resource for the conservation and management of MRRT.

The evolutionary pathways for local adaptation in mountain hares

Understanding the evolution of local adaptations is a central aim of evolutionary biology and key for the identification of unique populations and lineages of conservation relevance. By combining RAD sequencing and whole-genome sequencing, we identify genetic signatures of local adaptation in mountain hares (Lepus timidus) from isolated and distinctive habitats of its wide distribution: Ireland, the Alps and Fennoscandia. Demographic modelling suggested that the split of these mountain hares occurred around 20 thousand years ago, providing the opportunity to study adaptive evolution over a short timescale. Using genome-wide scans, we identified signatures of extreme differentiation among hares from distinct geographic areas that overlap with area-specific selective sweeps, suggesting targets for local adaptation. Several identified candidate genes are associated with traits related to the uniqueness of the different environments inhabited by the three groups of mountain hares, including coat colour, ability to live at high altitudes and variation in body size. In Irish mountain hares, a variant of ASIP, a gene previously implicated in introgression-driven winter coat colour variation in mountain and snowshoe hares (L. americanus), may underlie brown winter coats, reinforcing the repeated nature of evolution at ASIP moulding adaptive seasonal colouration. Comparative genomic analyses across several hare species suggested that mountain hares' adaptive variants appear predominantly species-specific. However, using coalescent simulations we also show instances where the candidate adaptive variants have been introduced via introgressive hybridization. Our work shows that standing adaptive variation, including that introgressed from other species, was a crucial component of the post-glacial dynamics of species.

Darwin's finches - an adaptive radiation constructed from ancestral genetic modules

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here we address this issue using Darwin's finches, which vary in size from an 8g warbler finch with a pointed beak to a 40g large ground finch with a massive blunt beak. We present a highly contiguous genome assembly for one of the species and investigate the genomic architecture underlying phenotypic diversity in the entire radiation. Admixture mapping for beak and body size in the small, medium and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins as old as the Darwin's finch phylogeny (1-2 million years). Genes expressed in the developing beak are overrepresented in these genomic regions. Frequencies of allelic variants at the 28 loci covary with phenotypic similarities in body and beak size across the Darwin's finch phylogeny. These ancestral haplotypes constitute genetic modules for selection, and act as key determinants of the exceptional phenotypic diversity of Darwin's finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.