Prevalent Introgression Underlies Convergent Evolution in the Diversification of Pungitius Sticklebacks

Contrast and Genomic Characterisation of Ancient and Recent Interspecific Introgression Between Deeply Diverged Moustache Toads (Leptobrachium)

Recent genomic analyses have provided new insights into the process of interspecific introgression and its consequences on species evolution. Most recent studies, however, focused on hybridization between recently radiated species, with few examining the genomic outcomes of ancient hybridization across deeply diverged species. Using whole genome data of moustache toads (Leptobrachium), we identified signals of three hybridization events among nine species that diverged at the Eocene. An ancient introgression from L. leishanense to the ancestral branch (C1) of L. liui introduced adaptive variants. The highly introgressed regions include genes with important functions in odorant detection and immune responses. These genes are preserved in all three descendent populations of L. liui_C1, and these regions likely have been positively selected over a long filtering process. A recent introgression occurred from L. huashen to L. tengchongense, with the introgressed regions being mostly neutral. Furthermore, one F1 hybrid individual was detected between sympatric L. ailaonicum and L. promustache. The signals of introgression largely disappeared after removing the hybrid individual, indicating an occasional hybridization but minimal introgression. Further examination of highly divergent but low introgressed genomic regions revealed both pre-mating isolation and genetic incompatibility as potential mechanisms of resisting introgression and maintaining species boundaries. Additionally, no large X-effect was found in these introgression events. Hybridization between deeply diverged amphibian species may be common, but detectable introgressions are likely less so, with recent introgression being mostly neutral and the rare ancient one potentially adaptive. Our findings complement recent genomic work, and together they provide a better understanding of the genomic characteristics of interspecific introgression and its significance in species adaptation and evolution.

Improved assembly of the Pungitius pungitius reference genome

The nine-spined stickleback (Pungitius pungitius) has been increasingly used as a model system in studies of local adaptation and sex chromosome evolution but its current reference genome assembly is far from perfect, lacking distinct sex chromosomes. We generated an improved assembly of the nine-spined stickleback reference genome (98.3% BUSCO completeness) with the aid of linked-read mapping. While the new assembly (v8) was of similar size as the earlier version (v7), we were able to assign 4.4 times more contigs to the linkage groups and improve the contiguity of the genome. Moreover, the new assembly contains a ∼22.8 Mb Y-linked scaffold (LG22) consisting mainly of previously assigned X-contigs, putative Y-contigs, putative centromere contigs, and highly repetitive elements. The male individual showed an even mapping depth on LG12 (pseudo X chromosome) and LG22 (Y-linked scaffold) in the segregating sites, suggesting near-pure X and Y representation in the v8 assembly. A total of 26,803 genes were annotated, and about 33% of the assembly was found to consist of repetitive elements. The high proportion of repetitive elements in LG22 (53.10%) suggests it can be difficult to assemble the complete sequence of the species' Y chromosome. Nevertheless, the new assembly is a significant improvement over the previous version and should provide a valuable resource for genomic studies of stickleback fishes.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.