Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier

Accelerated differentiation of neo-W nuclear-encoded mitochondrial genes between two climate-associated bird lineages signals potential co-evolution with mitogenomes

There is considerable evidence for mitochondrial-nuclear co-adaptation as a key evolutionary driver. Hypotheses regarding the roles of sex-linkage have emphasized Z-linked nuclear genes with mitochondrial function (N-mt genes), whereas it remains contentious whether the perfect co-inheritance of W genes with mitogenomes could hinder or facilitate co-adaptation. Young (neo-) sex chromosomes that possess relatively many N-mt genes compared to older chromosomes provide unprecedented hypothesis-testing opportunities. Eastern Yellow Robin (EYR) lineages in coastal and inland habitats with different climates are diverged in mitogenomes, and in a ~ 15.4 Mb nuclear region enriched with N-mt genes, in contrast with otherwise-similar nuclear genomes. This nuclear region maps to passerine chromosome 1A, previously found to be neo-sex in the inland EYR genome. To compare sex-linked Chr1A-derived genes between lineages, we assembled and annotated the coastal EYR genome. We found that: (i) the coastal lineage shares a similar neo-sex system with the inland lineage, (ii) neo-W and neo-Z N-mt genes are not more diverged between lineages than are comparable non-N-mt genes, and showed little evidence for broad positive selection, (iii) however, W-linked N-mt genes are more diverged between lineages than are their Z-linked gametologs. The latter effect was ~7 times stronger for N-mt than non-N-mt genes, suggesting that W-linked N-mt genes might have diverged between lineages under environmental selection through co-evolution with mitogenomes. Finally, we identify a candidate gene driver for divergent selection, NDUFA12. Our data represent a rare example suggesting a possible role for W-associated mitochondrial-nuclear interactions in climate-associated adaptation and lineage differentiation.

Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Multiple karyotype differences between populations of the Hoplias malabaricus (Teleostei; Characiformes), a species complex in the gray area of the speciation process

Neotropical fishes exhibit remarkable karyotype diversity, whose evolution is poorly understood. Here, we studied genetic differences in 60 individuals, from 11 localities of one species, the wolf fish Hoplias malabaricus, from populations that include six different "karyomorphs". These differ in Y-X chromosome differentiation, and, in several cases, by fusions with autosomes that have resulted in multiple sex chromosomes. Other differences are also observed in diploid chromosome numbers and morphologies. In an attempt to start understanding how this diversity was generated, we analyzed within- and between-population differences in a genome-wide sequence data set. We detect clear genotype differences between karyomorphs. Even in sympatry, samples with different karyomorphs differ more in sequence than samples from allopatric populations of the same karyomorph, suggesting that they represent populations that are to some degree reproductively isolated. However, sequence divergence between populations with different karyomorphs is remarkably low, suggesting that chromosome rearrangements may have evolved during a brief evolutionary time. We suggest that the karyotypic differences probably evolved in allopatry, in small populations that would have allowed rapid fixation of rearrangements, and that they became sympatric after their differentiation. Further studies are needed to test whether the karyotype differences contribute to reproductive isolation detected between some H. malabaricus karyomorphs.

Chromosomal conservatism vs chromosomal megaevolution: enigma of karyotypic evolution in Lepidoptera

In the evolution of many organisms, periods of slow genome reorganization (= chromosomal conservatism) are interrupted by bursts of numerous chromosomal changes (= chromosomal megaevolution). Using comparative analysis of chromosome-level genome assemblies, we investigated these processes in blue butterflies (Lycaenidae). We demonstrate that the phase of chromosome number conservatism is characterized by the stability of most autosomes and dynamic evolution of the sex chromosome Z, resulting in multiple variants of NeoZ chromosomes due to autosome-sex chromosome fusions. In contrast during the phase of rapid chromosomal evolution, the explosive increase in chromosome number occurs mainly through simple chromosomal fissions. We show that chromosomal megaevolution is a highly non-random canalized process, and in two phylogenetically independent Lysandra lineages, the drastic parallel increase in number of fragmented chromosomes was achieved, at least partially, through reuse of the same ancestral chromosomal breakpoints. In species showing chromosome number doubling, we found no blocks of duplicated sequences or duplicated chromosomes, thus refuting the hypothesis of polyploidy. In the studied taxa, long blocks of interstitial telomere sequences (ITSs) consist of (TTAGG)_n arrays interspersed with telomere-specific retrotransposons. ITSs are sporadically present in rapidly evolving Lysandra karyotypes, but not in the species with ancestral chromosome number. Therefore, we hypothesize that the transposition of telomeric sequences may be triggers of the rapid chromosome number increase. Finally, we discuss the hypothetical genomic and population mechanisms of chromosomal megaevolution and argue that the disproportionally high evolutionary role of the Z sex chromosome can be additionally reinforced by sex chromosome-autosome fusions and Z-chromosome inversions.

Chromosomal rearrangements and the first indication of an ♀X1 X1 X2 X2 /♂X1 X2 Y sex chromosome system in Rineloricaria fishes (Teleostei: Siluriformes)

Rineloricaria is the most diverse genus within the freshwater fish subfamily Loricariinae, and it is widely distributed in the Neotropical region. Despite limited cytogenetic data, records from southern and south-eastern Brazil suggest a high rate of chromosomal rearrangements in this genus, mirrored in remarkable inter- and intraspecific karyotype variability. In the present work, we investigated the karyotype features of Rineloricaria teffeana, an endemic representative from northern Brazil, using both conventional and molecular cytogenetic techniques. We revealed different diploid chromosome numbers (2n) between sexes (33♂/34♀), which suggests the presence of an ♀X₁ X₁ X₂ X₂ /♂X₁ X₂ Y multiple sex chromosome system. The male-limited Y chromosome was the largest and the only biarmed element in the karyotype, implying Y-autosome fusion as the most probable mechanism behind its origination. C-banding revealed low amounts of constitutive heterochromatin, mostly confined to the (peri)centromeric regions of most chromosomes (including the X₂ and the Y) but also occupying the distal regions of a few chromosomal pairs. The chromosomal localization of the 18S ribosomal DNA (rDNA) clusters revealed a single site on chromosome pair 4, which was adjacent to the 5S rDNA cluster. Additional 5S rDNA loci were present on the autosome pair 8, X₁ chromosome, and in the presumed fusion point on the Y chromosome. The probe for telomeric repeat motif (TTAGGG)_n revealed signals of variable intensities at the ends of all chromosomes except for the Y chromosome, where no detectable signals were evidenced. Male-to-female comparative genomic hybridization revealed no sex-specific or sex-biased repetitive DNA accumulations, suggesting a presumably low level of neo-Y chromosome differentiation. We provide evidence that rDNA sites might have played a role in the formation of this putative multiple sex chromosome system and that chromosome fusions originate through different mechanisms among different Rineloricaria species.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

Effects of the neo-X chromosome on genomic signatures of hybridization in Rumex hastatulus

Natural hybrid zones provide opportunities for studies of the evolution of reproductive isolation in wild populations. Although recent investigations have found that the formation of neo-sex chromosomes is associated with reproductive isolation, the mechanisms remain unclear in most cases. Here, we assess the contemporary structure of gene flow in the contact zone between largely allopatric cytotypes of the dioecious plant Rumex hastatulus, a species with evidence of sex chromosome turn-over. Males to the west of the Mississippi river, USA, have an X and a single Y chromosome, whereas populations to the east of the river have undergone a chromosomal rearrangement giving rise to a larger X and two Y chromosomes. Using reduced-representation sequencing, we provide evidence that hybrids form readily and survive multiple backcross generations in the field, demonstrating the potential for ongoing gene flow between the cytotypes. Cline analysis of each chromosome separately captured no signals of difference in cline shape between chromosomes. However, principal component regression revealed a significant increase in the contribution of individual SNPs to inter-cytotype differentiation on the neo-X chromosome, but no correlation with recombination rate. Cline analysis revealed that the only SNPs with significantly steeper clines than the genome average were located on the neo-X. Our data are consistent with a role for neo-sex chromosomes in reproductive isolation between R. hastatulus cytotypes. Our investigation highlights the importance of studying plant hybrid zones for understanding the evolution of sex chromosomes.