Dosage sensitivity and the evolution of gene families in yeast

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.

Genome-Wide Identification and Expression Profiling of 2OGD Superfamily Genes from Three Brassica Plants

The 2-oxoglutarate and Fe(II)-dependent dioxygenase (2OGD) superfamily is the second largest enzyme family in the plant genome, and its members are involved in various oxygenation and hydroxylation reactions. Due to their important biochemical significance in metabolism, a systematic analysis of the plant 2OGD genes family is necessary. Here, we identified 160, 179, and 337 putative 2OGDs from Brassica rapa, Brassica oleracea, and Brassica napus. According to their gene structure, domain, phylogenetic features, function, and previous studies, we also divided 676 2OGDs into three subfamilies: DOXA, DOXB, and DOXC. Additionally, homologous and phylogenetic comparisons of three subfamily genes provided valuable insight into the evolutionary characteristics of the 2OGD genes from Brassica plants. Expression profiles derived from the transcriptome and Genevestigator database exhibited distinct expression patterns of the At2OGD, Br2OGD, and Bo2OGD genes in different developmental stages, tissues, or anatomical parts. Some 2OGD genes showed high expression levels in various tissues, such as callus, seed, silique, and root tissues, while other 2OGD genes were expressed at very low levels in other tissues. Analysis of six Bo2OGD genes in different tissues by qRT-PCR indicated that these genes are involved in the metabolism of gibberellin, which in turn regulates plant growth and development. Our working system analysed 2OGD gene families of three Brassica plants and laid the foundation for further study of their functional characterization.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Multilayered regulation of proteome stoichiometry

Cellular systems depend on multiprotein complexes whose functionalities require defined stoichiometries of subunit proteins. Proper stoichiometry is achieved by controlling the amount of protein synthesis and degradation even in the presence of genetic perturbations caused by changes in gene dosage. As a consequence of increased gene copy number, excess subunits unassembled into the complex are synthesized and rapidly degraded by the ubiquitin–proteasome system. This mechanism, called protein-level dosage compensation, is widely observed not only under such perturbed conditions but also in unperturbed physiological cells. Recent studies have shown that recognition of unassembled subunits and their selective degradation are intricately regulated. This review summarizes the nature, strategies, and increasing complexity of protein-level dosage compensation and discusses possible mechanisms for controlling proteome stoichiometry in multiple layers of biological processes.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

Genes Encoding Teleost Orthologs of Human Haploinsufficient and Monoallelically Expressed Genes Remain in Duplicate More Frequently Than the Whole Genome

Gene dosage is an important issue both in cell and evolutionary biology. Most genes are present in two copies or alleles in diploid eukariotic cells. The most outstanding exception is monoallelic gene expression (MA) that concerns genes localized on the X chromosome or in regions undergoing parental imprinting in eutherians, and many other genes scattered throughout the genome. In diploids, haploinsufficiency (HI) implies that a single functional copy of a gene in a diploid organism is insufficient to ensure a normal biological function. One of the most important mechanisms ensuring functional innovation during evolution is whole genome duplication (WGD). In addition to the two WGDs that have occurred in vertebrate genomes, the teleost genomes underwent an additional WGD, after their divergence from tetrapods. In the present work, we have studied on 57 teleost species whether the orthologs of human MA or HI genes remain more frequently in duplicates or returned more frequently in singleton than the rest of the genome. Our results show that the teleost orthologs of HI human genes remained more frequently in duplicate than the rest of the genome in all of the teleost species studied. No signal was observed for the orthologs of genes mapping to the human X chromosome or subjected to parental imprinting. Surprisingly, the teleost orthologs of the other human MA genes remained in duplicate more frequently than the rest of the genome for most teleost species. These results suggest that the teleost orthologs of MA and HI human genes also undergo selective pressures either related to absolute protein amounts and/or of dosage balance issues. However, these constraints seem to be different for MA genes in teleost in comparison with human genomes.

Revisiting the Neurospora crassa mitochondrial genome

The mitochondrial genome of Neurospora crassa has been less studied than its nuclear counterpart, yet it holds great potential for understanding the diversity and evolution of this important fungus. Here we describe a new mitochondrial DNA (mtDNA) complete sequence of a N. crassa wild type strain. The genome with 64 839 bp revealed 21 protein-coding genes and several hypothetical open reading frames with no significant homology to any described gene. Five large repetitive regions were identified across the genome, including partial or complete genes. The largest repeated region holds a partial nd2 section that was also detected in Neurospora intermedia, suggesting a rearrangement that occurred before the N. crassa speciation. Interestingly, N. crassa has a palindrome adjacent to the partial nd2 repeated region possibly related to the genomic rearrangement, which is absent in N. intermedia. Finally, we compared the sequences of the three available N. crassa complete mtDNAs and found low levels of intraspecific variability. Most differences among strains were due to small indels in noncoding regions. The revisiting of the N. crassa mtDNA forms the basis for future studies on mitochondrial genome organization and variability.

Expression attenuation as a mechanism of robustness against gene duplication

Many studies have focused on the mechanisms of long-term retention of gene duplicates, such as the gain of functions or reciprocal losses. However, such changes are more likely to occur if the duplicates are maintained for a long period. This time span will be short if duplication is immediately deleterious. We measured the distribution of fitness effects of gene duplication for 899 genes in budding yeast. We find that gene duplication is more likely to be deleterious than beneficial. However, contrary to previous models, in general, gene duplication does not affect fitness by altering the organization of protein complexes. We show that expression attenuation may protect complexes from the effects of gene duplication.

Gene duplication is ubiquitous and a major driver of phenotypic diversity across the tree of life, but its immediate consequences are not fully understood. Deleterious effects would decrease the probability of retention of duplicates and prevent their contribution to long-term evolution. One possible detrimental effect of duplication is the perturbation of the stoichiometry of protein complexes. Here, we measured the fitness effects of the duplication of 899 essential genes in the budding yeast using high-resolution competition assays. At least 10% of genes caused a fitness disadvantage when duplicated. Intriguingly, the duplication of most protein complex subunits had small to nondetectable effects on fitness, with few exceptions. We selected four complexes with subunits that had an impact on fitness when duplicated and measured the impact of individual gene duplications on their protein–protein interactions. We found that very few duplications affect both fitness and interactions. Furthermore, large complexes such as the 26S proteasome are protected from gene duplication by attenuation of protein abundance. Regulatory mechanisms that maintain the stoichiometric balance of protein complexes may protect from the immediate effects of gene duplication. Our results show that a better understanding of protein regulation and assembly in complexes is required for the refinement of current models of gene duplication.

A genome-scale yeast library with inducible expression of individual genes

The ability to switch a gene from off to on and monitor dynamic changes provides a powerful approach for probing gene function and elucidating causal regulatory relationships. Here, we developed and characterized YETI (Yeast Estradiol strains with Titratable Induction), a collection in which > 5,600 yeast genes are engineered for transcriptional inducibility with single-gene precision at their native loci and without plasmids. Each strain contains SGA screening markers and a unique barcode, enabling high-throughput genetics. We characterized YETI using growth phenotyping and BAR-seq screens, and we used a YETI allele to identify the regulon of Rof1, showing that it acts to repress transcription. We observed that strains with inducible essential genes that have low native expression can often grow without inducer. Analysis of data from eukaryotic and prokaryotic systems shows that native expression is a variable that can bias promoter-perturbing screens, including CRISPRi. We engineered a second expression system, Z3 EB42, that gives lower expression than Z3 EV, a feature enabling conditional activation and repression of lowly expressed essential genes that grow without inducer in the YETI library.

HAG1 and SWI3A/B control of male germ line development in P. patens suggests conservation of epigenetic reproductive control across land plants

KEY MESSAGE

Bryophytes as models to study the male germ line: loss-of-function mutants of epigenetic regulators HAG1 and SWI3a/b demonstrate conserved function in sexual reproduction. With the water-to-land transition, land plants evolved a peculiar haplodiplontic life cycle in which both the haploid gametophyte and the diploid sporophyte are multicellular. The switch between these phases was coined alternation of generations. Several key regulators that control the bauplan of either generation are already known. Analyses of such regulators in flowering plants are difficult due to the highly reduced gametophytic generation, and the fact that loss of function of such genes often is embryo lethal in homozygous plants. Here we set out to determine gene function and conservation via studies in bryophytes. Bryophytes are sister to vascular plants and hence allow evolutionary inferences. Moreover, embryo lethal mutants can be grown and vegetatively propagated due to the dominance of the bryophyte gametophytic generation. We determined candidates by selecting single copy orthologs that are involved in transcriptional control, and of which flowering plant mutants show defects during sexual reproduction, with a focus on the under-studied male germ line. We selected two orthologs, SWI3a/b and HAG1, and analyzed loss-of-function mutants in the moss P. patens. In both mutants, due to lack of fertile spermatozoids, fertilization and hence the switch to the diploid generation do not occur. Pphag1 additionally shows arrested male and impaired female gametangia development. We analyzed HAG1 in the dioecious liverwort M. polymorpha and found that in Mphag1 the development of gametangiophores is impaired. Taken together, we find that involvement of both regulators in sexual reproduction is conserved since the earliest divergence of land plants.

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.

Paralogs in the PKA Regulon Traveled Different Evolutionary Routes to Divergent Expression in Budding Yeast

Functional divergence of duplicate genes, or paralogs, is an important driver of novelty in evolution. In the model yeast Saccharomyces cerevisiae, there are 547 paralog gene pairs that survive from an interspecies Whole Genome Hybridization (WGH) that occurred ~100MYA. In this work, we report that ~1/6th (110) of these WGH paralogs pairs (or ohnologs) are differentially expressed with a striking pattern upon Protein Kinase A (PKA) inhibition. One member of each pair in this group has low basal expression that increases upon PKA inhibition, while the other has moderate and unchanging expression. For these genes, expression of orthologs upon PKA inhibition in the non-WGH species Kluyveromyces lactis and for PKA-related stresses in other budding yeasts shows unchanging expression, suggesting that lack of responsiveness to PKA was likely the typical ancestral phenotype prior to duplication. Promoter sequence analysis across related budding yeast species further revealed that the subsequent emergence of PKA-dependence took different evolutionary routes. In some examples, regulation by PKA and differential expression appears to have arisen following the WGH, while in others, regulation by PKA appears to have arisen in one of the two parental lineages prior to the WGH. More broadly, our results illustrate the unique opportunities presented by a WGH event for generating functional divergence by bringing together two parental lineages with separately evolved regulation into one species. We propose that functional divergence of two ohnologs can be facilitated through such regulatory divergence.

Shared Neurodevelopmental Perturbations Can Lead to Intellectual Disability in Individuals with Distinct Rare Chromosome Duplications

Chromosomal duplications are associated with a large group of human diseases that arise mainly from dosage imbalance of genes within the rearrangements. Phenotypes range widely but are often associated with global development delay, intellectual disability, autism spectrum disorders, and multiple congenital abnormalities. How different contiguous genes from a duplicated genomic region interact and dynamically affect the expression of each other remains unclear in most cases. Here, we report a genomic comparative delineation of genes located in duplicated chromosomal regions 8q24.13q24.3, 18p11.32p11.21, and Xq22.3q27.2 in three patients followed up at our genetics service who has the intellectual disability (ID) as a common phenotype. We integrated several genomic data levels by identification of gene content within the duplications, protein-protein interactions, and functional analysis on specific tissues. We found functional relationships among genes from three different duplicated chromosomal regions, reflecting interactions of protein-coding genes and their involvement in common cellular subnetworks. Furthermore, the sharing of common significant biological processes associated with ID has been demonstrated between proteins from the different chromosomal regions. Finally, we elaborated a shared model of pathways directly or indirectly related to the central nervous system (CNS), which could perturb cognitive function and lead to ID in the three duplication conditions.

Functional and evolutionary analysis of the Arabidopsis 4R-MYB protein SNAPc4 as part of the SNAP complex

Transcription initiation of the genes coding for small nuclear RNA (snRNA) has been extensively analyzed in humans and fruit fly, but only a single ortholog of a snRNA-activating protein complex (SNAPc) subunit has so far been characterized in plants. The genome of the model plant Arabidopsis thaliana encodes orthologs of all three core SNAPc subunits, including A. thaliana SNAP complex 4 (AtSNAPc4)-a 4R-MYB-type protein with four-and-a-half adjacent MYB repeat units. We report the conserved role of AtSNAPc4 as subunit of a protein complex involved in snRNA gene transcription and present genetic evidence that AtSNAPc4 is an essential gene in gametophyte and zygote development. We present experimental evidence that the three A. thaliana SNAPc subunits assemble into a SNAP complex and demonstrate the binding of AtSNAPc4 to snRNA promoters. In addition, co-localization studies show a link between AtSNAPc4 accumulation and Cajal bodies, known to aggregate at snRNA gene loci in humans. Moreover, we show the strong evolutionary conservation of single-copy 4R-MYB/SNAPc4 genes in a broad range of eukaryotes and present additional shared protein features besides the MYB domain, suggesting a conservation of the snRNA transcription initiation machinery along the course of the eukaryotic evolution.

Long QT Syndrome KCNH2 Variant Induces hERG1a/1b Subunit Imbalance in Patient-Specific Induced Pluripotent Stem Cell-Derived Cardiomyocytes

[Figure: see text].

Reduced Translational Efficiency of Eukaryotic Genes after Duplication Events

Control of gene expression has been found to be predominantly determined at the level of protein translation. However, to date, reduced expression from duplicated genes in eukaryotes for dosage maintenance has only been linked to transcriptional control involving epigenetic mechanisms. Here, we hypothesize that dosage maintenance following gene duplication also involves regulation at the protein level. To test this hypothesis, we compared transcriptome and proteome data of yeast models, Saccharomyces cerevisiae and Schizosaccharomyces pombe, and worm models, Caenorhabditis elegans and Caenorhabditis briggsae, to investigate lineage-specifically duplicated genes. Duplicated genes in both eukaryotic models exhibited a reduced protein-to-mRNA abundance ratio. Moreover, dosage sensitive genes, represented by genes encoding protein complex subunits, reduced their protein-to-mRNA abundance ratios more significantly than the other genes after duplication events. An analysis of ribosome profiling (Ribo-Seq) data further showed that reduced translational efficiency was more prominent for dosage sensitive genes than for the other genes. Meanwhile, no difference in protein degradation rate was associated with duplication events. Translationally repressed duplicated genes were also more likely to be inhibited at the level of transcription. Taken together, these results suggest that translation-mediated dosage control is partially contributed by natural selection and it enhances transcriptional control in maintaining gene dosage after gene duplication events during eukaryotic genome evolution.

Constructive Neutral Evolution 20 Years Later

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

Dosage sensitivity of JDPs, a valuable tool for understanding their function: a case study on Caj1 overexpression-mediated filamentous growth in budding yeast

J-domain proteins (JDPs) partner with Hsp70s to oversee proper synthesis, folding, transport and turnover of proteins in the cell. In any subcellular compartment, often multiple JDPs collaborate with a single Hsp70 to perform a variety of functions. Being co-localized, JDPs may exhibit complex genetic and physical interactions with each other, their clients as well as the Hsp70 partners. Even though most JDPs are highly specialized, redundancy between them is possible, making their functional analysis challenging. In the absence of assayable deletion phenotypes, protein overexpression appears to be a powerful alternative strategy to study JDP function. Here, we show that high levels of Caj1, one of the cytosolic JDPs, cause filamentous growth and G2/M arrest in yeast cells. Mutation in the critical HPD motif in the J-domain of Caj1 completely abolished these phenotypes, suggesting that Hsp70 co-chaperone function is important for the dominant-negative phenotypes exhibited by Caj1 overexpression. In this paper, we discuss the possible underlying mechanisms responsible for the pleiotropic phenotypes displayed by Caj1 overexpression in the light of current models proposed for dosage-sensitive genes (DSGs). Finally, we present generalized mechanisms of JDP overexpression-mediated dominant-negative phenotypes in budding yeast.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.

RIPK1 Promotes Energy Sensing by the mTORC1 Pathway

The mechanisms of cellular energy sensing and AMPK-mediated mTORC1 inhibition are not fully delineated. Here, we discover that RIPK1 promotes mTORC1 inhibition during energetic stress. RIPK1 is involved in mediating the interaction between AMPK and TSC2 and facilitate TSC2 phosphorylation at Ser1387. RIPK1 loss results in a high basal mTORC1 activity that drives defective lysosomes in cells and mice, leading to accumulation of RIPK3 and CASP8 and sensitization to cell death. RIPK1-deficient cells are unable to cope with energetic stress and are vulnerable to low glucose levels and metformin. Inhibition of mTORC1 rescues the lysosomal defects and vulnerability to energetic stress and prolongs the survival of RIPK1-deficient neonatal mice. Thus, RIPK1 plays an important role in the cellular response to low energy levels and mediates AMPK-mTORC1 signaling. These findings shed light on the regulation of mTORC1 during energetic stress and unveil a point of crosstalk between pro-survival and pro-death pathways.

The evolution of sex chromosome dosage compensation in animals

The evolution of heteromorphic sex chromosomes shall lead to gene expression dosage problems, as in at least one of the sexes, the sex-linked gene dose has been reduced by half. It has been proposed that the transcriptional output of the whole X or Z chromosome should be doubled for complete dosage compensation in heterogametic sex. However, owing to the variability of the existing methods to determine the transcriptional differences between sex chromosomes and autosomes (S:A ratios) in different studies, we collected more than 500 public RNA-Seq data set from multiple tissues and species in major clades and proposed a unified computational framework for unbiased and comparable measurement of the S:A ratios of multiple species. We also tested the evolution of dosage compensation more directly by assessing changes in the expression levels of the current sex-linked genes relative to those of the ancestral sex-linked genes. We found that in mammals and birds, the S:A ratio is approximately 0.5, whereas in insects, fishes, and flatworms, the S:A ratio is approximately 1.0. Further analysis showed that the fraction of dosage-sensitive housekeeping genes on the X/Z chromosome is significantly correlated with the S:A ratio. In addition, the degree of degeneration of the Y chromosome may be responsible for the change in the S:A ratio in mammals without a dosage compensation mechanism. Our observations offer unequivocal support for the sex chromosome insensitivity hypothesis in animals and suggest that dosage sensitivity states of sex chromosomes are a major factor underlying different evolutionary strategies of dosage compensation.

Exploring the complexity of soybean (Glycine max) transcriptional regulation using global gene co-expression networks

MAIN CONCLUSION

We report a soybean gene co-expression network built with data from 1284 RNA-Seq experiments, which was used to identify important regulators, modules and to elucidate the fates of gene duplicates. Soybean (Glycine max (L.) Merr.) is one of the most important crops worldwide, constituting a major source of protein and edible oil. Gene co-expression networks (GCN) have been extensively used to study transcriptional regulation and evolution of genes and genomes. Here, we report a soybean GCN using 1284 publicly available RNA-Seq samples from 15 distinct tissues. We found modules that are differentially regulated in specific tissues, comprising processes such as photosynthesis, gluconeogenesis, lignin metabolism, and response to biotic stress. We identified transcription factors among intramodular hubs, which probably integrate different pathways and shape the transcriptional landscape in different conditions. The top hubs for each module tend to encode proteins with critical roles, such as succinate dehydrogenase and RNA polymerase subunits. Importantly, gene essentiality was strongly correlated with degree centrality and essential hubs were enriched in genes involved in nucleic acids metabolism and regulation of cell replication. Using a guilt-by-association approach, we predicted functions for 93 of 106 hubs without functional description in soybean. Most of the duplicated genes had different transcriptional profiles, supporting their functional divergence, although paralogs originating from whole-genome duplications (WGD) are more often preserved in the same module than those from other mechanisms. Together, our results highlight the importance of GCN analysis in unraveling key functional aspects of the soybean genome, in particular those associated with hub genes and WGD events.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Nuclear gene proximity and protein interactions shape transcript covariations in mammalian single cells

Single-cell RNA sequencing studies on gene co-expression patterns could yield important regulatory and functional insights, but have so far been limited by the confounding effects of differentiation and cell cycle. We apply a tailored experimental design that eliminates these confounders, and report thousands of intrinsically covarying gene pairs in mouse embryonic stem cells. These covariations form a network with biological properties, outlining known and novel gene interactions. We provide the first evidence that miRNAs naturally induce transcriptome-wide covariations and compare the relative importance of nuclear organization, transcriptional and post-transcriptional regulation in defining covariations. We find that nuclear organization has the greatest impact, and that genes encoding for physically interacting proteins specifically tend to covary, suggesting importance for protein complex formation. Our results lend support to the concept of post-transcriptional RNA operons, but we further present evidence that nuclear proximity of genes may provide substantial functional regulation in mammalian single cells.

Gene expression covariation can be studied by single-cell RNA sequencing. Here the authors analyze intrinsically covarying gene pairs by eliminating the confounding effects in single-cell experiments and observe covariation of proximal genes and miRNA-induced covariation of target mRNAs.

Exploring the Complexity of Protein-Level Dosage Compensation that Fine-Tunes Stoichiometry of Multiprotein Complexes

Proper control of gene expression levels upon various perturbations is a fundamental aspect of cellular robustness. Protein-level dosage compensation is one mechanism buffering perturbations to stoichiometry of multiprotein complexes through accelerated proteolysis of unassembled subunits. Although N-terminal acetylation- and ubiquitin-mediated proteasomal degradation by the Ac/N-end rule pathway enables selective compensation of excess subunits, it is unclear how widespread this pathway contributes to stoichiometry control. Here we report that dosage compensation depends only partially on the Ac/N-end rule pathway. Our analysis of genetic interactions between 18 subunits and 12 quality control factors in budding yeast demonstrated that multiple E3 ubiquitin ligases and N-acetyltransferases are involved in dosage compensation. We find that N-acetyltransferases-mediated compensation is not simply predictable from N-terminal sequence despite their sequence specificity for N-acetylation. We also find that the compensation of Pop3 and Bet4 is due in large part to a minor N-acetyltransferase NatD. Furthermore, canonical NatD substrates histone H2A/H4 were compensated even in its absence, suggesting N-acetylation-independent stoichiometry control. Our study reveals the complexity and robustness of the stoichiometry control system.

Quality control of multiprotein complexes is important for maintaining homeostasis in cellular systems that are based on functional complexes. Proper stoichiometry of multiprotein complexes is achieved by the balance between protein synthesis and degradation. Recent studies showed that translation efficiency tends to scale with stoichiometry of their subunits. On the other hand, although protein N-terminal acetylation- and ubiquitin-mediated proteolysis pathway is involved in selective degradation of excess subunits, it is unclear how widespread this pathway contributes to stoichiometry control due to the lack of a systematic investigation using endogenous proteins. To better understand the landscape of the stoichiometry control system, we examined genetic interactions between 18 subunits and 12 quality control factors (E3 ubiquitin ligases and N-acetyltransferases), in total 114 combinations. Our data suggest that N-acetyltransferases are partially responsible for stoichiometry control and that N-acetylation-independent pathway is also involved in selective degradation of excess subunits. Therefore, this study reveals the complexity and robustness of the stoichiometry control system. Further dissection of this complexity will help to understand the mechanisms buffering gene expression perturbations and shaping proteome stoichiometry.

Experimental determination of evolutionary barriers to horizontal gene transfer

BACKGROUND

Horizontal gene transfer, the acquisition of genes across species boundaries, is a major source of novel phenotypes that enables microbes to rapidly adapt to new environments. How the transferred gene alters the growth - fitness - of the new host affects the success of the horizontal gene transfer event and how rapidly the gene spreads in the population. Several selective barriers - factors that impact the fitness effect of the transferred gene - have been suggested to impede the likelihood of horizontal transmission, however experimental evidence is scarce. The objective of this study was to determine the fitness effects of orthologous genes transferred from Salmonella enterica serovar Typhimurium to Escherichia coli to identify the selective barriers using highly precise experimental measurements.

RESULTS

We found that most gene transfers result in strong fitness costs. Previously identified evolutionary barriers - gene function and the number of protein-protein interactions - did not predict the fitness effects of transferred genes. In contrast, dosage sensitivity, gene length, and the intrinsic protein disorder significantly impact the likelihood of a successful horizontal transfer.

CONCLUSION

While computational approaches have been successful in describing long-term barriers to horizontal gene transfer, our experimental results identified previously underappreciated barriers that determine the fitness effects of newly transferred genes, and hence their short-term eco-evolutionary dynamics.

Parallel Concerted Evolution of Ribosomal Protein Genes in Fungi and Its Adaptive Significance

Ribosomal protein (RP) genes encode structural components of ribosomes, the cellular machinery for protein synthesis. A single functional copy has been maintained in most of 78–80 RP families in animals due to evolutionary constraints imposed by gene dosage balance. Some fungal species have maintained duplicate copies in most RP families. The mechanisms by which the RP genes were duplicated and maintained and their functional significance are poorly understood. To address these questions, we identified all RP genes from 295 fungi and inferred the timing and nature of gene duplication events for all RP families. We found that massive duplications of RP genes have independently occurred by different mechanisms in three distantly related lineages: budding yeasts, fission yeasts, and Mucoromycota. The RP gene duplicates in budding yeasts and Mucoromycota were mainly created by whole genome duplication events. However, duplicate RP genes in fission yeasts were likely generated by retroposition, which is unexpected considering their dosage sensitivity. The sequences of most RP paralogs have been homogenized by repeated gene conversion in each species, demonstrating parallel concerted evolution, which might have facilitated the retention of their duplicates. Transcriptomic data suggest that the duplication and retention of RP genes increased their transcript abundance. Physiological data indicate that increased ribosome biogenesis allowed these organisms to rapidly consume sugars through fermentation while maintaining high growth rates, providing selective advantages to these species in sugar-rich environments.

Pooled analysis of radiation hybrids identifies loci for growth and drug action in mammalian cells

Genetic screens in mammalian cells commonly focus on loss-of-function approaches. To evaluate the phenotypic consequences of extra gene copies, we used bulk segregant analysis (BSA) of radiation hybrid (RH) cells. We constructed six pools of RH cells, each consisting of ∼2500 independent clones, and placed the pools under selection in media with or without paclitaxel. Low pass sequencing identified 859 growth loci, 38 paclitaxel loci, 62 interaction loci, and three loci for mitochondrial abundance at genome-wide significance. Resolution was measured as ∼30 kb, close to single-gene. Divergent properties were displayed by the RH-BSA growth genes compared to those from loss-of-function screens, refuting the balance hypothesis. In addition, enhanced retention of human centromeres in the RH pools suggests a new approach to functional dissection of these chromosomal elements. Pooled analysis of RH cells showed high power and resolution and should be a useful addition to the mammalian genetic toolkit.

Dosage-sensitive genes in autism spectrum disorders: From neurobiology to therapy

Autism spectrum disorders (ASDs) are a group of heterogenous neurodevelopmental disorders affecting 1 in 59 children. Syndromic ASDs are commonly associated with chromosomal rearrangements or dosage imbalance involving a single gene. Many of these genes are dosage-sensitive and regulate transcription, protein homeostasis, and synaptic function in the brain. Despite vastly different molecular perturbations, syndromic ASDs share core symptoms including social dysfunction and repetitive behavior. However, each ASD subtype has a unique pathogenic mechanism and combination of comorbidities that require individual attention. We have learned a great deal about how these dosage-sensitive genes control brain development and behaviors from genetically-engineered mice. Here we describe the clinical features of eight monogenic neurodevelopmental disorders caused by dosage imbalance of four genes, as well as recent advances in using genetic mouse models to understand their pathogenic mechanisms and develop intervention strategies. We propose that applying newly developed quantitative molecular and neuroscience technologies will advance our understanding of the unique neurobiology of each disorder and enable the development of personalized therapy.

Transcriptional Rewiring, Adaptation, and the Role of Gene Duplication in the Metabolism of Ethanol of Saccharomyces cerevisiae

Ethanol is the main by-product of yeast sugar fermentation that affects microbial growth parameters, being considered a dual molecule, a nutrient and a stressor. Previous works demonstrated that the budding yeast arose after an ancient hybridization process resulted in a tier of duplicated genes within its genome, many of them with implications in this ethanol "produce-accumulate-consume" strategy. The evolutionary link between ethanol production, consumption, and tolerance versus ploidy and stability of the hybrids is an ongoing debatable issue. The implication of ancestral duplicates in this metabolic rewiring, and how these duplicates differ transcriptionally, remains unsolved. Here, we study the transcriptomic adaptive signatures to ethanol as a nonfermentative carbon source to sustain clonal yeast growth by experimental evolution, emphasizing the role of duplicated genes in the adaptive process. As expected, ethanol was able to sustain growth but at a lower rate than glucose. Our results demonstrate that in asexual populations a complete transcriptomic rewiring was produced, strikingly by downregulation of duplicated genes, mainly whole-genome duplicates, whereas small-scale duplicates exhibited significant transcriptional divergence between copies. Overall, this study contributes to the understanding of evolution after gene duplication, linking transcriptional divergence with duplicates' fate in a multigene trait as ethanol tolerance.IMPORTANCE Gene duplication events have been related with increasing biological complexity through the tree of life, but also with illnesses, including cancer. Early evolutionary theories indicated that duplicated genes could explore alternative functions due to relaxation of selective constraints in one of the copies, as the other remains as ancestral-function backup. In unicellular eukaryotes like yeasts, it has been demonstrated that the fate and persistence of duplicates depend on duplication mechanism (whole-genome or small-scale events), shaping their actual genomes. Although it has been shown that small-scale duplicates tend to innovate and whole-genome duplicates specialize in ancestral functions, the implication of duplicates' transcriptional plasticity and transcriptional divergence on environmental and metabolic responses remains largely obscure. Here, by experimental adaptive evolution, we show that Saccharomyces cerevisiae is able to respond to metabolic stress (ethanol as nonfermentative carbon source) due to the persistence of duplicated genes. These duplicates respond by transcriptional rewiring, depending on their transcriptional background. Our results shed light on the mechanisms that determine the role of duplicates, and on their evolvability.

Exploring whole-genome duplicate gene retention with complex genetic interaction analysis

Whole-genome duplication has played a central role in the genome evolution of many organisms, including the human genome. Most duplicated genes are eliminated, and factors that influence the retention of persisting duplicates remain poorly understood. We describe a systematic complex genetic interaction analysis with yeast paralogs derived from the whole-genome duplication event. Mapping of digenic interactions for a deletion mutant of each paralog, and of trigenic interactions for the double mutant, provides insight into their roles and a quantitative measure of their functional redundancy. Trigenic interaction analysis distinguishes two classes of paralogs: a more functionally divergent subset and another that retained more functional overlap. Gene feature analysis and modeling suggest that evolutionary trajectories of duplicated genes are dictated by combined functional and structural entanglement factors.

A Physiological Basis for Nonheritable Antibiotic Resistance

Antibiotics constitute one of the cornerstones of modern medicine. However, individuals may succumb to a bacterial infection if a pathogen survives exposure to antibiotics. The ability of bacteria to survive bactericidal antibiotics results from genetic changes in the preexisting bacterial genome, from the acquisition of genes from other organisms, and from nonheritable phenomena that give rise to antibiotic tolerance. Nonheritable antibiotic tolerance can be exhibited by a large fraction of the bacterial population or by a small subpopulation referred to as persisters. Nonheritable resistance to antibiotics has been ascribed to the activity of toxins that are part of toxin-antitoxin modules, to the universal energy currency ATP, and to the signaling molecule guanosine (penta) tetraphosphate. However, these molecules are dispensable for nonheritable resistance to antibiotics in many organisms. By contrast, nutrient limitation, treatment with bacteriostatic antibiotics, or expression of genes that slow bacterial growth invariably promote nonheritable resistance. We posit that antibiotic persistence results from conditions promoting feedback inhibition among core cellular processes, resulting phenotypically in a slowdown or halt in bacterial growth.

Autophagy-mediated apoptosis eliminates aneuploid cells in a mouse model of chromosome mosaicism

The high incidence of aneuploidy in the embryo is considered the principal cause for low human fecundity. However, the prevalence of aneuploidy dramatically declines as pregnancy progresses, with the steepest drop occurring as the embryo completes implantation. Despite the fact that the plasticity of the embryo in dealing with aneuploidy is fundamental to normal development, the mechanisms responsible for eliminating aneuploid cells are unclear. Here, using a mouse model of chromosome mosaicism, we show that aneuploid cells are preferentially eliminated from the embryonic lineage in a p53-dependent process involving both autophagy and apoptosis before, during and after implantation. Moreover, we show that diploid cells in mosaic embryos undertake compensatory proliferation during the implantation stages to confer embryonic viability. Together, our results indicate a close link between aneuploidy, autophagy, and apoptosis to refine the embryonic cell population and ensure only chromosomally fit cells proceed through development of the fetus.

The mechanisms behind the plasticity of embryos and how they deal with aneuploid cells are unclear. Here, the authors show that aneuploid cells in a mouse embryo are preferentially eliminated during pre- and peri-implantation development in a p53-dependent process involving both autophagy and apoptosis.

Supplementation of Aspergillus glaucus with gfdB gene encoding a glycerol 3-phosphate dehydrogenase in Aspergillus nidulans

In Aspergillus nidulans, there are two putative glycerol 3-phosphate dehydrogenases encoded by the genes gfdA and gfdB, while the genome of the osmophilic Aspergillus glaucus harbors only the ortholog of the A. nidulans gfdA gene. Our aim was to insert the gfdB gene into the genome of A. glaucus, and we reached this goal with the adaptation of the Agrobacterium tumefaciens-mediated transformation method. We tested the growth of the gfdB-complemented A. glaucus strains on a medium containing 2 mol l^-1 sorbitol in the presence of oxidative stress generating agents such as tert-butyl hydroperoxide, H₂ O₂ , menadione sodium bisulfite, as well as the cell wall integrity stress-inducing agent Congo Red and the heavy metal stress eliciting CdCl₂ . The growth of the complemented strains was significantly higher than that of the wild-type strain on media supplemented with these stress generating agents. The A. nidulans ΔgfdB mutant was also examined under the same conditions and resulted in a considerably lower growth than that of the control strain in all stress exposure experiments. Our results shed light on the fact that the gfdB gene from A. nidulans was also involved in the stress responses of the complemented A. glaucus strains supporting our hypothesis on the antioxidant function of GfdB in the Aspergilli. Nevertheless, the osmotolerant nature of A. glaucus could not be explained by the lack of the gfdB gene in A. glaucus, as we hypothesized earlier.

Epigenetic and transcriptomic consequences of excess X-chromosome material in 47,XXX syndrome-A comparison with Turner syndrome and 46,XX females

47,XXX (triple X) and Turner syndrome (45,X) are sex chromosomal abnormalities with detrimental effects on health with increased mortality and morbidity. In karyotypical normal females, X-chromosome inactivation balances gene expression between sexes and upregulation of the X chromosome in both sexes maintain stoichiometry with the autosomes. In 47,XXX and Turner syndrome a gene dosage imbalance may ensue from increased or decreased expression from the genes that escape X inactivation, as well as from incomplete X chromosome inactivation in 47,XXX. We aim to study genome-wide DNA-methylation and RNA-expression changes can explain phenotypic traits in 47,XXX syndrome. We compare DNA-methylation and RNA-expression data derived from white blood cells of seven women with 47,XXX syndrome, with data from seven female controls, as well as with seven women with Turner syndrome (45,X). To address these questions, we explored genome-wide DNA-methylation and transcriptome data in blood from seven females with 47,XXX syndrome, seven females with Turner syndrome, and seven karyotypically normal females (46,XX). Based on promoter methylation, we describe a demethylation of six X-chromosomal genes (AMOT, HTR2C, IL1RAPL2, STAG2, TCEANC, ZNF673), increased methylation for GEMIN8, and four differentially methylated autosomal regions related to four genes (SPEG, MUC4, SP6, and ZNF492). We illustrate how these changes seem compensated at the transcriptome level although several genes show differential exon usage. In conclusion, our results suggest an impact of the supernumerary X chromosome in 47,XXX syndrome on the methylation status of selected genes despite an overall comparable expression profile.

Evolutionary strategies drive a balance of the interacting gene products for the CBL and CIPK gene families

Genes encoding interacting proteins tend to be co-retained after whole-genome duplication (WGD). The preferential retention after WGD has been explained by the gene balance hypothesis (GBH). However, small-scale duplications could independently occur in the connected gene families. Certain evolutionary strategies might keep the dosage balanced. Here, we examined the gene duplication, interaction and expression patterns of calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) gene families to understand the underlying principles. The ratio of the CBL and CIPK gene numbers evolved from 5 : 7 in Physcomitrella to 10 : 26 in Arabidopsis, and retrotransposition, tandem duplication, and WGDs contributed to the expansion. Two pairs of CBLs and six pairs of CIPKs were retained after the α WGD in Arabidopsis, in which specific interaction patterns were identified. In some cases, two retained CBLs (CIPKs) might compete to interact with a sole CIPK (CBL). Results of gene expression analyses indicated that the relatively over-retained duplicates tend to show asymmetric expression, thus avoiding competition. In conclusion, our results suggested that the highly specific interaction, together with the differential gene expression pattern, jointly maintained the balanced dosage for the interacting CBL and CIPK proteins.

Aligning functional network constraint to evolutionary outcomes

BACKGROUND

Functional constraint through genomic architecture is suggested to be an important dimension of genome evolution, but quantitative evidence for this idea is rare. In this contribution, existing evidence and discussions on genomic architecture as constraint for convergent evolution, rapid adaptation, and genic adaptation are summarized into alternative, testable hypotheses. Network architecture statistics from protein-protein interaction networks are then used to calculate differences in evolutionary outcomes on the example of genomic evolution in yeast, and the results are used to evaluate statistical support for these longstanding hypotheses.

RESULTS

A discriminant function analysis lent statistical support to classifying the yeast interactome into hub, intermediate and peripheral nodes based on network neighborhood connectivity, betweenness centrality, and average shortest path length. Quantitative support for the existence of genomic architecture as a mechanistic basis for evolutionary constraint is then revealed through utilizing these statistical parameters of the protein-protein interaction network in combination with estimators of protein evolution.

CONCLUSIONS

As functional genetic networks are becoming increasingly available, it will now be possible to evaluate functional genetic network constraint against variables describing complex phenotypes and environments, for better understanding of commonly observed deterministic patterns of evolution in non-model organisms. The hypothesis framework and methodological approach outlined herein may help to quantify the extrinsic versus intrinsic dimensions of evolutionary constraint, and result in a better understanding of how fast, effectively, or deterministically organisms adapt.

Complex Molecular Evolution and Expression of Expansin Gene Families in Three Basic Diploid Species of Brassica

Expansins are a kind of structural proteins of the plant cell wall, and they enlarge cells by loosening the cell walls. Therefore, expansins are involved in many growth and development processes. The complete genomic sequences of Brassica rapa, Brassica oleracea and Brassica nigra provide effective platforms for researchers to study expansin genes, and can be compared with analogues in Arabidopsis thaliana. This study identified and characterized expansin families in B. rapa, B. oleracea, and B. nigra. Through the comparative analysis of phylogeny, gene structure, and physicochemical properties, the expansin families were divided into four subfamilies, and then their expansion patterns and evolution details were explored accordingly. Results showed that after the three species underwent independent evolution following their separation from A. thaliana, the expansin families in the three species had increased similarities but fewer divergences. By searching divergences of promoters and coding sequences, significant positive correlations were revealed among orthologs in A. thaliana and the three basic species. Subsequently, differential expressions indicated extensive functional divergences in the expansin families of the three species, especially in reproductive development. Hence, these results support the molecular evolution of basic Brassica species, potential functions of these genes, and genetic improvement of related crops.

Identification of novel synaptonemal complex components in C. elegans

The synaptonemal complex (SC) is a tripartite protein scaffold that forms between homologous chromosomes during meiosis. Although the SC is essential for stable homologue pairing and crossover recombination in diverse eukaryotes, it is unknown how individual components assemble into the highly conserved SC structure. Here we report the biochemical identification of two new SC components, SYP-5 and SYP-6, in Caenorhabditis elegans. SYP-5 and SYP-6 are paralogous to each other and play redundant roles in synapsis, providing an explanation for why these genes have evaded previous genetic screens. Superresolution microscopy reveals that they localize between the chromosome axes and span the width of the SC in a head-to-head manner, similar to the orientation of other known transverse filament proteins. Using genetic redundancy and structure-function analyses to truncate C-terminal tails of SYP-5/6, we provide evidence supporting the role of SC in both limiting and promoting crossover formation.

Gene Balance Predicts Transcriptional Responses Immediately Following Ploidy Change in Arabidopsis thaliana

The gene balance hypothesis postulates that there is selection on gene copy number (gene dosage) to preserve the stoichiometric balance among interacting proteins. This presupposes that gene product abundance is governed by gene dosage and that gene dosage responses are consistent for interacting genes in a dosage-balance-sensitive network or complex. Gene dosage responses, however, have rarely been quantified, and the available data suggest that they are highly variable. We sequenced the transcriptomes of two synthetic autopolyploid accessions of Arabidopsis (Arabidopsis thaliana) and their diploid progenitors, as well as one natural tetraploid and its synthetic diploid produced via haploid induction, to estimate transcriptome size and dosage responses immediately following ploidy change. Similar to what has been observed in previous studies, overall transcriptome size does not exhibit a simple doubling in response to genome doubling, and individual gene dosage responses are highly variable in all three accessions, indicating that expression is not strictly coupled with gene dosage. Nonetheless, putatively dosage balance-sensitive gene groups (Gene Ontology terms, metabolic networks, gene families, and predicted interacting proteins) exhibit smaller and more coordinated dosage responses than do putatively dosage-insensitive gene groups, suggesting that constraints on dosage balance operate immediately following whole-genome duplication and that duplicate gene retention patterns are shaped by selection to preserve dosage balance.

The lasting after-effects of an ancient polyploidy on the genomes of teleosts

The ancestor of most teleost fishes underwent a whole-genome duplication event three hundred million years ago. Despite its antiquity, the effects of this event are evident both in the structure of teleost genomes and in how the surviving duplicated genes still operate to drive form and function. I inferred a set of shared syntenic regions that survive from the teleost genome duplication (TGD) using eight teleost genomes and the outgroup gar genome (which lacks the TGD). I then phylogenetically modeled the TGD's resolution via shared and independent gene losses and applied a new simulation-based statistical test for the presence of bias toward the preservation of genes from one parental subgenome. On the basis of that test, I argue that the TGD was likely an allopolyploidy. I find that duplicate genes surviving from this duplication in zebrafish are less likely to function in early embryo development than are genes that have returned to single copy at some point in this species' history. The tissues these ohnologs are expressed in, as well as their biological functions, lend support to recent suggestions that the TGD was the source of a morphological innovation in the structure of the teleost retina. Surviving duplicates also appear less likely to be essential than singletons, despite the fact that their single-copy orthologs in mouse are no less essential than other genes.

Sensitivity to gene dosage and gene expression affects genes with copy number variants observed among neuropsychiatric diseases

BACKGROUND

Copy number variants (CNVs) have been reported to be associated with diseases, traits, and evolution. However, it is hard to determine which gene should have priority as a target for further functional experiments if a CNV is rare or a singleton. In this study, we attempted to overcome this issue by using two approaches: by assessing the influences of gene dosage sensitivity and gene expression sensitivity. Dosage sensitive genes derived from two-round whole-genome duplication in previous studies. In addition, we proposed a cross-sectional omics approach that utilizes open data from GTEx to assess the effect of whole-genome CNVs on gene expression.

METHODS

Affymetrix Genome-Wide SNP Array 6.0 was used to detect CNVs by PennCNV and CNV Workshop. After quality controls for population stratification, family relationship and CNV detection, 287 patients with narcolepsy, 133 patients with essential hypersomnia, 380 patients with panic disorders, 164 patients with autism, 784 patients with Alzheimer disease and 1280 healthy individuals remained for the enrichment analysis.

RESULTS

Overall, significant enrichment of dosage sensitive genes was found across patients with narcolepsy, panic disorders and autism. Particularly, significant enrichment of dosage-sensitive genes in duplications was observed across all diseases except for Alzheimer disease. For deletions, less or no enrichment of dosage-sensitive genes with deletions was seen in the patients when compared to the healthy individuals. Interestingly, significant enrichments of genes with expression sensitivity in brain were observed in patients with panic disorder and autism. While duplications presented a higher burden, deletions did not cause significant differences when compared to the healthy individuals. When we assess the effect of sensitivity to genome dosage and gene expression at the same time, the highest ratio of enrichment was observed in the group including dosage-sensitive genes and genes with expression sensitivity only in brain. In addition, shared CNV regions among the five neuropsychiatric diseases were also investigated.

CONCLUSIONS

This study contributed the evidence that dosage-sensitive genes are associated with CNVs among neuropsychiatric diseases. In addition, we utilized open data from GTEx to assess the effect of whole-genome CNVs on gene expression. We also investigated shared CNV region among neuropsychiatric diseases.

Allele-specific single-cell RNA sequencing reveals different architectures of intrinsic and extrinsic gene expression noises

Gene expression noise refers to the variation of the expression level of a gene among isogenic cells in the same environment, and has two sources: extrinsic noise arising from the disparity of the cell state and intrinsic noise arising from the stochastic process of gene expression in the same cell state. Due to the low throughput of the existing method for measuring the two noise components, the architectures of intrinsic and extrinsic expression noises remain elusive. Using allele-specific single-cell RNA sequencing, we here estimate the two noise components of 3975 genes in mouse fibroblast cells. Our analyses verify predicted influences of several factors such as the TATA-box and microRNA targeting on intrinsic or extrinsic noises and reveal gene function-associated noise trends implicating the action of natural selection. These findings unravel differential regulations, optimizations, and biological consequences of intrinsic and extrinsic noises and can aid the construction of desired synthetic circuits.

Genome-Wide Genotype-Expression Relationships Reveal Both Copy Number and Single Nucleotide Differentiation Contribute to Differential Gene Expression between Stickleback Ecotypes

Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.