Tandem chimeric transcription as a means to increase protein diversity in the human genome

The “one-gene, one-protein” rule, coined by Beadle and Tatum, has been fundamental to molecular biology. The rule implies that the genetic complexity of an organism depends essentially on its gene number. The discovery, however, that alternative gene splicing and transcription are widespread phenomena dramatically altered our understanding of the genetic complexity of higher eukaryotic organisms; in these, a limited number of genes may potentially encode a much larger number of proteins. Here we investigate yet another phenomenon that may contribute to generate additional protein diversity. Indeed, by relying on both computational and experimental analysis, we estimate that at least 4%–5% of the tandem gene pairs in the human genome can be eventually transcribed into a single RNA sequence encoding a putative chimeric protein. While the functional significance of most of these chimeric transcripts remains to be determined, we provide strong evidence that this phenomenon does not correspond to mere technical artifacts and that it is a common mechanism with the potential of generating hundreds of additional proteins in the human genome.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

GENCODE: producing a reference annotation for ENCODE

Background: The GENCODE consortium was formed to identify and map all protein-coding genes within the ENCODE regions. This was achieved by a combination of initial manualannotation by the HAVANA team, experimental validation by the GENCODE consortium and a refinement of the annotation based on these experimental results.Results: The GENCODE gene features are divided into eight different categories of which onlythe first two (known and novel coding sequence) are confidently predicted to be protein-codinggenes. 5’ rapid amplification of cDNA ends (RACE) and RT-PCR were used to experimentallyverify the initial annotation. Of the 420 coding loci tested, 229 RACE products have beensequenced. They supported 5’ extensions of 30 loci and new splice variants in 50 loci. In addition,46 loci without evidence for a coding sequence were validated, consisting of 31 novel and 15putative transcripts. We assessed the comprehensiveness of the GENCODE annotation byattempting to validate all the predicted exon boundaries outside the GENCODE annotation. Outof 1,215 tested in a subset of the ENCODE regions, 14 novel exon pairs were validated, only twoof them in intergenic regions.Conclusions: In total, 487 loci, of which 434 are coding, have been annotated as part of theGENCODE reference set available from the UCSC browser. Comparison of GENCODEannotation with RefSeq and ENSEMBL show only 40% of GENCODE exons are contained withinthe two sets, which is a reflection of the high number of alternative splice forms with uniqueexons annotated. Over 50% of coding loci have been experimentally verified by 5’ RACE forEGASP and the GENCODE collaboration is continuing to refine its annotation of 1% humangenome with the aid of experimental validation.

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.

Differentiated evolutionary rates in alternative exons and the implications for splicing regulation

Background: Alternatively spliced exons play an important role in the diversification of gene function in most metazoans and are highly regulated by conserved motifs in exons and introns. Two contradicting properties have been associated to evolutionary conserved alternative exons: higher sequence conservation and higher rate of non-synonymous substitutions, relative to constitutive exons. In order to clarify this issue, we have performed an analysis of the evolution of alternative and constitutive exons, using a large set of protein coding exons conserved between human and mouse and taking into account the conservation of the transcript exonic structure. Further, we have also defined a measure of the variation of the arrangement of exonic splicing enhancers (ESE-conservation score) to study the evolution of splicing regulatory sequences. We have used this measure to correlate the changes in the arrangement of ESEs with the divergence of exon and intron sequences. Results: We find evidence for a relation between the lack of conservation of the exonic structure and the weakening of the sequence evolutionary constraints in alternative and constitutive exons. Exons in transcripts with non-conserved exonic structures have higher synonymous (dS) and non-synonymous (dN) substitution rates than exons in conserved structures. Moreover, alternative exons in transcripts with non-conserved exonic structure are the least constrained in sequence evolution, and at high EST-inclusion levels they are found to be very similar to constitutive exons, whereas alternative exons in transcripts with conserved exonic structure have a dS significantly lower than average at all EST-inclusion levels. We also find higher conservation in the arrangement of ESEs in constitutive exons compared to alternative ones. Additionally, the sequence conservation at flanking introns remains constant for constitutive exons at all ESE-conservation values, but increases for alternative exons at high ESE-conservation values. Conclusion: We conclude that most of the differences in dN observed between alternative and constitutive exons can be explained by the conservation of the transcript exonic structure. Low dS values are more characteristic of alternative exons with conserved exonic structure, but not of those with non-conserved exonic structure. Additionally, constitutive exons are characterized by a higher conservation in the arrangement of ESEs, and alternative exons with an ESE-conservation similar to that of constitutive exons are characterized by a conservation of the flanking intron sequences higher than average, indicating the presence of more intronic regulatory signals.

Evaluation of the chicken transcriptome by SAGE of B cells and the DT40 cell line

Background: The understanding of whole genome sequences in higher eukaryotes depends to a large degree on the reliable definition of transcription units including exon/intron structures, translated open reading frames (ORFs) and flanking untranslated regions. The best currently available chicken transcript catalog is the Ensembl build based on the mappings of a relatively small number of full length cDNAs and ESTs to the genome as well as genome sequence derived in silico gene predictions.Results: We use Long Serial Analysis of Gene Expression (LongSAGE) in bursal lymphocytes and the DT40 cell line to verify the quality and completeness of the annotated transcripts. 53.6% of the more than 38,000 unique SAGE tags (unitags) match to full length bursal cDNAs, the Ensembl transcript build or the genome sequence. The majority of all matching unitags show single matches to the genome, but no matches to the genome derived Ensembl transcript build. Nevertheless, most of these tags map close to the 3' boundaries of annotated Ensembl transcripts.Conclusions: These results suggests that rather few genes are missing in the current Ensembl chicken transcript build, but that the 3' ends of many transcripts may not have been accurately predicted. The tags with no match in the transcript sequences can now be used to improve gene predictions, pinpoint the genomic location of entirely missed transcripts and optimize the accuracy of gene finder software.

GENCODE: the reference human genome annotation for The ENCODE Project.

The GENCODE Consortium aims to identify all gene features in the human genome using a combination of computational analysis, manual annotation, and experimental validation. Since the first public release of this annotation data set, few new protein-coding loci have been added, yet the number of alternative splicing transcripts annotated has steadily increased. The GENCODE 7 release contains 20,687 protein-coding and 9640 long noncoding RNA loci and has 33,977 coding transcripts not represented in UCSC genes and RefSeq. It also has the most comprehensive annotation of long noncoding RNA (lncRNA) loci publicly available with the predominant transcript form consisting of two exons. We have examined the completeness of the transcript annotation and found that 35% of transcriptional start sites are supported by CAGE clusters and 62% of protein-coding genes have annotated polyA sites. Over one-third of GENCODE protein-coding genes are supported by peptide hits derived from mass spectrometry spectra submitted to Peptide Atlas. New models derived from the Illumina Body Map 2.0 RNA-seq data identify 3689 new loci not currently in GENCODE, of which 3127 consist of two exon models indicating that they are possibly unannotated long noncoding loci. GENCODE 7 is publicly available from gencodegenes.org and via the Ensembl and UCSC Genome Browsers.

Differential Gene Expression between Fall- and Spring-Run Chinook Salmon Assessed by Long Serial Analysis of Gene Expression

Of all Pacific salmonids, Chinook salmon Oncorhynchus tshawytscha display the greatest variability in return times to freshwater. The molecular mechanisms of these differential return times have not been well described. Current methods, such as long serial analysis of gene expression (LongSAGE) and microarrays, allow gene expression to be analyzed for thousands of genes simultaneously. To investigate whether differential gene expression is observed between fall- and spring-run Chinook salmon from California's Central Valley, LongSAGE libraries were constructed. Three libraries containing between 25,512 and 29,372 sequenced tags (21 base pairs/tag) were generated using messenger RNA from the brains of adult Chinook salmon returning in fall and spring and from one ocean-caught Chinook salmon. Tags were annotated to genes using complementary DNA libraries from Atlantic salmon Salmo salar and rainbow trout O. mykiss. Differentially expressed genes, as estimated by differences in the number of sequence tags, were found in all pairwise comparisons of libraries (freshwater versus saltwater = 40 genes; fall versus spring = 11 genes: and spawning versus nonspawning = 51 genes). The gene for ependymin, an extracellular glycoprotein involved in behavioral plasticity in fish, exhibited the most differential expression among the three groupings. Reverse transcription polymerase chain reaction analysis verified the differential expression of ependymin between the fall- and spring-run samples. These LongSAGE libraries, the first reported for Chinook salmon, provide a window of the transcriptional changes during Chinook salmon return migration to freshwater and spawning and increase the amount of expressed sequence data.

The Drosophila cdc25 homolog twine is required for meiosis.

We have identified a second cdc25 homolog in Drosophila. In contrast to string (the first homolog identified in Drosophila) this second homolog, twine, does not function in the mitotic cell cycle, but is specialized for meiosis. Expression of twine was observed exclusively in male and female gonads. twine transcripts are present in germ cells during meiosis, and appear only late during gametogenesis, well after the end of the mitotic germ cell divisions. The sterile Drosophila mutant, mat(2)synHB5, which had previously been isolated and mapped to the same genomic region as twine (35F), was found to carry a missense mutation in the twine gene. This missense mutation in twine abolished its ability to complement a mutation in Schizosaccharomyces pombe cdc25. Phenotypic analysis of mat(2)synHB5 mutant flies revealed a complete block of meiosis in males and severe meiotic defects in females.

Transcriptional and functional profiles of voltage-gated Na(+) channels in injured and non-injured DRG neurons in the SNI model of neuropathic pain.

Changes in expression and function of voltage-gated sodium channels (VGSC) in dorsal root ganglion (DRG) neurons may play a major role in the genesis of peripheral hyperexcitability that occurs in neuropathic pain. We present here the first description of changes induced by spared nerve injury (SNI) to Na(v)1 mRNA levels and tetrodotoxin-sensitive and -resistant (TTX-S/TTX-R) Na(+) currents in injured and adjacent non-injured small DRG neurons. VGSC transcripts were down-regulated in injured neurons except for Na(v)1.3, which increased, while they were either unchanged or increased in non-injured neurons. TTX-R current densities were reduced in injured neurons and the voltage dependence of steady-state inactivation for TTX-R was positively shifted in injured and non-injured neurons. TTX-S current densities were not affected by SNI, while the rate of recovery from inactivation was accelerated in injured neurons. Our results describe altered neuronal electrogenesis following SNI that is likely induced by a complex regulation of VGSCs.

CNOT3 is a modifier of PRPF31 mutations in retinitis pigmentosa with incomplete penetrance.

Heterozygous mutations in the PRPF31 gene cause autosomal dominant retinitis pigmentosa (adRP), a hereditary disorder leading to progressive blindness. In some cases, such mutations display incomplete penetrance, implying that certain carriers develop retinal degeneration while others have no symptoms at all. Asymptomatic carriers are protected from the disease by a higher than average expression of the PRPF31 allele that is not mutated, mainly through the action of an unknown modifier gene mapping to chromosome 19q13.4. We investigated a large family with adRP segregating an 11-bp deletion in PRPF31. The analysis of cell lines derived from asymptomatic and affected individuals revealed that the expression of only one gene among a number of candidates within the 19q13.4 interval significantly correlated with that of PRPF31, both at the mRNA and protein levels, and according to an inverse relationship. This gene was CNOT3, encoding a subunit of the Ccr4-not transcription complex. In cultured cells, siRNA-mediated silencing of CNOT3 provoked an increase in PRPF31 expression, confirming a repressive nature of CNOT3 on PRPF31. Furthermore, chromatin immunoprecipitation revealed that CNOT3 directly binds to a specific PRPF31 promoter sequence, while next-generation sequencing of the CNOT3 genomic region indicated that its variable expression is associated with a common intronic SNP. In conclusion, we identify CNOT3 as the main modifier gene determining penetrance of PRPF31 mutations, via a mechanism of transcriptional repression. In asymptomatic carriers CNOT3 is expressed at low levels, allowing higher amounts of wild-type PRPF31 transcripts to be produced and preventing manifestation of retinal degeneration.

Identification and characterization of the members of the pho1 gene family in moss and rice

Inorganic phosphate (Pi) is one of the main nutrients limiting plant growth anddevelopment in many agro-ecosystems. In plants, phosphate is acquired from the soil by theroots, and is then transferred to the shoot via the xylem. In the model plant Arabidopsisthaliana, PHO1 was previously identified as being involved in loading Pi into the xylem ofroots. AtPHO1, belongs to a multigenic family composed of 10 additional members, namelyAtPHO1;H1 to AtPHO1;10. In this study, we aimed at further investigating the role of thePHO1 gene family in Pi homeostasis in plants, and to this end we isolated and characterizedthe PHO1 members of two main model plants, the moss Physcomitrella patens and the riceOryza sativa.In the bryophyte P. patens, bioinformatic analyses revealed the presence of seven AtPHO1homologues, highly similar to AtPHO1. The seven moss PHO1 genes, namely PpPHO1;1 toPpPHO1;7 appeared to be differentially regulated, both at the tissue level and in response toPi status. However only PpPHO1;1 and PpPHO1;7 were specifically up-regulated upon Pistarvation, suggesting a potential role in Pi homeostasis. We also characterized the responseof P. patens to Pi starvation, showing that higher and lower plants share some commonstrategies to adapt to Pi-deficiency.In the second part, focusing on the monocotyledon rice, we showed the existence of threePHO1 homologues OsPHO1;1 to OsPHO1;3, with the unique particularity of each havingNatural Antisense Transcripts (NATs). Molecular analyses revealed that both the sense andthe antisense OsPHO1;2 transcripts were by far the most abundantly expressed transcripts ofthe family, preferentially expressed in the roots. The stable expression of OsPHO1;2 in allconditions tested, in opposition with the highly induced antisense transcript upon Pistarvation, suggest a putative role for the antisense in regulating the sense transcript.Moreover, mutant analyses revealed that OsPHO1;2 plays a key role in Pi homeostasis, intransferring Pi from the root to the shoot. Finally, complementing the pho1 mutant inArabidopsis, characterized by low Pi in the shoot and reduced growth, with the riceOsPHO1;2 gene revealed a new role for PHO1 in Pi signaling. Indeed, the complementedplants showed normal growth, with however low Pi content.

Developmental stage-specific expression of cyclic adenosine 3',5'-monophosphate response element-binding protein CREB during spermatogenesis involves alternative exon splicing.

Spermatogenesis is a temporally regulated developmental process by which the gonadotropin-responsive somatic Sertoli and Leydig cells act interdependently to direct the maturation of the germinal cells. The metabolism of Sertoli and Leydig cells is regulated by the pituitary gonadotropins FSH and LH, which, in turn, activate adenylate cyclase. Because the cAMP-second messenger pathway is activated by FSH and LH, we postulated that the cAMP-responsive element-binding protein (CREB) plays a physiological role in Sertoli and Leydig cells, respectively. Immunocytochemical analyses of rat testicular sections show a remarkably high expression of CREB in the haploid round spermatids and, to some extent, in pachytene spermatocytes and Sertoli cells. Although most of the CREB antigen is detected in the nuclei, some CREB antigen is also present in the cytoplasm. Remarkably, the cytoplasmic CREB results from the translation of a unique alternatively spliced transcript of the CREB gene that incorporates an exon containing multiple stop codons inserted immediately up-stream of the exons encoding the DNA-binding domain of CREB. Thus, the RNA containing the alternatively spliced exon encodes a truncated transcriptional transactivator protein lacking both the DNA-binding domain and nuclear translocation signal of CREB. Most of the CREB transcripts detected in the germinal cells contain the alternatively spliced exon, suggesting a function of the exon to modulate the synthesis of CREB. In the Sertoli cells we observed a striking cyclical (12-day periodicity) increase in the levels of CREB mRNA that coincides with the splicing out of the restrictive exon containing the stop codons. Because earlier studies established that FSH-stimulated cAMP levels in Sertoli cells are also cyclical, and the CREB gene promoter contains cAMP-responsive enhancers, we suggest that the alternative RNA splicing controls a positive autoregulation of CREB gene expression mediated by cAMP.

Velocity estimates for signal propagation leading to systemic jasmonic acid accumulation in wounded Arabidopsis.

The wound response prohormone jasmonic acid (JA) accumulates rapidly in tissues both proximal and distal to injury sites in plants. Using quantitative liquid chromatography-mass spectrometry after flash freezing of tissues, we found that JA accumulated within 30 s of injury in wounded Arabidopsis leaves (p = 3.5 e(-7)). JA augmentation distal to wounds was strongest in unwounded leaves with direct vascular connections to wounded leaves wherein JA levels increased significantly within 120 s of wounding (p = 0.00027). This gave conservative and statistically robust temporal boundaries for the average velocity of the long distance signal leading to distal JA accumulation in unwounded leaves of 3.4-4.5 cm min(-1). Like JA, transcripts of the JA synthesis gene LIPOXYGENASE2 (LOX2) and the jasmonate response gene JAZ10.3 also accumulated to higher levels in directly interconnected leaves than in indirectly connected leaves. JA accumulation in a lox2-1 mutant plant was initiated rapidly after wounding then slowed progressively compared with the wild type (WT). Despite this, JAZ10.3 expression in the two genotypes was similar. Free cyclopentenone jasmonate levels were similar in both resting WT and lox2-1. In contrast, bound cyclopentenone jasmonates (arabidopsides) were far lower in lox2-1 than in the WT. The major roles of LOX2 are to generate arabidopsides and the large levels of JA that accumulate proximal to the wound. LOX2 is not essential for some of the most rapid events elicited by wounding.

Early change in defence mechanisms and coping in short-term dynamic psychotherapy: relations with symptoms and alliance.

Several patient-related variables have already been investigated as predictors of change in psychodynamic psychotherapy. Defensive functioning is one of them. However, few studies have investigated adaptational processes, encompassing defence mechanisms and coping, from an integrative or comparative viewpoint. This study includes 32 patients, mainly diagnosed with adjustment disorder and undergoing time-limited psychodynamic psychotherapy lasting up to 40 sessions, and will focus on early change in defence and coping. Observer-rater methodology was applied to the transcripts of two sessions of the first part of the psychotherapeutic process. It is assumed that the contextual-relational variable of therapeutic alliance intervenes as moderator on change in adaptational processes. Results corroborated the hypothesis, but only for coping, whereas for defences, overall functioning remained stable over the first 20 sessions of psychotherapy. These results are discussed within the framework of disentangling processes underlying adaptation, i.e., related to issues on trait and state aspects, as well as the role of the therapeutic alliance.