Open Collections

Description

Abstract

Comparisons of genomic variation among closely related species often show more differentiation in mitochondrial DNA (mtDNA) and sex chromosomes than in autosomes, a pattern expected due to the differing effective population sizes and evolutionary dynamics of these genomic components. Yet, introgression can cause species pairs to deviate dramatically from general differentiation trends. The yellowhammer (Emberiza citrinella) and pine bunting (E. leucocephalos) are hybridizing avian sister species that differ greatly in appearance and moderately in nuclear DNA, but that show no mtDNA differentiation. This discordance is best explained by adaptive mtDNA introgression—a process that can select for co-introgression at nuclear genes with mitochondrial functions (mitonuclear genes). To better understand these discordant differentiation patterns and characterize nuclear differentiation in this system, we investigated genome-wide differentiation between allopatric yellowhammers and pine buntings and compared it to what was seen previously in mtDNA. We found significant nuclear differentiation that was highly heterogeneous across the genome, with a particularly wide differentiation peak on the sex chromosome Z. We further investigated mitonuclear gene co-introgression between yellowhammers and pine buntings and found support for this process in the direction of pine buntings into yellowhammers. Genomic signals indicative of co-introgression were common in mitonuclear genes coding for subunits of the mitoribosome and electron transport chain complexes. Such introgression of mitochondrial DNA and mitonuclear genes provides a possible explanation for the patterns of high genomic heterogeneity in genomic differentiation seen among some species groups.

Full methods can be found in the linked 2023 Heredity paper with additional details provided in the attached scripts and text files. Briefly, DNA was extracted from blood and tissue samples taken from avian individuals and then sequenced using a Genotyping-By-Sequencing approach. DNA reads were processed using the scripts provided in the "GBS_Reads_processing_to_VCF.txt" file and then analyzed using the scripts provided in both the "Heredity_Pub_PopulationGenetic_Analysis_Code.R" and "Heredity_Pub_Mitonuclear_Analysis_Code.R"  files. The former of these two R files contains general genetic analyses while the latter contains mitonuclear analyses.