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.

Probing the different chaperone activities of the bacterial HSP70-HSP40 system using a thermolabile luciferase substrate.

During mild heat-stress, a native thermolabile polypeptide may partially unfold and transiently expose water-avoiding hydrophobic segments that readily tend to associate into a stable misfolded species, rich in intra-molecular non-native beta-sheet structures. When the concentration of the heat-unfolded intermediates is elevated, the exposed hydrophobic segments tend to associate with other molecules into large stable insoluble complexes, also called "aggregates." In mammalian cells, stress- and mutation-induced protein misfolding and aggregation may cause degenerative diseases and aging. Young cells, however, effectively counteract toxic protein misfolding with a potent network of molecular chaperones that bind hydrophobic surfaces and actively unfold otherwise stable misfolded and aggregated polypeptides. Here, we followed the behavior of a purified, initially mostly native thermolabile luciferase mutant, in the presence or absence of the Escherichia coli DnaK-DnaJ-GrpE chaperones and/or of ATP, at 22 °C or under mild heat-stress. We concomitantly measured luciferase enzymatic activity, Thioflavin-T fluorescence, and light-scattering to assess the effects of temperature and chaperones on the formation, respectively, of native, unfolded, misfolded, and/or of aggregated species. During mild heat-denaturation, DnaK-DnaJ-GrpE+ATP best maintained, although transiently, high luciferase activity and best prevented heat-induced misfolding and aggregation. In contrast, the ATP-less DnaK and DnaJ did not maintain optimal luciferase activity and were less effective at preventing luciferase misfolding and aggregation. We present a model accounting for the experimental data, where native, unfolded, misfolded, and aggregated species spontaneously inter-convert, and in which DnaK-DnaJ-GrpE+ATP specifically convert stable misfolded species into unstable unfolded intermediates.

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.

Active solubilization and refolding of stable protein aggregates by cooperative unfolding action of individual hsp70 chaperones.

Hsp70 is a central molecular chaperone that passively prevents protein aggregation and uses the energy of ATP hydrolysis to solubilize, translocate, and mediate the proper refolding of proteins in the cell. Yet, the molecular mechanism by which the active Hsp70 chaperone functions are achieved remains unclear. Here, we show that the bacterial Hsp70 (DnaK) can actively unfold misfolded structures in aggregated polypeptides, leading to gradual disaggregation. We found that the specific unfolding and disaggregation activities of individual DnaK molecules were optimal for large aggregates but dramatically decreased for small aggregates. The active unfolding of the smallest aggregates, leading to proper global refolding, required the cooperative action of several DnaK molecules per misfolded polypeptide. This finding suggests that the unique ATP-fueled locking/unlocking mechanism of the Hsp70 chaperones can recruit random chaperone motions to locally unfold misfolded structures and gradually disentangle stable aggregates into refoldable proteins.

Hsp110 is a bona fide chaperone using ATP to unfold stable misfolded polypeptides and reciprocally collaborate with hsp70 to solubilize protein aggregates.

Structurally and sequence-wise, the Hsp110s belong to a subfamily of the Hsp70 chaperones. Like the classical Hsp70s, members of the Hsp110 subfamily can bind misfolding polypeptides and hydrolyze ATP. However, they apparently act as a mere subordinate nucleotide exchange factors, regulating the ability of Hsp70 to hydrolyze ATP and convert stable protein aggregates into native proteins. Using stably misfolded and aggregated polypeptides as substrates in optimized in vitro chaperone assays, we show that the human cytosolic Hsp110s (HSPH1 and HSPH2) are bona fide chaperones on their own that collaborate with Hsp40 (DNAJA1 and DNAJB1) to hydrolyze ATP and unfold and thus convert stable misfolded polypeptides into natively refolded proteins. Moreover, equimolar Hsp70 (HSPA1A) and Hsp110 (HSPH1) formed a powerful molecular machinery that optimally reactivated stable luciferase aggregates in an ATP- and DNAJA1-dependent manner, in a disaggregation mechanism whereby the two paralogous chaperones alternatively activate the release of bound unfolded polypeptide substrates from one another, leading to native protein refolding.

Reactivation of protein aggregates by mortalin and Tid1-the human mitochondrial Hsp70 chaperone system.

The mitochondrial 70-kDa heat shock protein (mtHsp70), also known in humans as mortalin, is a central component of the mitochondrial protein import motor and plays a key role in the folding of matrix-localized mitochondrial proteins. MtHsp70 is assisted by a member of the 40-kDa heat shock protein co-chaperone family named Tid1 and a nucleotide exchange factor. Whereas, yeast mtHsp70 has been extensively studied in the context of protein import in the mitochondria, and the bacterial 70-kDa heat shock protein was recently shown to act as an ATP-fuelled unfolding enzyme capable of detoxifying stably misfolded polypeptides into harmless natively refolded proteins, little is known about the molecular functions of the human mortalin in protein homeostasis. Here, we developed novel and efficient purification protocols for mortalin and the two spliced versions of Tid1, Tid1-S, and Tid1-L and showed that mortalin can mediate the in vitro ATP-dependent reactivation of stable-preformed heat-denatured model aggregates, with the assistance of Mge1 and either Tid1-L or Tid1-S co-chaperones or yeast Mdj1. Thus, in addition of being a central component of the protein import machinery, human mortalin together with Tid1, may serve as a protein disaggregating machine which, for lack of Hsp100/ClpB disaggregating co-chaperones, may carry alone the scavenging of toxic protein aggregates in stressed, diseased, or aging human mitochondria.

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.