Functional importance of conserved nucleotides at the histone RNA 3' processing site.

Histone pre-mRNA 3' processing is controlled by a hairpin element preceding the processing site that interacts with a hairpin-binding protein (HBP) and a downstream spacer element that serves as anchoring site for the U7 snRNP. In addition, the nucleotides following the hairpin and surrounding the processing site (ACCCA'CA) are conserved among vertebrate histone genes. Single to triple nucleotide mutations of this sequence were tested for their ability to be processed in nuclear extract from animal cells. Changing the first four nucleotides had no qualitative and little if any quantitative effects on histone RNA 3' processing in mouse K21 cell extract, where processing of this gene is virtually independent of the HBP. A gel mobility shift assay revealing HBP interactions and a processing assay in HeLa cell extract (where the contribution of HBP to efficient processing is more important) showed that only one of these mutations, predicted to extend the hairpin by one base pair, affected the interaction with HBP. Mutations in the next three nucleotides affected both the cleavage efficiency and the choice of processing sites. Analysis of these novel sites indicated a preference for the nucleotide 5' of the cleavage site in the order A > C > U > G. Moreover, a guanosine in the 3' position inhibited cleavage. The preference for an A is shared with the cleavage/polyadenylation reaction, but the preference order for the other nucleotides is different [Chen F, MacDonald CC, Wilusz J, 1995, Nucleic Acids Res 23:2614-2620].

Common target genes of palatal and gingival fibroblasts for EMD: the microarray approach.

BACKGROUND AND OBJECTIVE Connective tissue grafts are frequently applied, together with Emdogain(®) , for root coverage. However, it is unknown whether fibroblasts from the gingiva and from the palate respond similarly to Emdogain. The aim of this study was therefore to evaluate the effect of Emdogain(®) on fibroblasts from palatal and gingival connective tissue using a genome-wide microarray approach. MATERIAL AND METHODS Human palatal and gingival fibroblasts were exposed to Emdogain(®) and RNA was subjected to microarray analysis followed by gene ontology screening with Database for Annotation, Visualization and Integrated Discovery functional annotation clustering, Kyoto Encyclopedia of Genes and Genomes pathway analysis and the Search Tool for the Retrieval of Interacting Genes/Proteins functional protein association network. Microarray results were confirmed by quantitative RT-PCR analysis. RESULTS The transcription levels of 106 genes were up-/down-regulated by at least five-fold in both gingival and palatal fibroblasts upon exposure to Emdogain(®) . Gene ontology screening assigned the respective genes into 118 biological processes, six cellular components, eight molecular functions and five pathways. Among the striking patterns observed were the changing expression of ligands targeting the transforming growth factor-beta and gp130 receptor family as well as the transition of mesenchymal epithelial cells. Moreover, Emdogain(®) caused changes in expression of receptors for chemokines, lipids and hormones, and for transcription factors such as SMAD3, peroxisome proliferator-activated receptor gamma and those of the ETS family. CONCLUSION The present data suggest that Emdogain(®) causes substantial alterations in gene expression, with similar patterns observed in palatal and gingival fibroblasts.

An intergenic non-coding RNA targets and regulates the ribosome in H. volcanii

As translation is the final step in gene expression it is particularly important to understand the processes involved in translation regulation. It was shown in the last years that a class of RNA, the non-protein-coding RNAs (ncRNAs), is involved in regulation of gene expression via various mechanisms (e.g. gene silencing by microRNAs). Almost all of these ncRNA discovered so far target the mRNA in order to modulate protein biosynthesis, this is rather unexpected considering the crucial role of the ribosome during gene expression. However, recent data from our laboratory showed that there is a new class of ncRNAs, which target the ribosome itself [Gebetsberger et al., 2012/ Pircher et al, 2014]. These so called ribosome-associated ncRNAs (rancRNAs) have an impact on translation regulation, mainly by interfering / modulating the rate of protein biosynthesis. The main goal of this project is to identify and describe novel potential regulatory rancRNAs in H. volcanii with the focus on intergenic candidates. Northern blot analyses already revealed interactions with the ribosome and showed differential expression of rancRNAs during different growth phases or under specific stress conditions. To investigate the biological relevance of these rancRNAs, knock-outs were generated in H. volcanii which were used for phenotypic characterization studies. The rancRNA s194 showed association with the 50S ribosomal subunit in vitro and in vivo and was capable of inhibiting peptide bond formation and seems to inhibit translation in vitro. These preliminary data for the rancRNA s194 make it an interesting candidate for further functional studies to identify the molecular mechanisms by which rancRNAs can modulate protein biosynthesis. Characterization of further rancRNA candidates are also underway.

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-bound non-protein-coding RNA (ncRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [1].

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress- dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-protein-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

tRNA-derived fragments target the small ribosomal subunit to fine-tune translation

Post-transcriptional cleavage of RNA molecules to generate smaller fragments is a widespread mechanism that enlarges the structural and functional complexity of cellular RNomes. In particular, fragments deriving from both precursor and mature tRNAs represent one of the rapidly growing classes of post-transcriptional RNA pieces. Importantly, these tRNA-derived fragments (tRFs) possess distinct expression patterns, abundance, cellular localizations, or biological roles compared with their parental tRNA molecules [1]. Here we present evidence that tRFs from the halophilic archaeon Haloferax volcanii directly bind to ribosomes. In a previous genomic screen for ribosome-associated small RNAs we have identified a 26 residue long fragment originating from the 5’ part of valine tRNA (Val-tRF) to be by far the most abundant tRF in H. volcanii [2]. The Val-tRF is processed in a stress-dependent manner and was found to primarily target the small ribosomal subunit in vitro and in vivo. Translational activity was markedly reduced in the presence of Val-tRF, while control RNA fragments of similar length did not show inhibition of protein biosynthesis. Crosslinking experiments and subsequent primer extension analyses revealed the Val-tRF interaction site to surround the mRNA path in the 30S subunit. In support of this, binding experiments demonstrated that Val-tRF does compete with mRNAs for ribosome binding. Therefore this tRF represents a ribosome-associated non-coding RNA (rancRNA) capable of regulating gene expression in H. volcanii under environmental stress conditions probably by fine-tuning the rate of protein production [3].

Variable effects of the conserved RNA hairpin element upon 3' end processing of histone pre-mRNA in vitro.

We have studied the requirements for efficient histone-specific RNA 3' processing in nuclear extract from mammalian tissue culture cells. Processing is strongly impaired by mutations in the pre-mRNA spacer element that reduce the base-pairing potential with U7 RNA. Moreover, by exchanging the hairpin and spacer elements of two differently processed H4 genes, we find that this difference is exclusively due to the spacer element. Finally, processing is inhibited by the addition of competitor RNAs, if these contain a wild-type spacer sequence, but not if their spacer element is mutated. Conversely, the importance of the hairpin for histone RNA 3' processing is highly variable: A hairpin mutant of the H4-12 gene is processed with almost wild-type efficiency in extract from K21 mouse mastocytoma cells but is strongly affected in HeLa cell extract, whereas an identical hairpin mutant of the H4-1 gene is affected in both extracts. The hairpin defect of H4-12-specific RNA in HeLa cells can be overcome by a compensatory mutation that increases the base complementarity to U7 snRNA. Very similar results were also obtained in RNA competition experiments: processing of H4-12-specific RNA can be competed by RNA carrying a wild-type hairpin element in extract from HeLa, but not K21 cells, whereas processing of H4-1-specific RNA can be competed in both extracts. With two additional histone genes we obtained results that were in one case intermediate and in the other similar to those obtained with H4-1. These results suggest that hairpin binding factor(s) can cooperatively support the ability of U7 snRNPs to form an active processing complex, but is(are) not directly involved in the processing mechanism.

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

Correction of aberrant FGFR1 RNA splicing through blocking intronic regulatory elements

Analysis of the human genome has revealed that more than 74% of human genes undergo alternative RNA splicing. Aberrations in alternative RNA splicing have been associated with several human disorders, including cancer. ^ We studied the aberrant expression of alternative RNA splicing isoforms of the Fibroblast Growth Factor Receptor 1 (FGFR1) gene in a human glioblastoma cancer model. Normal glial cells express the FGFR1α, which contains three extracellular domains. In tumors the most abundant isoform is the FGFR1β, which lacks the first extracellular domain due to the skipping of a single exon, termed alpha. The skipping of the α-exon is regulated by two intronic silencing sequences within the precursor mRNA. Since we observed no mutations on these elements in tumor cells, we hypothesized that the over-expression of regulatory proteins that recognize these sequences is responsible for the aberrant expression of splicing isoforms. Hence, we blocked the formation of protein complexes on the ISS using antisense RNA oligonucleotides in vitro. We also evaluated the impact of the ISS antisense oligonucleotides on the endogenous FGFR1 splicing, in a glioblastoma cell model. By targeting intronic regulatory elements we were able to increase the level of alpha exon inclusion up to 90% in glioblastoma cells. The effect was dose dependent, sequence specific and reproducible in glioblastoma and other cancer cells, which also exhibit an alpha exon skipping phenotype. Targeting FGFR1 endogenous ISS1 and ISS2 sequences did not have an additive or synergistic effect, which suggest a regulatory splicing mechanism that requires the interaction of complexes formed on these elements. An increase in the levels of the FGFR1α isoform resulted in a reduction in cell invasiveness. Also, a significant increase in the levels of caspase 3/7 activities, which is indicative of an elevation in apoptosis levels, suggests that expression of FGFR1β might be relevant for tumor survival. These studies demonstrate that it is possible to prevent aberrant expression of exon skipping events through the targeting of intronic regulatory elements, providing an important new therapeutic tool for the correction of human disease caused by alternative RNA splicing. ^

RNA polymerase II large subunit degradation induced by ecteinascidin 743: Molecular characterization and subsequent rational investigation of antitumor mechanisms in cancers with clinical response

Ecteinascidin 743 (Et-743), which is a novel DNA minor groove alkylator with a unique spectrum of antitumor activity, is currently being evaluated in phase II/III clinical trials. Although the precise molecular mechanisms responsible for the observed antitumor activity are poorly understood, recent data suggests that post-translational modifications of RNA polymerase II Large Subunit (RNAPII LS) may play a central role in the cellular response to this promising anticancer agent. The stalling of an actively transcribing RNAPII LS at Et-743-DNA adducts is the initial cellular signal for transcription-coupled nucleotide excision repair (TC-NER). In this manner, Et-743 poisons TC-NER and produces DNA single strand breaks. Et-743 also inhibits the transcription and RNAPII LS-mediated expression of selected genes. Because the poisoning of TC-NER and transcription inhibition are critical components of the molecular response to Et-743 treatment, we have investigated if changes in RNAPII LS contribute to the disruption of these two cellular pathways. In addition, we have studied changes in RNAPII LS in two tumors for which clinical responses were reported in phase I/II clinical trials: renal cell carcinoma and Ewing's sarcoma. Our results demonstrate that Et-743 induces degradation of the RNAPII LS that is dependent on active transcription, a functional 26S proteasome, and requires functional TC-NER, but not global genome repair. Additionally, we have provided the first experimental data indicating that degradation of RNAPII LS might lead to the inhibition of activated gene transcription. A set of studies performed in isogenic renal carcinoma cells deficient in von Hippel-Lindau protein, which is a ubiquitin-E3-ligase for RNAPII LS, confirmed the central role of RNAPII LS degradation in the sensitivity to Et-743. Finally, we have shown that RNAPII LS is also degraded in Ewing's sarcoma tumors following Et-743 treatment and provide data to suggest that this event plays a role in decreased expression of the Ewing's sarcoma oncoprotein, EWS-Fli1. Altogether, these data implicate degradation of RNAPII LS as a critical event following Et-743 exposure and suggest that the clinical activity observed in renal carcinoma and Ewing's sarcoma may be mediated by disruption of molecular pathways requiring a fully functional RNAPII LS. ^

Regulation of RNA splicing by nonsense mutations

Translation termination as a result of premature nonsense codon-incorporation in a RNA transcript can lead to the production of aberrant proteins with gain-of-function or dominant negative properties that could have deletrious effects on the cell. T-cell Receptor (TCR) genes acquire premature termination codons two-thirds of the time as a result of the error-prone programmed rearrangement events that normally occur during T-cell development. My studies have focused on the fate of TCR precursor mRNAs in response to in-frame nonsense mutations. ^ Previous published studies from our laboratory have shown that TCR precursor mRNAs are subject to nonsense mediated upregulation of pre-mRNA (NMUP). In this dissertation, I performed substitution and deletion analysis to characterize specific regions of TCR which are required to elicit NMUP. I performed frame- and factor-dependence studies to determine its relationship with other nonsense codon induced responses using several approaches including (i) translation dependence studies (ii) deletion and mutational analysis, as well as (iii) siRNA mediated knockdown of proteins involved. I also addressed the underlying molecular mechanism for this pre-mRNA upregulation by (i) RNA half-life studies using a c-fos inducible promoter, and (ii) a variety of assays to determine pre-mRNA splicing efficiency. ^ Using these approaches, I have identified a region of TCR that is both necessary and sufficient to elicit (NMUP). I have also found that neither cytoplasmic translation machinery nor the protein UPF1 are involved in eliciting this nuclear event. I have shown that the NMUP can be induced not only by nonsense and frameshift mutations, but also missense mutations that disrupt a cis splicing element in the exon that contains the mutation. However, the effect of nonsense mutations on pre-mRNA is unique and distinguishable from that of missense mutations in that nonsense mutations can upregulate pre-mRNA in a frame-dependent manner. Lastly, I provide evidence that NMUP occurs by a mechanism in which nonsense mutations inhibit the splicing of introns. In summary, I have found that TCR precursor mRNAs are subject to multiple forces involving both RNA splicing and translation that can either increase or decrease the levels of these precursor mRNAs. ^

Micro RNA expression alterations in patients with glioblastoma

Glioblastoma, also known as glioblastoma multiform or GBM, is the most common and most malignant primary brain tumor. The clinical history of patients with glioblastoma is short, usually less than 3 months in more than 50% of cases after diagnosis. Currently, the methods of glioblastoma treatment are chemotherapy, radiotherapy and surgery. Even with the more effective treatment options, patients with glioblastoma most likely have a median survival time of 10 to 12 months. It is necessary to seek other treatment methods, including gene-targeted treatment. The success of gene-targeted treatment depends critically on the knowledge of genes that may be the cause of, or contribute to disease. To establish a correlate between glioblastoma survival timeline and micro RNA expression alteration, a study of 91 glioblastoma patients was conducted at the University of Texas M. D. Anderson Cancer Center. These 91 glioblastoma patients were newly diagnosed from 2002 to 2007. Statistical analysis was conducted to test the association of miRNA expression alteration between long-term survival and short-term survival glioblastoma. The completion of this proposed study will provide a better understanding of the regulatory role of miRNA in glioblastoma progression.^

Identification of germ cell tumor modifier genes from the 129.MOLF-Chr19 consomic mouse strain

Naturally occurring genetic variants confer susceptibility to disease in the human population, including in testicular germ cell tumor development. Disease susceptibility loci for testicular germ cell tumors have been identified by genetic mapping in humans and mice. However, the identity of many of the susceptibility genes remains unclear. My study utilized a chromosome substitution strain, the 129.MOLF-Chr 19 (or M19 strain), to identify candidate testicular germ cell tumor susceptibility genes. Males of this strain have a high incidence of germ cell tumors in the testes. By forward genetic approaches, five susceptibility loci were fine-mapped and the genetic interactions were dissected. In addition, I identified three protein-coding genes and one micro-RNA as testicular tumor susceptibility genes by genomic screening. Using reverse genetic approaches, I verified one of the candidates, Splicing factor 1, as a modifier of testicular tumor. Deficiency of SF1 significantly reduces the incidence of testicular tumors in mice. This study highlights the advantage of the 129.MOLF-Chr 19 consomic strain in disease gene identification and validation. It also sets the stage to elucidate the molecular mechanisms of tumorigenesis in the testis. ^

GENETIC INTERACTIONS OF FACTORS THAT REGULATE ALTERNATIVE RNA SPLICING IN THE MALE GERM LINE OF DROSOPHILA

Alternative RNA splicing is a critical process that contributes variety to protein functions, and further controls cell differentiation and normal development. Although it is known that most eukaryotic genes produce multiple transcripts in which splice site selection is regulated, how RNA binding proteins cooperate to activate and repress specific splice sites is still poorly understood. In addition how the regulation of alternative splicing affects germ cell development is also not well known. In this study, Drosophila Transformer 2 (Tra2) was used as a model to explore both the mechanism of its repressive function on its own pre-mRNA splicing, and the effect of the splicing regulation on spermatogenesis in testis. Half-pint (Hfp), a protein known as splicing activator, was identified in an S2 cell-based RNAi screen as a co-repressor that functions in combination with Tra2 in the splicing repression of the M1 intron. Its repressive splicing function is found to be sequence specific and is dependent on both the weak 3’ splice site and an intronic splicing silencer within the M1 intron. In addition we found that in vivo, two forms of Hfp are expressed in a cell type specific manner. These alternative forms differ at their amino terminus affecting the presence of a region with four RS dipeptides. Using assays in Drosophila S2 cells, we determined that the alternative N terminal domain is necessary in repression. This difference is probably due to differential localization of the two isoforms in the nucleus and cytoplasm. Our in vivo studies show that both Hfp and Tra2 are required for normal spermatogenesis and cooperate in repression of M1 splicing in spermatocytes. But interestingly, Tra2 and Hfp antagonize each other’s function in regulating germline specific alternative splicing of Taf1 (TBP associated factor 1). Genetic and cytological studies showed that mutants of Hfp and Taf1 both cause similar defects in meiosis and spermatogenesis. These results suggest Hfp regulates normal spermatogenesis partially through the regulation of taf1 splicing. These observations indicate that Hfp regulates tra2 and taf1 activity and play an important role in germ cell differentiation of male flies.

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.