Evidence for transcript networks composed of chimeric RNAs in human cells.

The classic organization of a gene structure has followed the Jacob and Monod bacterial gene model proposed more than 50 years ago. Since then, empirical determinations of the complexity of the transcriptomes found in yeast to human has blurred the definition and physical boundaries of genes. Using multiple analysis approaches we have characterized individual gene boundaries mapping on human chromosomes 21 and 22. Analyses of the locations of the 5' and 3' transcriptional termini of 492 protein coding genes revealed that for 85% of these genes the boundaries extend beyond the current annotated termini, most often connecting with exons of transcripts from other well annotated genes. The biological and evolutionary importance of these chimeric transcripts is underscored by (1) the non-random interconnections of genes involved, (2) the greater phylogenetic depth of the genes involved in many chimeric interactions, (3) the coordination of the expression of connected genes and (4) the close in vivo and three dimensional proximity of the genomic regions being transcribed and contributing to parts of the chimeric RNAs. The non-random nature of the connection of the genes involved suggest that chimeric transcripts should not be studied in isolation, but together, as an RNA network.

Versatile gene-specific sequence tags for Arabidopsis functional genomics: transcript profiling and reverse genetics applications.

Microarray transcript profiling and RNA interference are two new technologies crucial for large-scale gene function studies in multicellular eukaryotes. Both rely on sequence-specific hybridization between complementary nucleic acid strands, inciting us to create a collection of gene-specific sequence tags (GSTs) representing at least 21,500 Arabidopsis genes and which are compatible with both approaches. The GSTs were carefully selected to ensure that each of them shared no significant similarity with any other region in the Arabidopsis genome. They were synthesized by PCR amplification from genomic DNA. Spotted microarrays fabricated from the GSTs show good dynamic range, specificity, and sensitivity in transcript profiling experiments. The GSTs have also been transferred to bacterial plasmid vectors via recombinational cloning protocols. These cloned GSTs constitute the ideal starting point for a variety of functional approaches, including reverse genetics. We have subcloned GSTs on a large scale into vectors designed for gene silencing in plant cells. We show that in planta expression of GST hairpin RNA results in the expected phenotypes in silenced Arabidopsis lines. These versatile GST resources provide novel and powerful tools for functional genomics.

Novel testis germ cell-specific transcript of the CREB gene contains an alternatively spliced exon with multiple in-frame stop codons.

The gene encoding the cAMP-responsive transcription factor CREB consists of multiple small exons some of which undergo alternative RNA splicing. We describe the finding of a novel transcript of the CREB gene expressed at high levels in the germ cells of the rat testis. The transcript contains an alternatively spliced exon inserted within the sequence encoding the transcriptional transactivation domain of CREB and this exon contains multiple in-frame stop codons. Furthermore, the exon is conserved in both rat and human genes (75% nucleotide identity). Although the function(s) of this RNA or the truncated CREB protein predicted to result from the translation of this unusual transcript is unknown, the high level of expression in the testicular germ cells and remarkable conservation of sequences in rat and human suggests that it may have a unique biological function in these cells.

CATMA, a comprehensive genome-scale resource for silencing and transcript profiling of Arabidopsis genes.

BACKGROUND: The Complete Arabidopsis Transcript MicroArray (CATMA) initiative combines the efforts of laboratories in eight European countries 1 to deliver gene-specific sequence tags (GSTs) for the Arabidopsis research community. The CATMA initiative offers the power and flexibility to regularly update the GST collection according to evolving knowledge about the gene repertoire. These GST amplicons can easily be reamplified and shared, subsets can be picked at will to print dedicated arrays, and the GSTs can be cloned and used for other functional studies. This ongoing initiative has already produced approximately 24,000 GSTs that have been made publicly available for spotted microarray printing and RNA interference. RESULTS: GSTs from the CATMA version 2 repertoire (CATMAv2, created in 2002) were mapped onto the gene models from two independent Arabidopsis nuclear genome annotation efforts, TIGR5 and PSB-EuGène, to consolidate a list of genes that were targeted by previously designed CATMA tags. A total of 9,027 gene models were not tagged by any amplified CATMAv2 GST, and 2,533 amplified GSTs were no longer predicted to tag an updated gene model. To validate the efficacy of GST mapping criteria and design rules, the predicted and experimentally observed hybridization characteristics associated to GST features were correlated in transcript profiling datasets obtained with the CATMAv2 microarray, confirming the reliability of this platform. To complete the CATMA repertoire, all 9,027 gene models for which no GST had yet been designed were processed with an adjusted version of the Specific Primer and Amplicon Design Software (SPADS). A total of 5,756 novel GSTs were designed and amplified by PCR from genomic DNA. Together with the pre-existing GST collection, this new addition constitutes the CATMAv3 repertoire. It comprises 30,343 unique amplified sequences that tag 24,202 and 23,009 protein-encoding nuclear gene models in the TAIR6 and EuGène genome annotations, respectively. To cover the remaining untagged genes, we identified 543 additional GSTs using less stringent design criteria and designed 990 sequence tags matching multiple members of gene families (Gene Family Tags or GFTs) to cover any remaining untagged genes. These latter 1,533 features constitute the CATMAv4 addition. CONCLUSION: To update the CATMA GST repertoire, we designed 7,289 additional sequence tags, bringing the total number of tagged TAIR6-annotated Arabidopsis nuclear protein-coding genes to 26,173. This resource is used both for the production of spotted microarrays and the large-scale cloning of hairpin RNA silencing vectors. All information about the resulting updated CATMA repertoire is available through the CATMA database http://www.catma.org.

Affymetrix Whole-Transcript Human Gene 1.0 ST array is highly concordant with standard 3' expression arrays.

The recently released Affymetrix Human Gene 1.0 ST array has two major differences compared with standard 3' based arrays: (i) it interrogates the entire mRNA transcript, and (ii) it uses DNA targets. To assess the impact of these differences on array performance, we performed a series of comparative hybridizations between the Human Gene 1.0 ST and the Affymetrix HG-U133 Plus 2.0 and the Illumina HumanRef-8 BeadChip arrays. Additionally, both RNA and DNA targets were hybridized on HG-U133 Plus 2.0 arrays. The results show that the overall reproducibility of the Gene 1.0 ST array is best. When looking only at the high intensity probes, the reproducibility of the Gene 1.0 ST array and the Illumina BeadChip array is equally good. Concordance of array results was assessed using different inter-platform mappings. Agreements are best between the two labeling protocols using HG-U133 Plus 2.0 array. The Gene 1.0 ST array is most concordant with the HG-U133 array hybridized with cDNA targets. This may reflect the impact of the target type. Overall, the high degree of correspondence provides strong evidence for the reliability of the Gene 1.0 ST array.

Eeny meeny miny moe, catch a transcript by the toe, or how to enumerate eukaryotic transcripts.

In this issue of Genes & Development, Revyakin and colleagues (pp. 1691-1702) measure the relation between individual RNA polymerase II transcription events and transcription factor assembly by counting RNA transcripts retained on the template DNA using single-molecule fluorescence.

Analysis of epidermis- and mesophyll-specific transcript accumulation in powdery mildew-inoculated wheat leaves.

Powdery mildew is an important disease of wheat caused by the obligate biotrophic fungus Blumeria graminis f. sp. tritici. This pathogen invades exclusively epidermal cells after penetrating directly through the cell wall. Because powdery mildew colonizes exclusively epidermal cells, it is of importance not only to identify genes which are activated, but also to monitor tissue specificity of gene activation. Acquired resistance of wheat to powdery mildew can be induced by a previous inoculation with the non-host pathogen B. graminis f. sp. hordei, the causal agent of barley powdery mildew. The establishment of the resistant state is accompanied by the activation of genes. Here we report the tissue-specific cDNA-AFLP analysis and cloning of transcripts accumulating 6 and 24 h after the resistance-inducing inoculation with B. graminis f. sp. hordei. A total of 25,000 fragments estimated to represent about 17,000 transcripts were displayed. Out of these, 141 transcripts, were found to accumulate after Bgh inoculation using microarray hybridization analysis. Forty-four accumulated predominantly in the epidermis whereas 76 transcripts accumulated mostly in mesophyll tissue.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.

A novel vasopressin-induced transcript promotes MAP kinase activation and ENaC downregulation.

In the principal cell of the renal collecting duct, vasopressin regulates the expression of a gene network responsible for sodium and water reabsorption through the regulation of the water channel and the epithelial sodium channel (ENaC). We have recently identified a novel vasopressin-induced transcript (VIT32) that encodes for a 142 amino acid vasopressin-induced protein (VIP32), which has no homology with any protein of known function. The Xenopus oocyte expression system revealed two functions: (i) when injected alone, VIT32 cRNA rapidly induces oocyte meiotic maturation through the activation of the maturation promoting factor, the amphibian homolog of the universal M phase trigger Cdc2/cyclin; and (ii) when co-injected with the ENaC, VIT32 cRNA selectively downregulates channel activity, but not channel cell surface expression. In the kidney principal cell, VIP32 may be involved in the downregulation of transepithelial sodium transport observed within a few hours after vasopressin treatment. VIP32 belongs to a novel gene family ubiquitously expressed in oocyte and somatic cells that may be involved in G to M transition and cell cycling.

Transcript profiling suggests that differential metabolic adaptation of mice to a high fat diet is associated with changes in liver to muscle lipid fluxes.

Genetically homogenous C57Bl/6 mice display differential metabolic adaptation when fed a high fat diet for 9 months. Most become obese and diabetic, but a significant fraction remains lean and diabetic or lean and non-diabetic. Here, we performed microarray analysis of "metabolic" transcripts expressed in liver and hindlimb muscles to evaluate: (i) whether expressed transcript patterns could indicate changes in metabolic pathways associated with the different phenotypes, (ii) how these changes differed from the early metabolic adaptation to short term high fat feeding, and (iii) whether gene classifiers could be established that were characteristic of each metabolic phenotype. Our data indicate that obesity/diabetes was associated with preserved hepatic lipogenic gene expression and increased plasma levels of very low density lipoprotein and, in muscle, with an increase in lipoprotein lipase gene expression. This suggests increased muscle fatty acid uptake, which may favor insulin resistance. In contrast, the lean mice showed a strong reduction in the expression of hepatic lipogenic genes, in particular of Scd-1, a gene linked to sensitivity to diet-induced obesity; the lean and non-diabetic mice presented an additional increased expression of eNos in liver. After 1 week of high fat feeding the liver gene expression pattern was distinct from that seen at 9 months in any of the three mouse groups, thus indicating progressive establishment of the different phenotypes. Strikingly, development of the obese phenotype involved re-expression of Scd-1 and other lipogenic genes. Finally, gene classifiers could be established that were characteristic of each metabolic phenotype. Together, these data suggest that epigenetic mechanisms influence gene expression patterns and metabolic fates.

Loss-of-function mutations of an inhibitory upstream ORF in the human hairless transcript cause Marie Unna hereditary hypotrichosis.

Marie Unna hereditary hypotrichosis (MUHH) is an autosomal dominant form of genetic hair loss. In a large Chinese family carrying MUHH, we identified a pathogenic initiation codon mutation in U2HR, an inhibitory upstream ORF in the 5' UTR of the gene encoding the human hairless homolog (HR). U2HR is predicted to encode a 34-amino acid peptide that is highly conserved among mammals. In 18 more families from different ancestral groups, we identified a range of defects in U2HR, including loss of initiation, delayed termination codon and nonsense and missense mutations. Functional analysis showed that these classes of mutations all resulted in increased translation of the main HR physiological ORF. Our results establish the link between MUHH and U2HR, show that fine-tuning of HR protein levels is important in control of hair growth, and identify a potential mechanism for preventing hair loss or promoting hair removal.

Low E2F1 transcript levels are a strong determinant of favorable breast cancer outcome.

INTRODUCTION: We investigated whether mRNA levels of E2F1, a key transcription factor involved in proliferation, differentiation and apoptosis, could be used as a surrogate marker for the determination of breast cancer outcome. METHODS: E2F1 and other proliferation markers were measured by quantitative RT-PCR in 317 primary breast cancer patients from the Stiftung Tumorbank Basel. Correlations to one another as well as to the estrogen receptor and ERBB2 status and clinical outcome were investigated. Results were validated and further compared with expression-based prognostic profiles using The Netherlands Cancer Institute microarray data set reported by Fan and colleagues. RESULTS: E2F1 mRNA expression levels correlated strongly with the expression of other proliferation markers, and low values were mainly found in estrogen receptor-positive and ERBB2-negative phenotypes. Patients with low E2F1-expressing tumors were associated with favorable outcome (hazard ratio = 4.3 (95% confidence interval = 1.8-9.9), P = 0.001). These results were consistent in univariate and multivariate Cox analyses, and were successfully validated in The Netherlands Cancer Institute data set. Furthermore, E2F1 expression levels correlated well with the 70-gene signature displaying the ability of selecting a common subset of patients at good prognosis. Breast cancer patients' outcome was comparably predictable by E2F1 levels, by the 70-gene signature, by the intrinsic subtype gene classification, by the wound response signature and by the recurrence score. CONCLUSION: Assessment of E2F1 at the mRNA level in primary breast cancer is a strong determinant of breast cancer patient outcome. E2F1 expression identified patients at low risk of metastasis irrespective of the estrogen receptor and ERBB2 status, and demonstrated similar prognostic performance to different gene expression-based predictors.

Spatio-temporal transcript profiling of rice roots and shoots in response to phosphate starvation and recovery.

Using rice (Oryza sativa) as a model crop species, we performed an in-depth temporal transcriptome analysis, covering the early and late stages of Pi deprivation as well as Pi recovery in roots and shoots, using next-generation sequencing. Analyses of 126 paired-end RNA sequencing libraries, spanning nine time points, provided a comprehensive overview of the dynamic responses of rice to Pi stress. Differentially expressed genes were grouped into eight sets based on their responses to Pi starvation and recovery, enabling the complex signaling pathways involved in Pi homeostasis to be untangled. A reference annotation-based transcript assembly was also generated, identifying 438 unannotated loci that were differentially expressed under Pi starvation. Several genes also showed induction of unannotated splice isoforms under Pi starvation. Among these, PHOSPHATE2 (PHO2), a key regulator of Pi homeostasis, displayed a Pi starvation-induced isoform, which was associated with increased translation activity. In addition, microRNA (miRNA) expression profiles after long-term Pi starvation in roots and shoots were assessed, identifying 20 miRNA families that were not previously associated with Pi starvation, such as miR6250. In this article, we present a comprehensive spatio-temporal transcriptome analysis of plant responses to Pi stress, revealing a large number of potential key regulators of Pi homeostasis in plants.

A conserved transcript pattern in response to a specialist and a generalist herbivore.

Transcript patterns elicited in response to attack reveal, at the molecular level, how plants respond to aggressors. These patterns are fashioned both by inflicted physical damage as well as by biological components displayed or released by the attacker. Different types of attacking organisms might therefore be expected to elicit different transcription programs in the host. Using a large-scale DNA microarray, we characterized gene expression in damaged as well as in distal Arabidopsis thaliana leaves in response to the specialist insect, Pieris rapae. More than 100 insect-responsive genes potentially involved in defense were identified, including genes involved in pathogenesis, indole glucosinolate metabolism, detoxification and cell survival, and signal transduction. Of these 114 genes, 111 were induced in Pieris feeding, and only three were repressed. Expression patterns in distal leaves were markedly similar to those of local leaves. Analysis of wild-type and jasmonate mutant plants, coupled with jasmonate treatment, showed that between 67 and 84% of Pieris-regulated gene expression was controlled, totally or in part, by the jasmonate pathway. This was correlated with increased larval performance on the coronatine insensitive1 glabrous1 (coi1-1 gl1) mutant. Independent mutations in COI1 and GL1 led to a faster larval weight gain, but the gl1 mutation had relatively little effect on the expression of the insect-responsive genes examined. Finally, we compared transcript patterns in Arabidopis in response to larvae of the specialist P. rapae and to a generalist insect, Spodoptera littoralis. Surprisingly, given the complex nature of insect salivary components and reported differences between species, almost identical transcript profiles were observed. This study also provides a robustly characterized gene set for the further investigation of plant-insect interaction.

Novel targets of the CbrAB/Crc carbon catabolite control system revealed by transcript abundance in Pseudomonas aeruginosa.

The opportunistic human pathogen Pseudomonas aeruginosa is able to utilize a wide range of carbon and nitrogen compounds, allowing it to grow in vastly different environments. The uptake and catabolism of growth substrates are organized hierarchically by a mechanism termed catabolite repression control (Crc) whereby the Crc protein establishes translational repression of target mRNAs at CA (catabolite activity) motifs present in target mRNAs near ribosome binding sites. Poor carbon sources lead to activation of the CbrAB two-component system, which induces transcription of the small RNA (sRNA) CrcZ. This sRNA relieves Crc-mediated repression of target mRNAs. In this study, we have identified novel targets of the CbrAB/Crc system in P. aeruginosa using transcriptome analysis in combination with a search for CA motifs. We characterized four target genes involved in the uptake and utilization of less preferred carbon sources: estA (secreted esterase), acsA (acetyl-CoA synthetase), bkdR (regulator of branched-chain amino acid catabolism) and aroP2 (aromatic amino acid uptake protein). Evidence for regulation by CbrAB, CrcZ and Crc was obtained in vivo using appropriate reporter fusions, in which mutation of the CA motif resulted in loss of catabolite repression. CbrB and CrcZ were important for growth of P. aeruginosa in cystic fibrosis (CF) sputum medium, suggesting that the CbrAB/Crc system may act as an important regulator during chronic infection of the CF lung.