Dosage compensation and inverse effects in triple X metafemales of Drosophila

Drift drives the evolution of chromosome number II: The impact of range size on genome evolution in Carnivora

Chromosome number is a fundamental genomic trait that is often the first recorded characteristic of a genome. Across large clades, a common pattern emerges: many or even most lineages exhibit relative stasis, while a handful of lineages or species exhibit striking variation. Despite recent developments in comparative methods, most of this heterogeneity is still poorly understood. It is essential to understand why some lineages have rapid rates of chromosome number evolution, as it can impact a variety of other traits. Previous research suggests that biased female meiotic drive may shape rates of karyotype evolution in some mammals. However, Carnivora exhibits variation that this female meiotic drive model cannot explain. We hypothesize that variation in effective population size may underlie rate variation in Carnivora. To test this hypothesis, we estimated rates of fusions and fissions while accounting for range size, which we use as a proxy for effective population size. We reason fusions and fissions are deleterious or underdominant and that only in lineages with small range sizes will these changes be able to fix due to genetic drift. In this study, we find that the rates of fusions and fissions are elevated in taxa with small range sizes relative to those with large range sizes. Based on these findings, we conclude that 1) naturally occurring structural mutations that change chromosome number are underdominant or mildly deleterious, and 2) when population sizes are small, structural rearrangements may play an important role in speciation and reduction in gene flow among populations.

Inverse and Proportional Trans Modulation of Gene Expression in Human Aneuploidies

Genomic imbalance in aneuploidy is often detrimental to organisms. To gain insight into the molecular basis of aneuploidies in humans, we analyzed transcriptome data from several autosomal and sex chromosome aneuploidies. The results showed that in human aneuploid cells, genes located on unvaried chromosomes are inversely or proportionally trans-modulated, while a subset of genes on the varied chromosomes are compensated. Less genome-wide modulation is found for sex chromosome aneuploidy compared with autosomal aneuploidy due to X inactivation and the retention of dosage sensitive regulators on both sex chromosomes to limit the effective dosage change. We also found that lncRNA and mRNA can have different responses to aneuploidy. Furthermore, we analyzed the relationship between dosage-sensitive transcription factors and their targets, which illustrated the modulations and indicates genomic imbalance is related to stoichiometric changes in components of gene regulatory complexes.In summary, this study demonstrates the existence of trans-acting effects and compensation mechanisms in human aneuploidies and contributes to our understanding of gene expression regulation in unbalanced genomes and disease states.

Characteristics and expression of lncRNA and transposable elements in Drosophila aneuploidy

Aneuploidy can globally affect the expression of the whole genome, which is detrimental to organisms. Dosage-sensitive regulators usually have multiple intermolecular interactions, and changes in their stoichiometry are responsible for the dysregulation of the regulatory network. Currently, studies on noncoding genes in aneuploidy are relatively rare. We studied the characteristics and expression profiles of long noncoding RNAs (lncRNAs) and transposable elements (TEs) in aneuploid Drosophila. It is found that lncRNAs and TEs are affected by genomic imbalance and appear to be more sensitive to an inverse dosage effect than mRNAs. Several dosage-sensitive lncRNAs and TEs were detected for their expression patterns during embryogenesis, and their biological functions in the ovary and testes were investigated using tissue-specific RNAi. This study advances our understanding of the noncoding sequences in imbalanced genomes and provides a novel perspective for the study of aneuploidy-related human diseases such as cancer.

Genome-Wide Expression Analysis of Long Noncoding RNAs and Their Target Genes in Metafemale Drosophila

Aneuploidy is usually more detrimental than altered ploidy of the entire set of chromosomes. To explore the regulatory mechanism of gene expression in aneuploidy, we analyzed the transcriptome sequencing data of metafemale Drosophila. The results showed that most genes on the X chromosome undergo dosage compensation, while the genes on the autosomal chromosomes mainly present inverse dosage effects. Furthermore, long noncoding RNAs (lncRNAs) have been identified as key regulators of gene expression, and they are more sensitive to dosage changes than mRNAs. We analyzed differentially expressed mRNAs (DEGs) and differentially expressed lncRNAs (DELs) in metafemale Drosophila and performed functional enrichment analyses of DEGs and the target genes of DELs, and we found that they are involved in several important biological processes. By constructing lncRNA-mRNA interaction networks and calculating the maximal clique centrality (MCC) value of each node in the network, we also identified two key candidate lncRNAs (CR43940 and CR42765), and two of their target genes, Sin3A and MED1, were identified as inverse dosage modulators. These results suggest that lncRNAs play an important role in the regulation of genomic imbalances. This study may deepen the understanding of the gene expression regulatory mechanisms in aneuploidy from the perspective of lncRNAs.

Integrative epigenome profiling of 47XXY provides insights into whole genomic DNA hypermethylation and active chromatin accessibility

Klinefelter syndrome (KS, 47XXY) is a disorder characterized by sex chromosomal aneuploidy, which may lead to changes in epigenetic regulations of gene expression. To define epigenetic architectures in 47XXY, we annotated DNA methylation in euploid males (46XY) and females (46XX), and 47XXY individuals using whole genome bisulfite sequencing (WGBS) and integrated chromatin accessbilty, and detected abnormal hypermethylation in 47XXY. Furthermore, we detected altered chromatin accessibility in 47XXY, in particular in chromosome X, using Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) in cultured amniotic cells. Our results construct the whole genome-wide DNA methylation map in 47XXY, and provide new insights into the early epigenomic dysregulation resulting from an extra chromosome X in 47XXY.

Genome-Wide Effects on Gene Expression Between Parental and Filial Generations of Trisomy 11 and 12 of Rice

Aneuploid refers to the gene dosage imbalance due to copy number alterations. Aneuploidy is generally harmful to the growth, development and reproduction of organisms according to the numerous research. However, it has rarely been reported on whether aneuploid have a relevant pattern of genome expression between the parental and its offspring generations. In this study, mRNA sequencing analysis was performed on rice (Oryza sativa L.) primary trisomes 11 and 12, same primary trisomes and normal individuals in their filial generation. We systematically summarized the changes in gene expression patterns that occur on cis genes and on trans genes between parental and filial generations. In T11 and T12, the ratio of cis-gene expression showed intermediate type in parents and dosage compensation in filial generations, which maybe due to more genes being downregulated. The trans genes were also affected by aneuploidy and manifested as cis-related. The strains with normal chromosomes in filial generations, there are still aneuploid-sensitive genes differentially expressed in their genomes, indicating that the effect of aneuploidy is far-reaching and could not be easily eliminated. Meanwhile, among these differentially expressed genes, genes with low-expression level were more likely to be upregulated, while genes with medium- and high-expression level were easy to be downregulated. For the different types of rice aneuploid, upregulated genes were mainly associated with genomic imbalance while downregulated genes were mainly influenced by the specific added chromosome. In conclusion, our results provide new insights into the genetic characterization and evolution of biological aneuploidy genomes.

Four layer multi-omics reveals molecular responses to aneuploidy in Leishmania

Aneuploidy causes system-wide disruptions in the stochiometric balances of transcripts, proteins, and metabolites, often resulting in detrimental effects for the organism. The protozoan parasite Leishmania has an unusually high tolerance for aneuploidy, but the molecular and functional consequences for the pathogen remain poorly understood. Here, we addressed this question in vitro and present the first integrated analysis of the genome, transcriptome, proteome, and metabolome of highly aneuploid Leishmania donovani strains. Our analyses unambiguously establish that aneuploidy in Leishmania proportionally impacts the average transcript- and protein abundance levels of affected chromosomes, ultimately correlating with the degree of metabolic differences between closely related aneuploid strains. This proportionality was present in both proliferative and non-proliferative in vitro promastigotes. However, as in other Eukaryotes, we observed attenuation of dosage effects for protein complex subunits and in addition, non-cytoplasmic proteins. Differentially expressed transcripts and proteins between aneuploid Leishmania strains also originated from non-aneuploid chromosomes. At protein level, these were enriched for proteins involved in protein metabolism, such as chaperones and chaperonins, peptidases, and heat-shock proteins. In conclusion, our results further support the view that aneuploidy in Leishmania can be adaptive. Additionally, we believe that the high karyotype diversity in vitro and absence of classical transcriptional regulation make Leishmania an attractive model to study processes of protein homeostasis in the context of aneuploidy and beyond.

Dosage Compensation in Drosophila: Its Canonical and Non-Canonical Mechanisms

Dosage compensation equalizes gene expression in a single male X chromosome with that in the pairs of autosomes and female X chromosomes. In the fruit fly Drosophila, canonical dosage compensation is implemented by the male-specific lethal (MSL) complex functioning in all male somatic cells. This complex contains acetyl transferase males absent on the first (MOF), which performs H4K16 hyperacetylation specifically in the male X chromosome, thus facilitating transcription of the X-linked genes. However, accumulating evidence points to an existence of additional, non-canonical dosage compensation mechanisms operating in somatic and germline cells. In this review, we discuss current advances in the understanding of both canonical and non-canonical mechanisms of dosage compensation in Drosophila.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

Pseudoautosomal Region 1 Overdosage Affects the Global Transcriptome in iPSCs From Patients With Klinefelter Syndrome and High-Grade X Chromosome Aneuploidies

Klinefelter syndrome (KS) is the most prevalent aneuploidy in males and is characterized by a 47,XXY karyotype. Less frequently, higher grade sex chromosome aneuploidies (HGAs) can also occur. Here, using a paradigmatic cohort of KS and HGA induced pluripotent stem cells (iPSCs) carrying 49,XXXXY, 48,XXXY, and 47,XXY karyotypes, we identified the genes within the pseudoautosomal region 1 (PAR1) as the most susceptible to dosage-dependent transcriptional dysregulation and therefore potentially responsible for the progressively worsening phenotype in higher grade X aneuploidies. By contrast, the biallelically expressed non-PAR escape genes displayed high interclonal and interpatient variability in iPSCs and differentiated derivatives, suggesting that these genes could be associated with variable KS traits. By interrogating KS and HGA iPSCs at the single-cell resolution we showed that PAR1 and non-PAR escape genes are not only resilient to the X-inactive specific transcript (XIST)-mediated inactivation but also that their transcriptional regulation is disjointed from the absolute XIST expression level. Finally, we explored the transcriptional effects of X chromosome overdosage on autosomes and identified the nuclear respiratory factor 1 (NRF1) as a key regulator of the zinc finger protein X-linked (ZFX). Our study provides the first evidence of an X-dosage-sensitive autosomal transcription factor regulating an X-linked gene in low- and high-grade X aneuploidies.

Dissecting cricket genomes for the advancement of entomology and entomophagy

Significant advances in biophysical methods such as next-generation sequencing technologies have now opened the way to conduct evolutionary and applied research based on the genomic information of greatly diverse insects. Crickets belonging to Orthoptera (Insecta: Polyneoptera), one of the most flourishing groups of insects, have contributed to the development of multiple scientific fields including developmental biology and neuroscience and have been attractive targets in evolutionary ecology for their diverse ecological niches. In addition, crickets have recently gained recognition as food and feed. However, the genomic information underlying their biological basis and application research toward breeding is currently underrepresented. In this review, we summarize the progress of genomics of crickets. First, we outline the phylogenetic position of crickets in insects and then introduce recent studies on cricket genomics and transcriptomics in a variety of fields. Furthermore, we present findings from our analysis of polyneopteran genomes, with a particular focus on their large genome sizes, chromosome number, and repetitive sequences. Finally, how the cricket genome can be beneficial to the food industry is discussed. This review is expected to enhance greater recognition of how important the cricket genomes are to the multiple biological fields and how basic research based on cricket genome information can contribute to tackling global food security.

MYC dosage compensation is mediated by miRNA-transcription factor interactions in aneuploid cancer

We hypothesize that dosage compensation of critical genes arises from systems-level properties for cancer cells to withstand the negative effects of aneuploidy. We identified several candidate genes in cancer multiomics data and developed a biocomputational platform to construct a mathematical model of their interaction network with micro-RNAs and transcription factors, where the property of dosage compensation emerged for MYC and was dependent on the kinetic parameters of its feedback interactions with three micro-RNAs. These circuits were experimentally validated using a genetic tug-of-war technique to overexpress an exogenous MYC, leading to overexpression of the three microRNAs involved and downregulation of endogenous MYC. In addition, MYC overexpression or inhibition of its compensating miRNAs led to dosage-dependent cytotoxicity in MYC-amplified colon cancer cells. Finally, we identified negative correlation of MYC dosage compensation with patient survival in TCGA breast cancer patients, highlighting the potential of this mechanism to prevent aneuploid cancer progression.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

Modulation of Global Gene Expression by Aneuploidy and CNV of Dosage Sensitive Regulatory Genes

Aneuploidy, which disrupts the genetic balance due to partial genome dosage changes, is usually more detrimental than euploidy variation. To investigate the modulation of gene expression in aneuploidy, we analyzed the transcriptome sequencing data of autosomal and sex chromosome trisomy in Drosophila. The results showed that most genes on the varied chromosome (cis) present dosage compensation, while the remainder of the genome (trans) produce widespread inverse dosage effects. Some altered functions and pathways were identified as the common characteristics of aneuploidy, and several possible regulatory genes were screened for an inverse dosage effect. Furthermore, we demonstrated that dosage changes of inverse regulator Inr-a/pcf11 can produce a genome-wide inverse dosage effect. All these findings suggest that the mechanism of genomic imbalance is related to the changes in the stoichiometric relationships of macromolecular complex members that affect the overall function. These studies may deepen the understanding of gene expression regulatory mechanisms.

Interaction of Male Specific Lethal complex and genomic imbalance on global gene expression in Drosophila

The inverse dosage effect caused by chromosome number variations shows global consequences in genomic imbalance including sexual dimorphism and an X chromosome-specific response. To investigate the relationship of the MSL complex to genomic imbalance, we over-expressed MSL2 in autosomal and sex chromosomal aneuploids, and analyzed the different transcriptomes. Some candidate genes involved in regulatory mechanisms have also been tested during embryogenesis using TSA-FISH. Here we show that the de novo MSL complex assembled on the X chromosomes in females further reduced the global expression level on the basis of 2/3 down-regulation caused by the inverse dosage effect in trisomy through epigenetic modulations rather than induced dosage compensation. Plus, the sexual dimorphism effect in unbalanced genomes was further examined due to the pre-existing of the MSL complex in males. All these results demonstrate the dynamic functions of the MSL complex on global gene expression in different aneuploid genomes.

Lack of MECP2 gene transcription on the duplicated alleles of two related asymptomatic females with Xq28 duplications and opposite X-chromosome inactivation skewing

Xq28 duplication syndrome (MIM# 300815) is a severe neurodevelopmental disorder in males due to MeCP2 overexpression. Most females with MECP2 duplication are asymptomatic carriers, but there are phenotypic heterogeneities. Skewed X-chromosome inactivation (XCI) can protect females from exhibiting clinical phenotypes. Herein we reported two asymptomatic females (mother and grandmother) with interstitial Xq28 duplication. AR and RP2 assays showed that both had extremely skewed XCI, the Xq28 duplicated chromosome was inactivated in the mother, but was surprisingly activated in the grandmother. Interestingly, by combining RNA sequencing and whole-exome sequencing, we confirmed that XIST only expressed in the Xq28 duplication chromosomes of the two females, indicating that the Xq28 duplication chromosomes were inactive. Meanwhile, MECP2 and most XCI genes in the duplicated X-chromosomes were not transcriptionally expressed or upregulated, precluding major clinical phenotypes in the two females, especially the grandmother. We showed that XCI status detected using RNA sequencing was more relevant for establishing the clinical phenotype of MECP2 duplication in females. It suggested that there were other factors maintaining the XCI status in addition to DNA methylation, a possible additional inhibition mechanism occurred at the transcriptional level in the unmethylated X-chromosome, counter balancing the MECP2 duplication's detrimental phenotype effects.

Rapid Evolution of Autosomal Binding Sites of the Dosage Compensation Complex in Drosophila melanogaster and Its Association With Transcription Divergence

How pleiotropy influences evolution of protein sequence remains unclear. The male-specific lethal (MSL) complex in Drosophila mediates dosage compensation by 2-fold upregulation of the X chromosome in males. Nevertheless, several MSL proteins also bind autosomes and likely perform functions not related to dosage compensation. Here, we study the evolution of MOF, MSL1, and MSL2 biding sites in Drosophila melanogaster and its close relative Drosophila simulans. We found pervasive expansion of the MSL binding sites in D. melanogaster, particularly on autosomes. The majority of these newly-bound regions are unlikely to function in dosage compensation and associated with an increase in expression divergence between D. melanogaster and D. simulans. While dosage-compensation related sites show clear signatures of adaptive evolution, these signatures are even more marked among autosomal regions. Our study points to an intriguing avenue of investigation of pleiotropy as a mechanism promoting rapid protein sequence evolution.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

An empirical bayesian approach for testing gene expression fold change and its application in detecting global dosage effects

We are motivated by biological studies intended to understand global gene expression fold change. Biologists have generally adopted a fixed cutoff to determine the significance of fold changes in gene expression studies (e.g. by using an observed fold change equal to two as a fixed threshold). Scientists can also use a t-test or a modified differential expression test to assess the significance of fold changes. However, these methods either fail to take advantage of the high dimensionality of gene expression data or fail to test fold change directly. Our research develops a new empirical Bayesian approach to substantially improve the power and accuracy of fold-change detection. Specifically, we more accurately estimate gene-wise error variation in the log of fold change. We then adopt a t-test with adjusted degrees of freedom for significance assessment. We apply our method to a dosage study in Arabidopsis and a Down syndrome study in humans to illustrate the utility of our approach. We also present a simulation study based on real datasets to demonstrate the accuracy of our method relative to error variance estimation and power in fold-change detection. Our developed R package with a detailed user manual is publicly available on GitHub at https://github.com/cuiyingbeicheng/Foldseq.

Lineage-specific patterns of chromosome evolution are the rule not the exception in Polyneoptera insects

The structure of a genome can be described at its simplest by the number of chromosomes and the sex chromosome system it contains. Despite over a century of study, the evolution of genome structure on this scale remains recalcitrant to broad generalizations that can be applied across clades. To address this issue, we have assembled a dataset of 823 karyotypes from the insect group Polyneoptera. This group contains orders with a range of variations in chromosome number, and offer the opportunity to explore the possible causes of these differences. We have analysed these data using both phylogenetic and taxonomic approaches. Our analysis allows us to assess the importance of rates of evolution, phylogenetic history, sex chromosome systems, parthenogenesis and genome size on variation in chromosome number within clades. We find that fusions play a key role in the origin of new sex chromosomes, and that orders exhibit striking differences in rates of fusions, fissions and polyploidy. Our results suggest that the difficulty in finding consistent rules that govern evolution at this scale may be due to the presence of many interacting forces that can lead to variation among groups.

Evidence for Dosage Compensation in Coccinia grandis, a Plant with a Highly Heteromorphic XY System

About 15,000 angiosperms are dioecious, but the mechanisms of sex determination in plants remain poorly understood. In particular, how Y chromosomes evolve and degenerate, and whether dosage compensation evolves as a response, are matters of debate. Here, we focus on Coccinia grandis, a dioecious cucurbit with the highest level of X/Y heteromorphy recorded so far. We identified sex-linked genes using RNA sequences from a cross and a model-based method termed SEX-DETector. Parents and F1 individuals were genotyped, and the transmission patterns of SNPs were then analyzed. In the >1300 sex-linked genes studied, maximum X-Y divergence was 0.13-0.17, and substantial Y degeneration is implied by an average Y/X expression ratio of 0.63 and an inferred gene loss on the Y of ~40%. We also found reduced Y gene expression being compensated by elevated expression of corresponding genes on the X and an excess of sex-biased genes on the sex chromosomes. Molecular evolution of sex-linked genes in C. grandis is thus comparable to that in Silene latifolia, another dioecious plant with a strongly heteromorphic XY system, and cucurbits are the fourth plant family in which dosage compensation is described, suggesting it might be common in plants.

Magnitude of modulation of gene expression in aneuploid maize depends on the extent of genomic imbalance

Aneuploidy has profound effects on an organism, typically more so than polyploidy, and the basis of this contrast is not fully understood. A dosage series of the maize long arm of chromosome 1 (1L) was used to compare relative global gene expression in different types and degrees of aneuploidy to gain insights into how the magnitude of genomic imbalance as well as hypoploidy affects global gene expression. While previously available methods require a selective examination of specific genes, RNA sequencing provides a whole-genome view of gene expression in aneuploids. Most studies of global aneuploidy effects have concentrated on individual types of aneuploids because multiple dose aneuploidies of the same genomic region are difficult to produce in most model genetic organisms. The genetic toolkit of maize allows the examination of multiple ploidies and 1-4 doses of chromosome arms. Thus, a detailed examination of expression changes both on the varied chromosome arms and elsewhere in the genome is possible, in both hypoploids and hyperploids, compared with euploid controls. Previous studies observed the inverse trans effect, in which genes not varied in DNA dosage were expressed in a negative relationship to the varied chromosomal region. This response was also the major type of changes found globally in this study. Many genes varied in dosage showed proportional expression changes, though some were seen to be partly or fully dosage compensated. It was also found that the effects of aneuploidy were progressive, with more severe aneuploids producing effects of greater magnitude.

Aneuploidy and gene expression: is there dosage compensation?

Aneuploidy (i.e., abnormal chromosome number) is the leading cause of miscarriage and congenital defects in humans. Moreover, aneuploidy is ubiquitous in cancer. The deleterious phenotypes associated with aneuploidy are likely a result of the imbalance in the levels of gene products derived from the additional chromosome(s). Here, we summarize the current knowledge on how the presence of extra chromosomes impacts gene expression. We describe studies that have found a strict correlation between gene dosage and transcript levels as wells as studies that have found a less stringent correlation, hinting at the possible existence of dosage compensation mechanisms. We conclude by peering into the epigenetic changes found in aneuploid cells and outlining current knowledge gaps and potential areas of future investigation.

Immediate Dosage Compensation Is Triggered by the Deletion of Y-Linked Genes in Silene latifolia

The loss of functional genes from non-recombining sex-specific chromosomes [1, 2], such as the Y chromosomes in mammals [3] or W chromosomes in birds [4], should result in an imbalance of gene products for sex-linked genes [5]. Different chromosome-wide systems that rebalance gene expression are known to operate in organisms with relatively old sex chromosomes [6]; e.g., Drosophila overexpress X-linked genes in males [7], while mammals shut down one of the X chromosomes in females [8]. It is not known how long it takes for a chromosome-wide dosage compensation system to evolve. To shed light on the early evolution of dosage compensation, we constructed a high-density Y-deletion map and used deletion mutants to manipulate gene dose and analyze gene expression in white campion (Silene latifolia), which evolved dioecy and sex chromosomes only 11 million years ago [9]. We demonstrate that immediate dosage compensation can be triggered by deletions in a large portion of the p arm of the Y chromosome. Our results indicate that dosage compensation in S. latifolia does not have to evolve gene by gene because a system to upregulate gene expression is already operating on part of the X chromosome, which likely represents an intermediate step in the evolution of a chromosome-wide dosage compensation system in this species.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

Sex-chromosome dosage effects on gene expression in humans

A fundamental question in the biology of sex differences has eluded direct study in humans: How does sex-chromosome dosage (SCD) shape genome function? To address this, we developed a systematic map of SCD effects on gene function by analyzing genome-wide expression data in humans with diverse sex-chromosome aneuploidies (XO, XXX, XXY, XYY, and XXYY). For sex chromosomes, we demonstrate a pattern of obligate dosage sensitivity among evolutionarily preserved X-Y homologs and update prevailing theoretical models for SCD compensation by detecting X-linked genes that increase expression with decreasing X- and/or Y-chromosome dosage. We further show that SCD-sensitive sex-chromosome genes regulate specific coexpression networks of SCD-sensitive autosomal genes with critical cellular functions and a demonstrable potential to mediate previously documented SCD effects on disease. These gene coexpression results converge with analysis of transcription factor binding site enrichment and measures of gene expression in murine knockout models to spotlight the dosage-sensitive X-linked transcription factor ZFX as a key mediator of SCD effects on wider genome expression. Our findings characterize the effects of SCD broadly across the genome, with potential implications for human phenotypic variation.

Dosage Compensation of the X Chromosome: A Complex Epigenetic Assignment Involving Chromatin Regulators and Long Noncoding RNAs

X chromosome regulation represents a prime example of an epigenetic phenomenon where coordinated regulation of a whole chromosome is required. In flies, this is achieved by transcriptional upregulation of X chromosomal genes in males to equalize the gene dosage differences in females. Chromatin-bound proteins and long noncoding RNAs (lncRNAs) constituting a ribonucleoprotein complex known as the male-specific lethal (MSL) complex or the dosage compensation complex mediate this process. MSL complex members decorate the male X chromosome, and their absence leads to male lethality. The male X chromosome is also enriched with histone H4 lysine 16 acetylation (H4K16ac), indicating that the chromatin compaction status of the X chromosome also plays an important role in transcriptional activation. How the X chromosome is specifically targeted and how dosage compensation is mechanistically achieved are central questions for the field. Here, we review recent advances, which reveal a complex interplay among lncRNAs, the chromatin landscape, transcription, and chromosome conformation that fine-tune X chromosome gene expression.

Dosage-Dependent Expression Variation Suppressed on the Drosophila Male X Chromosome

DNA copy number variation is associated with many high phenotypic heterogeneity disorders. We systematically examined the impact of Drosophila melanogaster deletions on gene expression profiles to ask whether increased expression variability owing to reduced gene dose might underlie this phenotypic heterogeneity. Indeed, we found that one-dose genes have higher gene expression variability relative to two-dose genes. We then asked whether this increase in variability could be explained by intrinsic noise within cells due to stochastic biochemical events, or whether expression variability is due to extrinsic noise arising from more complex interactions. Our modeling showed that intrinsic gene expression noise averages at the organism level and thus cannot explain increased variation in one-dose gene expression. Interestingly, expression variability was related to the magnitude of expression compensation, suggesting that regulation, induced by gene dose reduction, is noisy. In a remarkable exception to this rule, the single X chromosome of males showed reduced expression variability, even compared with two-dose genes. Analysis of sex-transformed flies indicates that X expression variability is independent of the male differentiation program. Instead, we uncovered a correlation between occupancy of the chromatin-modifying protein encoded by males absent on the first (mof) and expression variability, linking noise suppression to the specialized X chromosome dosage compensation system. MOF occupancy on autosomes in both sexes also lowered transcriptional noise. Our results demonstrate that gene dose reduction can lead to heterogeneous responses, which are often noisy. This has implications for understanding gene network regulatory interactions and phenotypic heterogeneity. Additionally, chromatin modification appears to play a role in dampening transcriptional noise.

Global Analysis of Gene Expression in Response to Whole-Chromosome Aneuploidy in Hexaploid Wheat

Aneuploidy, a condition of unbalanced chromosome content, represents a large-effect mutation that bears significant relevance to human health and microbe adaptation. As such, extensive studies of aneuploidy have been conducted in unicellular model organisms and cancer cells. Aneuploidy also frequently is associated with plant polyploidization, but its impact on gene expression and its relevance to polyploid genome evolution/functional innovation remain largely unknown. Here, we used a panel of diverse types of whole-chromosome aneuploidy of hexaploid wheat (Triticum aestivum), all under the common genetic background of cv Chinese Spring, to systemically investigate the impact of aneuploidy on genome-, subgenome-, and chromosome-wide gene expression. Compared with prior findings in haploid or diploid aneuploid systems, we unravel additional and novel features of alteration in global gene expression resulting from the two major impacts of aneuploidy, cis- and trans-regulation, as well as dosage compensation. We show that the expression-altered genes map evenly along each chromosome, with no evidence for coregulating aggregated expression domains. However, chromosomes and subgenomes in hexaploid wheat are unequal in their responses to aneuploidy with respect to the number of genes being dysregulated. Strikingly, homeologous chromosomes do not differ from nonhomologous chromosomes in terms of aneuploidy-induced trans-acting effects, suggesting that the three constituent subgenomes of hexaploid wheat are largely uncoupled at the transcriptional level of gene regulation. Together, our findings shed new insights into the functional interplay between homeologous chromosomes and interactions between subgenomes in hexaploid wheat, which bear implications to further our understanding of allopolyploid genome evolution and efforts in breeding new allopolyploid crops.

Segmental Duplication of Chromosome 11 and its Implications for Cell Division and Genome-wide Expression in Rice

Segmental duplication is a major structural variation that occurs in chromosomes. Duplication leads to the production of gene copies with increased numbers of related repeat segments, causing the global genome to be in a state of imbalance. In addition, if the added segment contains a centromeric specific DNA, the duplicated chromosome will have structural multiple centromeres. We identified a segmental duplication containing structurally tricentric regions derived from the short arm of chromosome 11 (11L∙ + 11L∙ + 11S∙11S∙11S∙11S, “∙” represents the centromeric DNA repeat loci), and analyzed its implications for cell division and genome-wide expression. In the variant, only the middle centromere of 11S∙11S∙11S∙11S is functionally active. As a result, the structurally tricentric chromosome was stable in mitosis, because it is actually a functional monocentric chromosome. However, the structurally tricentric chromosome, which usually formed a bivalent, was either arranged on the equatorial plane or was lagging, which affected its separation during meiosis. Furthermore, RNA-seq and RT-qPCR analysis showed that the segmental duplication affected genome-wide expression patterns. 34.60% of genes in repeat region showed positive dosage effect. Thus, the genes on chromosome arm 11S-2 didn’t exhibit obviously dosage compensation, as illustrated by no peak around a ratio of 1.00. However, the gene dosage effect will reduce after sexual reproduction of a generation.

Sexual Dimorphism of Body Size Is Controlled by Dosage of the X-Chromosomal Gene Myc and by the Sex-Determining Gene tra in Drosophila

Drosophila females are larger than males. In this article, we describe how X-chromosome dosage drives sexual dimorphism of body size through two means: first, through unbalanced expression of a key X-linked growth-regulating gene, and second, through female-specific activation of the sex-determination pathway. X-chromosome dosage determines phenotypic sex by regulating the genes of the sex-determining pathway. In the presence of two sets of X-chromosome signal elements (XSEs), Sex-lethal (Sxl) is activated in female (XX) but not male (XY) animals. Sxl activates transformer (tra), a gene that encodes a splicing factor essential for female-specific development. It has previously been shown that null mutations in the tra gene result in only a partial reduction of body size of XX animals, which shows that other factors must contribute to size determination. We tested whether X dosage directly affects animal size by analyzing males with duplications of X-chromosomal segments. Upon tiling across the X chromosome, we found four duplications that increase male size by >9%. Within these, we identified several genes that promote growth as a result of duplication. Only one of these, Myc, was found not to be dosage compensated. Together, our results indicate that both Myc dosage and tra expression play crucial roles in determining sex-specific size in Drosophila larvae and adult tissue. Since Myc also acts as an XSE that contributes to tra activation in early development, a double dose of Myc in females serves at least twice in development to promote sexual size dimorphism.

Parallel Universes for Models of X Chromosome Dosage Compensation in Drosophila: A Review

Dosage compensation in Drosophila involves an approximately 2-fold increase in expression of the single X chromosome in males compared to the per gene expression in females with 2 X chromosomes. Two models have been considered for an explanation. One proposes that the male-specific lethal (MSL) complex that is associated with the male X chromosome brings histone modifiers to the sex chromosome to increase its expression. The other proposes that the inverse effect which results from genomic imbalance would tend to upregulate the genome approximately 2-fold, but the MSL complex sequesters histone modifiers from the autosomes to the X to mute this autosomal male-biased expression. On the X, the MSL complex must override the high level of resulting histone modifications to prevent overcompensation of the X chromosome. Each model is evaluated in terms of fitting classical genetic and recent molecular data. Potential paths toward resolving the models are suggested.

Genes and Small RNA Transcripts Exhibit Dosage-Dependent Expression Pattern in Maize Copy-Number Alterations

Copy-number alterations are widespread in animal and plant genomes, but their immediate impact on gene expression is still unclear. In animals, copy-number alterations usually exhibit dosage effects, except for sex chromosomes which tend to be dosage compensated. In plants, genes within small duplications (<100 kb) often exhibit dosage-dependent expression, whereas large duplications (>50 Mb) are more often dosage compensated. However, little or nothing is known about expression in moderately-sized (1-50 Mb) segmental duplications, and about the response of small RNAs to dosage change. Here, we compared maize (Zea mays) plants with two, three, and four doses of a 14.6-Mb segment of chromosome 1 that contains ∼300 genes. Plants containing the duplicated segment exhibit dosage-dependent effects on ear length and flowering time. Transcriptome analyses using GeneChip and RNA-sequencing methods indicate that most expressed genes and unique small RNAs within the duplicated segments exhibit dosage-dependent transcript levels. We conclude that dosage effect is the predominant regulatory response for both genes and unique small RNA transcripts in the segmental dosage series we tested. To our knowledge this is the first analysis of small RNA expression in plant gene dosage variants. Because segmental duplications comprise a significant proportion of eukaryotic genomes, these findings provide important new insight into the regulation of genes and small RNAs in response to dosage changes.

Dosage Compensation in the African Malaria Mosquito Anopheles gambiae

Dosage compensation is the fundamental process by which gene expression from the male monosomic X chromosome and from the diploid set of autosomes is equalized. Various molecular mechanisms have evolved in different organisms to achieve this task. In Drosophila, genes on the male X chromosome are upregulated to the levels of expression from the two X chromosomes in females. To test whether a similar mechanism is operating in immature stages of Anopheles mosquitoes, we analyzed global gene expression in the Anopheles gambiae fourth instar larvae and pupae using high-coverage RNA-seq data. In pupae of both sexes, the median expression ratios of X-linked to autosomal genes (X:A) were close to 1.0, and within the ranges of expression ratios between the autosomal pairs, consistent with complete compensation. Gene-by-gene comparisons of expression in males and females revealed mild female bias, likely attributable to a deficit of male-biased X-linked genes. In larvae, male to female ratios of the X chromosome expression levels were more female biased than in pupae, suggesting that compensation may not be complete. No compensation mechanism appears to operate in male germline of early pupae. Confirmation of the existence of dosage compensation in A. gambiae lays the foundation for research into the components of dosage compensation machinery in this important vector species.

Messenger RNA profile analysis deciphers new Esrrb responsive genes in prostate cancer cells

Background

Orphan nuclear receptor estrogen related receptor β (Esrrb or ERRβ) is well known in stem cells and early embryonic development. However, little is known about its function in cancer.

Method

We investigated the mRNA profile alterations induced by Esrrb expression and its synthetic ligand DY131 in human prostate cancer DU145 cells via RNA-Seq analysis.

Results

We distinguished 67 mRNAs differentially expressed by Esrrb alone. Although DY131 alone did not change any gene, treatment of DY131 in the presence of Esrrb altered 1161 mRNAs. These observations indicated Esrrb had both ligand-independent and ligand-dependent activity. When Esrrb was expressed, DY131 treatment further regulated 15 Esrrb-altered mRNAs. DY131 acted as an antagonist for 11 of 15 mRNAs (wdr52, f13a1, pxdn, spns2, loc100506599, tagln, loc441454, tkel1, sema3f, zcwpw2, sdc2) and as an agonist for 4 of the 15 mRNAs (rarres3, oasl, padi2, ddx60). Gene ontology analyses showed altered genes are related to transcription and translation regulation, cell proliferation and apoptosis regulation, and cellular metabolism.

Conclusion

Our results characterized mRNA profiles in DU145 prostate cancer cells driven by Esrrb expression and Esrrb ligand DY131, and provided multiple markers to characterize Esrrb’s function in Esrrb research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-015-0049-1) contains supplementary material, which is available to authorized users.

Primary Sex Determination in Drosophila melanogaster Does Not Rely on the Male-Specific Lethal Complex

It has been proposed that the Male Specific Lethal (MSL) complex is active in Drosophila melanogaster embryos of both sexes prior to the maternal-to-zygotic transition. Elevated gene expression from the two X chromosomes of female embryos is proposed to facilitate the stable establishment of Sex-lethal (Sxl) expression, which determines sex and represses further activity of the MSL complex, leaving it active only in males. Important supporting data included female-lethal genetic interactions between the seven msl genes and either Sxl or scute and sisterlessA, two of the X-signal elements (XSE) that regulate early Sxl expression. Here I report contrary findings that there are no female-lethal genetic interactions between the msl genes and Sxl or its XSE regulators. Fly stocks containing the msl3(1) allele were found to exhibit a maternal-effect interaction with Sxl, scute, and sisterlessA mutations, but genetic complementation experiments showed that msl3 is neither necessary nor sufficient for the female-lethal interactions, which appear to be due to an unidentified maternal regulator of Sxl. Published data cited as evidence for an early function of the MSL complex in females, including a maternal effect of msl2, have been reevaluated and found not to support a maternal, or other effect, of the MSL complex in sex determination. These findings suggest that the MSL complex is not involved in primary sex determination or in X chromosome dosage compensation prior to the maternal-to-zygotic transition.

From Gigabyte to Kilobyte: A Bioinformatics Protocol for Mining Large RNA-Seq Transcriptomics Data

RNA-Seq techniques generate hundreds of millions of short RNA reads using next-generation sequencing (NGS). These RNA reads can be mapped to reference genomes to investigate changes of gene expression but improved procedures for mining large RNA-Seq datasets to extract valuable biological knowledge are needed. RNAMiner--a multi-level bioinformatics protocol and pipeline--has been developed for such datasets. It includes five steps: Mapping RNA-Seq reads to a reference genome, calculating gene expression values, identifying differentially expressed genes, predicting gene functions, and constructing gene regulatory networks. To demonstrate its utility, we applied RNAMiner to datasets generated from Human, Mouse, Arabidopsis thaliana, and Drosophila melanogaster cells, and successfully identified differentially expressed genes, clustered them into cohesive functional groups, and constructed novel gene regulatory networks. The RNAMiner web service is available at http://calla.rnet.missouri.edu/rnaminer/index.html.

X chromosome inactivation and active X upregulation in therian mammals: facts, questions, and hypotheses

X chromosome inactivation is a mechanism that modulates the expression of X-linked genes in eutherian females (XX). Ohno proposed that to achieve a proper balance between X-linked and autosomal genes, those on the active X should also undergo a 2-fold upregulation. Although some support for Ohno's hypothesis has been provided through the years, recent genomic studies testing this hypothesis have brought contradictory results and fueled debate. Thus far, there are as many results in favor as against Ohno's hypothesis, depending on the nature of the datasets and the various assumptions and thresholds involved in the analyses. However, they have confirmed the importance of dosage balance between X-linked and autosomal genes involved in stoichiometric relationships. These facts as well as questions and hypotheses are discussed below.

Transcriptomic analysis identifies gene networks regulated by estrogen receptor α (ERα) and ERβ that control distinct effects of different botanical estrogens

The estrogen receptors (ERs) ERα and ERβ mediate the actions of endogenous estrogens as well as those of botanical estrogens (BEs) present in plants. BEs are ingested in the diet and also widely consumed by postmenopausal women as dietary supplements, often as a substitute for the loss of endogenous estrogens at menopause. However, their activities and efficacies, and similarities and differences in gene expression programs with respect to endogenous estrogens such as estradiol (E2) are not fully understood. Because gene expression patterns underlie and control the broad physiological effects of estrogens, we have investigated and compared the gene networks that are regulated by different BEs and by E2. Our aim was to determine if the soy and licorice BEs control similar or different gene expression programs and to compare their gene regulations with that of E2. Gene expression was examined by RNA-Seq in human breast cancer (MCF7) cells treated with control vehicle, BE or E2. These cells contained three different complements of ERs, ERα only, ERα+ERβ, or ERβ only, reflecting the different ratios of these two receptors in different human breast cancers and in different estrogen target cells. Using principal component, hierarchical clustering, and gene ontology and interactome analyses, we found that BEs regulated many of the same genes as did E2. The genes regulated by each BE, however, were somewhat different from one another, with some genes being regulated uniquely by each compound. The overlap with E2 in regulated genes was greatest for the soy isoflavones genistein and S-equol, while the greatest difference from E2 in gene expression pattern was observed for the licorice root BE liquiritigenin. The gene expression pattern of each ligand depended greatly on the cell background of ERs present. Despite similarities in gene expression pattern with E2, the BEs were generally less stimulatory of genes promoting proliferation and were more pro-apoptotic in their gene regulations than E2. The distinctive patterns of gene regulation by the individual BEs and E2 may underlie differences in the activities of these soy and licorice-derived BEs in estrogen target cells containing different levels of the two ERs.

Genetic degeneration of old and young Y chromosomes in the flowering plant Rumex hastatulus

Heteromorphic sex chromosomes have originated independently in many species, and a common feature of their evolution is the degeneration of the Y chromosome, characterized by a loss of gene content and function. Despite being of broad significance to our understanding of sex chromosome evolution, the genetic changes that occur during the early stages of Y-chromosome degeneration are poorly understood, especially in plants. Here, we investigate sex chromosome evolution in the dioecious plant Rumex hastatulus, in which X and Y chromosomes have evolved relatively recently and occur in two distinct systems: an ancestral XX/XY system and a derived XX/XY1Y2 system. This polymorphism provides a unique opportunity to investigate the effect of sex chromosome age on patterns of divergence and gene degeneration within a species. Despite recent suppression of recombination and low X-Y divergence in both systems, we find evidence that Y-linked genes have started to undergo gene loss, causing ∼ 28% and ∼ 8% hemizygosity of the ancestral and derived X chromosomes, respectively. Furthermore, genes remaining on Y chromosomes have accumulated more amino acid replacements, contain more unpreferred changes in codon use, and exhibit significantly reduced gene expression compared with their X-linked alleles, with the magnitude of these effects being greatest for older sex-linked genes. Our results provide evidence for reduced selection efficiency and ongoing Y-chromosome degeneration in a flowering plant, and indicate that Y degeneration can occur soon after recombination suppression between sex chromosomes.

Differential effect of aneuploidy on the X chromosome and genes with sex-biased expression in Drosophila

Global analysis of gene expression via RNA sequencing was conducted for trisomics for the left arm of chromosome 2 (2L) and compared with the normal genotype. The predominant response of genes on 2L was dosage compensation in that similar expression occurred in the trisomic compared with the diploid control. However, the male and female trisomic/normal expression ratio distributions for 2L genes differed in that females also showed a strong peak of genes with increased expression and males showed a peak of reduced expression relative to the opposite sex. For genes in other autosomal regions, the predominant response to trisomy was reduced expression to the inverse of the altered chromosomal dosage (2/3), but a minor peak of increased expression in females and further reduced expression in males were also found, illustrating a sexual dimorphism for the response to aneuploidy. Moreover, genes with sex-biased expression as revealed by comparing amounts in normal males and females showed responses of greater magnitude to trisomy 2L, suggesting that the genes involved in dosage-sensitive aneuploid effects also influence sex-biased expression. Each autosomal chromosome arm responded to 2L trisomy similarly, but the ratio distributions for X-linked genes were distinct in both sexes, illustrating an X chromosome-specific response to aneuploidy.