Ribosomal DNA deletions modulate genome-wide gene expression: "rDNA-sensitive" genes and natural variation

The ribosomal rDNA gene array is an epigenetically-regulated repeated gene locus. While rDNA copy number varies widely between and within species, the functional consequences of subtle copy number polymorphisms have been largely unknown. Deletions in the Drosophila Y-linked rDNA modifies heterochromatin-induced position effect variegation (PEV), but it has been unknown if the euchromatic component of the genome is affected by rDNA copy number. Polymorphisms of naturally occurring Y chromosomes affect both euchromatin and heterochromatin, although the elements responsible for these effects are unknown. Here we show that copy number of the Y-linked rDNA array is a source of genome-wide variation in gene expression. Induced deletions in the rDNA affect the expression of hundreds to thousands of euchromatic genes throughout the genome of males and females. Although the affected genes are not physically clustered, we observed functional enrichments for genes whose protein products are located in the mitochondria and are involved in electron transport. The affected genes significantly overlap with genes affected by natural polymorphisms on Y chromosomes, suggesting that polymorphic rDNA copy number is an important determinant of gene expression diversity in natural populations. Altogether, our results indicate that subtle changes to rDNA copy number between individuals may contribute to biologically relevant phenotypic variation.     

DOES LOCAL ADAPTATION CAUSE HIGH POPULATION DIFFERENTIATION OF SILENE LATIFOLIA Y CHROMOSOMES?

Natural selection can reduce the effective population size of the nonrecombining Y chromosome, whereas local adaptation of Y-linked genes can increase the population divergence and overall intra-species polymorphism of Y-linked sequences. The plant Silene latifolia evolved a Y chromosome relatively recently, and most known X-linked genes have functional Y homologues, making the species interesting for comparisons of X- and Y-linked diversity and subdivision. Y-linked genes show higher population differentiation, compared to X-linked genes, and this might be maintained by local adaptation in Y-linked genes (or low sequence diversity). Here we attempt to test between these causes by investigating DNA polymorphism and population differentiation using a larger set of Y-linked and X-linked S. latifolia genes (than used previously), and show that net sequence divergence for Y-linked sequences (measured by D(a) , also known as δ) is low, and not consistently higher than X-linked genes. This does not support local adaptation, instead, the higher values of differentiation measures for the Y-linked genes probably result largely from reduced total variation on the Y chromosome, which in turn reflect deterministic processes lowering effective population sizes of evolving Y-chromosomes.     

Indirect evidence from DNA sequence diversity for genetic degeneration of the Y-chromosome in dioecious species of the plant Silene: the SlY4/SlX4 and DD44-X/DD44-Y gene pairs

The action of natural selection is expected to reduce the effective population size of a nonrecombining chromosome, and this is thought to be the chief factor leading to genetic degeneration of Y-chromosomes, which cease recombining during their evolution from ordinary chromosomes. Low effective population size of Y chromosomes can be tested by studying DNA sequence diversity of Y-linked genes. In the dioecious plant, Silene latifolia, which has sex chromosomes, one comparison (SlX1 vs. SlY1) indeed finds lower Y diversity compared with the homologous X-linked gene, and one Y-linked gene with no X-linked homologue has lower species-wide diversity than a homologous autosomal copy (SlAp3Y vs. SlAp3A). To test whether this is a general pattern for Y-linked genes, we studied two further recently described X and Y homologous gene pairs in samples from several populations of S. latifolia and S. dioica. Diversity is reduced for both Y-linked genes, compared with their X-linked homologues. Our new data are analysed to show that the low Y effective size cannot be explained by different levels of gene flow for the X vs. the Y chromosomes, either between populations or between these closely related species. Thus, all four Y-linked genes that have now been studied in these plants (the two studied here, and two previously studied genes, have low diversity). This supports other evidence for an ongoing degeneration process in these species.     

Unequal rates of Y chromosome gene divergence during speciation of the family Ursidae

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.     

Evolution of the Ribosomal DNA Spacers of Drosophila melanogaster : Different Patterns of Variation on X and Y Chromosomes

Length variation of the ribosomal gene spacers of Drosophila melanogaster was studied. Analysis of 47 X chromosomal and 47 Y chromosomal linked rDNA arrays collected from five continents indicates that the arrays on the two chromosomes differ qualitatively. The Y-linked arrays from around the world share little or no similarity for either their overall length or the organization of their spacers. Most of the X-linked arrays do, however, share a major length spacer of 5.1 kb. In addition, those X-linked arrays that have a major 5.1-kb band have similar spacer organization as demonstrated by genomic DNA digestions with several restriction enzymes. These data strongly support the hypothesis that spacer length patterns on only X-linked genes are maintained primarily by natural selection.     

Characterization of a Small Family (CAIII) of Microsatellite-Containing Sequences with X-Y Homology

Four X-linked loci showing homology with a previously described Y-linked polymorphic locus (DYS413) were identified and characterized. By fluorescent in situ hybridization (FISH), somatic cell hybrids, and YAC screening, the X-linked members of this small family of sequences (CAIII) all map in Xp22, while the Y members map in Yq11. These loci contribute to the overall similarity of the two genomic regions. All of the CAIII loci contain an internal microsatellite of the (CA)_n type. The microsatellites display extensive length polymorphism in two of the X-linked members as well as in the Y members. In addition, common sequence variants are found in the portions flanking the microsatellites in two of the X-linked members. Our results indicate that, during the evolution of this family, length variation on the Y chromosome was accumulated at a rate not slower than that on the X chromosome. Finally, these sequences represent a model system with which to analyze human populations for similar X- and Y-linked polymorphisms. Received: 29 July 1996 / Accepted: 15 January 1997     

Evidence for degeneration of the Y chromosome in the dioecious plant Silene latifolia

The human Y--probably because of its nonrecombining nature--has lost 97% of its genes since X and Y chromosomes started to diverge [1, 2]. There are clear signs of degeneration in the Drosophila miranda neoY chromosome (an autosome fused to the Y chromosome), with neoY genes showing faster protein evolution [3-6], accumulation of unpreferred codons [6], more insertions of transposable elements [5, 7], and lower levels of expression [8] than neoX genes. In the many other taxa with sex chromosomes, Y degeneration has hardly been studied. In plants, many genes are expressed in pollen [9], and strong pollen selection may oppose the degeneration of plant Y chromosomes [10]. Silene latifolia is a dioecious plant with young heteromorphic sex chromosomes [11, 12]. Here we test whether the S. latifolia Y chromosome is undergoing genetic degeneration by analyzing seven sex-linked genes. S. latifolia Y-linked genes tend to evolve faster at the protein level than their X-linked homologs, and they have lower expression levels. Several Y gene introns have increased in length, with evidence for transposable-element accumulation. We detect signs of degeneration in most of the Y-linked gene sequences analyzed, similar to those of animal Y-linked and neo-Y chromosome genes.     

Nucleotide diversity in Silene latifolia autosomal and sex-linked genes

The plant Silene latifolia has separate sexes and sex chromosomes, and is of interest for studying the early stages of sex chromosome evolution, especially the evolution of non-recombining regions on the Y chromosome. Hitch-hiking processes associated with ongoing genetic degeneration of the non-recombining Y chromosome are predicted to reduce Y-linked genes'' effective population sizes, and S. latifolia Y-linked genes indeed have lower diversity than X-linked ones. We tested whether this represents a true diversity reduction on the Y, versus the alternative possibility, elevated diversity at X-linked genes, by collecting new data on nucleotide diversity for autosomal genes, which had previously been little studied. We find clear evidence that Y-linked genes have reduced diversity. However, another alternative explanation for a low Y effective size is a high variance in male reproductive success. Autosomal genes should then also have lower diversity than expected, relative to the X, but this is not found in our loci. Taking into account the higher mutation rate of Y-linked genes, their low sequence diversity indicates a strong effect of within-population hitch-hiking on the Y chromosome.     

Analysis and evolution of two functional Y-linked loci in a plant sex chromosome system.

White campion (Silene latifolia) is one of the few examples of plants with separate sexes and with X and Y sex chromosomes. The presence or absence of the Y chromosome determines which type of reproductive organs--male or female--will develop. Recently, we characterized the first active gene located on a plant Y chromosome, SlY1, and its X-linked homolog, SlX1. These genes encode WD-repeat proteins likely to be involved in cell proliferation. Here, we report the characterization of a novel Y-linked gene, SlY4, which also has a homolog on the X chromosome, SlX4. Both SlY4 and SlX4 potentially encode fructose-2,6-bisphosphatases. A comparative molecular analysis of the two sex-linked loci (SlY1/SlX1 and SlY4/SlX4) suggests selective constraint on both X- and Y-linked genes and thus that both X- and Y-linked copies are functional. Divergence between SlY4 and SlX4 is much greater than that between the SlY1 and SlX1 genes. These results suggest that, as for human XY-linked genes, the sex-linked plant loci ceased recombining at different times and reveal distinct events in the evolutionary history of the sex chromosomes.     

Molecular characterization of ribosomal gene variation within and among NORs segregating in specialized populations of chicken.

The molecular organization of the 18S, 5.8S, and 28S ribosomal RNA gene repeat units, located at the single nucleolus organizer region (NOR) locus in the chicken, was investigated in genetically distinct populations of research and commercial chickens. Substantial gene repeat variation within and among NORs was documented. Intact ribosomal gene repeat size ranged from 11 kb to over 50 kb. Unique combinations of ribosomal genes, of different size, were specific to particular populations. It was determined that the basis for the ribosomal gene repeat size variation was intergenic spacer (IGS) length heterogeneity. Interestingly, in different populations, the location of the variation that contributes to length heterogeneity was specific to particular IGS subregions. In addition to IGS variation, an inbred line of Red Jungle Fowl exhibited coding region variation. Ribosomal gene copy number variation was also studied, and line averages ranged from 279 to 368. Average rDNA array size (a function of copy number and gene repeat length) was calculated for each of the populations and found to vary over a range of two megabases, from 5 to 7 Mb.     

Positive and negative selection on mammalian Y chromosomes

Y chromosomes are genetically degenerate in most organisms studied. The loss of genes from Y chromosomes is thought to be due to the inefficiency of purifying selection in nonrecombining regions, which leads to the accumulation of deleterious mutations via the processes of hitchhiking, background selection, and Muller's ratchet. As the severity of these processes depends on the number of functional genes linked together on the nonrecombining Y, it is not clear whether these processes are still at work on the old, gene-poor mammalian Y chromosomes. If purifying selection is indeed less efficient in the Y-linked, compared to the X-linked genes, deleterious nonsynonymous substitutions are expected to accumulate faster on the Y chromosome. However, positive selection on Y-linked genes could also increase the rate of amino acid-changing substitutions. Thus, the previous reports of an elevated nonsynonymous substitution rate in Y-linked genes are still open to interpretation. Here, we report evidence for positive selection in two out of three studied mammalian Y-linked genes, suggesting that adaptive Darwinian evolution may be common on mammalian Y chromosomes. Taking positive selection into account, we demonstrate that purifying selection is less efficient in mammalian Y-linked genes compared to their X-linked homologues, suggesting that these genes continue to degenerate.     

Plant sex chromosomes: molecular structure and function

Recent molecular and genomic studies carried out in a number of model dioecious plant species, including Asparagus officinalis, Carica papaya, Silene latifolia, Rumex acetosa and Marchantia polymorpha, have shed light on the molecular structure of both homomorphic and heteromorphic sex chromosomes, and also on the gene functions they have maintained since their evolution from a pair of autosomes. The molecular structure of sex chromosomes in species from different plant families represents the evolutionary pathway followed by sex chromosomes during their evolution. The degree of Y chromosome degeneration that accompanies the suppression of recombination between the Xs and Ys differs among species. The primitive Ys of A. officinalis and C. papaya have only diverged from their homomorphic Xs in a short male-specific and non-recombining region (MSY), while the heteromorphic Ys of S. latifolia, R. acetosa and M. polymorpha have diverged from their respective Xs. As in the Y chromosomes of mammals and Drosophila, the accumulation of repetitive DNA, including both transposable elements and satellite DNA, has played an important role in the divergence and size enlargement of plant Ys, and consequently in reducing gene density. Nevertheless, the degeneration process in plants does not appear to have reached the Y-linked genes. Although a low gene density has been found in the sequenced Y chromosome of M. polymorpha, most of its genes are essential and are expressed in the vegetative and reproductive organs in both male and females. Similarly, most of the Y-linked genes that have been isolated and characterized up to now in S. latifolia are housekeeping genes that have X-linked homologues, and are therefore expressed in both males and females. Only one of them seems to be degenerate with respect to its homologous region in the X. Sequence analysis of larger regions in the homomorphic X and Y chromosomes of papaya and asparagus, and also in the heteromorphic sex chromosomes of S. latifolia and R. acetosa, will reveal the degenerative changes that the Y-linked gene functions have experienced during sex chromosome evolution.     

DNA diversity in sex-linked and autosomal genes of the plant species Silene latifolia and Silene dioica

The relatively recent origin of sex chromosomes in the plant genus Silene provides an opportunity to study the early stages of sex chromosome evolution and, potentially, to test between the different population genetic processes likely to operate in nonrecombining chromosomes such as Y chromosomes. We previously reported much lower nucleotide polymorphism in a Y-linked gene (SlY1) of the plant Silene latifolia than in the homologous X-linked gene (SlX1). Here, we report a more extensive study of nucleotide diversity in these sex-linked genes, including a larger S. latifolia sample and a sample from the closely related species Silene dioica, and we also study the diversity of an autosomal gene, CCLS37.1. We demonstrate that nucleotide diversity in the Y-linked genes of both S. latifolia and S. dioica is very low compared with that of the X-linked gene. However, the autosomal gene also has low DNA polymorphism, which may be due to a selective sweep. We use a single individual of the related hermaphrodite species Silene conica, as an outgroup to show that the low SlY1 diversity is not due to a lower mutation rate than that for the X-linked gene. We also investigate several other possibilities for the low SlY1 diversity, including differential gene flow between the two species for Y-linked, X-linked, and autosomal genes. The frequency spectrum of nucleotide polymorphism on the Y chromosome deviates significantly from that expected under a selective-sweep model. However, we detect population subdivision in both S. latifolia and S. dioica, so it is not simple to test for selective sweeps. We also discuss the possibility that Y-linked diversity is reduced due to highly variable male reproductive success, and we conclude that this explanation is unlikely.     

Localization of Murine X and Autosomal Sequences Homologous to the Human Y Located Testis-Determining Region

Recently a candidate gene for the primary testis-determining factor (TDF) encoding a zinc finger protein (ZFY) has been cloned from the human Y chromosome. A highly homologous X-linked copy has also been identified. Using this human sequence it is possible to identify two Y loci, an X and an autosomal locus in the mouse (Zfy-1, Zfy-2, Zfx and Zfa, respectively). Suprisingly ZFY is more homologous to the mouse X and autosomal sequences than it is to either of the Y-linked loci. Both Zfy-1 and Zfy-2 are present in the Sxr region of the Y but Zfy-2 is absent in the Sxr deletion variant Sxrb (or Sxr") suggesting it is not necessary for male determination. Extensive backcross analyses map Zfa to mouse chromosome 10 and Zfx to a 5-cM interval between anonymous X probe MDXS120 and the tabby locus (Ta). We also show that the mouse androgen receptor locus (m-AR) believed to underlie the testicular feminization mutation (Tfm) shows complete linkage to Zfx. Comparative mapping indicates that in man these genes lie in separate conserved DNA segments.     

The Molecular through Ecological Genetics of Abnormal Abdomen in DROSOPHILA MERCATORUM. I. Basic Genetics

The abnormal abdomen syndrome (aa) in Drosophila mercatorum is characterized by the persistence of juvenilized cuticle on the adult abdomen. The aa phenotype is shown to depend on at least two X-linked genetic elements that are about one map unit apart near the centromeric end of the X chromosome. These two genetic elements are necessary for aa expression; one behaves as a dominant element and the other as a recessive. Overlaying these genetic studies upon molecular work reported elsewhere, it is argued that the dominant element is the presence of a 5 kb insertion in a majority of the X-linked repeats coding for the 28S ribosomal RNA. The recessive element appears to be a locus controlling differential replication of noninserted over inserted 28S genes during polytenization. The aa syndrome requires both the presence of the inserted repeats and the failure to preferentially amplify noninserted repeats. Given the necessary X-linked elements for aa, a variety of modifiers are revealed. First, aa expression in males is Y-linked, apparently corresponding to a deletion of the 18S/28S rDNA gene cluster normally found on the Y. Moreover, all major autosomes can modify the penetrance of aa.     

Chromosomal Copy Number Variation in Saccharomyces pastorianus Is Evidence for Extensive Genome Dynamics in Industrial Lager Brewing Strains

Lager brewing strains of Saccharomyces pastorianus are natural interspecific hybrids originating from the spontaneous hybridization of Saccharomyces cerevisiae and Saccharomyces eubayanus. Over the past 500 years, S. pastorianus has been domesticated to become one of the most important industrial microorganisms. Production of lager-type beers requires a set of essential phenotypes, including the ability to ferment maltose and maltotriose at low temperature, the production of flavors and aromas, and the ability to flocculate. Understanding of the molecular basis of complex brewing-related phenotypic traits is a prerequisite for rational strain improvement. While genome sequences have been reported, the variability and dynamics of S. pastorianus genomes have not been investigated in detail. Here, using deep sequencing and chromosome copy number analysis, we showed that S. pastorianus strain CBS1483 exhibited extensive aneuploidy. This was confirmed by quantitative PCR and by flow cytometry. As a direct consequence of this aneuploidy, a massive number of sequence variants was identified, leading to at least 1,800 additional protein variants in S. pastorianus CBS1483. Analysis of eight additional S. pastorianus strains revealed that the previously defined group I strains showed comparable karyotypes, while group II strains showed large interstrain karyotypic variability. Comparison of three strains with nearly identical genome sequences revealed substantial chromosome copy number variation, which may contribute to strain-specific phenotypic traits. The observed variability of lager yeast genomes demonstrates that systematic linking of genotype to phenotype requires a three-dimensional genome analysis encompassing physical chromosomal structures, the copy number of individual chromosomes or chromosomal regions, and the allelic variation of copies of individual genes.     

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.