GC level and expression of human coding sequences

qPrimerDB: A Powerful and User-Friendly Database for qPCR Primer Design

Real-time quantitative polymerase chain reaction (qPCR) is a powerful tool for analyzing and quantifying gene expression, and its primer design is the first and most important step. In order to improve the efficiency and effectiveness of primer design, we designed a database qPrimerDB, based on thermodynamics gene-specific for multispecies qPCR primers design.In this chapter, we explained the working principle of the database and detailed the step-by-step practical steps with examples. The valuable and time-saving qPrimerDB database is publicly accessible at http://biodb.swu.edu.cn/qprimerdb and will be routinely updated.

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

BACKGROUND

The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq).

RESULTS

Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability.

CONCLUSION

We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

A comparative analysis of the mobility of 45 proteins in the synaptic bouton

Many proteins involved in synaptic transmission are well known, and their features, as their abundance or spatial distribution, have been analyzed in systematic studies. This has not been the case, however, for their mobility. To solve this, we analyzed the motion of 45 GFP-tagged synaptic proteins expressed in cultured hippocampal neurons, using fluorescence recovery after photobleaching, particle tracking, and modeling. We compared synaptic vesicle proteins, endo- and exocytosis cofactors, cytoskeleton components, and trafficking proteins. We found that movement was influenced by the protein association with synaptic vesicles, especially for membrane proteins. Surprisingly, protein mobility also correlated significantly with parameters as the protein lifetimes, or the nucleotide composition of their mRNAs. We then analyzed protein movement thoroughly, taking into account the spatial characteristics of the system. This resulted in a first visualization of overall protein motion in the synapse, which should enable future modeling studies of synaptic physiology.

Evidence of distinct gene functional patterns in GC-poor and GC-rich isochores in Bos taurus

Vertebrate genomes are mosaics of megabase-size DNA segments with a fairly homogeneous base composition, called isochores. They are divided into five families characterized by different guanine-cytosine (GC) levels and linked to several functional and structural properties. The increased availability of fully sequenced genomes allows the investigation of isochores in several species, assessing their level of conservation across vertebrate genomes. In this work, we characterized the isochores in Bos taurus using the ARS-UCD1.2 genome version. The comparison of our results with the well-studied human isochores and those of other mammals revealed a large conservation in isochore families, in number, average GC levels and gene density. Exceptions to the established increase in gene density with the increase in isochores (GC%) were observed for the following gene biotypes: tRNA, small nuclear RNA, small nucleolar RNA and pseudogenes that have their maximum number in H2 and H1 isochores. Subsequently, we assessed the ontology of all gene biotypes looking for functional classes that are statistically over- or under-represented in each isochore. Receptor activity and sensory perception pathways were significantly over-represented in L1 and L2 (GC-poor) isochores. This was also validated for the horse genome. Our analysis of housekeeping genes confirmed a preferential localization in GC-rich isochores, as reported in other species. Finally, we assessed the SNP distribution of a bovine high-density SNP chip across the isochores, finding a higher density in the GC-rich families, reflecting a potential bias in the chip, widely used for genetic selection and biodiversity studies.

qPrimerDB: a thermodynamics-based gene-specific qPCR primer database for 147 organisms

Real-time quantitative polymerase chain reaction (qPCR) is one of the most important methods for analyzing the expression patterns of target genes. However, successful qPCR experiments rely heavily on the use of high-quality primers. Various qPCR primer databases have been developed to address this issue, but these databases target only a few important organisms. Here, we developed the qPrimerDB database, founded on an automatic gene-specific qPCR primer design and thermodynamics-based validation workflow. The qPrimerDB database is the most comprehensive qPCR primer database available to date, with a web front-end providing gene-specific and pre-computed primer pairs across 147 important organisms, including human, mouse, zebrafish, yeast, thale cress, rice and maize. In this database, we provide 3331426 of the best primer pairs for each gene, based on primer pair coverage, as well as 47760359 alternative gene-specific primer pairs, which can be conveniently batch downloaded. The specificity and efficiency was validated for qPCR primer pairs for 66 randomly selected genes, in six different organisms, through qPCR assays and gel electrophoresis. The qPrimerDB database represents a valuable, timesaving resource for gene expression analysis. This resource, which will be routinely updated, is publically accessible at http://biodb.swu.edu.cn/qprimerdb.

The codon sequences predict protein lifetimes and other parameters of the protein life cycle in the mouse brain

The homeostasis of the proteome depends on the tight regulation of the mRNA and protein abundances, of the translation rates, and of the protein lifetimes. Results from several studies on prokaryotes or eukaryotic cell cultures have suggested that protein homeostasis is connected to, and perhaps regulated by, the protein and the codon sequences. However, this has been little investigated for mammals in vivo. Moreover, the link between the coding sequences and one critical parameter, the protein lifetime, has remained largely unexplored, both in vivo and in vitro. We tested this in the mouse brain, and found that the percentages of amino acids and codons in the sequences could predict all of the homeostasis parameters with a precision approaching experimental measurements. A key predictive element was the wobble nucleotide. G-/C-ending codons correlated with higher protein lifetimes, protein abundances, mRNA abundances and translation rates than A-/U-ending codons. Modifying the proportions of G-/C-ending codons could tune these parameters in cell cultures, in a proof-of-principle experiment. We suggest that the coding sequences are strongly linked to protein homeostasis in vivo, albeit it still remains to be determined whether this relation is causal in nature.

Cellular Genomic Sites of Hepatitis B Virus DNA Integration

Infection with the Hepatitis B Virus (HBV) is one of the strongest risk-factors for liver cancer (hepatocellular carcinoma, HCC). One of the reported drivers of HCC is the integration of HBV DNA into the host cell genome, which may induce pro-carcinogenic pathways. These reported pathways include: induction of chromosomal instability; generation of insertional mutagenesis in key cancer-associated genes; transcription of downstream cancer-associated cellular genes; and/or formation of a persistent source of viral protein expression (particularly HBV surface and X proteins). The contribution of each of these specific mechanisms towards carcinogenesis is currently unclear. Here, we review the current knowledge of specific sites of HBV DNA integration into the host genome, which sheds light on these mechanisms. We give an overview of previously-used methods to detect HBV DNA integration and the enrichment of integration events in specific functional and structural cellular genomic sites. Finally, we posit a theoretical model of HBV DNA integration during disease progression and highlight open questions in the field.

Transcriptome activity of isochores during preimplantation process in human and mouse

This work investigates the role of isochores during preimplantation process. Using RNA-seq data from human and mouse preimplantation stages, we created the spatio-temporal transcriptional profiles of the isochores during preimplantation. We found that from early to late stages, GC-rich isochores increase their expression while GC-poor ones decrease it. Network analysis revealed that modules with few coexpressed isochores are GC-poorer than medium-large ones, characterized by an opposite expression as preimplantation advances, decreasing and increasing respectively. Our results reveal a functional contribution of the isochores, supporting the presence of structural-functional interactions during maturation and early-embryonic development.

All-or-(N)One - an epistemological characterization of the human tumorigenic neuronal paralogous FAM72 gene loci

FAM72 is a novel neuronal progenitor cell (NPC) self-renewal supporting protein expressed under physiological conditions at low levels in other tissues. Accumulating data indicate the potential pivotal tumourigenic effects of FAM72. Our in silico human genome-wide analysis (GWA) revealed that the FAM72 gene family consists of four human-specific paralogous members, all of which are located on chromosome (chr) 1. Unique asymmetric FAM72 segmental gene duplications are most likely to have occurred in conjunction with the paired genomic neighbour SRGAP2 (SLIT-ROBO Rho GTPase activating protein), as both genes have four paralogues in humans but only one vertebra-emerging orthologue in all other species. No species with two or three FAM72/SRGAP2 gene pairs could be identified, and the four exclusively human-defining ohnologues, with different mutation patterns in Homo neanderthalensis and Denisova hominin, may remain under epigenetic control through long non-coding (lnc) RNAs.

A comparative study and a phylogenetic exploration of the compositional architectures of mammalian nuclear genomes

For the past four decades the compositional organization of the mammalian genome posed a formidable challenge to molecular evolutionists attempting to explain it from an evolutionary perspective. Unfortunately, most of the explanations adhered to the "isochore theory," which has long been rebutted. Recently, an alternative compositional domain model was proposed depicting the human and cow genomes as composed mostly of short compositionally homogeneous and nonhomogeneous domains and a few long ones. We test the validity of this model through a rigorous sequence-based analysis of eleven completely sequenced mammalian and avian genomes. Seven attributes of compositional domains are used in the analyses: (1) the number of compositional domains, (2) compositional domain-length distribution, (3) density of compositional domains, (4) genome coverage by the different domain types, (5) degree of fit to a power-law distribution, (6) compositional domain GC content, and (7) the joint distribution of GC content and length of the different domain types. We discuss the evolution of these attributes in light of two competing phylogenetic hypotheses that differ from each other in the validity of clade Euarchontoglires. If valid, the murid genome compositional organization would be a derived state and exhibit a high similarity to that of other mammals. If invalid, the murid genome compositional organization would be closer to an ancestral state. We demonstrate that the compositional organization of the murid genome differs from those of primates and laurasiatherians, a phenomenon previously termed the "murid shift," and in many ways resembles the genome of opossum. We find no support to the "isochore theory." Instead, our findings depict the mammalian genome as a tapestry of mostly short homogeneous and nonhomogeneous domains and few long ones thus providing strong evidence in favor of the compositional domain model and seem to invalidate clade Euarchontoglires.

Impact of GC content on gene expression pattern in chicken

Background

GC content varies greatly between different genomic regions in many eukaryotes. In order to determine whether this organization named isochore organization influences gene expression patterns, the relationship between GC content and gene expression has been investigated in man and mouse. However, to date, this question is still a matter for debate. Among the avian species, chicken (Gallus gallus) is the best studied representative with a complete genome sequence. The distinctive features and organization of its sequence make it a good model to explore important issues in genome structure and evolution.

Methods

Only nuclear genes with complete information on protein-coding sequence with no evidence of multiple-splicing forms were included in this study. Chicken protein coding sequences, complete mRNA sequences (or full length cDNA sequences), and 5^′ untranslated region sequences (5^′ UTR) were downloaded from Ensembl and chicken expression data originated from a previous work. Three indices i.e. expression level, expression breadth and maximum expression level were used to measure the expression pattern of a given gene. CpG islands were identified using hgTables of the UCSC Genome Browser. Correlation analysis between variables was performed by SAS Proprietary Software Release 8.1.

Results

In chicken, the GC content of 5^′ UTR is significantly and positively correlated with expression level, expression breadth, and maximum expression level, whereas that of coding sequences and introns and at the third coding position are negatively correlated with expression level and expression breadth, and not correlated with maximum expression level. These significant trends are independent of recombination rate, chromosome size and gene density. Furthermore, multiple linear regression analysis indicated that GC content in genes could explain approximately 10% of the variation in gene expression.

Conclusions

GC content is significantly associated with gene expression pattern and could be one of the important regulation factors in the chicken genome.

Transcriptome map of mouse isochores in embryonic and neonatal cortex

Several studies on adult tissues agree on the presence of a positive effect of the genomic and genic base composition on mammalian gene expression. Recent literature supports the idea that during developmental processes GC-poor genomic regions are preferentially implicated. We investigate the relationship between the compositional properties of the isochores and of the genes with their respective expression activity during developmental processes. Using RNA-seq data from two distinct developmental stages of the mouse cortex, embryonic day 18 (E18) and postnatal day 7 (P7), we established for the first time a developmental-related transcriptome map of the mouse isochores. Additionally, for each stage we estimated the correlation between isochores' GC level and their expression activity, and the genes' expression patterns for each isochore family. Our analyses add evidence supporting the idea that during development GC-poor isochores are preferentially implicated, and confirm the positive effect of genes' GC level on their expression activity.

An unbiased approach to identify genes involved in development in a turtle with temperature-dependent sex determination

BACKGROUND

Many reptiles exhibit temperature-dependent sex determination (TSD). The initial cue in TSD is incubation temperature, unlike genotypic sex determination (GSD) where it is determined by the presence of specific alleles (or genetic loci). We used patterns of gene expression to identify candidates for genes with a role in TSD and other developmental processes without making a priori assumptions about the identity of these genes (ortholog-based approach). We identified genes with sexually dimorphic mRNA accumulation during the temperature sensitive period of development in the Red-eared slider turtle (Trachemys scripta), a turtle with TSD. Genes with differential mRNA accumulation in response to estrogen (estradiol-17β; E(2)) exposure and developmental stages were also identified.

RESULTS

Sequencing 767 clones from three suppression-subtractive hybridization libraries yielded a total of 581 unique sequences. Screening a macroarray with a subset of those sequences revealed a total of 26 genes that exhibited differential mRNA accumulation: 16 female biased and 10 male biased. Additional analyses revealed that C16ORF62 (an unknown gene) and MALAT1 (a long noncoding RNA) exhibited increased mRNA accumulation at the male producing temperature relative to the female producing temperature during embryonic sexual development. Finally, we identified four genes (C16ORF62, CCT3, MMP2, and NFIB) that exhibited a stage effect and five genes (C16ORF62, CCT3, MMP2, NFIB and NOTCH2) showed a response to E(2) exposure.

CONCLUSIONS

Here we report a survey of genes identified using patterns of mRNA accumulation during embryonic development in a turtle with TSD. Many previous studies have focused on examining the turtle orthologs of genes involved in mammalian development. Although valuable, the limitations of this approach are exemplified by our identification of two genes (MALAT1 and C16ORF62) that are sexually dimorphic during embryonic development. MALAT1 is a noncoding RNA that has not been implicated in sexual differentiation in other vertebrates and C16ORF62 has an unknown function. Our results revealed genes that are candidates for having roles in turtle embryonic development, including TSD, and highlight the need to expand our search parameters beyond protein-coding genes.

Transcriptome map of mouse isochores

Background

The availability of fully sequenced genomes and the implementation of transcriptome technologies have increased the studies investigating the expression profiles for a variety of tissues, conditions, and species. In this study, using RNA-seq data for three distinct tissues (brain, liver, and muscle), we investigate how base composition affects mammalian gene expression, an issue of prime practical and evolutionary interest.

Results

We present the transcriptome map of the mouse isochores (DNA segments with a fairly homogeneous base composition) for the three different tissues and the effects of isochores' base composition on their expression activity. Our analyses also cover the relations between the genes' expression activity and their localization in the isochore families.

Conclusions

This study is the first where next-generation sequencing data are used to associate the effects of both genomic and genic compositional properties to their corresponding expression activity. Our findings confirm previous results, and further support the existence of a relationship between isochores and gene expression. This relationship corroborates that isochores are primarily a product of evolutionary adaptation rather than a simple by-product of neutral evolutionary processes.

Pack-MULEs: Recycling and reshaping genes through GC-biased acquisition

The availability of genomic sequences provided new opportunities to decipher how plant genomes evolve. One recent discovery about plant genomes is the abundance of Pack-MULEs, a special group of transposable elements that duplicate, amplify and recombine gene fragments in many species at a very large scale. Despite the widespread occurrence of Pack-MULEs, their function remains an enigma. Our analysis using maize, rice and Arabidopsis genomic sequences indicates that the acquisition of genic sequences by Pack-MULEs is not random. Pack-MULEs in grasses specifically acquire and amplify GC-rich gene fragments. The resulting GC-rich elements have the ability to form independent transcripts with negative GC gradient, which refers to the decline of GC content along the orientation of transcription of genes. In other cases, Pack-MULEs insert near the 5' region of "normal" genes, and consequently form additional 5' exons or replace the original 5' exon of genes. In this manner, Pack-MULEs raise the GC content of the 5' termini of genes, modify the gene structure and contribute to the increased number of genes with negative GC gradient in grasses. The possible consequence of such activity is discussed.

Transcriptome sequencing of a highly salt tolerant mangrove species Sonneratia alba using Illumina platform

Mangroves are critical and threatened marine resources, yet few transcriptomic and genomic data are available in public databases. The transcriptome of a highly salt tolerant mangrove species, Sonneratia alba, was sequenced using the Illumina Genome Analyzer in this study. Over 15 million 75-bp paired-end reads were assembled into 30,628 unique sequences with an average length of 581 bp. Of them, 2358 SSRs were detected, with di-nucleotide repeats (59.2%) and tri-nucleotide repeats (37.7%) being the most common. Analysis of codon usage bias based on 20,945 coding sequences indicated that genes of S. alba were less biased than those of some microorganisms and Drosophila and that codon usage variation in S. alba was due primarily to compositional mutation bias, while translational selection has a relatively weak effect. Genome-wide gene ontology (GO) assignments showed that S. alba shared a similar GO slim classification with Arabidopsis thaliana. High percentages of sequences assigned to GO slim category 'mitochondrion' and four KEGG pathways, such as carbohydrates and secondary metabolites metabolism, may contribute to salt adaptation of S. alba. In addition, 1266 unique sequences matched to 273 known salt responsive genes (gene families) in other species were screened as candidates for salt tolerance of S. alba, and some of these genes showed fairly high coverage depth. At last, we identified four genes with signals of strong diversifying selection (K(a)/K(s)>1) by comparing the transcriptome sequences of S. alba with 249 known ESTs from its congener S. caseolaris. This study demonstrated a successful application of the Illumina platform to de novo assembly of the transcriptome of a non-model organism. Abundant SSR markers, salt responsive genes and four genes with signature of natural selection obtained from S. alba provide abundant sequence sources for future genetic diversity, salt adaptation and speciation studies.

GC3 biology in corn, rice, sorghum and other grasses

BACKGROUND

The third, or wobble, position in a codon provides a high degree of possible degeneracy and is an elegant fault-tolerance mechanism. Nucleotide biases between organisms at the wobble position have been documented and correlated with the abundances of the complementary tRNAs. We and others have noticed a bias for cytosine and guanine at the third position in a subset of transcripts within a single organism. The bias is present in some plant species and warm-blooded vertebrates but not in all plants, or in invertebrates or cold-blooded vertebrates.

RESULTS

Here we demonstrate that in certain organisms the amount of GC at the wobble position (GC3) can be used to distinguish two classes of genes. We highlight the following features of genes with high GC3 content: they (1) provide more targets for methylation, (2) exhibit more variable expression, (3) more frequently possess upstream TATA boxes, (4) are predominant in certain classes of genes (e.g., stress responsive genes) and (5) have a GC3 content that increases from 5'to 3'. These observations led us to formulate a hypothesis to explain GC3 bimodality in grasses.

CONCLUSIONS

Our findings suggest that high levels of GC3 typify a class of genes whose expression is regulated through DNA methylation or are a legacy of accelerated evolution through gene conversion. We discuss the three most probable explanations for GC3 bimodality: biased gene conversion, transcriptional and translational advantage and gene methylation.

Compositional genome contexts affect gene expression control in sea urchin embryo

Gene expression is widely perceived as exclusively controlled by the information contained in cis-regulatory regions. These are built in a modular way, each module being a cluster of binding sites for the transcription factors that control the level, the location and the time at which gene transcription takes place. On the other hand, results from our laboratory have shown that gene expression is affected by the compositional properties (GC levels) of the isochores in which genes are embedded, i.e. the genome context. To clarify how compositional genomic properties affect the way cis-regulatory information is utilized, we have changed the genome context of a GFP-reporter gene containing the complete cis-regulatory region of the gene spdeadringer (spdri), expressed during sea urchin embryogenesis. We have observed that GC levels higher or lower than those found in the natural genome context can alter the reporter expression pattern. We explain this as the result of an interference with the functionality of specific modules in the gene's cis-regulatory region. From these observations we derive the notion that the compositional properties of the genome context can affect cis-regulatory control of gene expression. Therefore although the way a gene works depends on the information contained in its cis-regulatory region, availability of such information depends on the compositional properties of the genomic context.