A third approach to gene prediction suggests thousands of additional human transcribed regions

CASC15 polymorphisms are correlated with cervical cancer susceptibility in Chinese women

Background

Cervical cancer is a frequent, common cancer in women, and causes high cancer‐related deaths among women in our world. Accumulating studies provided an important evidence for long noncoding RNA (lncRNA) polymorphisms in the susceptibility of various cancer. Here, we recruited 494 cervical cancer cases and 504 unrelated controls to assess the relationship between CASC15 (OMIM# 616610) polymorphisms and cervical cancer susceptibility.

Methods

Agena MassARRAY platform was conducted to genotype CASC15 polymorphisms. Odds ratios (ORs) and 95% confidence intervals (CIs) were analyzed through logistic regression to adjust for confounding factors, such as age and gender.

Results

Our study suggested that rs12212674 (NC_000006.12:g.22086845T>A) “A” allele was significantly associated with an increased risk of cervical cancer (OR = 1.31, 95% CI = 1.01–1.69, p = .041). The result was demonstrated in the log‐additive model (OR = 1.32, 95% CI = 1.02–1.72, p = .037). After age stratification, we also found that the “TT” genotype of rs4712653 (NC_000006.11:g.22125964T>C) in CASC15 was interaction with a higher cervical cancer risk in subjects aged ≤51 years in the co‐dominant model (OR = 2.08, 95% CI = 1.02–4.25, p = .044) and the recessive model (OR = 2.11, 95% CI = 1.05–4.24, p = .036). Whereas no significant correlation was found among other SNPs of CASC15 polymorphisms and the risk of cervical cancer. MDR analysis illustrated that the interaction between rs7740084 (NC_000006.11:g.21727531G>A), rs1555529 (NC_000006.11:g.21691704A>G), and rs12212674 had a certain effect on the progress of cervical cancer.

Conclusion

Our results revealed a potential interaction between CASC15 polymorphisms and cervical cancer susceptibility. The results provided important insights into CASC15 function in the development of cervical cancer.

A Case/Control Study: AGBL1 Polymorphism Related to Lung Cancer Risk in Chinese Nonsmoking Females

ATP/GTP binding protein like 1 (AGBL1) plays a role in controlling the length of polyglutamate side chains. Polymorphism rs4513061 in AGBL1 is suspected to influence the risk of lung cancer. A case/control study was performed involving 556 cases and 563 controls from a hospital participating in donation. The relationship between rs4513061 and the risk of lung cancer and the interaction between rs4513061 and environmental exposure were determined by the chi-square tests, logistic regression analysis, and crossover analysis. The survival analysis was conducted by Cox proportional hazard regression. The results showed that rs4513061 polymorphism is associated with the risk of lung cancer. The stratified analysis suggested the protective effect of rs4513061 to different histological types of lung cancer, including lung adenocarcinoma (AA vs. GG: odds ratio [OR] = 0.505, 95% confidence interval [CI] = 0.337-0.756, p < 0.001), squamous cell lung cancer (AG vs. GG: OR = 0.488, 95% CI = 0.269-0.883, p = 0.018), and small-cell lung cancer (AA vs. GG: OR = 0.421, 95% CI = 0.216-0.819, p = 0.011). Nevertheless, there was no significant interaction between rs4513061 and cooking oil fume. Significant impact was not observed between the rs4513061 polymorphism and survival time of lung cancer. Our study indicated that rs4513061 in AGBL1 decreases the risk of lung cancer in nonsmoking females from northeast China.

The Regulatory Role of Non-coding RNAs on Programmed Cell Death Four in Inflammation and Cancer

Programmed cell death 4 (PDCD4) is a tumor suppressor gene implicated in many cellular functions, including transcription, translation, apoptosis, and the modulation of different signal transduction pathways. The downstream mechanisms of PDCD4 have been well-discussed, but its upstream regulators have not been systematically summarized. Noncoding RNAs (ncRNAs) are gene transcripts with no protein-coding potential but play a pivotal role in the regulation of the pathogenesis of solid tumors, cardiac injury, and inflamed tissue. In recent studies, many ncRNAs, especially microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), were found to interact with PDCD4 to manipulate its expression through transcriptional regulation and function as oncogenes or tumor suppressors. For example, miR-21, as a classic oncogene, was identified as the key regulator of PDCD4 by targeting its 3′-untranslated region (UTR) to promote tumor proliferation, migration, and invasion in colon, breast, and bladder carcinoma. Therefore, we reviewed the recently emerging pleiotropic regulation of PDCD4 by ncRNAs in cancer and inflammatory disorders and aimed to shed light on the mechanisms of associated diseases, which could be conducive to the development of novel treatment strategies for PDCD4-induced diseases.

Long non-coding RNA CASC15 promotes melanoma progression by epigenetically regulating PDCD4

Background

Long non-coding RNAs (LncRNAs) have been identified as critical regulators in a variety of cancer types. Cancer susceptibility candidate 15 (CASC15), a lncRNA located at chromosome 6p22.3, has been discovered to participate in melanoma progression and phenotype switching. Nevertheless, the roles and molecular mechanisms of CASC15 in melanoma are far from being understood.

Results

We found that CASC15 expression was up-regulated in melanoma tissues and associated with advanced pathological stages. Function experiments displayed that CASC15 knockdown hindered proliferation, facilitated apoptosis and suppressed invasion, while CASC15 overexpression facilitated proliferation and invasion in melanoma cells. Further mechanistic analysis showed that CASC15 epigenetically silenced the expression of programmed cell death 4 (PDCD4) by recruiting EZH2 and increasing H3K27me3 level at the promoter region of PDCD4. Additionally, PDCD4 overexpression inhibited proliferation, enhanced apoptosis and decreased invasion of melanoma cells. Moreover, CASC15-knockdown-induced anti-cancer effects were abated by PDCD4 down-regulation. Furthermore, depletion of CASC15 blocked tumor growth of melanoma by up-regulating PDCD4 in vivo.

Conclusions

CASC15 acts as an oncogene by negatively regulating PDCD4 expression via recruiting EZH2 and subsequently increasing H3K27me3 level. Together, our study indicates that CASC15/EZH2/PDCD4 may serve as a promising therapeutic target for melanoma intervention.

CASC15 gene polymorphisms reduce neuroblastoma risk in Chinese children

In this case-control study, we analyzed the association between three single nucleotide polymorphisms (SNPs) in the CASC15 gene (rs6939340 A>G, rs4712653 T>C, and rs9295536 C>A) and neuroblastoma susceptibility in the Guangdong and Henan populations of China. We genotyped and analyzed 118 cases and 281 control subjects from Henan province and combined them with previously published data from the Guangdong population. In the Henan population, only the rs6939340 G>A variant homozygote AA was associated with decreased neuroblastoma risk [AA vs. GG: adjusted odds ratio (OR) = 0.47, 95% confidence interval (CI) = 0.23-0.98; P=0.045]. All three polymorphisms, individually and in combination, were associated with decreased neuroblastoma susceptibility in the Guangdong population. Moreover, subjects carrying 1-3 of these protective genotypes had lower neuroblastoma susceptibility than non-carriers (adjusted OR=0.65, 95% CI=0.51-0.84, P=0.0007). These results show that all three genetic variants of CASC15 identified in a genome-wide association study (GWAS) decrease neuroblastoma risk in two distinct Chinese populations.

The CASC15 Long Intergenic Noncoding RNA Locus Is Involved in Melanoma Progression and Phenotype Switching

In recent years, considerable advances have been made in the characterization of protein-coding alterations involved in the pathogenesis of melanoma. However, despite their growing implication in cancer, little is known about the role of long non-coding RNAs in melanoma progression. We hypothesized that copy number alterations of intergenic non-protein coding domains could help identify long intergenic non-coding RNAs (lincRNAs) associated with metastatic cutaneous melanoma. Among several candidates, our approach uncovered the chromosome 6p22.3 CASC15 lincRNA locus as a frequently gained genomic segment in metastatic melanoma tumors and cell lines. The locus was actively transcribed in metastatic melanoma cells, and up-regulation of CASC15 expression was associated with metastatic progression to brain metastasis in a mouse xenograft model. In clinical specimens, CASC15 levels increased during melanoma progression and were independent predictors of disease recurrence in a cohort of 141 patients with AJCC stage III lymph node metastasis. Moreover, siRNA knockdown experiments revealed that CASC15 regulates melanoma cell phenotype switching between proliferative and invasive states. Accordingly, CASC15 levels correlated with known gene signatures corresponding to melanoma proliferative and invasive phenotypes. These findings support a key role for CASC15 in metastatic melanoma.

Realistic artificial DNA sequences as negative controls for computational genomics

A common practice in computational genomic analysis is to use a set of ‘background’ sequences as negative controls for evaluating the false-positive rates of prediction tools, such as gene identification programs and algorithms for detection of cis-regulatory elements. Such ‘background’ sequences are generally taken from regions of the genome presumed to be intergenic, or generated synthetically by ‘shuffling’ real sequences. This last method can lead to underestimation of false-positive rates. We developed a new method for generating artificial sequences that are modeled after real intergenic sequences in terms of composition, complexity and interspersed repeat content. These artificial sequences can serve as an inexhaustible source of high-quality negative controls. We used artificial sequences to evaluate the false-positive rates of a set of programs for detecting interspersed repeats, ab initio prediction of coding genes, transcribed regions and non-coding genes. We found that RepeatMasker is more accurate than PClouds, Augustus has the lowest false-positive rate of the coding gene prediction programs tested, and Infernal has a low false-positive rate for non-coding gene detection. A web service, source code and the models for human and many other species are freely available at http://repeatmasker.org/garlic/.

Systematic asymmetric nucleotide exchanges produce human mitochondrial RNAs cryptically encoding for overlapping protein coding genes

GenBank's EST database includes RNAs matching exactly human mitochondrial sequences assuming systematic asymmetric nucleotide exchange-transcription along exchange rules: A→G→C→U/T→A (12 ESTs), A→U/T→C→G→A (4 ESTs), C→G→U/T→C (3 ESTs), and A→C→G→U/T→A (1 EST), no RNAs correspond to other potential asymmetric exchange rules. Hypothetical polypeptides translated from nucleotide-exchanged human mitochondrial protein coding genes align with numerous GenBank proteins, predicted secondary structures resemble their putative GenBank homologue's. Two independent methods designed to detect overlapping genes (one based on nucleotide contents analyses in relation to replicative deamination gradients at third codon positions, and circular code analyses of codon contents based on frame redundancy), confirm nucleotide-exchange-encrypted overlapping genes. Methods converge on which genes are most probably active, and which not, and this for the various exchange rules. Mean EST lengths produced by different nucleotide exchanges are proportional to (a) extents that various bioinformatics analyses confirm the protein coding status of putative overlapping genes; (b) known kinetic chemistry parameters of the corresponding nucleotide substitutions by the human mitochondrial DNA polymerase gamma (nucleotide DNA misinsertion rates); (c) stop codon densities in predicted overlapping genes (stop codon readthrough and exchanging polymerization regulate gene expression by counterbalancing each other). Numerous rarely expressed proteins seem encoded within regular mitochondrial genes through asymmetric nucleotide exchange, avoiding lengthening genomes. Intersecting evidence between several independent approaches confirms the working hypothesis status of gene encryption by systematic nucleotide exchanges.

DNA-energetics-based analyses suggest additional genes in prokaryotes

We present here a novel methodology for predicting new genes in prokaryotic genomes on the basis of inherent energetics of DNA. Regions of higher thermodynamic stability were identified, which were filtered based on already known annotations to yield a set of potentially new genes. These were then processed for their compatibility with the stereo-chemical properties of proteins and tripeptide frequencies of proteins in Swissprot data, which results in a reliable set of new genes in a genome. Quite surprisingly, the methodology identifies new genes even in well-annotated genomes. Also, the methodology can handle genomes of any GC-content, size and number of annotated genes.

Coding constraints modulate chemically spontaneous mutational replication gradients in mitochondrial genomes

Distances from heavy and light strand replication origins determine duration mitochondrial DNA remains singlestranded during replication. Hydrolytic deaminations from A->G and C->T occur more on single- than doublestranded DNA. Corresponding replicational nucleotide gradients exist across mitochondrial genomes, most at 3rd, least 2^nd codon positions. DNA singlestrandedness during RNA transcription causes gradients mainly in long-lived species with relatively slow metabolism (high transcription/replication ratios). Third codon nucleotide contents, evolutionary results of mutation cumulation, follow replicational, not transcriptional gradients in Homo; observed human mutations follow transcriptional gradients. Synonymous third codon position transitions potentially alter adaptive off frame information. No mutational gradients occur at synonymous positions forming off frame stops (these adaptively stop early accidental frameshifted protein synthesis), nor in regions coding for putative overlapping genes according to an overlapping genetic code reassigning stop codons to amino acids. Deviation of 3rd codon nucleotide contents from deamination gradients increases with coding importance of main frame 3rd codon positions in overlapping genes (greatest if these are 2^nd position in overlapping genes). Third codon position deamination gradients calculated separately for each codon family are strongest where synonymous transitions are rarely pathogenic; weakest where transitions are frequently pathogenic. Synonymous mutations affect translational accuracy, such as error compensation of misloaded tRNAs by codon-anticodon mismatches (prevents amino acid misinsertion despite tRNA misacylation), a potential cause of pathogenic mutations at synonymous codon positions. Indeed, codon-family-specific gradients are inversely proportional to error compensation associated with gradient-promoted transitions. Deamination gradients reflect spontaneous chemical reactions in singlestranded DNA, but functional coding constraints modulate gradients.

Genomic signatures of germline gene expression

Transcribed regions in the human genome differ from adjacent intergenic regions in transposable element density, crossover rates, and asymmetric substitution and sequence composition patterns. We tested whether these differences reflect selection or are instead a byproduct of germline transcription, using publicly available gene expression data from a variety of germline and somatic tissues. Crossover rate shows a strong negative correlation with gene expression in meiotic tissues, suggesting that crossover is inhibited by transcription. Strand-biased composition (G+T content) and A → G versus T → C substitution asymmetry are both positively correlated with germline gene expression. We find no evidence for a strand bias in allele frequency data, implying that the substitution asymmetry reflects a mutation rather than a fixation bias. The density of transposable elements is positively correlated with germline expression, suggesting that such elements preferentially insert into regions that are actively transcribed. For each of the features examined, our analyses favor a nonselective explanation for the observed trends and point to the role of germline gene expression in shaping the mammalian genome.

Transposon identification using profile HMMs

BACKGROUND

Transposons are "jumping genes" that account for large quantities of repetitive content in genomes. They are known to affect transcriptional regulation in several different ways, and are implicated in many human diseases. Transposons are related to microRNAs and viruses, and many genes, pseudogenes, and gene promoters are derived from transposons or have origins in transposon-induced duplication. Modeling transposon-derived genomic content is difficult because they are poorly conserved. Profile hidden Markov models (profile HMMs), widely used for protein sequence family modeling, are rarely used for modeling DNA sequence families. The algorithm commonly used to estimate the parameters of profile HMMs, Baum-Welch, is prone to prematurely converge to local optima. The DNA domain is especially problematic for the Baum-Welch algorithm, since it has only four letters as opposed to the twenty residues of the amino acid alphabet.

RESULTS

We demonstrate with a simulation study and with an application to modeling the MIR family of transposons that two recently introduced methods, Conditional Baum-Welch and Dynamic Model Surgery, achieve better estimates of the parameters of profile HMMs across a range of conditions.

CONCLUSIONS

We argue that these new algorithms expand the range of potential applications of profile HMMs to many important DNA sequence family modeling problems, including that of searching for and modeling the virus-like transposons that are found in all known genomes.

The completion of the Mammalian Gene Collection (MGC)

Since its start, the Mammalian Gene Collection (MGC) has sought to provide at least one full-protein-coding sequence cDNA clone for every human and mouse gene with a RefSeq transcript, and at least 6200 rat genes. The MGC cloning effort initially relied on random expressed sequence tag screening of cDNA libraries. Here, we summarize our recent progress using directed RT-PCR cloning and DNA synthesis. The MGC now contains clones with the entire protein-coding sequence for 92% of human and 89% of mouse genes with curated RefSeq (NM-accession) transcripts, and for 97% of human and 96% of mouse genes with curated RefSeq transcripts that have one or more PubMed publications, in addition to clones for more than 6300 rat genes. These high-quality MGC clones and their sequences are accessible without restriction to researchers worldwide.

Conserved introns reveal novel transcripts in Drosophila melanogaster

Noncoding RNAs that are-like mRNAs-spliced, capped, and polyadenylated have important functions in cellular processes. The inventory of these mRNA-like noncoding RNAs (mlncRNAs), however, is incomplete even in well-studied organisms, and so far, no computational methods exist to predict such RNAs from genomic sequences only. The subclass of these transcripts that is evolutionarily conserved usually has conserved intron positions. We demonstrate here that a genome-wide comparative genomics approach searching for short conserved introns is capable of identifying conserved transcripts with a high specificity. Our approach requires neither an open reading frame nor substantial sequence or secondary structure conservation in the surrounding exons. Thus it identifies spliced transcripts in an unbiased way. After applying our approach to insect genomes, we predict 369 introns outside annotated coding transcripts, of which 131 are confirmed by expressed sequence tags (ESTs) and/or noncoding FlyBase transcripts. Of the remaining 238 novel introns, about half are associated with protein-coding genes-either extending coding or untranslated regions or likely belonging to unannotated coding genes. The remaining 129 introns belong to novel mlncRNAs that are largely unstructured. Using RT-PCR, we verified seven of 12 tested introns in novel mlncRNAs and 11 of 17 introns in novel coding genes. The expression level of all verified mlncRNA transcripts is low but varies during development, which suggests regulation. As conserved introns indicate both purifying selection on the exon-intron structure and conserved expression of the transcript in related species, the novel mlncRNAs are good candidates for functional transcripts.

Small molecule activators of pre-mRNA 3' cleavage

3' Cleavage and polyadenylation are obligatory steps in the biogenesis of most mammalian pre-mRNAs. In vitro reconstitution of the 3' cleavage reaction from human cleavage factors requires high concentrations of creatine phosphate (CP), though how CP activates cleavage is not known. Previously, we proposed that CP might work by competitively inhibiting a cleavage-suppressing serine/threonine (S/T) phosphatase. Here we show that fluoride/EDTA, a general S/T phosphatase inhibitor, activates in vitro cleavage in place of CP. Subsequent testing of inhibitors specific for different S/T phosphatases showed that inhibitors of the PPM family of S/T phosphatases, which includes PP2C, but not the PPP family, which includes PP1, PP2A, and PP2B, activated 3' cleavage in vitro. In particular, NCI 83633, an inhibitor of PP2C, activated extensive 3' cleavage at a concentration 50-fold below that required by fluoride or CP. The testing of structural analogs led to the identification of a more potent compound that activated 3' cleavage at 200 microM. While testing CP analogs to understand the origin of its cleavage activation effect, we found phosphocholine to be a more effective activator than CP. The minimal structural determinants of 3' cleavage activation by phosphocholine were identified. Our results describe a much improved small molecule activator of in vitro pre-mRNA cleavage, identify the molecular determinants of cleavage activation by phosphoamines such as phosphocholine, and suggest that a PPM family phosphatase is involved in the negative regulation of mammalian pre-mRNA 3' cleavage.

Transcription induces strand-specific mutations at the 5' end of human genes

A regional analysis of nucleotide substitution rates along human genes and their flanking regions allows us to quantify the effect of mutational mechanisms associated with transcription in germ line cells. Our analysis reveals three distinct patterns of substitution rates. First, a sharp decline in the deamination rate of methylated CpG dinucleotides, which is observed in the vicinity of the 5' end of genes. Second, a strand asymmetry in complementary substitution rates, which extends from the 5' end to 1 kbp downstream from the 3' end, associated with transcription-coupled repair. Finally, a localized strand asymmetry, an excess of C-->T over G-->A substitution in the nontemplate strand confined to the first 1-2 kbp downstream of the 5' end of genes. We hypothesize that higher exposure of the nontemplate strand near the 5' end of genes leads to a higher cytosine deamination rate. Up to now, only the somatic hypermutation (SHM) pathway has been known to mediate localized and strand-specific mutagenic processes associated with transcription in mammalia. The mutational patterns in SHM are induced by cytosine deaminase, which just targets single-stranded DNA. This DNA conformation is induced by R-loops, which preferentially occur at the 5' ends of genes. We predict that R-loops are extensively formed in the beginning of transcribed regions in germ line cells.

Finishing touches: post-translational modification of protein factors involved in mammalian pre-mRNA 3' end formation

In eukaryotes, a pre-messenger RNA (pre-mRNA) must undergo several processing reactions before it is exported to the cytoplasm for translation. One of these reactions, endonucleolytic 3' cleavage at the polyadenylation site, prepares the pre-mRNA for addition of the poly(A) tail and defines the 3' untranslated region (UTR), which typically contains important gene expression regulatory sequences. While the protein factors responsible for the endonucleolytic cleavage have been largely identified, the means by which their action is limited to the 3' end of the transcription unit and coordinated with other co-transcriptional events remains unclear. In this review, we summarize and review recent findings revealing that the mammalian 3' cleavage factors undergo extensive post-translational modification. These modifications include: arginine methylation, lysine sumoylation, lysine acetylation, and the phosphorylation of serine, threonine and tyrosine residues. Every cleavage factor, though not every subunit, is affected. Human Fip1 and the 59 kDa subunit of cleavage factor I emerge as the most frequently modified core cleavage factor subunits. We outline and compare the various proteomic methods that have uncovered these modifications, and review emerging hypotheses concerning their function. The roles of these covalent but reversible modifications in other systems suggest that 3' end formation in mammals relies upon post-translational modification for proper function and regulation.

Mammalian non-LTR retrotransposons: for better or worse, in sickness and in health

Transposable elements (TEs) have shared an exceptionally long coexistence with their host organisms and have come to occupy a significant fraction of eukaryotic genomes. The bulk of the expansion occurring within mammalian genomes has arisen from the activity of type I retrotransposons, which amplify in a "copy-and-paste" fashion through an RNA intermediate. For better or worse, the sequences of these retrotransposons are now wedded to the genomes of their mammalian hosts. Although there are several reported instances of the positive contribution of mobile elements to their host genomes, these discoveries have occurred alongside growing evidence of the role of TEs in human disease and genetic instability. Here we examine, with a particular emphasis on human retrotransposon activity, several newly discovered aspects of mammalian retrotransposon biology. We consider their potential impact on host biology as well as their ultimate implications for the nature of the TE-host relationship.

Genomic DNA from animals shows contrasting strand bias in large and small subsequences

Background

For eukaryotes, there is almost no strand bias with regard to base composition, with exceptions for origins of replication and transcription start sites and transcribed regions. This paper revisits the question for subsequences of DNA taken at random from the genome.

Results

For a typical mammal, for example mouse or human, there is a small strand bias throughout the genomic DNA: there is a correlation between (G - C) and (A - T) on the same strand, (that is between the difference in the number of guanine and cytosine bases and the difference in the number of adenine and thymine bases). For small subsequences – up to 1 kb – this correlation is weak but positive; but for large windows – around 50 kb to 2 Mb – the correlation is strong and negative. This effect is largely independent of GC%. Transcribed and untranscribed regions give similar correlations both for small and large subsequences, but there is a difference in these regions for intermediate sized subsequences. An analysis of the human genome showed that position within the isochore structure did not affect these correlations. An analysis of available genomes of different species shows that this contrast between large and small windows is a general feature of mammals and birds. Further down the evolutionary tree, other organisms show a similar but smaller effect. Except for the nematode, all the animals analysed showed at least a small effect.

Conclusion

The correlations on the large scale may be explained by DNA replication. Transcription may be a modifier of these effects but is not the fundamental cause. These results cast light on how DNA mutations affect the genome over evolutionary time. At least for vertebrates, there is a broad relationship between body temperature and the size of the correlation. The genome of mammals and birds has a structure marked by strand bias segments.

Strand bias structure in mouse DNA gives a glimpse of how chromatin structure affects gene expression

Background

On a single strand of genomic DNA the number of As is usually about equal to the number of Ts (and similarly for Gs and Cs), but deviations have been noted for transcribed regions and origins of replication.

Results

The mouse genome is shown to have a segmented structure defined by strand bias. Transcription is known to cause a strand bias and numerous analyses are presented to show that the strand bias in question is not caused by transcription. However, these strand bias segments influence the position of genes and their unspliced length. The position of genes within the strand bias structure affects the probability that a gene is switched on and its expression level. Transcription has a highly directional flow within this structure and the peak volume of transcription is around 20 kb from the A-rich/T-rich segment boundary on the T-rich side, directed away from the boundary. The A-rich/T-rich boundaries are SATB1 binding regions, whereas the T-rich/A-rich boundary regions are not.

Conclusion

The direct cause of the strand bias structure may be DNA replication. The strand bias segments represent a further biological feature, the chromatin structure, which in turn influences the ease of transcription.

A personal journey of discovery: developing technology and changing biology

This autobiographical article describes my experiences in developing chemically based, biological technologies for deciphering biological information: DNA, RNA, proteins, interactions, and networks. The instruments developed include protein and DNA sequencers and synthesizers, as well as ink-jet technology for synthesizing DNA chips. Diverse new strategies for doing biology also arose from novel applications of these instruments. The functioning of these instruments can be integrated to generate powerful new approaches to cloning and characterizing genes from a small amount of protein sequence or to using gene sequences to synthesize peptide fragments so as to characterize various properties of the proteins. I also discuss the five paradigm changes in which I have participated: the development and integration of biological instrumentation; the human genome project; cross-disciplinary biology; systems biology; and predictive, personalized, preventive, and participatory (P4) medicine. Finally, I discuss the origins, the philosophy, some accomplishments, and the future trajectories of the Institute for Systems Biology.

Transcription factor expression in lipopolysaccharide-activated peripheral-blood-derived mononuclear cells

Transcription factors play a key role in integrating and modulating biological information. In this study, we comprehensively measured the changing abundances of mRNAs over a time course of activation of human peripheral-blood-derived mononuclear cells ("macrophages") with lipopolysaccharide. Global and dynamic analysis of transcription factors in response to a physiological stimulus has yet to be achieved in a human system, and our efforts significantly advanced this goal. We used multiple global high-throughput technologies for measuring mRNA levels, including massively parallel signature sequencing and GeneChip microarrays. We identified 92 of 1,288 known human transcription factors as having significantly measurable changes during our 24-h time course. At least 42 of these changes were previously unidentified in this system. Our data demonstrate that some transcription factors operate in a functional range below 10 transcripts per cell, whereas others operate in a range three orders of magnitude greater. The highly reproducible response of many mRNAs indicates feedback control. A broad range of activation kinetics was observed; thus, combinatorial regulation by small subsets of transcription factors would permit almost any timing input to cis-regulatory elements controlling gene transcription.

Structured RNAs in the ENCODE selected regions of the human genome

Functional RNA structures play an important role both in the context of noncoding RNA transcripts as well as regulatory elements in mRNAs. Here we present a computational study to detect functional RNA structures within the ENCODE regions of the human genome. Since structural RNAs in general lack characteristic signals in primary sequence, comparative approaches evaluating evolutionary conservation of structures are most promising. We have used three recently introduced programs based on either phylogenetic-stochastic context-free grammar (EvoFold) or energy directed folding (RNAz and AlifoldZ), yielding several thousand candidate structures (corresponding to approximately 2.7% of the ENCODE regions). EvoFold has its highest sensitivity in highly conserved and relatively AU-rich regions, while RNAz favors slightly GC-rich regions, resulting in a relatively small overlap between methods. Comparison with the GENCODE annotation points to functional RNAs in all genomic contexts, with a slightly increased density in 3'-UTRs. While we estimate a significant false discovery rate of approximately 50%-70% many of the predictions can be further substantiated by additional criteria: 248 loci are predicted by both RNAz and EvoFold, and an additional 239 RNAz or EvoFold predictions are supported by the (more stringent) AlifoldZ algorithm. Five hundred seventy RNAz structure predictions fall into regions that show signs of selection pressure also on the sequence level (i.e., conserved elements). More than 700 predictions overlap with noncoding transcripts detected by oligonucleotide tiling arrays. One hundred seventy-five selected candidates were tested by RT-PCR in six tissues, and expression could be verified in 43 cases (24.6%).

Widespread mRNA polyadenylation events in introns indicate dynamic interplay between polyadenylation and splicing

mRNA polyadenylation and pre-mRNA splicing are two essential steps for the maturation of most human mRNAs. Studies have shown that some genes generate mRNA variants involving both alternative polyadenylation and alternative splicing. Polyadenylation in introns can lead to conversion of an internal exon to a 3' terminal exon, which is termed composite terminal exon, or usage of a 3' terminal exon that is otherwise skipped, which is termed skipped terminal exon. Using cDNA/EST and genome sequences, we identified polyadenylation sites in introns for all currently known human genes. We found that approximately 20% human genes have at least one intronic polyadenylation event that can potentially lead to mRNA variants, most of which encode different protein products. The conservation of human intronic poly(A) sites in mouse and rat genomes is lower than that of poly(A) sites in 3'-most exons. Quantitative analysis of a number of mRNA variants generated by intronic poly(A) sites suggests that the intronic polyadenylation activity can vary under different cellular conditions for most genes. Furthermore, we found that weak 5' splice site and large intron size are the determining factors controlling the usage of composite terminal exon poly(A) sites, whereas skipped terminal exon poly(A) sites tend to be associated with strong polyadenylation signals. Thus, our data indicate that dynamic interplay between polyadenylation and splicing leads to widespread polyadenylation in introns and contributes to the complexity of transcriptome in the cell.

Genetic mapping at 3-kilobase resolution reveals inositol 1,4,5-triphosphate receptor 3 as a risk factor for type 1 diabetes in Sweden

We mapped the genetic influences for type 1 diabetes (T1D), using 2,360 single-nucleotide polymorphism (SNP) markers in the 4.4-Mb human major histocompatibility complex (MHC) locus and the adjacent 493 kb centromeric to the MHC, initially in a survey of 363 Swedish T1D cases and controls. We confirmed prior studies showing association with T1D in the MHC, most significantly near HLA-DR/DQ. In the region centromeric to the MHC, we identified a peak of association within the inositol 1,4,5-triphosphate receptor 3 gene (ITPR3; formerly IP3R3). The most significant single SNP in this region was at the center of the ITPR3 peak of association (P=1.7 x 10(-4) for the survey study). For validation, we typed an additional 761 Swedish individuals. The P value for association computed from all 1,124 individuals was 1.30 x 10(-6) (recessive odds ratio 2.5; 95% confidence interval [CI] 1.7-3.9). The estimated population-attributable risk of 21.6% (95% CI 10.0%-31.0%) suggests that variation within ITPR3 reflects an important contribution to T1D in Sweden. Two-locus regression analysis supports an influence of ITPR3 variation on T1D that is distinct from that of any MHC class II gene.