Estimation of protein coding density in a corpus of DNA sequence data

A number of experimental methods have been reported for estimating the number of genes in a genome, or the closely related coding density of a genome, defined as the fraction of base pairs in codons. Recently, DNA sequence data representative of the genome as a whole have become available for several organisms, making the problem of estimating coding density amenable to sequence analytic methods. Estimates of coding density for a single genome vary widely, so that methods with characterized error bounds have become increasingly desirable. We present a method to estimate the protein coding density in a corpus of DNA sequence data, in which a ‘coding statistic’ is calculated for a large number of windows of the sequence under study, and the distribution of the statistic is decomposed into two normal distributions, assumed to be the distributions of the coding statistic in the coding and noncoding fractions of the sequence windows. The accuracy of the method is evaluated using known data and application is made to the yeast chromosome III sequence and to C.elegans cosmid sequences. It can also be applied to fragmentary data, for example a collection of short sequences determined in the course of STS mapping.

A short motif in Drosophila SECIS Binding Protein 2 provides differential binding affinity to SECIS RNA hairpins

Selenoproteins contain the amino acid selenocysteine which is encoded by a UGA Sec codon. Recoding UGA Sec requires a complex mechanism, comprising the cis-acting SECIS RNA hairpin in the 3′UTR of selenoprotein mRNAs, and trans-acting factors. Among these, the SECIS Binding Protein 2 (SBP2) is central to the mechanism. SBP2 has been so far functionally characterized only in rats and humans. In this work, we report the characterization of the Drosophila melanogaster SBP2 (dSBP2). Despite its shorter length, it retained the same selenoprotein synthesis-promoting capabilities as the mammalian counterpart. However, a major difference resides in the SECIS recognition pattern: while human SBP2 (hSBP2) binds the distinct form 1 and 2 SECIS RNAs with similar affinities, dSBP2 exhibits high affinity toward form 2 only. In addition, we report the identification of a K (lysine)-rich domain in all SBP2s, essential for SECIS and 60S ribosomal subunit binding, differing from the well-characterized L7Ae RNA-binding domain. Swapping only five amino acids between dSBP2 and hSBP2 in the K-rich domain conferred reversed SECIS-binding properties to the proteins, thus unveiling an important sequence for form 1 binding.

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.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Genome organization and gene expression shape the transposable element distribution in the Drosophila melanogaster euchromatin.

The distribution of transposable elements (TEs) in a genome reflects a balance between insertion rate and selection against new insertions. Understanding the distribution of TEs therefore provides insights into the forces shaping the organization of genomes. Past research has shown that TEs tend to accumulate in genomic regions with low gene density and low recombination rate. However, little is known about the factors modulating insertion rates across the genome and their evolutionary significance. One candidate factor is gene expression, which has been suggested to increase local insertion rate by rendering DNA more accessible. We test this hypothesis by comparing the TE density around germline- and soma-expressed genes in the euchromatin of Drosophila melanogaster. Because only insertions that occur in the germline are transmitted to the next generation, we predicted a higher density of TEs around germline-expressed genes than soma-expressed genes. We show that the rate of TE insertions is greater near germline- than soma-expressed genes. However, this effect is partly offset by stronger selection for genome compactness (against excess noncoding DNA) on germline-expressed genes. We also demonstrate that the local genome organization in clusters of coexpressed genes plays a fundamental role in the genomic distribution of TEs. Our analysis shows that-in addition to recombination rate-the distribution of TEs is shaped by the interaction of gene expression and genome organization. The important role of selection for compactness sheds a new light on the role of TEs in genome evolution. Instead of making genomes grow passively, TEs are controlled by the forces shaping genome compactness, most likely linked to the efficiency of gene expression or its complexity and possibly their interaction with mechanisms of TE silencing.

The A-kinase anchoring protein (AKAP)-Lbc-signaling complex mediates alpha1 adrenergic receptor-induced cardiomyocyte hypertrophy.

In response to various pathological stresses, the heart undergoes a pathological remodeling process that is associated with cardiomyocyte hypertrophy. Because cardiac hypertrophy can progress to heart failure, a major cause of lethality worldwide, the intracellular signaling pathways that control cardiomyocyte growth have been the subject of intensive investigation. It has been known for more than a decade that the small molecular weight GTPase RhoA is involved in the signaling pathways leading to cardiomyocyte hypertrophy. Although some of the hypertrophic pathways activated by RhoA have now been identified, the identity of the exchange factors that modulate its activity in cardiomyocytes is currently unknown. In this study, we show that AKAP-Lbc, an A-kinase anchoring protein (AKAP) with an intrinsic Rho-specific guanine nucleotide exchange factor activity, is critical for activating RhoA and transducing hypertrophic signals downstream of alpha1-adrenergic receptors (ARs). In particular, our results indicate that suppression of AKAP-Lbc expression by infecting rat neonatal ventricular cardiomyocytes with lentiviruses encoding AKAP-Lbc-specific short hairpin RNAs strongly reduces both alpha1-AR-mediated RhoA activation and hypertrophic responses. Interestingly, alpha1-ARs promote AKAP-Lbc activation via a pathway that requires the alpha subunit of the heterotrimeric G protein G12. These findings identify AKAP-Lbc as the first Rho-guanine nucleotide exchange factor (GEF) involved in the signaling pathways leading to cardiomyocytes hypertrophy.

Circulating microRNAs as biomarkers for detection of autologous blood transfusion.

MicroRNAs (miRNAs) are small non-coding RNAs that regulate various biological processes. Cell-free miRNAs measured in blood plasma have emerged as specific and sensitive markers of physiological processes and disease. In this study, we investigated whether circulating miRNAs can serve as biomarkers for the detection of autologous blood transfusion, a major doping technique that is still undetectable. Plasma miRNA levels were analyzed using high-throughput quantitative real-time PCR. Plasma samples were obtained before and at several time points after autologous blood transfusion (blood bag storage time 42 days) in 10 healthy subjects and 10 controls without transfusion. Other serum markers of erythropoiesis were determined in the same samples. Our results revealed a distinct change in the pattern of circulating miRNAs. Ten miRNAs were upregulated in transfusion samples compared with control samples. Among these, miR-30b, miR-30c, and miR-26b increased significantly and showed a 3.9-, 4.0-, and 3.0-fold change, respectively. The origin of these miRNAs was related to pulmonary and liver tissues. Erythropoietin (EPO) concentration decreased after blood reinfusion. A combination of miRNAs and EPO measurement in a mathematical model enhanced the efficiency of autologous transfusion detection through miRNA analysis. Therefore, our results lay the foundation for the development of miRNAs as novel blood-based biomarkers to detect autologous transfusion.

Hierarchical management of carbon sources is regulated similarly by the CbrA/B systems in Pseudomonas aeruginosa and Pseudomonas putida.

The CbrA/B system in pseudomonads is involved in the utilization of carbon sources and carbon catabolite repression (CCR) through the activation of the small RNAs crcZ in Pseudomonas aeruginosa, and crcZ and crcY in Pseudomonas putida. Interestingly, previous works reported that the CbrA/B system activity in P. aeruginosa PAO1 and P. putida KT2442 responded differently to the presence of different carbon sources, thus raising the question of the exact nature of the signal(s) detected by CbrA. Here, we demonstrated that the CbrA/B/CrcZ(Y) signal transduction pathway is similarly activated in the two Pseudomonas species. We show that the CbrA sensor kinase is fully interchangeable between the two species and, moreover, responds similarly to the presence of different carbon sources. In addition, a metabolomics analysis supported the hypothesis that CCR responds to the internal energy status of the cell, as the internal carbon/nitrogen ratio seems to determine CCR and non-CCR conditions. The strong difference found in the 2-oxoglutarate/glutamine ratio between CCR and non-CCR conditions points to the close relationship between carbon and nitrogen availability, or the relationship between the CbrA/B and NtrB/C systems, suggesting that both regulatory systems sense the same sort or interrelated signal.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

GLUTAMATE RECEPTOR-LIKE genes mediate leaf-to-leaf wound signalling.

Wounded leaves communicate their damage status to one another through a poorly understood process of long-distance signalling. This stimulates the distal production of jasmonates, potent regulators of defence responses. Using non-invasive electrodes we mapped surface potential changes in Arabidopsis thaliana after wounding leaf eight and found that membrane depolarizations correlated with jasmonate signalling domains in undamaged leaves. Furthermore, current injection elicited jasmonoyl-isoleucine accumulation, resulting in a transcriptome enriched in RNAs encoding key jasmonate signalling regulators. From among 34 screened membrane protein mutant lines, mutations in several clade 3 GLUTAMATE RECEPTOR-LIKE genes (GLRs 3.2, 3.3 and 3.6) attenuated wound-induced surface potential changes. Jasmonate-response gene expression in leaves distal to wounds was reduced in a glr3.3 glr3.6 double mutant. This work provides a genetic basis for investigating mechanisms of long-distance wound signalling in plants and indicates that plant genes related to those important for synaptic activity in animals function in organ-to-organ wound signalling.

RsmY, a small regulatory RNA, is required in concert with RsmZ for GacA-dependent expression of biocontrol traits in Pseudomonas fluorescens CHA0.

In the plant-beneficial soil bacterium and biocontrol model organism Pseudomonas fluorescens CHA0, the GacS/GacA two-component system upregulates the production of biocontrol factors, i.e. antifungal secondary metabolites and extracellular enzymes, under conditions of slow, non-exponential growth. When activated, the GacS/GacA system promotes the transcription of a small regulatory RNA (RsmZ), which sequesters the small RNA-binding protein RsmA, a translational regulator of genes involved in biocontrol. The gene for a second GacA-regulated small RNA (RsmY) was detected in silico in various pseudomonads, and was cloned from strain CHA0. RsmY, like RsmZ, contains several characteristic GGA motifs. The rsmY gene was expressed in strain CHA0 as a 118 nt transcript which was most abundant in stationary phase, as revealed by Northern blot and transcriptional fusion analysis. Transcription of rsmY was enhanced by the addition of the strain's own supernatant extract containing a quorum-sensing signal and was abolished in gacS or gacA mutants. An rsmA mutation led to reduced rsmY expression, via a gacA-independent mechanism. Overexpression of rsmY restored the expression of target genes (hcnA, aprA) to gacS or gacA mutants. Whereas mutants deleted for either the rsmY or the rsmZ structural gene were not significantly altered in the synthesis of extracellular products (hydrogen cyanide, 2,4-diacetylphloroglucinol, exoprotease), an rsmY rsmZ double mutant was strongly impaired in this production and in its biocontrol properties in a cucumber-Pythium ultimum microcosm. Mobility shift assays demonstrated that multiple molecules of RsmA bound specifically to RsmY and RsmZ RNAs. In conclusion, two small, untranslated RNAs, RsmY and RsmZ, are key factors that relieve RsmA-mediated regulation of secondary metabolism and biocontrol traits in the GacS/GacA cascade of strain CHA0.

Diabetes mellitus, a microRNA-related disease?

Diabetes mellitus is a complex disease resulting in altered glucose homeostasis. In both type 1 and type 2 diabetes mellitus, pancreatic β cells cannot secrete appropriate amounts of insulin to regulate blood glucose level. Moreover, in type 2 diabetes mellitus, altered insulin secretion is combined with a resistance of insulin-target tissues, mainly liver, adipose tissue, and skeletal muscle. Both environmental and genetic factors are known to contribute to the development of the disease. Growing evidence indicates that microRNAs (miRNAs), a class of small noncoding RNA molecules, are involved in the pathogenesis of diabetes. miRNAs function as translational repressors and are emerging as important regulators of key biological processes. Here, we review recent studies reporting changes in miRNA expression in tissues isolated from different diabetic animal models. We also describe the role of several miRNAs in pancreatic β cells and insulin-target tissues. Finally, we discuss the possible use of miRNAs as blood biomarkers to prevent diabetes development and as tools for gene-based therapy to treat both type 1 and type 2 diabetes mellitus.

Mapping of the genes coding for the two major vaccinia virus core polypeptides.

We have mapped the genes coding for two major structural polypeptides of the vaccinia virus core by hybrid selection and transcriptional mapping. First, RNA was selected by hybridization to restriction fragments of the vaccinia virus genome, translated in vitro and the products were immunoprecipitated with antibodies against the two polypeptides. This approach allowed us to map the genes to the left hand end of the largest Hind III restriction fragment of 50 kilobase pairs. Second, transcriptional mapping of this region of the genome revealed the presence of the two expected RNAs. Both RNAs are transcribed from the leftward reading strand and the 5'-ends of the genes are separated by about 7.5 kilobase pairs of DNA. Thus, two genes encoding structural polypeptides with a similar location in the vaccinia virus particle are clustered at approximately 105 kilobase pairs from the left hand end of the 180 kilobase pair vaccinia virus genome.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

Cross-species GacA-controlled induction of antibiosis in pseudomonads.

Signal extracts prepared from culture supernatants of Pseudomonas fluorescens CHA0 and Pseudomonas aeruginosa PAO stimulated GacA-dependent expression of small RNAs and hence of antibiotic compounds in both hosts. Pseudomonas corrugata LMG2172 and P. fluorescens SBW25 also produced signal molecules stimulating GacA-controlled antibiotic synthesis in strain CHA0, illustrating a novel, N-acyl-homoserine lactone-independent type of interspecies communication.