FoxM1B regulates NEDD4-1 expression, leading to cellular transformation and full malignant phenotype in immortalized human astrocytes.

Our recent studies have shown that the FoxM1B transcription factor is overexpressed in human glioma tissues and that the level of its expression correlates directly with glioma grade. However, whether FoxM1B plays a role in the early development of glioma (i.e., in transformation) is unknown. In this study, we found that the FoxM1B molecule causes cellular transformation and tumor formation in normal human astrocytes (NHA) immortalized by p53 and pRB inhibition. Moreover, brain tumors that arose from intracranial injection of FoxM1B-expressing immortalized NHAs displayed glioblastoma multiforme (GBM) phenotypes, suggesting that FoxM1B overexpression in immortalized NHAs not only transforms the cells but also leads to GBM formation. Mechanistically, our results showed that overexpression of FoxM1B upregulated NEDD4-1, an E3 ligase that mediates the degradation and downregulation of phosphatase and tensin homologue (PTEN) in multiple cell lines. Decreased PTEN in turn resulted in the hyperactivation of Akt, which led to phosphorylation and cytoplasmic retention of FoxO3a. Blocking Akt activation with phosphoinositide 3-kinase/Akt inhibitors inhibited the FoxM1B-induced transformation of immortalized NHAs. Furthermore, overexpression of FoxM1B in immortalized NHAs increased the expression of survivin, cyclin D1, and cyclin E, which are important molecules for tumor growth. Collectively, these results indicate that overexpression of FoxM1B, in cooperation with p53 and pRB inhibition in NHA cells, promotes astrocyte transformation and GBM formation through multiple mechanisms.

Conversion of myoblasts to physiologically active neuronal phenotype.

Repressor element 1 (RE1)-silencing transcription factor (REST)/neuron-restrictive silencer factor (NRSF) can repress several terminal neuronal differentiation genes by binding to a specific DNA sequence (RE1/neuron-restrictive silencer element [NRSE]) present in their regulatory regions. REST-VP16 binds to the same RE1/NRSE, but activates these REST/NRSF target genes. However, it is unclear whether REST-VP16 expression is sufficient to cause formation of functional neurons either from neural stem cells or from heterologous stem cells. Here we show that the expression of REST-VP16 in myoblasts grown under muscle differentiation conditions blocked entry into the muscle differentiation pathway, countered endogenous REST/NRSF-dependent repression, activated the REST/NRSF target genes, and, surprisingly, activated other neuronal differentiation genes and converted the myoblasts to a physiologically active neuronal phenotype. Furthermore, in vitro differentiated neurons produced by REST-VP16-expressing myoblasts, when injected into mouse brain, survived, incorporated into the normal brain, and did not form tumors. This is the first instance in which myoblasts were converted to a neuronal phenotype. Our results suggest that direct activation of REST/NRSF target genes with a single transgene, REST-VP16, is sufficient to activate other terminal neuronal differentiation genes and to override the muscle differentiation pathways, and they suggest that this approach provides an efficient way of triggering neuronal differentiation in myoblasts and possibly other stem cells.

MCAM/MUC18 REGULATES MELANOMA PROGRESSION BY MODULATING THE EXPRESSION OF ID-1

The acquisition of the metastatic melanoma phenotype is associated with increased expression of the melanoma cell adhesion molecule MCAM/MUC18 (CD146). However, the mechanism by which MUC18 contributes to melanoma metastasis remains unclear. Herein, we stably silenced MUC18 expression utilizing lentivirus-incorporated small hairpin RNA, in two metastatic melanoma cell lines, A375SM and C8161, and conducted cDNA microarray analysis. We identified and validated that the transcriptional regulator, Inhibitor of DNA Binding-1 (Id-1), previously shown to function as an oncogene in several malignancies, was downregulated by 5.6-fold following MUC18 silencing. Additionally, we found that MUC18 regulated Id-1 expression at the transcriptional level via ATF-3. Interestingly, ATF-3 was upregulated by 6.9 fold in our cDNA microarray analysis following MUC18 silencing. ChIP analysis showed increased binding of ATF-3 to the Id-1 promoter after MUC18 silencing, while mutation of the ATF-3 binding site on the Id-1 promoter increased Id-1 promoter activity in MUC18-silenced cells. These Data suggest that MUC18 silencing promotes inhibition of Id-1 expression by increasing ATF-3 expression and binding to the Id-1 promoter. Rescue of MUC18 reverted the expression of Id-1 and ATF-3, thus validating that they are not off-target effects of MUC18. To further assess the role of Id-1 in melanoma invasion and metastasis, we overexpressed Id-1 in MUC18-silenced cells. Overexpression of Id-1 in MUC18-silenced cells resulted in increased cell invasion, as well as increased expression and activity of MMP-2. Our data further reveal that Id-1 regulates MMP-2 at the transcriptional level through Sp1 and Ets-1. This is the first report to demonstrate that MUC18 does not act exclusively in cell adherence, but is also involved in cell signaling that regulates the expression of genes, such as Id-1 and ATF-3, thus contributing to the metastatic melanoma phenotype.

ROLE OF VITAMIN-D3 AND RETINOIC ACID IN A HUMAN THP-1 MACROPHAGE MODEL OF MYCOBACTERIUM TUBERCULOSIS INFECTION

Mycobacterium tuberculosis (Mtb) replicates within the human macrophages and we investigated the activating effects of retinoic acid (RA) and vitamin D3 (VD) on macrophages in relation to the viability of Mtb. A combination of these vitamins (RAVD) enhanced the receptors on THP-1 macrophage (Mannose receptor and DC-SIGN) that increased mycobacterial uptake but inhibited thesubsequent intracellular growth of Mtb by inducing reactive oxygen species (ROS) and autophagy. RAVD also enhanced antigen presenting and homing receptors in THPs that suggested an activated phenotype for THPs following RAVD treatment. RAVD mediated activation was also associated with a marked phenotypic change in Mtb infected THPs that fused with adjacent cells to formmultinucleate giant cells (MNGCs). Typically MNGCs occurred over 30 days of in vitro culture and contained non-replicating persisting Mtb for as long as 60 days in culture. We propose that the RAVD mediated inhibition of replicating Mtb leading to persistence of non-replicating Mtb within THPs may provide a novel human macrophage model simulating formation of MNGCs in humanlungs.

Skeletal abnormalities in mice lacking extracellular matrix proteins, thrombospondin-1, thrombospondin-3, thrombospondin-5, and type IX collagen.

Thrombospondin-5 (TSP5) is a large extracellular matrix glycoprotein found in musculoskeletal tissues. TSP5 mutations cause two skeletal dysplasias, pseudoachondroplasia and multiple epiphyseal dysplasia; both show a characteristic growth plate phenotype with retention of TSP5, type IX collagen (Col9), and matrillin-3 in the rough endoplasmic reticulum. Whereas most studies focus on defining the disease process, few functional studies have been performed. TSP5 knockout mice have no obvious skeletal abnormalities, suggesting that TSP5 is not essential in the growth plate and/or that other TSPs may compensate. In contrast, Col9 knockout mice have diminished matrillin-3 levels in the extracellular matrix and early-onset osteoarthritis. To define the roles of TSP1, TSP3, TSP5, and Col9 in the growth plate, all knockout and combinatorial strains were analyzed using histomorphometric techniques. While significant alterations in growth plate organization were found in certain single knockout mouse strains, skeletal growth was only mildly disturbed. In contrast, dramatic changes in growth plate organization in TSP3/5/Col9 knockout mice resulted in a 20% reduction in limb length, corresponding to similar short stature in humans. These studies show that type IX collagen may regulate growth plate width; TSP3, TSP5, and Col9 appear to contribute to growth plate organization; and TSP1 may help define the timing of growth plate closure when other extracellular proteins are absent.

Genetic and hormonal factors modulate spreading depression and transient hemiparesis in mouse models of familial hemiplegic migraine type 1.

Familial hemiplegic migraine type 1 (FHM1) is an autosomal dominant subtype of migraine with aura that is associated with hemiparesis. As with other types of migraine, it affects women more frequently than men. FHM1 is caused by mutations in the CACNA1A gene, which encodes the alpha1A subunit of Cav2.1 channels; the R192Q mutation in CACNA1A causes a mild form of FHM1, whereas the S218L mutation causes a severe, often lethal phenotype. Spreading depression (SD), a slowly propagating neuronal and glial cell depolarization that leads to depression of neuronal activity, is the most likely cause of migraine aura. Here, we have shown that transgenic mice expressing R192Q or S218L FHM1 mutations have increased SD frequency and propagation speed; enhanced corticostriatal propagation; and, similar to the human FHM1 phenotype, more severe and prolonged post-SD neurological deficits. The susceptibility to SD and neurological deficits is affected by allele dosage and is higher in S218L than R192Q mutants. Further, female S218L and R192Q mutant mice were more susceptible to SD and neurological deficits than males. This sex difference was abrogated by ovariectomy and senescence and was partially restored by estrogen replacement, implicating ovarian hormones in the observed sex differences in humans with FHM1. These findings demonstrate that genetic and hormonal factors modulate susceptibility to SD and neurological deficits in FHM1 mutant mice, providing a potential mechanism for the phenotypic diversity of human migraine and aura.

Defective pituitary function and gamete interactions in $\beta$1,4-galactosyltransferase null mice

Despite much attention, the function of oligosaccharide chains of glycoproteins remains largely unknown. Our understanding of oligosaccharide function in vivo has been limited to the use of reagents and targeted mutations that eliminate entire oligosaccharide chains. However, most, if not all biological functions for oligosaccharides have been attributed to specific terminal sequences on these oligosaccharides, yet there have been few studies to examine the consequences of modifying terminal oligosaccharide structures in vivo. To address this issue, mice were created bearing a targeted mutation in $\beta$1,4-galactosyltransferase, an enzyme responsible for elaboration of many of the proposed biologically-active carbohydrate epitopes. Most galactosyltransferase-null mice died within the first few weeks after birth and were characterized by stunted growth, thin skin, sparse hair, and dehydration. In addition, the adrenal cortices were poorly stratified and spermatogenesis was delayed. The few surviving adults had puffy skin (myxedema), difficulty delivering pups at birth (dystocia), and failed to lactate (agalactosis). All of these defects are consistant with endocrine insufficiency, which was confirmed by markedly decreased levels of serum thyroxine. The anterior pituitary gland appeared functionally delayed in newborn mutant mice, since the constituent cells were quiescent and nonsecretory, unlike that of control littermates. However, the anterior pituitary acquired a normal secretory phenotype during neonatal development, although it remained abnormally small and its glycoprotein hormones were devoid of $\beta$1,4-galactosyl residues. These results support in vitro studies suggesting that incomplete glycosylation of pituitary hormones leads to the creation of hormone antagonists that down regulate subsequent endocrine function producing polyglandular endocrine insufficiency. More surprisingly, the fact that some mice survive this neonatal period indicates the presence of a previously unrecognized compensatory pathway for glycoprotein hormone glycosylation and/or action.^ In addition to its well-studied biosynthetic function in the Golgi complex, a GalTase isoform is also expressed on the sperm surface where it functions as a gamete receptor during fertilization by binding to its oligosaccharide ligand on the egg coat glycoprotein, ZP3. Aggregation of GalTase by multivalent ZP3 oligosaccharides activates a G-protein cascade leading to the acrosome reaction. Although GalTase-null males are fertile, the mutant sperm bind less ZP3 than wild-type sperm, and are unable to undergo the acrosome reaction in response to either zona pellucida glycoproteins or to anti-GalTase anti-serum, as do wild-type sperm. However, mutant and wild-type sperm undergo the acrosome reaction normally in response to calcium ionophore which bypasses the requirement for ZP3 binding. Interestingly, the phenotype of the GalTase-null sperm is reciprocal to that of sperm that overexpress surface GalTAse and which bind more ZP3 leading to precocious acrosome reactions. These results confirm that GalTase functions as at least one of the sperm receptors for ZP3, and that GalTase participates in the ZP3-induced signal transduction pathway during zona pellucida-induced acrosome reactions. ^

Effects of galectin-1 and galectin-3 expression on prostate cancer cell adhesion, differentiation, proliferation, and tumorigenicity

The vertebrate $\beta$-galactoside-binding lectins galectin-1 and galectin-3 have been proposed to function in diverse cellular processes such as adhesion, proliferation, differentiation, and tumorigenesis. Experiments were initiated to further study the functional properties of these molecules. A prostate cancer cell line, LNCaP, was identified which expressed neither galectin. This line was stably transfected with cDNA for either galectin-1 or galectin-3. The resultant clones were used to study effects on critical cell processes. LNCaP cells expressing galectin-1 on the surface were found to bind more rapidly than control lines to the human extracellular matrix proteins laminin and fibronectin, although overall binding was not increased. To analyze effects on differentiation, LNCaP cells were studied which had either been transfected with galectin-1 or which had been induced to express endogenous galectin-1 by treatment with the differentiation agent sodium butyrate. In both cases, cells displayed a slower rate of growth and increased rate of apoptosis. A transient decrease in expression of prostate specific antigen was seen in the butyrate treated cells but not in the transfected cells. To investigate the role of galectins in the process of malignant transformation and progression, immunohistochemical analysis was performed on formalin-fixed, paraffin-embedded sections of human prostate tissue, the premalignant lesion prostatic intraepithelial neoplasia, primary adenocarcinoma of the prostate, and foci of metastatic prostate cancer. Galectin-1 expression was relatively constant throughout in contrast to galectin-3 which demonstrated significantly less expression in primary and metastatic tumors. LNCaP cells transfected with galectin-3 cDNA displayed lower proliferation rates, increased spontaneous apoptosis, and G1 growth phase arrest compared to controls. Four of six galectin-3 lines tested were less tumorigenic in nude mice than controls. The following conclusions are drawn regarding the role of galectin-1 and galectin-3 expression in the context of prostate cancer: (1) galectin-1 may participate in the early stages of cancer cell adhesion to extracellular matrix proteins; (2) galectin-1 expression results in a differentiated phenotype and may contribute to differentiation induction by butyrate; (3) galectin-3 expression correlates inversely with prostate cell tumorigenesis and prostate cancer metastasis. ^

The role of conserved region 1 of Escherichia coli sigma-70 in transcription initiation

The sigma (σ) subunit of eubacterial RNA polymerase is essential for initiation of transcription at promoter sites. σ factor directs the RNA polymerase core subunits ( a2bb′ ) to the promoter consensus elements and thereby confers selectivity for transcription initiation. The N-terminal domain (region 1.1) of Escherichia coli σ70 has been shown to inhibit DNA binding by the C-terminal DNA recognition domains when σ is separated from the core subunits. Since DNA recognition by RNA polymerase is the first step in transcription, it seemed plausible that region 1 might also influence initiation processes subsesquent to DNA binding. This study explores the functional roles of regions 1.1 and 1.2 of σ70 in transcription initiation. Analysis in vitro of the transcriptional properties of a series of N-terminally truncated σ70 derivates revealed a critical role for region 1.1 at several key stages of initiation. Deletion of the first 75 to 100 amino acids of σ70 (region 1.1) resulted in both a slow rate of transition from a closed promoter complex to a DNA-strand-separated open complex, as well as a reduced efficiency of transition from the open complex to a transcriptionally active open complex. These effects were partially reversed by addition of a polypeptide containing region 1.1 in trans. Therefore, region 1.1 not only modulates DNA binding but is important for efficient transcription initiation, once a closed complex has formed. A deletion of the first 133 amino acids which removes both regions 1.1 and 1.2 resulted in arrest of initiation at the earliest closed complex, suggesting that region 1.2 is required for open complex formation. Mutagenesis of region 1.1 uncovered a mechanistically important role for isoleucine at position 53 (I53). Substitution of I53 with alanine created a σ factor that associated with the core subunits to form holoenzyme, but the holoenzyme was severely deficient for promoter binding. The I53A phenotype was suppressed in vivo by truncation of five amino acids from the C-terminus of σ 70. These observations are consistent with a model in which σ 70I53A fails to undergo a critical conformational change upon association with the core subunits, which is needed to expose the DNA-binding domains and confer promoter recognition capability upon holoenzyme. To understand the basis of the autoinhibitory properties of the σ70 N-terminal domain, in the absence of core RNA polymerase, a preliminary physical assessment of the interdomain interactions within the σ70 subunit was launched. Results support a model in which N-terminal amino acids are in close proximity to residues in the C-terminus of the σ 70 polypeptide. ^

Effects of butyrate on the colon cancer cell phenotype-chemosensitization and upregulation of proteins implicated in cancer progression

Colon cancer is the second leading cause of cancer mortality in the U.S. Surgery is the only truly effective human colon cancer (HCC) therapy due to marked intrinsic drug resistance. The inefficacy of therapies developed for metastatic HCC suggests that advances in colon cancer chemoprevention and chemotherapy will be needed to reduce HCC mortality. The dietary fiber metabolite butyrate (NaB) is a candidate cancer chemopreventive agent that inhibits growth, promotes differentiation and stimulates apoptosis of HCC cells. Epidemiological and experimental studies suggest that dietary fiber protects against the development of HCC, however, recent large prospective trials have not found significant protection. ^ The first central hypothesis of this dissertation project is that the diversity of phenotypic changes induced by NaB in HCC cells includes molecular alterations that oppose its chemopreventive action and thereby limit its efficacy. We investigated the effect of NaB on the expression/activity of epidermal growth factor receptor (EGFR) and cyclooxygenase-2 (COX-2) in HCC HT29 cells. NaB treatment induced a 13-fold increase in EGFR expression in concert with its chemopreventive action in vitro, i.e., induction of growth suppression and G1 arrest, apoptosis and a differentiated phenotype. NaB-induced EGFR was active based on multiple lines of evidence. The EGFR was: (1) heavily phosphorylated at Tyrosine (P-Tyr); (2) associated with the cytoskeleton; (3) localized at the cell surface, and activated in response to EGF; and (4) NaB treatment of the cells induced activation of the EGFR effector Erk1/2. NaB treatment also induced a 7-fold increase in COX-2 expression. The NaB-induced COX-2 was active based on significantly increased PGE2 production. ^ The second central hypothesis is that NaB treatment would render HCC cells more chemosensitive to chemotherapy agents based on the increased apoptotic index induced by NaB. NaB treatment chemosensitized HT29 cells to 5-FU and doxorubicin, despite increases in the expression of P-glycoprotein and a related drug resistance protein (MRP). ^ These results raise the intriguing possibility that the chemopreventive effects of fiber may require concomitant treatment with EGFR and/or COX-2 inhibitors. Similarly, NaB may be a rational drug to combine with existing chemotherapeutic agents for the management of advanced HCC patients. ^

TUMOR PROGRESSION: ANALYSIS OF THE INSTABILITY OF THE METASTATIC PHENOTYPE, SENSITIVITY TO RADIATION AND CHEMOTHERAPY (PHENOTYPIC DRIFT, BREAST CANCER)

The major complications for tumor therapy are (i) tumor spread (metastasis); (ii) the mixed nature of tumors (heterogeneity); and (iii) the capacity of tumors to evolve (progress). To study these tumor characteristics, the rat 13762NF mammary adenocarcinoma was cloned and studied for metastatic properties and sensitivities to therapy (chemotherapy, radiation and hyperthermia). The cell clones were heterogeneous and no correlation between metastatic potential and therapeutic sensitivities was observed. Further, these phenotypes were unstable during passage in vitro; yet, the changes were clone dependent and reproducible using different cryoprotected cell stocks. To understand the phenotypic instability, subclones were isolated from low and high passage cell clones. Each subclone possessed a unique composite phenotype. Again, no apparent correlation was seen between metastatic potential and sensitivity to therapy. The results demonstrated that (1) tumor cells are heterogeneous for multiple phenotypes; (2) tumor cells are unstable for multiple phenotypes; (3) the magnitude, direction and time of occurrence of phenotypic drift is clone dependent; (4) the sensitivity of cell clones to ionizing radiation (gamma or heat) and chemotherapy agents is independent of their metastatic potential; (5) shifts in metastatic potential and sensitivity to therapy may occur simultaneously but are not linked; and (6) tumor cells independently diverge to form several subpopulations with unique phenotypic profiles. ^

Loss of activator protein-2alpha results in overexpression of the thrombin receptor and correlates with the malignant phenotype of human melanoma

Increasing evidence demonstrates that the thrombin receptor (protease activated receptor-1, PAR-1) plays a major role in tumor invasion and contributes to the metastatic phenotype of human melanoma. We demonstrate that the metastatic potential of human melanoma cells correlates with overexpression of PAR-1. The promoter of the PAR-1 gene contains multiple putative AP-2 and Sp1 consensus elements. We provide evidence that an inverse correlation exists between the expression of AP-2 and the expression of PAR-1 in human melanoma cells. Re-expression of AP-2 in WM266-4 melanoma cells (AP-2 negative) resulted in decreased mRNA and protein expression of PAR-1 and significantly reduced the tumor potential in nude mice. ChIP analysis of the PAR-1 promoter regions bp −365 to −329 (complex 1) and bp −206 to −180 (complex 2) demonstrates that in metastatic cells Sp1 is predominantly binding to the PAR-1 promoter, while in nonmetastatic cells AP-2 is bound. In vitro analysis of complex 1 demonstrates that AP-2 and Sp1 bind to this region in a mutually exclusive manner. Transfection experiments with full-length and progressive deletions of the PAR-1 promoter luciferase constructs demonstrated that metastatic cells had increased promoter activity compared to low and nonmetastatic melanoma cells. Our data shows that exogenous AP-2 expression decreased promoter activity, while transient expression of Sp1 further activated expression of the reporter gene. Mutational analysis of complex 1 within PAR-1 luciferase constructs further demonstrates that the regulation of PAR-1 is mediated through interactions with AP-2 and Sp1. Moreover, loss of AP-2 in metastatic cells alters the AP-2 to Sp1 ratio and DNA-binding activity resulting in overexpression of PAR-1. In addition, we evaluated the expression of AP-2 and PAR-1 utilizing a tissue microarray of 93 melanocytic lesions spanning from benign nevi to melanoma metastasis. We report loss of AP-2 expression in malignant tumors compared to benign tissue while PAR-1 was expressed more often in metastatic melanoma cells than in benign melanocytes. We propose that loss of AP-2 results in increased expression of PAR-1, which in turn results in upregulation of gene products that contribute to the metastatic phenotype of melanoma. ^

Novel Imaging-Based Techniques Reveal a Role for PD-1/PD-L1 in Tumor Immune Surveillance in the Lung

The binding of immune inhibitory receptor Programmed Death 1 (PD-1) on T cells to its ligand PD-L1 has been implicated as a major contributor to tumor induced immune suppression. Clinical trials of PD-L1 blockade have proven effective in unleashing therapeutic anti-tumor immune responses in a subset of patients with advanced melanoma, yet current response rates are low for reasons that remain unclear. Hypothesizing that the PD-1/PD-L1 pathway regulates T cell surveillance within the tumor microenvironment, we employed intravital microscopy to investigate the in vivo impact of PD-L1 blocking antibody upon tumor-associated immune cell migration. However, current analytical methods of intravital dynamic microscopy data lack the ability to identify cellular targets of T cell interactions in vivo, a crucial means for discovering which interactions are modulated by therapeutic intervention. By developing novel imaging techniques that allowed us to better analyze tumor progression and T cell dynamics in the microenvironment; we were able to explore the impact of PD-L1 blockade upon the migratory properties of tumor-associated immune cells, including T cells and antigen presenting cells, in lung tumor progression. Our results demonstrate that early changes in tumor morphology may be indicative of responsiveness to anti-PD-L1 therapy. We show that immune cells in the tumor microenvironment as well as tumors themselves express PD-L1, but immune phenotype alone is not a predictive marker of effective anti-tumor responses. Through a novel method in which we quantify T cell interactions, we show that T cells are largely engaged in interactions with dendritic cells in the tumor microenvironment. Additionally, we show that during PD-L1 blockade, non-activated T cells are recruited in greater numbers into the tumor microenvironment and engage more preferentially with dendritic cells. We further show that during PD-L1 blockade, activated T cells engage in more confined, immune synapse-like interactions with dendritic cells, as opposed to more dynamic, kinapse-like interactions with dendritic cells when PD-L1 is free to bind its receptor. By advancing the contextual analysis of anti-tumor immune surveillance in vivo, this study implicates the interaction between T cells and tumor-associated dendritic cells as a possible modulator in targeting PD-L1 for anti-tumor immunotherapy.

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.