Evolution of 5' untranslated region length and gene expression reprogramming in yeasts

The DEAD-Box RNA Helicase Ded1 Is Associated with Translating Ribosomes

DEAD-box RNA helicases are ATP-dependent RNA binding proteins and RNA-dependent ATPases that possess weak, nonprocessive unwinding activity in vitro, but they can form long-lived complexes on RNAs when the ATPase activity is inhibited. Ded1 is a yeast DEAD-box protein, the functional ortholog of mammalian DDX3, that is considered important for the scanning efficiency of the 48S pre-initiation complex ribosomes to the AUG start codon. We used a modified PAR-CLIP technique, which we call quicktime PAR-CLIP (qtPAR-CLIP), to crosslink Ded1 to 4-thiouridine-incorporated RNAs in vivo using UV light centered at 365 nm. The irradiation conditions are largely benign to the yeast cells and to Ded1, and we are able to obtain a high efficiency of crosslinking under physiological conditions. We find that Ded1 forms crosslinks on the open reading frames of many different mRNAs, but it forms the most extensive interactions on relatively few mRNAs, and particularly on mRNAs encoding certain ribosomal proteins and translation factors. Under glucose-depletion conditions, the crosslinking pattern shifts to mRNAs encoding metabolic and stress-related proteins, which reflects the altered translation. These data are consistent with Ded1 functioning in the regulation of translation elongation, perhaps by pausing or stabilizing the ribosomes through its ATP-dependent binding.

A highly contiguous genome assembly reveals sources of genomic novelty in the symbiotic fungus Rhizophagus irregularis

The root systems of most plant species are aided by the soil-foraging capacities of symbiotic arbuscular mycorrhizal (AM) fungi of the Glomeromycotina subphylum. Despite recent advances in our knowledge of the ecology and molecular biology of this mutualistic symbiosis, our understanding of the AM fungi genome biology is just emerging. Presented here is a close to T2T genome assembly of the model AM fungus Rhizophagus irregularis DAOM197198, achieved through Nanopore long-read DNA sequencing and Hi-C data. This haploid genome assembly of R. irregularis, alongside short- and long-read RNA-Sequencing data, was used to produce a comprehensive annotation catalog of gene models, repetitive elements, small RNA loci, and DNA cytosine methylome. A phylostratigraphic gene age inference framework revealed that the birth of genes associated with nutrient transporter activity and transmembrane ion transport systems predates the emergence of Glomeromycotina. While nutrient cycling in AM fungi relies on genes that existed in ancestor lineages, a burst of Glomeromycotina-restricted genetic innovation is also detected. Analysis of the chromosomal distribution of genetic and epigenetic features highlights evolutionarily young genomic regions that produce abundant small RNAs, suggesting active RNA-based monitoring of genetic sequences surrounding recently evolved genes. This chromosome-scale view of the genome of an AM fungus genome reveals previously unexplored sources of genomic novelty in an organism evolving under an obligate symbiotic life cycle.

A New Expression System Based on Psychrotolerant Debaryomyces macquariensis Yeast and Its Application to the Production of Cold-Active β-d-Galactosidase from Paracoccus sp. 32d

Yeasts provide attractive host/vector systems for heterologous gene expression. The currently used yeast-based expression platforms include mesophilic and thermotolerant species. A eukaryotic expression system working at low temperatures could be particularly useful for the production of thermolabile proteins and proteins that tend to form insoluble aggregates. For this purpose, an expression system based on an Antarctic psychrotolerant yeast Debaryomyces macquariensis strain D50 that is capable of growing at temperatures ranging from 0 to 30 °C has been developed. The optimal physical culture conditions for D. macquariensis D50 in a fermenter are as follows: temperature 20 °C, pH 5.5, aeration rate of 1.5 vvm, and a stirring speed of 300 rpm. Four integrative plasmid vectors equipped with an expression cassette containing the constitutive GAP promoter and CYC1 transcriptional terminator from D. macquariensis D50 were constructed and used to clone and express a gene-encoding cold-active β-d-galactosidase of Paracoccus sp. 32d. The yield was 1150 U/L of recombinant yeast culture. Recombinant D. macquariensis D50 strains were mitotically stable under both selective and non-selective conditions. The D. macquariensis D50 host/vector system has been successfully utilized for the synthesis of heterologous thermolabile protein, and it can be an alternative to other microbial expression systems.

Genetic Diversity of Transcription Factor Genes in Triticum and Mining for Promising Haplotypes for Beneficial Agronomic Traits

Transcription factor (TF) is a class of the sequence-specific DNA-binding proteins that modulate the transcription of target genes, and thus regulate their expressions. Variations in TF are the crucial determinants for phenotypic traits. Although much progress has been made in the functions of TF genes in wheat, one of the most important staple crops globally, the diversity of TF genes in wheat and its progenitors are not well understood, especially the agronomically promising haplotypes have not yet been characterized. Here, we identified a total of 6,023 TF genes from hexaploid wheat through a genome-search method and classified them into 59 gene families based on the conserved domain. The characteristics and dN/dS values of these genes showed evidently selective effects. Based on re-sequencing data, we found a strong genetic bottleneck among these TF genes on A and D subgenomes while no found in B subgenome during wheat domestication. Combined with selective signals and known QTLs on the whole genome, 21 TF genes were preliminarily found to be associated with yield-related traits. The haplotype frequency of these TF genes was further investigated in bread wheat and its progenitors and 13 major haplotypes were the casual loci related to key traits. Finally, the tissue-specific TF genes were also identified using RNA-seq analysis. This study provided insights into the diversity and evolution of TF genes and the identified TF genes and excellent haplotypes associating with traits will contribute to wheat genetic improvement.

Data-driven identification of inherent features of eukaryotic stress-responsive genes

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

New Genomic Signals Underlying the Emergence of Human Proto-Genes

De novo genes are novel genes which emerge from non-coding DNA. Until now, little is known about de novo genes’ properties, correlated to their age and mechanisms of emergence. In this study, we investigate four related properties: introns, upstream regulatory motifs, 5′ Untranslated regions (UTRs) and protein domains, in 23,135 human proto-genes. We found that proto-genes contain introns, whose number and position correlates with the genomic position of proto-gene emergence. The origin of these introns is debated, as our results suggest that 41% of proto-genes might have captured existing introns, and 13.7% of them do not splice the ORF. We show that proto-genes which emerged via overprinting tend to be more enriched in core promotor motifs, while intergenic and intronic genes are more enriched in enhancers, even if the TATA motif is most commonly found upstream in these genes. Intergenic and intronic 5′ UTRs of proto-genes have a lower potential to stabilise mRNA structures than exonic proto-genes and established human genes. Finally, we confirm that proteins expressed by proto-genes gain new putative domains with age. Overall, we find that regulatory motifs inducing transcription and translation of previously non-coding sequences may facilitate proto-gene emergence. Our study demonstrates that introns, 5′ UTRs, and domains have specific properties in proto-genes. We also emphasize that the genomic positions of de novo genes strongly impacts these properties.

Modeling the ribosomal small subunit dynamic in Saccharomyces cerevisiae based on TCP-seq data

Translation Complex Profile Sequencing (TCP-seq), a protocol that was developed and implemented on Saccharomyces cerevisiae, provides the footprints of the small subunit (SSU) of the ribosome (with additional factors) across the entire transcriptome of the analyzed organism. In this study, based on the TCP-seq data, we developed for the first-time a predictive model of the SSU density and analyzed the effect of transcript features on the dynamics of the SSU scan in the 5′UTR. Among others, our model is based on novel tools for detecting complex statistical relations tailored to TCP-seq. We quantitatively estimated the effect of several important features, including the context of the upstream AUG, the upstream ORF length and the mRNA folding strength. Specifically, we suggest that around 50% of the variance related to the read counts (RC) distribution near a start codon can be attributed to the AUG context score. We provide the first large scale direct quantitative evidence that shows that indeed AUG context affects the small sub-unit movement. In addition, we suggest that strong folding may cause the detachment of the SSU from the mRNA. We also identified a number of novel sequence motifs that can affect the SSU scan; some of these motifs affect transcription factors and RNA binding proteins. The results presented in this study provide a better understanding of the biophysical aspects related to the SSU scan along the 5′UTR and of translation initiation in S. cerevisiae, a fundamental step toward a comprehensive modeling of initiation.

Vps34 and TOR Kinases Coordinate HAC1 mRNA Translation in the Presence or Absence of Ire1-Dependent Splicing

In the budding yeast Saccharomyces cerevisiae, an mRNA, called HAC1, exists in a translationally repressed form in the cytoplasm. Under conditions of cellular stress, such as when unfolded proteins accumulate inside the endoplasmic reticulum (ER), an RNase Ire1 removes an intervening sequence (intron) from the HAC1 mRNA by nonconventional cytosolic splicing. Removal of the intron results in translational derepression of HAC1 mRNA and production of a transcription factor that activates expression of many enzymes and chaperones to increase the protein-folding capacity of the cell. Here, we show that Ire1-mediated RNA cleavage requires Watson-Crick base pairs in two RNA hairpins, which are located at the HAC1 mRNA exon-intron junctions. Then, we show that the translational derepression of HAC1 mRNA can occur independent of cytosolic splicing. These results are obtained from HAC1 variants that translated an active Hac1 protein from the unspliced mRNA. Additionally, we show that the phosphatidylinositol-3-kinase Vps34 and the nutrient-sensing kinases TOR and GCN2 are key regulators of HAC1 mRNA translation and consequently the ER stress responses. Collectively, our data suggest that the cytosolic splicing and the translational derepression of HAC1 mRNA are coordinated by unique and parallel networks of signaling pathways.

Genome-wide analysis of European sea bass provides insights into the evolution and functions of single-exon genes

Several studies have attempted to understand the origin and evolution of single-exon genes (SEGs) in eukaryotic organisms, including fishes, but few have examined the functional and evolutionary relationships between SEGs and multiple-exon gene (MEG) paralogs, in particular the conservation of promoter regions. Given that SEGs originate via the reverse transcription of mRNA from a "parental" MEGs, such comparisons may enable identifying evolutionarily-related SEG/MEG paralogs, which might fulfill equivalent physiological functions. Here, the relationship of SEG proportion with MEG count, gene density, intron count, and chromosome size was assessed for the genome of the European sea bass, Dicentrarchus labrax. Then, SEGs with an MEG parent were identified, and promoter sequences of SEG/MEG paralogs were compared, to identify highly conserved functional motifs. The results revealed a total count of 1,585 (8.3% of total genes) SEGs in the European sea bass genome, which was correlated with MEG count but not with gene density. The significant correlation of SEG content with the number of MEGs suggests that SEGs were continuously and independently generated over evolutionary time following species divergence through retrotranscription events, followed by tandem duplications. Functional annotation showed that the majority of SEGs are functional, as is evident from their expression in RNA-seq data used to support homology-based genome annotation. Differences in 5'UTR and 3'UTR lengths between SEG/MEG paralogs observed in this study may contribute to gene expression divergence between them and therefore lead to the emergence of new SEG functions. The comparison of nonsynonymous to synonymous changes (Ka/Ks) between SEG/MEG parents showed that 74 of them are under positive selection (Ka/Ks > 1; p = .0447). An additional fifteen SEGs with an MEG parent have a common promoter, which implies that they are under the influence of common regulatory networks.

5' Untranslated Region Elements Show High Abundance and Great Variability in Homologous ABCA Subfamily Genes

The 12 members of the ABCA subfamily in humans are known for their ability to transport cholesterol and its derivatives, vitamins, and xenobiotics across biomembranes. Several ABCA genes are causatively linked to inborn diseases, and the role in cancer progression and metastasis is studied intensively. The regulation of translation initiation is implicated as the major mechanism in the processes of post-transcriptional modifications determining final protein levels. In the current bioinformatics study, we mapped the features of the 5' untranslated regions (5'UTR) known to have the potential to regulate translation, such as the length of 5'UTRs, upstream ATG codons, upstream open-reading frames, introns, RNA G-quadruplex-forming sequences, stem loops, and Kozak consensus motifs, in the DNA sequences of all members of the subfamily. Subsequently, the conservation of the features, correlations among them, ribosome profiling data as well as protein levels in normal human tissues were examined. The 5'UTRs of ABCA genes contain above-average numbers of upstream ATGs, open-reading frames and introns, as well as conserved ones, and these elements probably play important biological roles in this subfamily, unlike RG4s. Although we found significant correlations among the features, we did not find any correlation between the numbers of 5'UTR features and protein tissue distribution and expression scores. We showed the existence of single nucleotide variants in relation to the 5'UTR features experimentally in a cohort of 105 breast cancer patients. 5'UTR features presumably prepare a complex playground, in which the other elements such as RNA binding proteins and non-coding RNAs play the major role in the fine-tuning of protein expression.

The origin and evolution of a distinct mechanism of transcription initiation in yeasts

The molecular process of transcription by RNA Polymerase II is highly conserved among eukaryotes (“classic model”). A distinct way of locating transcription start sites (TSSs) has been identified in a budding yeast Saccharomyces cerevisiae (“scanning model”). Herein, we applied genomic approaches to elucidate the origin of the scanning model and its underlying genetic mechanisms. We first identified TSSs at single-nucleotide resolution for 12 yeast species using the nAnT-iCAGE technique, which significantly improved the annotations of these genomes by providing accurate 5′ boundaries for protein-coding genes. We then inferred the initiation mechanism of each species based on its TSS maps and genome sequences. We discovered that the scanning model likely originated after the split of Yarrowia lipolytica and the other budding yeasts. Species that use the scanning model showed an adenine-rich region immediately upstream of the TSS that might facilitate TSS selection. Both initiation mechanisms share a strong preference for pyrimidine–purine dinucleotides surrounding the TSS. Our results suggest that the purine is required to accurately recruit the first nucleotide, thereby increasing the chances of a messenger RNA of being capped during mRNA maturation, which is critical for efficient translation initiation during protein biosynthesis. Based on our findings, we propose a model for TSS selection in the scanning-model species, as well as a model for the stepwise process responsible for the origin and evolution of the scanning model.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Translational control of one-carbon metabolism underpins ribosomal protein phenotypes in cell division and longevity

A long-standing problem is how cells that lack one of the highly similar ribosomal proteins (RPs) often display distinct phenotypes. Yeast and other organisms live longer when they lack specific ribosomal proteins, especially of the large 60S subunit of the ribosome. However, longevity is neither associated with the generation time of RP deletion mutants nor with bulk inhibition of protein synthesis. Here, we queried actively dividing RP mutants through the cell cycle. Our data link transcriptional, translational, and metabolic changes to phenotypes associated with the loss of paralogous RPs. We uncovered translational control of transcripts encoding enzymes of methionine and serine metabolism, which are part of one-carbon (1C) pathways. Cells lacking Rpl22Ap, which are long-lived, have lower levels of metabolites associated with 1C metabolism. Loss of 1C enzymes increased the longevity of wild type cells. 1C pathways exist in all organisms and targeting the relevant enzymes could represent longevity interventions.

5'-UTR and ORF elements, as well as the 3'-UTR regulate the translation of Cyclin

Mammalian oocyte maturation is wholly dependent on the translation of accumulated maternal transcripts. Therefore, measuring the translation of specific genes, especially Ccnb1 and Ccnb2, which are key regulators of the oocyte cell cycle in mice, is essential to monitor oocyte cell cycle progression. For this purpose, almost all previous research has used a reporter construct containing the 3'-untranslated region (UTR) of Ccnb. It is based on the concept that the 3'-UTR is the main modulator of translation. Here, we investigated the expression pattern of Renilla luciferase (RL) reporters combining the 5'-UTR and/or open reading frame (ORF) as well as the 3'-UTR (RL-3', 5'-RL-3', RL-ORF-3', and 5'-RL-ORF-3') of Ccnb1 and Ccnb2 in somatic cells and mouse oocytes. The addition of the 5'-UTR and/or ORF of Ccnb altered the expression of the RL-3' reporter in HEK293T cells and mouse oocytes. The ORF tended to suppress RL expression, whereas the 5'-UTR enhanced the expression in most cases. The increased rate in expression was the highest when only the 3'-UTR of Ccnb1 (RL-3') was used, whereas the 5'-RL-ORF-3' reporter showed a relatively lower increase during oocyte maturation. For Ccnb2, the RL-ORF-3' reporter showed the largest increase, and other reporters exhibited a similar increase in expression during oocyte maturation. Results show that the expression of these genes is modulated not only by the 3'-UTR but also by the 5'-UTR and ORF. Therefore, special caution should be taken when using only the 3'-UTR to monitor the expression of specific genes.

Influence of KISS1 gene polymorphisms on the risk of polycystic ovary syndrome and its associated variables, in Saudi women

BACKGROUND

Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder, affecting millions of women worldwide. The role of genetic polymorphisms of the KISS1 gene on the development of PCOS is still obscure. This study was designed to investigate the probable influence of KISS1 gene polymorphisms on PCOS and its associated variables: BMI, waist-hip ratio, kisspeptin, LH, FSH, and LH-FSH ratio.

METHODS

The study comprised 104 PCOS women and 109 controls, with age ranging from 19 to 36 years. BMI, waist-hip ratio, and circulating levels of kisspeptin, LH, and FSH were measured. DNA was extracted, and genotyping of the KISS1 gene was carried out by nucleotide sequencing. The PCOS-associated variables were analyzed in different genotypes of single nucleotide polymorphisms (SNPs) of the KISS1 gene.

RESULTS

The values of waist-hip ratio (WHR), LH, and LH-FSH ratio were significantly higher in PCOS women than controls. BMI, kisspeptin, and FSH levels exhibited no significant difference between the groups. Six novel SNPs of KISS1 gene were identified. Three: rs372790354G > A, rs12998G > A, and rs35431622A > T were investigated. Among these SNPs, the genotype and allele frequencies of rs372790354 showed significant association with PCOS (GA: p = 0.018, AA: p = 0.022, mutant allele-A: p = 0.021) and the G allele was protective. The values of LH, kisspeptin, and WHR of PCOS women were significantly influenced (p < 0.05) by the AA genotype of rs372790354. The other two SNPs rs12998G > A and rs35431622A > T revealed no significant influence on PCOS and associated variables. Haplotypes were constructed, but there was no significant difference between the patients and controls.

CONCLUSION

In conclusion, this is the first study, which reports a significant influence of KISS1 gene polymorphism (rs372790354G > A) on PCOS and its associated variables. However, more extensive research is necessary to confirm these findings.

Molecular Characterization of Auxin Efflux Carrier- ABCB1 in hexaploid wheat

Auxin is an important phytohormone that regulates response, differentiation, and development of plant cell, tissue, and organs. Along with its local production, long-distance transport coordinated by the efflux/influx membrane transporters is instrumental in plant development and architecture. In the present study, we cloned and characterized a wheat (Triticum aestivum) auxin efflux carrier ABCB1. The TaABCB1 was physically localized to the proximal 15% of the short arm of wheat homoeologous group 7 chromosomes. Size of the Chinese spring (CS) homoeologs genomic copies ranged from 5.3–6.2 kb with the 7A copy being the largest due to novel insertions in its third intron. The three homoeologous copies share 95–97% sequence similarity at the nucleotide, 98–99% amino acid, and overall Q-score of 0.98 at 3-D structure level. Though detected in all analyzed tissues, TaABCB1 predominantly expressed in the meristematic tissues likely due to the presence of meristem-specific activation regulatory element identified in the promoter region. RNAi plants of TaABCB1 gene resulted in reduced plant height and increased seed width. Promoter analysis revealed several responsive elements detected in the promoter region including that for different hormones as auxin, gibberellic acid, jasmonic acid and abscisic acid, light, and circadian regulated elements.

Pervasive and dynamic transcription initiation in Saccharomyces cerevisiae

Transcription initiation is finely regulated to ensure proper expression and function of genes. The regulated transcription initiation in response to various environmental stimuli in a classic model organism Saccharomyces cerevisiae has not been systematically investigated. In this study, we generated quantitative maps of transcription start sites (TSSs) at a single-nucleotide resolution for S. cerevisiae grown in nine different conditions using no-amplification nontagging Cap analysis of gene expression (nAnT-iCAGE) sequencing. We mapped ∼1 million well-supported TSSs, suggesting highly pervasive transcription initiation in the compact genome of the budding yeast. The comprehensive TSS maps allowed us to identify core promoters for ∼96% verified protein-coding genes. We corrected misannotation of translation start codon for 122 genes and suggested an alternative start codon for 57 genes. We found that 56% of yeast genes are controlled by multiple core promoters, and alternative core promoter usage by a gene is widespread in response to changing environments. Most core promoter shifts are coupled with altered gene expression, indicating that alternative core promoter usage might play an important role in controlling gene transcriptional activities. Based on their activities in responding to environmental cues, we divided core promoters into constitutive class (55%) and inducible class (45%). The two classes of core promoters display distinctive patterns in transcriptional abundance, chromatin structure, promoter shape, and sequence context. In summary, our study improved the annotation of the yeast genome and demonstrated a much more pervasive and dynamic nature of transcription initiation in yeast than previously recognized.

Molecular characterization of laccase genes from the basidiomycete Trametes hirsuta Bm-2 and analysis of the 5' untranslated region (5'UTR)

The aim of this study was to identify and characterize laccase genes produced by Trametes hirsuta Bm-2 in a liquid medium, both with and without induction. The amplification of 5'and 3'regions of laccase sequences was obtained by the RACE-PCR method, and these were assembled to obtain a cDNA of total length. Two new laccase genes were isolated from basal medium (lac-B) and lignocellulosic grapefruit substrate (lac-T), both encoding open reading frames of 2566 bp. Both laccase-predicted proteins consisted of 521 amino acids, four copper-binding regions, a signal peptide, and five potential glycosilation sites (Asn-Xaa-Ser/Tre). Moreover, the deduced amino acid sequences share about 76-85% identity with other laccases of WRF. Sequence comparison showed 47 synonymous point mutations between lac-B and lac-T. In addition, 5' untranslated regions (UTR) of laccase genes lac-B and lac-T showed differences in length and number of regulatory elements that may affect transcriptional or translational expression of these genes.

Variation of gene expression in plants is influenced by gene architecture and structural properties of promoters

In higher eukaryotes, gene architecture and structural properties of promoters have emerged as significant factors influencing variation in number of transcripts (expression level) and specificity of gene expression in a tissue (expression breadth), which eventually shape the phenotype. In this study, transcriptome data of different tissue types at various developmental stages of A. thaliana, O. sativa, S. bicolor and Z. mays have been used to understand the relationship between properties of gene components and its expression. Our findings indicate that in plants, among all gene architecture and structural properties of promoters, compactness of genes in terms of intron content is significantly linked to gene expression level and breadth, whereas in human an exactly opposite scenario is seen. In plants, for the first time we have carried out a quantitative estimation of effect of a particular trait on expression level and breadth, by using multiple regression analysis and it confirms that intron content of primary transcript (as %) is a powerful determinant of expression breadth. Similarly, further regression analysis revealed that among structural properties of the promoters, stability is negatively linked to expression breadth, while DNase1 sensitivity strongly governs gene expression breadth in monocots and gene expression level in dicots. In addition, promoter regions of tissue specific genes are found to be enriched with TATA box and Y-patch motifs. Finally, multi copy orthologous genes in plants are found to be longer, highly regulated and tissue specific.

Genetic Variations of DAOA (rs947267 and rs3918342) and COMT Genes (rs165599 and rs4680) in Schizophrenia and Bipolar I Disorder

Introduction:

Genetic and environmental factors are involved in the incidence of schizophrenia and bipolar disorder. Many reports confirm that several common genes are connected with these two psychotic disorders. Several neurotransmitters may be involved in the molecular mechanisms of schizophrenia and bipolar disorder. We aimed to estimate the role of two talent genes: DAOA in neurotransmission of glutamate and COMT in neurotransmission of dopamine to guide the treatment of schizophrenia and bipolar disorder.

Methods:

Blood samples (n=100 for schizophrenia, n=100 for bipolar I disorder and n=127 for case control) were collected from individuals unrelated in the southwest of Iran. The SNPs (rs947267 and rs3918342 for DAOA gene/rs165599 and rs4680 for COMT gene) were genotyped using the PCR-RFLP method. Our finding was studied by logistic regression and Mantel-Haenszel Chi-square tests.

Results:

We observed an association in rs3918342, rs165599 and rs4680 single nucleotide polymorphisms and schizophrenia and bipolar I disorder. In addition, our data demonstrated that the rs947267 was related to bipolar I disorder but there was no association between this SNP and schizophrenia.

Conclusion:

In conclusion, this result supports the hypothesis that variations in DAOA and COMT genes may play a role in schizophrenia and bipolar disorder.

Fluorescent Guide RNAs Facilitate Development of Layered Pol II-Driven CRISPR Circuits

Efficient clustered regularly interspaced short palindromic repeat (CRISPR) guide RNA (gRNA) expression from RNA Polymerase II (Pol II) promoters will aid in construction of complex CRISPR-based synthetic gene networks. Yet, we require tools to properly visualize gRNA directly to quantitatively study the corresponding network behavior. To address this need, we employed a fluorescent gRNA (fgRNA) to visualize synthetic CRISPR network dynamics without affecting gRNA functionality. We show that studying gRNA dynamics directly enables circuit modification and improvement of network function in Pol II-driven CRISPR circuits. This approach generates information necessary for optimizing the overall function of these networks and provides insight into the hurdles remaining in Pol II-regulated gRNA expression.

Toward Predictable 5'UTRs in Saccharomyces cerevisiae: Development of a yUTR Calculator

Fine-tuning biosynthetic pathways is crucial for the development of economic feasible microbial cell factories. Therefore, the use of computational models able to predictably design regulatory sequences for pathway engineering proves to be a valuable tool, especially for modifying genes at the translational level. In this study we developed a computational approach for the de novo design of 5'-untranslated regions (5'UTRs) in Saccharomyces cerevisiae with a predictive outcome on translation initiation rate. On the basis of existing data, a partial least-squares (PLS) regression model was trained and showed good performance on predicting protein abundances of an independent test set. This model was further used for the construction of a "yUTR calculator" that can design 5'UTR sequences with a diverse range of desired translation efficiencies. The predictive power of our yUTR calculator was confirmed in vivo by different representative case studies. As such, these results show the great potential of data driven approaches for reliable pathway engineering in S. cerevisiae.

Estimating the parameters of background selection and selective sweeps in Drosophila in the presence of gene conversion

We used whole-genome resequencing data from a population of Drosophila melanogaster to investigate the causes of the negative correlation between the within-population synonymous nucleotide site diversity (π_S ) of a gene and its degree of divergence from related species at nonsynonymous nucleotide sites (K_A ). By using the estimated distributions of mutational effects on fitness at nonsynonymous and UTR sites, we predicted the effects of background selection at sites within a gene on π_S and found that these could account for only part of the observed correlation between π_S and K_A We developed a model of the effects of selective sweeps that included gene conversion as well as crossing over. We used this model to estimate the average strength of selection on positively selected mutations in coding sequences and in UTRs, as well as the proportions of new mutations that are selectively advantageous. Genes with high levels of selective constraint on nonsynonymous sites were found to have lower strengths of positive selection and lower proportions of advantageous mutations than genes with low levels of constraint. Overall, background selection and selective sweeps within a typical gene reduce its synonymous diversity to ∼75% of its value in the absence of selection, with larger reductions for genes with high K_A Gene conversion has a major effect on the estimates of the parameters of positive selection, such that the estimated strength of selection on favorable mutations is greatly reduced if it is ignored.

Identification of Arabidopsis genic and non-genic promoters by paired-end sequencing of TSS tags

Information about transcription start sites (TSSs) provides baseline data for the analysis of promoter architecture. In this paper we used paired- and single-end deep sequencing to analyze Arabidopsis TSS tags from several libraries prepared from roots, shoots, flowers and etiolated seedlings. The clustering of approximately 33 million mapped TSS tags led to the identification of 324 461 promoters that covered 79.7% (21 672/27 206) of protein-coding genes in the Arabidopsis genome. In addition we identified intragenic, antisense and orphan promoters that were not associated with any gene models. Of these, intragenic promoters exhibited unique characteristics regarding dinucleotide sequences at TSSs and core promoter element composition, suggesting that these promoters use different mechanisms of transcriptional initiation. An analysis of base composition with regard to promoter position revealed a low GC content throughout the promoter region and several local strand biases that were evident for TATA-type promoters, but not for Coreless-type promoters. Most observed strand biases coincided with strand biases of single nucleotide polymorphism rate. Our analysis also revealed that transcription of a gene is supported by an average of 2.7 genic promoters, among which one specific promoter, designated as a top promoter, substantially determines the expression level of the gene.

Comparation of the effects of different 5'-untranslated regions (UTRs) on gene expression in HEK293 cells

OBJECTIVES

To evaluate four 5'-UTRs on GFP expression in HEK293T cells.

RESULTS

The recombinant plasmids were constructed by restriction enzyme digestion, digestion and DNA sequencing. Quantitative real-time PCR and western blotting results showed that the transcription and translation level of PPRV-GFP mRNA was significantly lower than that of the other reporters. The transcription and translation level of ChEF1-GFP was the highest in HEK293T cells.

CONCLUSIONS

Different UTRs can significantly affect protein expression. Additionally, the findings also will be useful in biological applications that require tuning of gene expression and system optimization.

Evolution of translation initiation sequences using in vitro yeast ribosome display

We report a novel in vitro yeast ribosome display method based on cell-free protein synthesis (CFPS) using linear DNA templates. We demonstrate that our platform can enrich a target gene from a model library by 100-fold per round of selection. We demonstrate the utility of our approach by evolving cap-independent translation initiation (CITI) sequences, which result in a 13-fold increase in CFPS yields after four rounds of selection, and a threefold further increase by placing the beneficial short sequences in tandem. We also show that 12 of the selected CITI sequences permit precise control of gene expression in vitro over a range of up to 80-fold by enhancing translation (and not as cryptic promoters). These 12 sequences are then shown to tune protein expression in vivo, though likely due to a different mechanism. Looking forward, yeast ribosome display holds promise for evolving libraries of proteins and DNA regulatory parts for protein engineering and synthetic biology. Biotechnol. Bioeng. 2016;113: 1777-1786. © 2016 Wiley Periodicals, Inc.

Comparisons between Arabidopsis thaliana and Drosophila melanogaster in relation to Coding and Noncoding Sequence Length and Gene Expression

There is a continuing interest in the analysis of gene architecture and gene expression to determine the relationship that may exist. Advances in high-quality sequencing technologies and large-scale resource datasets have increased the understanding of relationships and cross-referencing of expression data to the large genome data. Although a negative correlation between expression level and gene (especially transcript) length has been generally accepted, there have been some conflicting results arising from the literature concerning the impacts of different regions of genes, and the underlying reason is not well understood. The research aims to apply quantile regression techniques for statistical analysis of coding and noncoding sequence length and gene expression data in the plant, Arabidopsis thaliana, and fruit fly, Drosophila melanogaster, to determine if a relationship exists and if there is any variation or similarities between these species. The quantile regression analysis found that the coding sequence length and gene expression correlations varied, and similarities emerged for the noncoding sequence length (5′ and 3′ UTRs) between animal and plant species. In conclusion, the information described in this study provides the basis for further exploration into gene regulation with regard to coding and noncoding sequence length.

Control of relative timing and stoichiometry by a master regulator

Developmental processes in cells require a series of complex steps. Often only a single master regulator activates genes in these different steps. This poses several challenges: some targets need to be ordered temporally, while co-functional targets may need to be synchronized in both time and expression level. Here we study in single cells the dynamic activation patterns of early meiosis genes in budding yeast, targets of the meiosis master regulator Ime1. We quantify the individual roles of the promoter and protein levels in expression pattern control, as well as the roles of individual promoter elements. We find a consistent expression pattern difference between a non-cofunctional pair of genes, and a highly synchronized activation of a co-functional pair. We show that dynamic control leading to these patterns is distributed between promoter, gene and external regions. Through specific reciprocal changes to the promoters of pairs of genes, we show that different genes can use different promoter elements to reach near identical activation patterns.

Allelic variation within the S-adenosyl-L-homocysteine hydrolase gene family is associated with wood properties in Chinese white poplar (Populus tomentosa)

Background

S-adenosyl-l-homocysteine hydrolase (SAHH) is the only eukaryotic enzyme capable of S-adenosyl-l-homocysteine (SAH) catabolism for the maintenance of cellular transmethylation potential. Recently, biochemical and genetic studies in herbaceous species have obtained important discoveries in the function of SAHH, and an extensive characterization of SAHH family in even one tree species is essential, but currently lacking.

Results

Here, we first identified the SAHH family from Populus tomentosa using molecular cloning method. Phylogenetic analyses of 28 SAHH proteins from dicotyledons, monocotyledons, and lower plants revealed that the sequences formed two monophyletic groups: the PtrSAHHA with PtoSAHHA and PtrSAHHB with PtoSAHHB. Examination of tissue-specific expression profiles of the PtoSAHH family revealed similar expression patterns; high levels of expression in xylem were found. Nucleotide diversity and linkage disequilibrium (LD) in the PtoSAHH family, sampled from P. tomentosa natural distribution, revealed that PtoSAHH harbors high single-nucleotide polymorphism (SNP) diversity (π=0.01059±0.00122 and 0.00930±0.00079,respectively) and low LD (r² > 0.1, within 800 bp and 2,200 bp, respectively). Using an LD-linkage analysis approach, two noncoding SNPs (PtoSAHHB_1065 and PtoSAHHA_2203) and the corresponding haplotypes were found to significantly associate with α-cellulose content, and a nonsynonymous SNP (PtoSAHHB_410) within the SAHH signature motifs showed significant association with fiber length, with an average of 3.14% of the phenotypic variance explained.

Conclusions

The present study demonstrates that PtoSAHHs were split off prior to the divergence of interspecies in Populus, and SAHHs may play a key role promoting transmethylation reactions in the secondary cell walls biosynthesis in trees. Hence, our findings provide insights into SAHH function and evolution in woody species and also offer a theoretical basis for marker-aided selection breeding to improve the wood quality of Populus.

The promoter and the 5'-untranslated region of rice metallothionein OsMT2b gene are capable of directing high-level gene expression in germinated rice embryos

Critical regions within the rice metallothionein OsMT2b gene promoter are identified and the 5'-untranslated region (5'-UTR) is found essential for the high-level promoter activity in germinated transgenic rice embryos. Many metallothionein (MT) genes are highly expressed in plant tissues. A rice subfamily p2 (type 2) MT gene, OsMT2b, has been shown previously to exhibit the most abundant gene expression in young rice seedling. In the present study, transient expression assays and a transgenic approach were employed to characterize the expression of the OsMT2b gene in rice. We found that the OsMT2b gene is strongly and differentially expressed in germinated rice embryos during seed germination and seedling development. Histochemical staining analysis of transgenic rice carrying OsMT2b::GUS chimeric gene showed that high-level GUS activity was detected in germinated embryos and at the meristematic part of other tissues during germination. Deletion analysis of the OsMT2b promoter revealed that the 5'-flanking region of the OsMT2b between nucleotides -351 and -121 relative to the transcriptional initiation site is important for promoter activity in rice embryos, and this region contains the consensus sequences of G box and TA box. Our study demonstrates that the 5'-untranslated region (5'-UTR) of OsMT2b gene is not only necessary for the OsMT2b promoter activity, but also sufficient to augment the activity of a minimal promoter in both transformed cell cultures and germinated transgenic embryos in rice. We also found that addition of the maize Ubi intron 1 significantly enhanced the OsMT2b promoter activity in rice embryos. Our studies reveal that OsMT2b351-ubi(In) promoter can be applied in plant transformation and represents potential for driving high-level production of foreign proteins in transgenic rice.

Identifying cis-regulatory changes involved in the evolution of aerobic fermentation in yeasts

Gene regulation change has long been recognized as an important mechanism for phenotypic evolution. We used the evolution of yeast aerobic fermentation as a model to explore how gene regulation has evolved and how this process has contributed to phenotypic evolution and adaptation. Most eukaryotes fully oxidize glucose to CO₂ and H₂O in mitochondria to maximize energy yield, whereas some yeasts, such as Saccharomyces cerevisiae and its relatives, predominantly ferment glucose into ethanol even in the presence of oxygen, a phenomenon known as aerobic fermentation. We examined the genome-wide gene expression levels among 12 different yeasts and found that a group of genes involved in the mitochondrial respiration process showed the largest reduction in gene expression level during the evolution of aerobic fermentation. Our analysis revealed that the downregulation of these genes was significantly associated with massive loss of binding motifs of Cbf1p in the fermentative yeasts. Our experimental assays confirmed the binding of Cbf1p to the predicted motif and the activator role of Cbf1p. In summary, our study laid a foundation to unravel the long-time mystery about the genetic basis of evolution of aerobic fermentation, providing new insights into understanding the role of cis-regulatory changes in phenotypic evolution.

Fungal genes in context: genome architecture reflects regulatory complexity and function

Gene context determines gene expression, with local chromosomal environment most influential. Comparative genomic analysis is often limited in scope to conserved or divergent gene and protein families, and fungi are well suited to this approach with low functional redundancy and relatively streamlined genomes. We show here that one aspect of gene context, the amount of potential upstream regulatory sequence maintained through evolution, is highly predictive of both molecular function and biological process in diverse fungi. Orthologs with large upstream intergenic regions (UIRs) are strongly enriched in information processing functions, such as signal transduction and sequence-specific DNA binding, and, in the genus Aspergillus, include the majority of experimentally studied, high-level developmental and metabolic transcriptional regulators. Many uncharacterized genes are also present in this class and, by implication, may be of similar importance. Large intergenic regions also share two novel sequence characteristics, currently of unknown significance: they are enriched for plus-strand polypyrimidine tracts and an information-rich, putative regulatory motif that was present in the last common ancestor of the Pezizomycotina. Systematic consideration of gene UIR in comparative genomics, particularly for poorly characterized species, could help reveal organisms’ regulatory priorities.

Allelic variation in a cellulose synthase gene (PtoCesA4) associated with growth and wood properties in Populus tomentosa

Lignocellulosic biomass from trees provides a renewable feedstock for biofuels, lumber, pulp, paper, and other uses. Dissecting the mechanism underlying natural variation of the complex traits controlling growth and lignocellulose biosynthesis in trees can enable marker-assisted breeding to improve wood quality and yield. Here, we combined linkage disequilibrium (LD)-based association analysis with traditional linkage analysis to detect the genetic effect of a Populus tomentosa cellulose synthase gene, PtoCesA4. PtoCesA4 is strongly expressed in developing xylem and leaves. Nucleotide diversity and LD in PtoCesA4, sampled from the P. tomentosa natural distribution, revealed that PtoCesA4 harbors high single nucleotide polymorphism (SNP) diversity (πT = 0.0080 and θw = 0.0098) and low LD (r(2) ≥ 0.1, within 1400 bp), demonstrating that the potential of a candidate-gene-based LD approach in understanding the molecular basis underlying quantitative variation in this species. By combining single SNP, multi-SNP, and haplotype-based associations in an association population of 460 individuals with single SNP linkage analysis in a family-based linkage populations (1200 individuals), we identified three strong associations (false discovery rate Q < 0.05) in both populations. These include two nonsynonymous markers (SNP49 associated with α-cellulose content and SNP59 associated with fiber width) and a noncoding marker (SNP18 associated with α-cellulose content). Variation in RNA transcript abundance among genotypic classes of SNP49 was confirmed in these two populations. Therefore, combining different methods allowed us to examine functional PtoCesA4 allelic variation underlying natural variation in complex quantitative traits related to growth and lignocellulosic biosynthesis.

Rate-limiting steps in yeast protein translation

Deep sequencing now provides detailed snapshots of ribosome occupancy on mRNAs. We leverage these data to parameterize a computational model of translation, keeping track of every ribosome, tRNA, and mRNA molecule in a yeast cell. We determine the parameter regimes in which fast initiation or high codon bias in a transgene increases protein yield and infer the initiation rates of endogenous Saccharomyces cerevisiae genes, which vary by several orders of magnitude and correlate with 5′ mRNA folding energies. Our model recapitulates the previously reported 5′-to-3′ ramp of decreasing ribosome densities, although our analysis shows that this ramp is caused by rapid initiation of short genes rather than slow codons at the start of transcripts. We conclude that protein production in healthy yeast cells is typically limited by the availability of free ribosomes, whereas protein production under periods of stress can sometimes be rescued by reducing initiation or elongation rates.

•

Computational model of translation tracks all ribosomes, tRNAs, and mRNAs in a cell

•

Translation is generally limited by initiation, not elongation

•

Model allows inference of initiation rates for all yeast genes

•

Ramp of 5′ ribosomes is caused by rapid initiation of short genes

Sequence features of yeast and human core promoters that are predictive of maximal promoter activity

The core promoter is the region in which RNA polymerase II is recruited to the DNA and acts to initiate transcription, but the extent to which the core promoter sequence determines promoter activity levels is largely unknown. Here, we identified several base content and k-mer sequence features of the yeast core promoter sequence that are highly predictive of maximal promoter activity. These features are mainly located in the region 75 bp upstream and 50 bp downstream of the main transcription start site, and their associations hold for both constitutively active promoters and promoters that are induced or repressed in specific conditions. Our results unravel several architectural features of yeast core promoters and suggest that the yeast core promoter sequence downstream of the TATA box (or of similar sequences involved in recruitment of the pre-initiation complex) is a major determinant of maximal promoter activity. We further show that human core promoters also contain features that are indicative of maximal promoter activity; thus, our results emphasize the important role of the core promoter sequence in transcriptional regulation.

The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics

In woody crop plants, the oligosaccharide components of the cell wall are essential for important traits such as bioenergy content, growth, and structural wood properties. UDP-glucuronate decarboxylase (UXS) is a key enzyme in the synthesis of UDP-xylose for the formation of xylans during cell wall biosynthesis. Here, we isolated a multigene family of seven members (PtUXS1-7) encoding UXS from Populus tomentosa, the first investigation of UXSs in a tree species. Analysis of gene structure and phylogeny showed that the PtUXS family could be divided into three groups (PtUXS1/4, PtUXS2/5, and PtUXS3/6/7), consistent with the tissue-specific expression patterns of each PtUXS. We further evaluated the functional consequences of nucleotide polymorphisms in PtUXS1. In total, 243 single-nucleotide polymorphisms (SNPs) were identified, with a high frequency of SNPs (1/18 bp) and nucleotide diversity (π_T = 0.01033, θ_w =0.01280). Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r²<0.1, P<0.001, within 700 bp). SNP- and haplotype-based association analysis showed that nine SNPs (Q <0.10) and 12 haplotypes (P<0.05) were significantly associated with growth and wood property traits in the association population (426 individuals), with 2.70% to 12.37% of the phenotypic variation explained. Four significant single-marker associations (Q <0.10) were validated in a linkage mapping population of 1200 individuals. Also, RNA transcript accumulation varies among genotypic classes of SNP10 was further confirmed in the association population. This is the first comprehensive study of the UXS gene family in woody plants, and lays the foundation for genetic improvements of wood properties and growth in trees using genetic engineering or marker-assisted breeding.

Weak 5'-mRNA secondary structures in short eukaryotic genes

Experimental studies of translation have found that short genes tend to exhibit greater densities of ribosomes than long genes in eukaryotic species. It remains an open question whether the elevated ribosome density on short genes is due to faster initiation or slower elongation dynamics. Here, we address this question computationally using 5′-mRNA folding energy as a proxy for translation initiation rates and codon bias as a proxy for elongation rates. We report a significant trend toward reduced 5′-secondary structure in shorter coding sequences, suggesting that short genes initiate faster during translation. We also find a trend toward higher 5′-codon bias in short genes, suggesting that short genes elongate faster than long genes. Both of these trends hold across a diverse set of eukaryotic taxa. Thus, the elevated ribosome density on short eukaryotic genes is likely caused by differential rates of initiation, rather than differential rates of elongation.