Gender modulates the expression of calcium-regulating proteins in pediatric atrial myocardium

A differential expression of sarcoplasmic reticulum calcium-ATPase (SERCA2a) and phospholamban (PLB) characterizes the remodeling process in heart failure and atrial arrhythmias in adult patients. Gender is known to modulate the course and prognosis of different forms of heart disease. We hypothesized that gender plays a role in molecular changes of myocardial calcium regulating components already in childhood. Moreover, we studied the influence of volume overloaded (VO) on SERCA2a and PLB in pediatric patients. Quantitative reverse transcription-polymerase chain reaction was used to measure mRNA expression of SERCA2a and PLB in atrial myocardium from 30 pediatric patients (12 girls, 18 boys). Eighteen patients had VO right atria, and 12 patients had not-overloaded atria (NO). Protein expression was studied by Western blot. In the entire population, SERCA2a and PLB expression was not different between girls and boys. If hemodynamic overload was taken into account, SERCA2a mRNA expression was significantly reduced in the VO group compared with the NO group (P = 0.021). The VO versus NO difference was restricted to boys, which corresponds to a highly significant interaction of gender versus VO status (P = 0.002). The PLB to SERCA2a protein ratio was significantly lower in girls (P = 0.028). The decrease in SERCA2a mRNA expression in VO atrial myocardium and the PLB to SERCA2a ratio of protein expression was modulated by gender in this pediatric population. To our knowledge, this study is the first to show the impact of gender on the differential expression of calcium-regulating components in pediatric cardiac patients.

Distinct microRNA Expression Profile in Prostate Cancer Patients with Early Clinical Failure and the Impact of let-7 as Prognostic Marker in High-Risk Prostate Cancer

Background The identification of additional prognostic markers to improve risk stratification and to avoid overtreatment is one of the most urgent clinical needs in prostate cancer (PCa). MicroRNAs, being important regulators of gene expression, are promising biomarkers in various cancer entities, though the impact as prognostic predictors in PCa is poorly understood. The aim of this study was to identify specific miRNAs as potential prognostic markers in high-risk PCa and to validate their clinical impact. Methodology and Principal Findings We performed miRNA-microarray analysis in a high-risk PCa study group selected by their clinical outcome (clinical progression free survival (CPFS) vs. clinical failure (CF)). We identified seven candidate miRNAs (let-7a/b/c, miR-515-3p/5p, -181b, -146b, and -361) that showed differential expression between both groups. Further qRT-PCR analysis revealed down-regulation of members of the let-7 family in the majority of a large, well-characterized high-risk PCa cohort (n = 98). Expression of let-7a/b/and -c was correlated to clinical outcome parameters of this group. While let-7a showed no association or correlation with clinical relevant data, let-7b and let-7c were associated with CF in PCa patients and functioned partially as independent prognostic marker. Validation of the data using an independent high-risk study cohort revealed that let-7b, but not let-7c, has impact as an independent prognostic marker for BCR and CF. Furthermore, we identified HMGA1, a non-histone protein, as a new target of let-7b and found correlation of let-7b down-regulation with HMGA1 over-expression in primary PCa samples. Conclusion Our findings define a distinct miRNA expression profile in PCa cases with early CF and identified let-7b as prognostic biomarker in high-risk PCa. This study highlights the importance of let-7b as tumor suppressor miRNA in high-risk PCa and presents a basis to improve individual therapy for high-risk PCa patients.

Plakoglobin as a regulator of desmocollin gene expression.

Desmosomes are cell adhesion junctions required for the normal development and maintenance of mammalian tissues and organs such as the skin, skin appendages, and the heart. The goal of this study was to investigate how desmocollins (DSCs), transmembrane components of desmosomes, are regulated at the transcriptional level. We hypothesized that differential expression of the Dsc2 and Dsc3 genes is a prerequisite for normal development of skin appendages. We demonstrate that plakoglobin (Pg) in conjunction with lymphoid enhancer-binding factor 1 (Lef-1) differentially regulates the proximal promoters of these two genes. Specifically, we found that Lef-1 acts as a switch activating Dsc2 and repressing Dsc3 in the presence of Pg. Interestingly, we also determined that NF-κB pathway components, the downstream effectors of the ectodysplasin-A (EDA)/ ectodysplasin-A receptor (EDAR)/NF-κB signaling cascade, can activate Dsc2 expression. We hypothesize that Lef-1 and EDA/EDAR/NF-κB signaling contribute to a shift in Dsc isoform expression from Dsc3 to Dsc2 in placode keratinocytes. It is tempting to speculate that this shift is required for the invasive growth of placode keratinocytes into the dermis, a crucial step in skin appendage formation.

Bicarbonate enhances expression of the endocarditis and biofilm associated pilus locus, ebpR-ebpABC, in Enterococcus faecalis.

BACKGROUND: We previously identified ebpR, encoding a potential member of the AtxA/Mga transcriptional regulator family, and showed that it is important for transcriptional activation of the Enterococcus faecalis endocarditis and biofilm associated pilus operon, ebpABC. Although ebpR is not absolutely essential for ebpABC expression (100-fold reduction), its deletion led to phenotypes similar to those of an ebpABC mutant such as absence of pili at the cell surface and, consequently, reduced biofilm formation. A non-piliated ebpABC mutant has been shown to be attenuated in a rat model of endocarditis and in a murine urinary tract infection model, indicating an important participation of the ebpR-ebpABC locus in virulence. However, there is no report relating to the environmental conditions that affect expression of the ebpR-ebpABC locus. RESULTS: In this study, we examined the effect of CO2/HCO3(-), pH, and the Fsr system on the ebpR-ebpABC locus expression. The presence of 5% CO2/0.1 M HCO3(-) increased ebpR-ebpABC expression, while the Fsr system was confirmed to be a weak repressor of this locus. The mechanism by which the Fsr system repressed the ebpR-ebpABC locus expression appears independent of the effects of CO2(-) bicarbonate. Furthermore, by using an ebpA::lacZ fusion as a reporter, we showed that addition of 0.1 M sodium bicarbonate to TSBG (buffered at pH 7.5), but not the presence of 5% CO2, induced ebpA expression in TSBG broth. In addition, using microarray analysis, we found 73 genes affected by the presence of sodium bicarbonate (abs(fold) > 2, P < 0.05), the majority of which belong to the PTS system and ABC transporter families. Finally, pilus production correlated with ebpA mRNA levels under the conditions tested. CONCLUSIONS: This study reports that the ebp locus expression is enhanced by the presence of bicarbonate with a consequential increase in the number of cells producing pili. Although the molecular basis of the bicarbonate effect remains unclear, the pathway is independent of the Fsr system. In conclusion, E. faecalis joins the growing family of pathogens that regulates virulence gene expression in response to bicarbonate and/or CO2.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Mutations upstream of fabI in triclosan resistant Staphylococcus aureus strains are associated with elevated fabI gene expression

Background The enoyl-acyl carrier protein (ACP) reductase enzyme (FabI) is the target for a series of antimicrobial agents including novel compounds in clinical trial and the biocide triclosan. Mutations in fabI and heterodiploidy for fabI have been shown to confer resistance in S. aureus strains in a previous study. Here we further determined the fabI upstream sequence of a selection of these strains and the gene expression levels in strains with promoter region mutations. Results Mutations in the fabI promoter were found in 18% of triclosan resistant clinical isolates, regardless the previously identified molecular mechanism conferring resistance. Although not significant, a higher rate of promoter mutations were found in strains without previously described mechanisms of resistance. Some of the mutations identified in the clinical isolates were also detected in a series of laboratory mutants. Microarray analysis of selected laboratory mutants with fabI promoter region mutations, grown in the absence of triclosan, revealed increased fabI expression in three out of four tested strains. In two of these strains, only few genes other than fabI were upregulated. Consistently with these data, whole genome sequencing of in vitro selected mutants identified only few mutations except the upstream and coding regions of fabI, with the promoter mutation as the most probable cause of fabI overexpression. Importantly the gene expression profiling of clinical isolates containing similar mutations in the fabI promoter also showed, when compared to unrelated non-mutated isolates, a significant up-regulation of fabI. Conclusions In conclusion, we have demonstrated the presence of C34T, T109G, and A101C mutations in the fabI promoter region of strains with fabI up-regulation, both in clinical isolates and/or laboratory mutants. These data provide further observations linking mutations upstream fabI with up-regulated expression of the fabI gene.

Bone-Conditioned Medium Changes Gene Expression in Bone-Derived Fibroblasts.

PURPOSE Autologous bone is used for augmentation in the course of oral implant placement. Bone grafts release paracrine signals that can modulate mesenchymal cell differentiation in vitro. The detailed genetic response of the bone-derived fibroblasts to these paracrine signals has remained elusive. Paracrine signals accumulate in bone-conditioned medium (BCM) prepared from porcine cortical bone chips. MATERIALS AND METHODS In this study, bone-derived fibroblasts were exposed to BCM followed by a whole genome expression profiling and downstream quantitative reverse transciptase polymerase chain reaction of the most strongly regulated genes. RESULTS The data show that ADM, IL11, IL33, NOX4, PRG4, and PTX3 were differentially expressed in response to BCM in bone-derived fibroblasts. The transforming growth factor beta (TGF-β) receptor 1 antagonist SB431542 blocked the effect of BCM on the expression of the gene panel, except for IL33. CONCLUSION These in vitro results extend existing evidence that cortical bone chips release paracrine signals that provoke a robust genetic response in mesenchymal cells that is not exclusively mediated via the TGF-β receptor. The present data provide further insights into the process of graft consolidation.

Profiling invasiveness in head and neck cancer: recent contributions of genomic and transcriptomic approaches.

High-throughput molecular profiling approaches have emerged as precious research tools in the field of head and neck translational oncology. Such approaches have identified and/or confirmed the role of several genes or pathways in the acquisition/maintenance of an invasive phenotype and the execution of cellular programs related to cell invasion. Recently published new-generation sequencing studies in head and neck squamous cell carcinoma (HNSCC) have unveiled prominent roles in carcinogenesis and cell invasion of mutations involving NOTCH1 and PI3K-patwhay components. Gene-expression profiling studies combined with systems biology approaches have allowed identifying and gaining further mechanistic understanding into pathways commonly enriched in invasive HNSCC. These pathways include antigen-presenting and leucocyte adhesion molecules, as well as genes involved in cell-extracellular matrix interactions. Here we review the major insights into invasiveness in head and neck cancer provided by high-throughput molecular profiling approaches.

Clinical usefulness of gene-expression profile to rule out acute rejection after heart transplantation: CARGO II.

AIMS A non-invasive gene-expression profiling (GEP) test for rejection surveillance of heart transplant recipients originated in the USA. A European-based study, Cardiac Allograft Rejection Gene Expression Observational II Study (CARGO II), was conducted to further clinically validate the GEP test performance. METHODS AND RESULTS Blood samples for GEP testing (AlloMap(®), CareDx, Brisbane, CA, USA) were collected during post-transplant surveillance. The reference standard for rejection status was based on histopathology grading of tissue from endomyocardial biopsy. The area under the receiver operating characteristic curve (AUC-ROC), negative (NPVs), and positive predictive values (PPVs) for the GEP scores (range 0-39) were computed. Considering the GEP score of 34 as a cut-off (>6 months post-transplantation), 95.5% (381/399) of GEP tests were true negatives, 4.5% (18/399) were false negatives, 10.2% (6/59) were true positives, and 89.8% (53/59) were false positives. Based on 938 paired biopsies, the GEP test score AUC-ROC for distinguishing ≥3A rejection was 0.70 and 0.69 for ≥2-6 and >6 months post-transplantation, respectively. Depending on the chosen threshold score, the NPV and PPV range from 98.1 to 100% and 2.0 to 4.7%, respectively. CONCLUSION For ≥2-6 and >6 months post-transplantation, CARGO II GEP score performance (AUC-ROC = 0.70 and 0.69) is similar to the CARGO study results (AUC-ROC = 0.71 and 0.67). The low prevalence of ACR contributes to the high NPV and limited PPV of GEP testing. The choice of threshold score for practical use of GEP testing should consider overall clinical assessment of the patient's baseline risk for rejection.

Glioma-specific alternative RNA splicing: A role for polypyrimidine tract binding protein-1

Alternative RNA splicing plays an integral role in cell fate determination and function, especially in the cells of the brain. Errors in RNA processing contribute to diseases such as cancer, where it leads to the production of oncogenic proteins or the loss of tumor suppressors. In silica mining suggests that hundreds of splice isoforms are misexpressed in the glial cell-derived glioma. However, there is little experimental evidence of the prevalence and contribution of these changes and whether they contribute to the formation and progression of this devastating malignancy. To determine the frequency of these aberrant events, global profiling of alternative RNA splice patterns in glioma and nontumor brain was conducted using an exon array. Most splicing changes were less than 5-fold in magnitude and 14 cassette exon events were validated, including 7 previously published events. To determine the possible causes of missplicing, the differential expression levels of splicing factors in these two tissues were also analyzed. Six RNA splicing factors had greater than 2-fold changes in expression. The highest differentially expressed factor was polypyrimidine tract binding protein-1 (PTB). Evaluation by immunohistochemistry determined that this factor was elevated in both early and late stages of glioma. Glial cell-specific PTB expression in the adult brain led me to examine the role of PTB in gliomagenesis. Downregulation of PTB slowed glioma cell proliferation and migration and enhanced cell adhesion to fibronectin and vitronectin. To determine whether PTB was affecting these processes through splicing, genome-wide exon expression levels were correlated with PTB levels. Surprisingly, previously reported PTB target transcripts were insensitive to changes in PTB levels in both patient samples and PTB-depleted glioma cells. Only one validated glioma-specific splice target, RTN4/Nogo, had a significant PTB-mediated splicing change. Downregulation of PTB enhanced inclusion of its alternative exon 3, which encodes an auxiliary domain within a neurite inhibitor protein. Overexpression of this splice isoform in glioma cells slowed proliferation in a manner similar to that observed in PTB knockdown cells. In summary, aberrant expression of splicing factors such as PTB in glioma may elicit changes in splicing patterns that enhance tumorigenesis. ^

Transcriptional and translational regulators of gene expression in germ cell development

Germ cell development is a highly coordinated process driven, in part, by regulatory mechanisms that control gene expression. Not only transcription, but also translation, is under regulatory control to direct proper germ cell development. In this dissertation, I have focused on two regulators of germ cell development. One is the homeobox protein RHOX10, which has the potential to be both a transcriptional and translational regulator in mouse male germ cell development. The other is the RNA-binding protein, Hermes, which functions as a translational regulator in Xenopus laevis female germ cell development. ^ Rhox10 is a member of reproductive homeobox gene X-(linked (Rhox) gene cluster, of which expression is developmentally regulated in developing mouse testes. To identify the cell types and developmental stages in which Rhox10 might function, I characterized its temporal and spatial expression pattern in mouse embryonic, neonatal, and adult tissues. Among other things, this analysis revealed that both the level and the subcellular localization of RHOX10 are regulated during germ cell development. To understand the role of Rhox10 in germ cell development, I generated transgenic mice expressing an artificial microRNA (miRNA) targeting Rhox10. While this artificial miRNA robustly downregulated RHOX10 protein expression in vitro, it did not significantly reduce RHOX10 expression in vivo. So I next elected to knockdown RHOX10 levels in spermatogonial stem cells (SSCs), which I found highly express both Rhox10 mRNA and RHOX10 protein. Using a recently developed in vitro culture system for SSCs combined with a short-hairpin RNA (shRNA) approach, I strongly depleted RHOX10 expression in SSCs. These RHOX10-depleted cells exhibited a defect in the ability to form stem cell clusters in vitro. Expression profiling analysis revealed many genes regulated by Rhox10, including many meiotic genes, which could be downstream of Rhox10 in a molecular pathway that controls SSC differentiation. ^ RNA recognition motif (RRM) containing protein, Hermes is localized in germ plasm, where dormant mRNAs are also located, of Xenopus oocytes, which implicates its role in translational regulator. To understand the function of Hermes in oocyte meiosis, I used a morpholino oligonucleotide (MO) based knockdown approach. Microinjection of Hermes MO into fully grown oocytes, which are arrested in meiotic prophase, caused acceleration of oocytes reentry into meiosis (i.e., maturation) upon progesterone induction. Using a candidate approach, I identified at least three targets of Hermes: Ringo/Spy, Xcat2, and Mos. Ringo/Spy and Mos are known to have functions in oocyte maturation, while Ringo/Spy, Xcat2 mRNA are localized in the germ plasm of oocytes, which drives germ cell specification after fertilization. This led me to propose that Hermes functions in both oocyte maturation and germ cell development through its ability to regulate 3 crucial target mRNAs. ^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

Differential stress response and susceptibility to oxidative injury in alveolar macrophages from C57BL/6J and C3H/HeJ mice following oxidant exposure

Studies have demonstrated a variable response to ozone among individuals and animal species and strains. For instance, C57BL/6J mice have a greater inflammatory response to ozone exposure than C3H/HeJ mice. In these studies, I utilized these strain differences in an effort to derive a mechanistic explanation to the variable strain sensitivity to ozone exposure. Therefore, alveolar macrophages (AM) from C57BL/6J and C3H/HeJ mice were exposed in vitro to hydrogen peroxide ($\rm H\sb2O\sb2$), heat and acetyl ceramide or in vivo to ozone. Necrosis and DNA fragmentation in macrophages from the two murine strains were determined to assess cytotoxicity following these treatments. In addition, synthesis and expression of the stress proteins, stress protein 72 (SP72) and heme oxygenase (HO-1), were examined following treatments. The in vitro experiments were conducted to eliminate the possibility of in vivo confounders (i.e., differences in breathing rates in the two strains) and thus directly implicate some inherent difference between cells from the two murine strains. $\rm H\sb2O\sb2$ and heat caused greater cytotoxicity in AM from C57BL/6J than C3H/HeJ mice and DNA fragmentation was a particularly sensitive indicator of cell injury. Similarly, AM from C57BL/6J mice were more sensitive to ozone exposure than cells from C3H/HeJ mice. Exposure to either 1 or 0.4 ppm ozone caused greater cytotoxicity in macrophages from C57BL/6J mice compared to macrophages from C3H/HeJ mice. The increased sensitivity of AM to injury was associated with decreased synthesis and expression of stress proteins. AM from C57BL/6J mice synthesized and expressed significantly less stress proteins in response to heat and ozone than AM from C3H/HeJ mice. Heat treatment resulted in greater synthesis and expression of SP72. In addition, macrophages from C57BL/6J mice expressed lower amounts of HO-1 than macrophages from C3H/HeJ mice following 0.4 ppm ozone exposure. Therefore, AM from C57BL/6J mice are more susceptible to oxidative injury than AM from C3H/HeJ mice which might be due to differential expression of stress proteins in these cells. ^

MicroRNA Regulation of Prostate Cancer Stem/Progenitor Cells and Prostate Cancer Development

Most human tumors contain a population of cells with stem cell properties, called cancer stem cells (CSCs), which are believed to be responsible for tumor establishment, metastasis, and resistance to clinical therapy. It’s crucial to understand the regulatory mechanisms unique to CSCs, so that we may design CSC-specific therapeutics. Recent discoveries of microRNA (miRNA) have provided a new avenue in understanding the regulatory mechanisms of cancer. However, how miRNAs may regulate CSCs is still poorly understood. Here, we present miRNA expression profiling in six populations of prostate cancer (PCa) stem/progenitor cells that possess distinct tumorigenic properties. Six miRNAs were identified to be commonly and differentially expressed, namely, four miRNAs (miR-34a, let-7b, miR-106a and miR-141) were under-expressed, and two miRNAs (miR-301 and miR-452) were over-expressed in the tumorigenic subsets compared to the corresponding marker-negative subpopulations. Among them, the expression patterns of miR-34, let-7b, miR-141 and miR-301 were further confirmed in the CD44+ human primary prostate cancer (HPCa) samples. We then showed that miR-34a functioned as a critical negative regulator in prostate CSCs and PCa development and metastasis. Over-expression of miR-34a in either bulk or CD44+ PCa cells significantly suppressed clonal expansion, tumor development and metastasis. Systemic delivery of miR-34a in tumor-bearing mice demonstrated a potent therapeutic effect again tumor progression and metastasis, leading to extended animal survival. Of great interest, we identified CD44 itself as a direct and relevant downstream target of miR-34a in mediating its tumor-inhibitory effects. Like miR-34a, let-7 manifests similar tumor suppressive effects in PCa cells. In addition, we observed differential mechanisms between let-7 and miR-34a on cell cycle, with miR-34a mainly inducing G1 cell-cycle arrest followed by cell senescence and let-7 inducing G2/M arrest. MiR-301, on the other hand, exerted a cell type dependent effect in regulating prostate CSC properties and PCa development. In summary, our work reveals that the prostate CSC populations display unique miRNA expression signatures and different miRNAs distinctively and coordinately regulate various aspects of CSC properties. Altogether, our results lay a scientific foundation for developing miRNA-based anti-cancer therapy.

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.