Sex chromosome evolution: molecular aspects of Y-chromosome degeneration in Drosophila

Rewinding the Ratchet: Rare Recombination Locally Rescues Neo-W Degeneration and Generates Plateaus of Sex-Chromosome Divergence

Natural selection is less efficient in the absence of recombination. As a result, nonrecombining sequences, such as sex chromosomes, tend to degenerate over time. Although the outcomes of recombination arrest are typically observed after many millions of generations, recent neo-sex chromosomes can give insight into the early stages of this process. Here, we investigate the evolution of neo-sex chromosomes in the Spanish marbled white butterfly, Melanargia ines, where a Z-autosome fusion has turned the homologous autosome into a nonrecombining neo-W chromosome. We show that these neo-sex chromosomes are likely limited to the Iberian population of M. ines, and that they arose around the time when this population split from North-African populations, around 1.5 million years ago. Recombination arrest of the neo-W chromosome has led to an excess of premature stop-codons and frame-shift mutations, and reduced gene expression compared to the neo-Z chromosome. Surprisingly, we identified two regions of ∼1 Mb at one end of the neo-W that are both less diverged from the neo-Z and less degraded than the rest of the chromosome, suggesting a history of rare but repeated genetic exchange between the two neo-sex chromosomes. These plateaus of neo-sex chromosome divergence suggest that neo-W degradation can be locally reversed by rare recombination between neo-W and neo-Z chromosomes.

Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome

Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

Accumulation of retrotransposons contributes to W chromosome differentiation in the willow beauty Peribatodes rhomboidaria (Lepidoptera: Geometridae)

The W chromosome of Lepidoptera is typically gene-poor, repeat-rich and composed of heterochromatin. Pioneering studies investigating this chromosome reported an abundance of mobile elements. However, the actual composition of the W chromosome varies greatly between species, as repeatedly demonstrated by comparative genomic hybridization (CGH) or fluorescence in situ hybridization (FISH). Here we present an analysis of repeats on the W chromosome in the willow beauty, Peribatodes rhomboidaria (Geometridae), a species in which CGH predicted an abundance of W-enriched or W-specific sequences. Indeed, comparative analysis of male and female genomes using RepeatExplorer identified ten putative W chromosome-enriched repeats, most of which are LTR or LINE mobile elements. We analysed the two most abundant: PRW LINE-like and PRW Bel-Pao. The results of FISH mapping and bioinformatic analysis confirmed their enrichment on the W chromosome, supporting the hypothesis that mobile elements are the driving force of W chromosome differentiation in Lepidoptera. As the W chromosome is highly underrepresented in chromosome-level genome assemblies of Lepidoptera, this recently introduced approach, combining bioinformatic comparative genome analysis with molecular cytogenetics, provides an elegant tool for studying this elusive and rapidly evolving part of the genome.

Origin and chromatin remodeling of young X/Y sex chromosomes in catfish with sexual plasticity

Assembly of a complete Y chromosome is a significant challenge in animals with an XX/XY sex-determination system. Recently, we created YY-supermale yellow catfish by crossing XY males with sex-reversed XY females, providing a valuable model for Y-chromosome assembly and evolution. Here, we assembled highly homomorphic Y and X chromosomes by sequencing genomes of the YY supermale and XX female in yellow catfish, revealing their nucleotide divergences with only less than 1% and with the same gene compositions. The sex-determining region (SDR) was identified to locate within a physical distance of 0.3 Mb by F_ST scanning. Strikingly, the incipient sex chromosomes were revealed to originate via autosome-autosome fusion and were characterized by a highly rearranged region with an SDR downstream of the fusion site. We found that the Y chromosome was at a very early stage of differentiation, as no clear evidence of evolutionary strata and classical structure features of recombination suppression for a rather late stage of Y-chromosome evolution were observed. Significantly, a number of sex-antagonistic mutations and the accumulation of repetitive elements were discovered in the SDR, which might be the main driver of the initial establishment of recombination suppression between young X and Y chromosomes. Moreover, distinct three-dimensional chromatin organizations of the Y and X chromosomes were identified in the YY supermales and XX females, as the X chromosome exhibited denser chromatin structure than the Y chromosome, while they respectively have significantly spatial interactions with female- and male-related genes compared with other autosomes. The chromatin configuration of the sex chromosomes as well as the nucleus spatial organization of the XX neomale were remodeled after sex reversal and similar to those in YY supermales, and a male-specific loop containing the SDR was found in the open chromatin region. Our results elucidate the origin of young sex chromosomes and the chromatin remodeling configuration in the catfish sexual plasticity.

The Rattlesnake W Chromosome: A GC-Rich Retroelement Refugium with Retained Gene Function Across Ancient Evolutionary Strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a "refugium" for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

Tempo of degeneration across independently evolved non-recombining regions

Recombination is beneficial over the long term, allowing more effective selection. Despite long-term advantages of recombination, local recombination suppression can evolve and lead to genomic degeneration, in particular on sex chromosomes. Here, we investigated the tempo of degeneration in nonrecombining regions, that is, the function curve for the accumulation of deleterious mutations over time, leveraging on 22 independent events of recombination suppression identified on mating-type chromosomes of anther-smut fungi, including newly identified ones. Using previously available and newly generated high-quality genome assemblies of alternative mating types of 13 Microbotryum species, we estimated degeneration levels in terms of accumulation of nonoptimal codons and nonsynonymous substitutions in nonrecombining regions. We found a reduced frequency of optimal codons in the nonrecombining regions compared with autosomes, that was not due to less frequent GC-biased gene conversion or lower ancestral expression levels compared with recombining regions. The frequency of optimal codons rapidly decreased following recombination suppression and reached an asymptote after ca. 3 Ma. The strength of purifying selection remained virtually constant at d_N/d_S = 0.55, that is, at an intermediate level between purifying selection and neutral evolution. Accordingly, nonsynonymous differences between mating-type chromosomes increased linearly with stratum age, at a rate of 0.015 per My. We thus develop a method for disentangling effects of reduced selection efficacy from GC-biased gene conversion in the evolution of codon usage and we quantify the tempo of degeneration in nonrecombining regions, which is important for our knowledge on genomic evolution and on the maintenance of regions without recombination.

Genomic deciphering of sex determination and unique immune system of a potential model species rare minnow (Gobiocypris rarus)

Gobiocypris rarus is sensitive to environmental pollution, especially to heavy metal and grass carp reovirus (GCRV). Hence, it has potential utility as a biological monitor. Genetic deciphering of its unique immune system will advance our understanding of its unique adaptive strategies, which provide cues for its better application. A de novo genome of rare minnow was obtained, and its sex determination mechanism is ZZ/ZW. We identified several specific mutation genes and specific lost genes of rare minnow, and these might be related to the sensitivity of rare minnow to environmental stimuli. We also analyzed the gene expression level of different organs/tissues and found that several IFIT genes may play key roles in GCRV resistance. In addition, knockout of the gene PCDH10L indicates that PCDH10L affects Pb²⁺-induced mortality in rare minnow. Rare minnow is ready for genetic manipulation and shows potential as an emerging experimental model.

Size Variation of the Nonrecombining Region on the Mating-Type Chromosomes in the Fungal Podospora anserina Species Complex

Sex chromosomes often carry large nonrecombining regions that can extend progressively over time, generating evolutionary strata of sequence divergence. However, some sex chromosomes display an incomplete suppression of recombination. Large genomic regions without recombination and evolutionary strata have also been documented around fungal mating-type loci, but have been studied in only a few fungal systems. In the model fungus Podospora anserina (Ascomycota, Sordariomycetes), the reference S strain lacks recombination across a 0.8-Mb region around the mating-type locus. The lack of recombination in this region ensures that nuclei of opposite mating types are packaged into a single ascospore (pseudohomothallic lifecycle). We found evidence for a lack of recombination around the mating-type locus in the genomes of ten P. anserina strains and six closely related pseudohomothallic Podospora species. Importantly, the size of the nonrecombining region differed between strains and species, as indicated by the heterozygosity levels around the mating-type locus and experimental selfing. The nonrecombining region is probably labile and polymorphic, differing in size and precise location within and between species, resulting in occasional, but infrequent, recombination at a given base pair. This view is also supported by the low divergence between mating types, and the lack of strong linkage disequilibrium, chromosomal rearrangements, transspecific polymorphism and genomic degeneration. We found a pattern suggestive of evolutionary strata in P. pseudocomata. The observed heterozygosity levels indicate low but nonnull outcrossing rates in nature in these pseudohomothallic fungi. This study adds to our understanding of mating-type chromosome evolution and its relationship to mating systems.

Epigenetics drive the evolution of sex chromosomes in animals and plants

We review how epigenetics affect sex chromosome evolution in animals and plants. In a few species, sex is determined epigenetically through the action of Y-encoded small RNAs. Epigenetics is also responsible for changing the sex of individuals through time, even in species that carry sex chromosomes, and could favour species adaptation through breeding system plasticity. The Y chromosome accumulates repeats that become epigenetically silenced which leads to an epigenetic conflict with the expression of Y genes and could accelerate Y degeneration. Y heterochromatin can be lost through ageing, which activates transposable elements and lowers male longevity. Y chromosome degeneration has led to the evolution of meiotic sex chromosome inactivation in eutherians (placentals) and marsupials, and dosage compensation mechanisms in animals and plants. X-inactivation convergently evolved in eutherians and marsupials via two independently evolved non-coding RNAs. In Drosophila, male X upregulation by the male specific lethal (MSL) complex can spread to neo-X chromosomes through the transposition of transposable elements that carry an MSL-binding motif. We discuss similarities and possible differences between plants and animals and suggest future directions for this dynamic field of research. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

Identification of the First Oomycete Mating-type Locus Sequence in the Grapevine Downy Mildew Pathogen, Plasmopara viticola

Mating types are self-incompatibility systems that promote outcrossing in plants, fungi, and oomycetes. Mating-type genes have been widely studied in plants and fungi but have yet to be identified in oomycetes, eukaryotic organisms closely related to brown algae that cause many destructive animal and plant diseases. We identified the mating-type locus of Plasmopara viticola, the oomycete responsible for grapevine downy mildew, one of the most damaging grapevine diseases worldwide. Using a genome-wide association approach, we identified a 570-kb repeat-rich non-recombining region controlling mating types, with two highly divergent alleles. We showed that one mating type was homozygous, whereas the other was heterozygous at this locus. The mating-type locus encompassed 40 genes, including one encoding a putative hormone receptor. Functional studies will, however, be required to validate the function of these genes and find the actual determinants of mating type. Our findings have fundamental implications for our understanding of the evolution of mating types, as they reveal a unique determinism involving an asymmetry of heterozygosity, as in sex chromosomes and unlike other mating-type systems. This identification of the mating-type locus in such an economically important crop pathogen also has applied implications, as outcrossing facilitates rapid evolution and resistance to harsh environmental conditions.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies—Leptidea juvernica, L. sinapis and L. reali—have each a unique set of multiple sex-chromosomes with 3–4 W and 3–4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

Differential Gene Expression between Fungal Mating Types Is Associated with Sequence Degeneration

Degenerative mutations in non-recombining regions, such as in sex chromosomes, may lead to differential expression between alleles if mutations occur stochastically in one or the other allele. Reduced allelic expression due to degeneration has indeed been suggested to occur in various sex-chromosome systems. However, whether an association occurs between specific signatures of degeneration and differential expression between alleles has not been extensively tested, and sexual antagonism can also cause differential expression on sex chromosomes. The anther-smut fungus Microbotryum lychnidis-dioicae is ideal for testing associations between specific degenerative signatures and differential expression because 1) there are multiple evolutionary strata on the mating-type chromosomes, reflecting successive recombination suppression linked to mating-type loci; 2) separate haploid cultures of opposite mating types help identify differential expression between alleles; and 3) there is no sexual antagonism as a confounding factor accounting for differential expression. We found that differentially expressed genes were enriched in the four oldest evolutionary strata compared with other genomic compartments, and that, within compartments, several signatures of sequence degeneration were greater for differentially expressed than non-differentially expressed genes. Two particular degenerative signatures were significantly associated with lower expression levels within differentially expressed allele pairs: upstream insertion of transposable elements and mutations truncating the protein length. Other degenerative mutations associated with differential expression included nonsynonymous substitutions and altered intron or GC content. The association between differential expression and allele degeneration is relevant for a broad range of taxa where mating compatibility or sex is determined by genes located in large regions where recombination is suppressed.

Distribution of CR1-like transposable element in woodpeckers (Aves Piciformes): Z sex chromosomes can act as a refuge for transposable elements

Birds have relatively few repetitive sequences compared to other groups of vertebrates; however, the members of order Piciformes (woodpeckers) have more of these sequences, composed mainly of transposable elements (TE). The TE most often found in birds is a retrotransposon chicken repeat 1 (CR1). Piciformes lineages were subjected to an expansion of the CR1 elements, carrying a larger fraction of transposable elements. This study compared patterns of chromosome distribution among five bird species, through chromosome mapping of the CR1 sequence and reconstructed their phylogenetic tree. We analyzed several members of Piciformes (Colaptes campestris, Colaptes melanochloros, Melanerpes candidus, and Veniliornis spilogaster), as well as Galliformes (Gallus gallus). Gallus gallus is the species with which most genomic and hence cytogenetic studies have been performed. The results showed that CR1 sequences are a monophyletic group and do not depend on their hosts. All species analyzed showed a hybridization signal by fluorescence in situ hybridization (FISH). In all species, the chromosomal distribution of CR1 was not restricted to heterochromatin regions in the macrochromosomes, principally pair 1 and the Z sex chromosome. Accumulation in the Z sex chromosomes can serve as a refuge for transposable elements. These results highlight the importance of transposable elements in host genomes and karyotype evolution.

Non-canonical Drosophila X chromosome dosage compensation and repressive topologically associated domains

Background

In animals with XY sex chromosomes, X-linked genes from a single X chromosome in males are imbalanced relative to autosomal genes. To minimize the impact of genic imbalance in male Drosophila, there is a dosage compensation complex (MSL) that equilibrates X-linked gene expression with the autosomes. There are other potential contributions to dosage compensation. Hemizygous autosomal genes located in repressive chromatin domains are often derepressed. If this homolog-dependent repression occurs on the X, which has no pairing partner, then derepression could contribute to male dosage compensation.

Results

We asked whether different chromatin states or topological associations correlate with X chromosome dosage compensation, especially in regions with little MSL occupancy. Our analyses demonstrated that male X chromosome genes that are located in repressive chromatin states are depleted of MSL occupancy; however, they show dosage compensation. The genes in these repressive regions were also less sensitive to knockdown of MSL components.

Conclusions

Our results suggest that this non-canonical dosage compensation is due to the same transacting derepression that occurs on autosomes. This mechanism would facilitate immediate compensation during the evolution of sex chromosomes from autosomes. This mechanism is similar to that of C. elegans, where enhanced recruitment of X chromosomes to the nuclear lamina dampens X chromosome expression as part of the dosage compensation response in XX individuals.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0232-y) contains supplementary material, which is available to authorized users.

Regulation of the X Chromosome in the Germline and Soma of Drosophila melanogaster Males

During the evolution of heteromorphic sex chromosomes, the sex-specific Y chromosome degenerates, while the X chromosome evolves new mechanisms of regulation. Using bioinformatic and experimental approaches, we investigate the expression of the X chromosome in Drosophila melanogaster. We observe nearly complete X chromosome dosage compensation in male somatic tissues, but not in testis. The X chromosome contains disproportionately fewer genes with high expression in testis than the autosomes, even after accounting for the lack of dosage compensation, which suggests that another mechanism suppresses their expression in the male germline. This is consistent with studies of reporter genes and transposed genes, which find that the same gene has higher expression when autosomal than when X-linked. Using a new reporter gene that is expressed in both testis and somatic tissues, we find that the suppression of X-linked gene expression is limited to genes with high expression in testis and that the extent of the suppression is positively correlated with expression level.

Evolution and conservation of Characidium sex chromosomes

Fish species exhibit substantial variation in the degree of genetic differentiation between sex chromosome pairs, and therefore offer the opportunity to study the full range of sex chromosome evolution. We used restriction-site associated DNA sequencing (RAD-seq) to study the sex chromosomes of Characidium gomesi, a species with conspicuous heteromorphic ZW/ZZ sex chromosomes. We screened 9863 single-nucleotide polymorphisms (SNPs), corresponding to ~1 marker/100 kb distributed across the genome for sex-linked variation. With this data set, we identified 26 female-specific RAD loci, putatively located on the W chromosome, as well as 148 sex-associated SNPs showing significant differentiation (average F_ST=0.144) between males and females, and therefore in regions of more recent divergence between the Z and W chromosomes. In addition, we detected 25 RAD loci showing extreme heterozygote deficiency in females but which were in Hardy-Weinberg equilibrium in males, consistent with degeneration of the W chromosome and therefore female hemizygosity. We validated seven female-specific and two sex-associated markers in a larger sample of C. gomesi, of which three localised to the W chromosome, thereby providing useful markers for sexing wild samples. Validated markers were evaluated in other populations and species of the genus Characidium, this exploration suggesting a rapid turnover of W-specific repetitive elements. Together, our analyses point to a complex origin for the sex chromosome of C. gomesi and highlight the utility of RAD-seq for studying the composition and evolution of sex chromosomes systems in wild populations.

Achiasmy: Male Fruit Flies Are Not Ready to Mix

Maintenance of the chromosomal copy number over generations and recombination between homologous chromosomes are hallmarks of meiotic cell division. This genetic exchange that take place during gamete formation leads to genetic diversity, the main driving force behind natural selection. Formation of chiasmata, the physical link between homologous chromosomes during meiosis, is a requisite for recombination. In addition, chiasmata also aid in proper segregation of homologous chromosomes and has a major impact on reproductive fitness. Given these facts it is intriguing that many insect species have forgone the need for genetic exchange between homologous chromosomes during meiosis. Geneticists for several decades knew that meiotic crossover and recombination is absent in Drosophila males and some female lepidopterans, a condition termed achiasmy. However, a good understanding of the mechanisms that cause achiasmy and the evolutionary benefits of achiasmy is currently lacking. In this article we will discuss possible genetic and molecular basis of achiasmy in male Drosophila.

Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes

The present review discusses the roles of repetitive sequences played in plant sex chromosome evolution, and highlights epigenetic modification as potential mechanism of repetitive sequences involved in sex chromosome evolution. Sex determination in plants is mostly based on sex chromosomes. Classic theory proposes that sex chromosomes evolve from a specific pair of autosomes with emergence of a sex-determining gene(s). Subsequently, the newly formed sex chromosomes stop recombination in a small region around the sex-determining locus, and over time, the non-recombining region expands to almost all parts of the sex chromosomes. Accumulation of repetitive sequences, mostly transposable elements and tandem repeats, is a conspicuous feature of the non-recombining region of the Y chromosome, even in primitive one. Repetitive sequences may play multiple roles in sex chromosome evolution, such as triggering heterochromatization and causing recombination suppression, leading to structural and morphological differentiation of sex chromosomes, and promoting Y chromosome degeneration and X chromosome dosage compensation. In this article, we review the current status of this field, and based on preliminary evidence, we posit that repetitive sequences are involved in sex chromosome evolution probably via epigenetic modification, such as DNA and histone methylation, with small interfering RNAs as the mediator.

Mapping of the Retrotransposable Elements Rex1 and Rex3 in Chromosomes of Eigenmannia (Teleostei, Gymnotiformes, Sternopygidae)

Transposable elements constitute a remarkable fraction of the eukaryote genome and show particular capacity to move and insert in specific regions of the genome. This study identified the retrotransposable elements Rex1 and Rex3 in the genomes of 6 cytotypes of Eigenmannia. The sequences were isolated by PCR, sequenced and physically mapped in the chromosomes of these cytotypes, aiming to investigate the organization and distribution of these elements in this fish group, mainly in the sex chromosomes. The FISH physical mapping revealed that both Rex1 and Rex3 elements are dispersed in small clusters throughout the chromosomes of all cytotypes analyzed. However, conspicuous blocks occur in several samples, including an accentuated accumulation of the Rex3 element in X1 and X2 chromosomes of Eigenmannia sp. 2 and in the X chromosome of E. virescens. The accumulations are coincident with heterochromatin-rich regions, suggesting that Rex3 played a role in the differentiation process of the sex chromosomes.

Degeneration of the nonrecombining regions in the mating-type chromosomes of the anther-smut fungi

Dimorphic mating-type chromosomes in fungi are excellent models for understanding the genomic consequences of recombination suppression. Their suppressed recombination and reduced effective population size are expected to limit the efficacy of natural selection, leading to genomic degeneration. Our aim was to identify the sequences of the mating-type chromosomes (a1 and a2) of the anther-smut fungi and to investigate degeneration in their nonrecombining regions. We used the haploid a1 Microbotryum lychnidis-dioicae reference genome sequence. The a1 and a2 mating-type chromosomes were both isolated electrophoretically and sequenced. Integration with restriction-digest optical maps identified regions of recombination and nonrecombination in the mating-type chromosomes. Genome sequence data were also obtained for 12 other Microbotryum species. We found strong evidence of degeneration across the genus in the nonrecombining regions of the mating-type chromosomes, with significantly higher rates of nonsynonymous substitution (dN/dS) than in nonmating-type chromosomes or in recombining regions of the mating-type chromosomes. The nonrecombining regions of the mating-type chromosomes also showed high transposable element content, weak gene expression, and gene losses. The levels of degeneration did not differ between the a1 and a2 mating-type chromosomes, consistent with the lack of homogametic/heterogametic asymmetry between them, and contrasting with X/Y or Z/W sex chromosomes.

Positive and purifying selection on the Drosophila Y chromosome

Y chromosomes, with their reduced effective population size, lack of recombination, and male-limited transmission, present a unique collection of constraints for the operation of natural selection. Male-limited transmission may greatly increase the efficacy of selection for male-beneficial mutations, but the reduced effective size also inflates the role of random genetic drift. Together, these defining features of the Y chromosome are expected to influence rates and patterns of molecular evolution on the Y as compared with X-linked or autosomal loci. Here, we use sequence data from 11 genes in 9 Drosophila species to gain insight into the efficacy of natural selection on the Drosophila Y relative to the rest of the genome. Drosophila is an ideal system for assessing the consequences of Y-linkage for molecular evolution in part because the gene content of Drosophila Y chromosomes is highly dynamic, with orthologous genes being Y-linked in some species whereas autosomal in others. Our results confirm the expectation that the efficacy of natural selection at weakly selected sites is reduced on the Y chromosome. In contrast, purifying selection on the Y chromosome for strongly deleterious mutations does not appear to be compromised. Finally, we find evidence of recurrent positive selection for 4 of the 11 genes studied here. Our results thus highlight the variable nature of the mode and impact of natural selection on the Drosophila Y chromosome.

The chromatin landscape of Drosophila: comparisons between species, sexes, and chromosomes

The chromatin landscape is key for gene regulation, but little is known about how it differs between sexes or between species. Here, we study the sex-specific chromatin landscape of Drosophila miranda, a species with young sex chromosomes, and compare it with Drosophila melanogaster. We analyze six histone modifications in male and female larvae of D. miranda (H3K4me1, H3K4me3, H3K36me3, H4K16ac, H3K27me3, and H3K9me2), and define seven biologically meaningful chromatin states that show different enrichments for transcribed and silent genes, repetitive elements, housekeeping, and tissue-specific genes. The genome-wide distribution of both active and repressive chromatin states differs between males and females. In males, active chromatin is enriched on the X, relative to females, due to dosage compensation of the hemizygous X. Furthermore, a smaller fraction of the euchromatic portion of the genome is in a repressive chromatin state in males relative to females. However, sex-specific chromatin states appear not to explain sex-biased expression of genes. Overall, conservation of chromatin states between male and female D. miranda is comparable to conservation between D. miranda and D. melanogaster, which diverged >30 MY ago. Active chromatin states are more highly conserved across species, while heterochromatin shows very low levels of conservation. Divergence in chromatin profiles contributes to expression divergence between species, with ∼26% of genes in different chromatin states in the two species showing species-specific or species-biased expression, an enrichment of approximately threefold over null expectation. Our data suggest that heteromorphic sex chromosomes in males (that is, a hypertranscribed X and an inactivated Y) may contribute to global redistribution of active and repressive chromatin marks between chromosomes and sexes.

How do y-chromosomes modulate genome-wide epigenetic States: genome folding, chromatin sinks, and gene expression

The Y chromosomes of Drosophila melanogaster and D. simulans contain only a handful of protein-coding genes, which are related to sperm mobility and reproductive fitness. Despite low or absent protein coding polymorphism, the Drosophila Y chromosome has been associated with natural phenotypic variation, including variation in the expression of hundreds to thousands of genes located on autosomes and on the X chromosome. Polymorphisms present in the large blocks of heterochromatin and consisting of differences in the amounts and kinds of sequences for satellite DNA and transposable elements may be the source of this modulation. Here we review the evidence and discuss mechanisms for global epigenetic regulation by repetitious elements in the Y chromosome. We also discuss how the discovery of this new function impacts the current knowledge about Y chromosome origin, its current dynamics, and future fate.

Rapid degeneration of noncoding DNA regions surrounding SlAP3X/Y after recombination suppression in the dioecious plant Silene latifolia

Silene latifolia is a dioecious plant with heteromorphic XY sex chromosomes. Previous studies of sex chromosome–linked genes have suggested a gradual divergence between the X-linked and the Y-linked genes in proportion to the distance from the pseudoautosomal region. However, such a comparison has yet to be made for the noncoding regions. To better characterize the nonrecombining region of the X and Y chromosomes, we sequenced bacterial artificial chromosome clones containing the sex chromosome–linked paralogs SlAP3X and SlAP3Y, including 115 kb and 73 kb of sequences, respectively, flanking these genes. The synonymous nucleotide divergence between SlAP3X and SlAP3Y indicated that recombination stopped approximately 3.4 million years ago. Sequence homology analysis revealed the presence of six long terminal repeat retrotransposon-like elements. Using the nucleotide divergence calculated between left and right long terminal repeat sequences, insertion dates were estimated to be 0.083–1.6 million years ago, implying that all elements detected were inserted after recombination stopped. A reciprocal sequence homology search facilitated the identification of four homologous noncoding DNA regions between the X and Y chromosomes, spanning 6.7% and 10.6% of the X chromosome–derived and Y chromosome–derived sequences, respectively, investigated. Genomic Southern blotting and fluorescence in situ hybridization showed that the noncoding DNA flanking SlAP3X/Y has homology to many regions throughout the genome, regardless of whether they were homologous between the X and Y chromosomes. This finding suggests that most noncoding DNA regions rapidly lose their counterparts because of the introduction of transposable elements and indels (insertion–deletions) after recombination has stopped.

Chromosomal organization and evolutionary history of Mariner transposable elements in Scarabaeinae coleopterans

Background

With the aim to increase the knowledge on the evolution of coleopteran genomes, we investigated through cytogenetics and nucleotide sequence analysis Mariner transposons in three Scarabaeinae species (Coprophanaeus cyanescens, C. ensifer and Diabroctis mimas).

Results

The cytogenetic mapping revealed an accumulation of Mariner transposon in the pericentromeric repetitive regions characterized as rich in heterochromatin and C₀t-1 DNA fraction (DNA enriched with high and moderately repeated sequences). Nucleotide sequence analysis of Mariner revealed the presence of two major groups of Mariner copies in the three investigated coleoptera species.

Conclusions

The Mariner is accumulated in the centromeric area of the coleopteran chromosomes probably as a consequence of the absence of recombination in the heterochromatic regions. Our analysis detected high diversification of Mariner sequences during the evolutionary history of the group. Furthermore, comparisons between the coleopterans sequences with other insects and mammals, suggest that the horizontal transfer (HT) could have acted in the spreading of the Mariner in diverse non-related animal groups.

Meta-analysis reveals that genes regulated by the Y chromosome in Drosophila melanogaster are preferentially localized to repressive chromatin

The Drosophila Y chromosome is a degenerated, heterochromatic chromosome with few functional genes. Despite this, natural variation on the Y chromosome in D. melanogaster has substantial trans-acting effects on the regulation of X-linked and autosomal genes. It is not clear, however, whether these genes simply represent a random subset of the genome or whether specific functional properties are associated with susceptibility to regulation by Y-linked variation. Here, we present a meta-analysis of four previously published microarray studies of Y-linked regulatory variation (YRV) in D. melanogaster. We show that YRV genes are far from a random subset of the genome: They are more likely to be in repressive chromatin contexts, be expressed tissue specifically, and vary in expression within and between species than non-YRV genes. Furthermore, YRV genes are more likely to be associated with the nuclear lamina than non-YRV genes and are generally more likely to be close to each other in the nucleus (although not along chromosomes). Taken together, these results suggest that variation on the Y chromosome plays a role in modifying how the genome is distributed across chromatin compartments, either via changes in the distribution of DNA-binding proteins or via changes in the spatial arrangement of the genome in the nucleus.

Breakage-fusion-bridge cycles and large insertions contribute to the rapid evolution of accessory chromosomes in a fungal pathogen

Chromosomal rearrangements are a major driver of eukaryotic genome evolution, affecting speciation, pathogenicity and cancer progression. Changes in chromosome structure are often initiated by mis-repair of double-strand breaks in the DNA. Mis-repair is particularly likely when telomeres are lost or when dispersed repeats misalign during crossing-over. Fungi carry highly polymorphic chromosomal complements showing substantial variation in chromosome length and number. The mechanisms driving chromosome polymorphism in fungi are poorly understood. We aimed to identify mechanisms of chromosomal rearrangements in the fungal wheat pathogen Zymoseptoria tritici. We combined population genomic resequencing and chromosomal segment PCR assays with electrophoretic karyotyping and resequencing of parents and offspring from experimental crosses to show that this pathogen harbors a highly diverse complement of accessory chromosomes that exhibits strong global geographic differentiation in numbers and lengths of chromosomes. Homologous chromosomes carried highly differentiated gene contents due to numerous insertions and deletions. The largest accessory chromosome recently doubled in length through insertions totaling 380 kb. Based on comparative genomics, we identified the precise breakpoint locations of these insertions. Nondisjunction during meiosis led to chromosome losses in progeny of three different crosses. We showed that a new accessory chromosome emerged in two viable offspring through a fusion between sister chromatids. Such chromosome fusion is likely to initiate a breakage-fusion-bridge (BFB) cycle that can rapidly degenerate chromosomal structure. We suggest that the accessory chromosomes of Z. tritici originated mainly from ancient core chromosomes through a degeneration process that included BFB cycles, nondisjunction and mutational decay of duplicated sequences. The rapidly evolving accessory chromosome complement may serve as a cradle for adaptive evolution in this and other fungal pathogens.

Recent gene-capture on the UV sex chromosomes of the moss Ceratodon purpureus

Sex chromosomes evolve from ordinary autosomes through the expansion and subsequent degeneration of a region of suppressed recombination that is inherited through one sex. Here we investigate the relative timing of these processes in the UV sex chromosomes of the moss Ceratodon purpureus using molecular population genetic analyses of eight newly discovered sex-linked loci. In this system, recombination is suppressed on both the female-transmitted (U) sex chromosome and the male-transmitted (V) chromosome. Genes on both chromosomes therefore should show the deleterious effects of suppressed recombination and sex-limited transmission, while purifying selection should maintain homologs of genes essential for both sexes on both sex chromosomes. Based on analyses of eight sex-linked loci, we show that the nonrecombining portions of the U and V chromosomes expanded in at least two events (~0.6-1.3 MYA and ~2.8-3.5 MYA), after the divergence of C. purpureus from its dioecious sister species, Trichodon cylindricus and Cheilothela chloropus. Both U- and V-linked copies showed reduced nucleotide diversity and limited population structure, compared to autosomal loci, suggesting that the sex chromosomes experienced more recent selective sweeps that the autosomes. Collectively these results highlight the dynamic nature of gene composition and molecular evolution on nonrecombining portions of the U and V sex chromosomes.

Hidden magicians of genome evolution

Transposable elements (TEs) represent genome's dynamic component, causing mutations and genetic variations. Transposable elements can invade eukaryotic genomes in a short span; these are silenced by homology-dependent gene silencing and some functional parts of silenced elements are utilized to perform novel cellular functions. However, during the past two decades, major interest has been focused on the positive contribution of these elements in the evolution of genomes. The interaction between mobile DNAs and their host genomes are quite diverse, ranging from modifications of gene structure to alterations in general genome architecture and can be regarded as hidden magicians in shaping evolution of genomes. Some of the prominent examples that impressively demonstrate the beneficial impact of TEs on host biology over evolutionary time include their role in structure and functions of eukaryotic genomes.

Rapid turnover of the W chromosome in geographical populations of wild silkmoths, Samia cynthia ssp

Our previous studies revealed a considerably high level of chromosomal polymorphism in wild silkmoths, Samia cynthia ssp. (Lepidoptera: Saturniidae). Geographical populations of this species complex differ in chromosome numbers and show derived sex chromosome systems including Z0/ZZ in S. cynthia ricini (2n = 27/28; Vietnam), neo-Wneo-Z/neo-Zneo-Z in S. cynthia walkeri (2n = 26/26; Sapporo, Hokkaido) and neo-WZ1Z2/Z1Z1Z2Z2 in S. cynthia subsp. indet. (2n = 25/26; Nagano, Honshu). In this study, we collected specimens of S. cynthia pryeri in Japanese islands Kyushu, Shikoku and Honshu, with an ancestral-like karyotype of 2n = 28 in both sexes and a WZ/ZZ sex chromosome system, except for one population, in which females have lost the W chromosome. However, the S. cynthia pryeri W chromosome showed a very unusual morphology: It was composed of a highly heterochromatic body, which remained condensed throughout the whole cell cycle and of a euchromatin-like "tail." We examined molecular composition of the W and neo-W chromosomes in S. cynthia subspecies by comparative genomic hybridisation and fluorescence in situ hybridisation with W chromosome painting probes prepared from laser-microdissected W chromatin of S. cynthia pryeri. These methods revealed that the molecular composition of highly heterochromatic part of the S. cynthia pryeri W chromosome is very different and lacks homology in the genomes of other subspecies, whereas the euchromatin-like part of the W chromosome corresponds to a heterochromatic part of the neo-W chromosomes in S. cynthia walkeri and S. cynthia subsp. indet. Our findings suggest that the curious WZ system of S. cynthia pryeri may represent an ancestral state of the Samia species complex but do not exclude an alternative hypothesis of its derived origin.

Y-chromosome evolution: emerging insights into processes of Y-chromosome degeneration

The human Y chromosome is intriguing not only because it harbours the master-switch gene that determines gender but also because of its unusual evolutionary history. The Y chromosome evolved from an autosome, and its evolution has been characterized by massive gene decay. Recent whole-genome and transcriptome analyses of Y chromosomes in humans and other primates, in Drosophila species and in plants have shed light on the current gene content of the Y chromosome, its origins and its long-term fate. Furthermore, comparative analysis of young and old Y chromosomes has given further insights into the evolutionary and molecular forces triggering Y-chromosome degeneration and into the evolutionary destiny of the Y chromosome.

Genomic variation in natural populations of Drosophila melanogaster

This report of independent genome sequences of two natural populations of Drosophila melanogaster (37 from North America and 6 from Africa) provides unique insight into forces shaping genomic polymorphism and divergence. Evidence of interactions between natural selection and genetic linkage is abundant not only in centromere- and telomere-proximal regions, but also throughout the euchromatic arms. Linkage disequilibrium, which decays within 1 kbp, exhibits a strong bias toward coupling of the more frequent alleles and provides a high-resolution map of recombination rate. The juxtaposition of population genetics statistics in small genomic windows with gene structures and chromatin states yields a rich, high-resolution annotation, including the following: (1) 5'- and 3'-UTRs are enriched for regions of reduced polymorphism relative to lineage-specific divergence; (2) exons overlap with windows of excess relative polymorphism; (3) epigenetic marks associated with active transcription initiation sites overlap with regions of reduced relative polymorphism and relatively reduced estimates of the rate of recombination; (4) the rate of adaptive nonsynonymous fixation increases with the rate of crossing over per base pair; and (5) both duplications and deletions are enriched near origins of replication and their density correlates negatively with the rate of crossing over. Available demographic models of X and autosome descent cannot account for the increased divergence on the X and loss of diversity associated with the out-of-Africa migration. Comparison of the variation among these genomes to variation among genomes from D. simulans suggests that many targets of directional selection are shared between these species.

Large-scale introgression shapes the evolution of the mating-type chromosomes of the filamentous ascomycete Neurospora tetrasperma

The significance of introgression as an evolutionary force shaping natural populations is well established, especially in animal and plant systems. However, the abundance and size of introgression tracts, and to what degree interspecific gene flow is the result of adaptive processes, are largely unknown. In this study, we present medium coverage genomic data from species of the filamentous ascomycete Neurospora, and we use comparative genomics to investigate the introgression landscape at the genomic level in this model genus. We revealed one large introgression tract in each of the three investigated phylogenetic lineages of Neurospora tetrasperma (sizes of 5.6 Mbp, 5.2 Mbp, and 4.1 Mbp, respectively). The tract is located on the chromosome containing the locus conferring sexual identity, the mating-type (mat) chromosome. The region of introgression is confined to the region of suppressed recombination and is found on one of the two mat chromosomes (mat a). We used Bayesian concordance analyses to exclude incomplete lineage sorting as the cause for the observed pattern, and multilocus genealogies from additional species of Neurospora show that the introgression likely originates from two closely related, freely recombining, heterothallic species (N. hispaniola and N. crassa/N. perkinsii). Finally, we investigated patterns of molecular evolution of the mat chromosome in Neurospora, and we show that introgression is correlated with reduced level of molecular degeneration, consistent with a shorter time of recombination suppression. The chromosome specific (mat) and allele specific (mat a) introgression reported herein comprise the largest introgression tracts reported to date from natural populations. Furthermore, our data contradicts theoretical predictions that introgression should be less likely on sex-determining chromosomes. Taken together, the data presented herein advance our general understanding of introgression as a force shaping eukaryotic genomes.

Introgression is a process by which genetic material from one species becomes infiltrated into another, genetically distinct species. Introgression usually occurs via sexual reproduction: individuals of two species mate and produce a hybrid offspring, then the offspring repeatedly backcross with one of the parental species. Introgression has long been recognized as a key process in evolution, as it may contribute to speciation, diversification, and adaptation to new environments. The importance and prevalence of introgression has been well established in plant and animal systems, and in this study we use a fungal model system, Neurospora, to study the introgression at the genomic level. We gathered genomic data from six genomes, and by comparative genomics we revealed genetic transfer of DNA regions of unprecedentedly large sizes, covering over 50% of the mating-type chromosomes, and used phylogenetic analyses to reveal the origin and direction of the transfer. Introgression was found solely on the mating-type chromosomes, which contradicts theoretical predictions for sex-determining chromosomes. We argue that this unexpected pattern is due to the fact that fungi do not have differentiated sexes (female/male) and thereby are free from sex-biased evolutionary forces. Instead, we suggest that introgression between fungal species may result in reinvigoration of genomic regions exposed to suppressed recombination.

Patterns of molecular evolution of an avian neo-sex chromosome

Newer parts of sex chromosomes, neo-sex chromosomes, offer unique possibilities for studying gene degeneration and sequence evolution in response to loss of recombination and population size decrease. We have recently described a neo-sex chromosome system in Sylvioidea passerines that has resulted from a fusion between the first half (10 Mb) of chromosome 4a and the ancestral sex chromosomes. In this study, we report the results of molecular analyses of neo-Z and neo-W gametologs and intronic parts of neo-Z and autosomal genes on the second half of chromosome 4a in three species within different Sylvioidea lineages (Acrocephalidea, Timaliidae, and Alaudidae). In line with hypotheses of neo-sex chromosome evolution, we observe 1) lower genetic diversity of neo-Z genes compared with autosomal genes, 2) moderate synonymous and weak nonsynonymous sequence divergence between neo-Z and neo-W gametologs, and 3) lower GC content on neo-W than neo-Z gametologs. Phylogenetic reconstruction of eight neo-Z and neo-W gametologs suggests that recombination continued after the split of Alaudidae from the rest of the Sylvioidea lineages (i.e., after ~42.2 Ma) and with some exceptions also after the split of Acrocephalidea and Timaliidae (i.e., after ~39.4 Ma). The Sylvioidea neo-sex chromosome shares classical evolutionary features with the ancestral sex chromosomes but, as expected from its more recent origin, shows weaker divergence between gametologs.

The genetics of sex chromosomes: evolution and implications for hybrid incompatibility

Heteromorphic sex chromosomes, where one sex has two different types of sex chromosomes, face very different evolutionary consequences than do autosomes. Two important features of sex chromosomes arise from being present in only one copy in one of the sexes: dosage compensation and the meiotic silencing of sex chromosomes. Other differences arise because sex chromosomes spend unequal amounts of time in each sex. Thus, the impact of evolutionary processes (mutation, selection, genetic drift, and meiotic drive) differs substantially between each sex chromosome, and between the sex chromosomes and the autosomes. Sex chromosomes also play a disproportionate role in Haldane's rule and other important patterns related to hybrid incompatibility, and thus speciation. We review the consequences of sex chromosomes on hybrid incompatibility. A theme running through this review is that epigenetic processes, notably those related to chromatin, may be more important to the evolution of sex chromosomes and the evolution of hybrid incompatibility than previously recognized.

Contrasted patterns of molecular evolution in dominant and recessive self-incompatibility haplotypes in Arabidopsis

Self-incompatibility has been considered by geneticists a model system for reproductive biology and balancing selection, but our understanding of the genetic basis and evolution of this molecular lock-and-key system has remained limited by the extreme level of sequence divergence among haplotypes, resulting in a lack of appropriate genomic sequences. In this study, we report and analyze the full sequence of eleven distinct haplotypes of the self-incompatibility locus (S-locus) in two closely related Arabidopsis species, obtained from individual BAC libraries. We use this extensive dataset to highlight sharply contrasted patterns of molecular evolution of each of the two genes controlling self-incompatibility themselves, as well as of the genomic region surrounding them. We find strong collinearity of the flanking regions among haplotypes on each side of the S-locus together with high levels of sequence similarity. In contrast, the S-locus region itself shows spectacularly deep gene genealogies, high variability in size and gene organization, as well as complete absence of sequence similarity in intergenic sequences and striking accumulation of transposable elements. Of particular interest, we demonstrate that dominant and recessive S-haplotypes experience sharply contrasted patterns of molecular evolution. Indeed, dominant haplotypes exhibit larger size and a much higher density of transposable elements, being matched only by that in the centromere. Overall, these properties highlight that the S-locus presents many striking similarities with other regions involved in the determination of mating-types, such as sex chromosomes in animals or in plants, or the mating-type locus in fungi and green algae.

Self-incompatibility is a common genetic system preventing selfing through recognition and rejection of self-pollen in hermaphroditic flowering plants. In the Brassicaceae family, this system is controlled by a single genomic region, called the S-locus, where many distinct specificities segregate in natural populations. In this study, we obtained genomic sequences comprising the S-locus in two closely related Brassicaceae species, Arabidopsis lyrata and A. halleri, and analyzed their diversity and patterns of molecular evolution. We report compelling evidence that the S-locus presents many similar properties with other genomic regions involved in the determination of mating-types in mammals, insects, plants, or fungi. In particular, in spite of their diversity, these genomic regions all show absence of similarity in intergenic sequences, large depth of genealogies, highly divergent organization, and accumulation of transposable elements. Moreover, some of these features were found to vary according to dominance of the S-locus specificities, suggesting that dominance/recessivity interactions are key drivers of the evolution of this genomic region.

Preservation of the Y transcriptome in a 10-million-year-old plant sex chromosome system

Classical genetic studies discovered loss of genes from the ancient sex chromosome systems of several animals (genetic degeneration), and complete genome sequencing confirms that the heterogametic sex is hemizygous for most sex-linked genes. Genetic degeneration is thought to result from the absence of recombination between the sex chromosome pair (reviewed by [1]) and is very rapid after sex chromosome-autosome fusions in Drosophila [2-4]. Plant sex chromosome systems allow study of the time course of degeneration, because they evolved from a state wholly without sex chromosomes (rather than after a large genome region fused to a preexisting sex chromosome), and, in several taxa, recombination stopped very recently. However, despite increasing genetic and physical mapping of plant nonrecombining sex-determining regions [5-8], it remains very difficult to discover sex-linked genes, and it is unclear whether Y-linked genes are losing full function. We therefore developed a high-throughput method using RNA-Seq to identify sex linkage in Silene latifolia. Recombination suppression between this plant's XY sex chromosome pair started only about 10 million years ago [9]. Our approach identifies several hundred new sex-linked genes, and we show that this young Y chromosome retains many genes, yet these already have slightly reduced gene expression and are accumulating changes likely to reduce protein functions.

Molecular evolution in nonrecombining regions of the Drosophila melanogaster genome

We study the evolutionary effects of reduced recombination on the Drosophila melanogaster genome, analyzing more than 200 new genes that lack crossing-over and employing a novel orthology search among species of the melanogaster subgroup. These genes are located in the heterochromatin of chromosomes other than the dot (fourth) chromosome. Noncrossover regions of the genome all exhibited an elevated level of evolutionary divergence from D. yakuba at nonsynonymous sites, lower codon usage bias, lower GC content in coding and noncoding regions, and longer introns. Levels of gene expression are similar for genes in regions with and without crossing-over, which rules out the possibility that the reduced level of adaptation that we detect is caused by relaxed selection due to lower levels of gene expression in the heterochromatin. The patterns observed are consistent with a reduction in the efficacy of selection in all regions of the genome of D. melanogaster that lack crossing-over, as a result of the effects of enhanced Hill-Robertson interference. However, we also detected differences among nonrecombining locations: The X chromosome seems to exhibit the weakest effects, whereas the fourth chromosome and the heterochromatic genes on the autosomes located most proximal to the centromere showed the largest effects. However, signatures of selection on both nonsynonymous mutations and on codon usage persist in all heterochromatic regions.