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.

Incidence and functional consequences of hMLH1 promoter hypermethylation in colorectal carcinoma

Inactivation of the genes involved in DNA mismatch repair is associated with microsatellite instability (MSI) in colorectal cancer. We report that hypermethylation of the 5′ CpG island of hMLH1 is found in the majority of sporadic primary colorectal cancers with MSI, and that this methylation was often, but not invariably, associated with loss of hMLH1 protein expression. Such methylation also occurred, but was less common, in MSI− tumors, as well as in MSI+ tumors with known mutations of a mismatch repair gene (MMR). No hypermethylation of hMSH2 was found. Hypermethylation of colorectal cancer cell lines with MSI also was frequently observed, and in such cases, reversal of the methylation with 5-aza-2′-deoxycytidine not only resulted in reexpression of hMLH1 protein, but also in restoration of the MMR capacity in MMR-deficient cell lines. Our results suggest that microsatellite instability in sporadic colorectal cancer often results from epigenetic inactivation of hMLH1 in association with DNA methylation.

A second promoter provides an alternative target for therapeutic up-regulation of utrophin in Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is an inherited muscle-wasting disease caused by the absence of a muscle cytoskeletal protein, dystrophin. We have previously shown that utrophin, the autosomal homologue of dystrophin, is able to compensate for the absence of dystrophin in a mouse model of DMD; we have therefore undertaken a detailed study of the transcriptional regulation of utrophin to identify means of effecting its up-regulation in DMD muscle. We have previously isolated a promoter element lying within the CpG island at the 5′ end of the gene and have shown it to be synaptically regulated in vivo. In this paper, we show that there is an alternative promoter lying within the large second intron of the utrophin gene, 50 kb 3′ to exon 2. The promoter is highly regulated and drives transcription of a widely expressed unique first exon that splices into a common full-length mRNA at exon 3. The two utrophin promoters are independently regulated, and we predict that they respond to discrete sets of cellular signals. These findings significantly contribute to understanding the molecular physiology of utrophin expression and are important because the promoter reported here provides an alternative target for transcriptional activation of utrophin in DMD muscle. This promoter does not contain synaptic regulatory elements and might, therefore, be a more suitable target for pharmacological manipulation than the previously described promoter.

A model system to study genomic imprinting of human genes

Somatic-cell hybrids have been shown to maintain the correct epigenetic chromatin states to study developmental globin gene expression as well as gene expression on the active and inactive X chromosomes. This suggests the potential use of somatic-cell hybrids containing either a maternal or a paternal human chromosome as a model system to study known imprinted genes and to identify as-yet-unknown imprinted genes. Testing gene expression by using reverse transcription followed by PCR, we show that functional imprints are maintained at four previously characterized 15q11–q13 loci in hybrids containing a single human chromosome 15 and at two chromosome 11p15 loci in hybrids containing a single chromosome 11. In contrast, three γ-aminobutyric acid type A receptor subunit genes in 15q12–q13 are nonimprinted. Furthermore, we have found that differential DNA methylation imprints at the SNRPN promoter and at a CpG island in 11p15 are also maintained in somatic-cell hybrids. Somatic-cell hybrids therefore are a valid and powerful system for studying known imprinted genes as well as for rapidly identifying new imprinted genes.

Preselection of retrovirally transduced bone marrow avoids subsequent stem cell gene silencing and age-dependent extinction of expression of human β-globin in engrafted mice

Transcriptional silencing of genes transferred into hematopoietic stem cells poses one of the most significant challenges to the success of gene therapy. If the transferred gene is not completely silenced, a progressive decline in gene expression as the mice age often is encountered. These phenomena were observed to various degrees in mouse transplant experiments using retroviral vectors containing a human β-globin gene, even when cis-linked to locus control region derivatives. Here, we have investigated whether ex vivo preselection of retrovirally transduced stem cells on the basis of expression of the green fluorescent protein driven by the CpG island phosphoglycerate kinase promoter can ensure subsequent long-term expression of a cis-linked β-globin gene in the erythroid lineage of transplanted mice. We observed that 100% of mice (n = 7) engrafted with preselected cells concurrently expressed human β-globin and the green fluorescent protein in 20–95% of their RBC for up to 9.5 mo posttransplantation, the longest time point assessed. This expression pattern was successfully transferred to secondary transplant recipients. In the presence of β-locus control region hypersensitive site 2 alone, human β-globin mRNA expression levels ranged from 0.15% to 20% with human β-globin chains detected by HPLC. Neither the proportion of positive blood cells nor the average expression levels declined with time in transplanted recipients. Although suboptimal expression levels and heterocellular position effects persisted, in vivo stem cell gene silencing and age-dependent extinction of expression were avoided. These findings support the further investigation of this type of vector for the gene therapy of human hemoglobinopathies.

Switch from monoallelic to biallelic human IGF2 promoter methylation during aging and carcinogenesis.

We have previously linked aging, carcinogenesis, and de novo methylation within the promoter of the estrogen receptor (ER) gene in human colon. We now examine the dynamics of this process for the imprinted gene for insulin-like growth factor II (IGF2). In young individuals, the P2-4 promoters of IGF2 are methylated exclusively on the silenced maternal allele. During aging, this promoter methylation becomes more extensive and involves the originally unmethylated allele. Most adult human tumors, including colon, breast, lung, and leukemias, exhibit increased methylation at the P2-4 IGF2 promoters, suggesting further spreading during the neoplastic process. In tumors, this methylation is associated with diminished or absent IGF2 expression from the methylated P3 promoter but maintained expression from P1, an upstream promoter that is not contained within the IGF2 CpG island. Our results demonstrate a remarkable evolution of methylation patterns in the imprinted promoter of the IGF2 gene during aging and carcinogenesis, and provide further evidence for a potential link between aberrant methylation and diseases of aging.

A compositional map of human chromosome band Xq28.

The molar fractions of guanine plus cytosine (GC) in DNA were determined for 36 yeast artificial chromosomes (YACs) which almost completely cover human chromosome band Xq28, a terminal reverse band, corresponding to about 8 Mb of DNA. This allowed the construction of the most complete compositional map to date of a chromosomal band; three regions were observed: (i) a proximal 3.5-Mb region formed by GC-poor L and GC-rich H1 isochores; (ii) a middle 2,2-Mb region essentially formed by a GC-rich H2 isochore and a very GC-rich H3 isochore separated by a GC-poor L isochore, YACs from this region being characterized by a striking compositional heterogeneity and instability; and (iii) a distal 1.3-Mb region exclusively formed by GC-poor L isochores. Gene and CpG island concentrations increased with the GC levels of the isochores, as expected. Xq28 exemplifies a subset of reverse bands which are different from the two other subsets, namely from telomeric bands, which are characterized by specific cytogenetic properties and by the predominance of H2 and H3 isochores, and from the majority of reverse bands, which do not contain H2 and H3 isochores.

Construction of a 2.8-megabase yeast artificial chromosome contig and cloning of the human methylthioadenosine phosphorylase gene from the tumor suppressor region on 9p21.

Many human malignant cells lack methylthioadenosine phosphorylase (MTAP) enzyme activity. The gene (MTAP) encoding this enzyme was previously mapped to the short arm of chromosome 9, band p21-22, a region that is frequently deleted in multiple tumor types. To clone candidate tumor suppressor genes from the deleted region on 9p21-22, we have constructed a long-range physical map of 2.8 megabases for 9p21 by using overlapping yeast artificial chromosome and cosmid clones. This map includes the type IIFN gene cluster, the recently identified candidate tumor suppressor genes CDKN2 (p16INK4A) and CDKN2B (p15INK4B), and several CpG islands. In addition, we have identified other transcription units within the yeast artificial chromosome contig. Sequence analysis of a 2.5-kb cDNA clone isolated from a CpG island that maps between the IFN genes and CDKN2 reveals a predicted open reading frame of 283 amino acids followed by 1302 nucleotides of 3' untranslated sequence. This gene is evolutionarily conserved and shows significant amino acid homologies to mouse and human purine nucleoside phosphorylases and to a hypothetical 25.8-kDa protein in the pet gene (coding for cytochrome bc1 complex) region of Rhodospirillum rubrum. The location, expression pattern, and nucleotide sequence of this gene suggest that it codes for the MTAP enzyme.

Hairpins are formed by the single DNA strands of the fragile X triplet repeats: structure and biological implications.

Inordinate expansion and hypermethylation of the fragile X DNA triplet repeat, (GGC)n.(GCC)n, are correlated with the ability of the individual G- and C-rich single strands to form hairpin structures. Two-dimensional NMR and gel electrophoresis studies show that both the G- and C-rich single strands form hairpins under physiological conditions. This propensity of hairpin formation is more pronounced for the C-rich strand than for the G-rich strand. This observation suggests that the C-rich strand is more likely to form hairpin or "slippage" structure and show asymmetric strand expansion during replication. NMR data also show that the hairpins formed by the C-rich strands fold in such a way that the cytosine at the CpG step of the stem is C.C paired. The presence of a C.C mismatch at the CpG site generates local flexibility, thereby providing analogs of the transition to the methyltransferase. In other words, the hairpins of the C-rich strand act as better substrates for the human methyltransferase than the Watson-Crick duplex or the G-rich strand. Therefore, hairpin formation could account for the specific methylation of the CpG island in the fragile X repeat that occurs during inactivation of the FMR1 gene during the onset of the disease.

Complex genetic findings in a female patient with pyruvate dehydrogenase complex deficiency: Null mutations in the PDHX gene associated with unusual expression of the testis-specific PDHA2 gene in her somatic cells

Human pyruvate dehydrogenase complex (PDC) catalyzes a key step in the generation of cellular energy and is composed by three catalytic elements (E1, E2, E3), one structural subunit (E3-binding protein), and specific regulatory elements, phosphatases and kinases (PDKs, PDPs). The E1α subunit exists as two isoforms encoded by different genes: PDHA1 located on Xp22.1 and expressed in somatic tissues, and the intronless PDHA2 located on chromosome 4 and only detected in human spermatocytes and spermatids. We report on a young adult female patient who has PDC deficiency associated with a compound heterozygosity in PDHX encoding the E3-binding protein. Additionally, in the patient and in all members of her immediate family, a full-length testis-specific PDHA2 mRNA and a 5′UTR-truncated PDHA1 mRNA were detected in circulating lymphocytes and cultured fibroblasts, being bothmRNAs translated into full-length PDHA2 and PDHA1 proteins, resulting in the co-existence of both PDHA isoforms in somatic cells.Moreover, we observed that DNA hypomethylation of a CpG island in the coding region of PDHA2 gene is associatedwith the somatic activation of this gene transcription in these individuals. This study represents the first natural model of the de-repression of the testis-specific PDHA2 gene in human somatic cells, and raises some questions related to the somatic activation of this gene as a potential therapeutic approach for most forms of PDC deficiency.

Increased DNA Methylation at the AXIN1 Gene in a Monozygotic Twin from a Pair Discordant for a Caudal Duplication Anomaly

The AXIN1 gene has been implicated in caudal duplication anomalies. Its coding region was sequenced in both members of a monozygotic ( MZ) twin pair discordant for a caudal duplication anomaly, but no mutation was found. Using bisulfite sequencing, we examined methylation at the promoter region of the AXIN1 gene in these twins and in twin and age-matched singleton controls. Methylation of the promoter region in peripheral blood mononucleated cells was variable among individuals, including MZ pairs. In the MZ pair discordant for the caudal duplication, this region of the affected twin was significantly more methylated than that of the unaffected twin (), which was significantly more P < .0001 methylated than those of the controls (). We have confirmed that this CpG island does function as a promoter P = .02 in vitro and that its activity is inversely proportional to the extent of methylation. This finding raises the possibility that hypermethylation of the AXIN1 promoter, by mechanisms as yet undetermined, is associated with the malformation. This case may be paradigmatic for some cases of MZ discordance.

The isolation and characterization of tescalcin, a novel gene differentially expressed in the mouse testis during the early stages of gonadal differentiation

Gonadal development is an ideal model to study organogenesis because a variety of developmental processes can be studied during the differentiation of the bipotential primordium into testis or ovary. To better understand this process, Representational Difference Analysis of cDNA was used to identify genes that are differentially expressed in mouse gonads at 13.5 days post-coitus. The analysis led to the identification of three testis specific genes and a sequence that was only expressed in the ovary. The male genes identified: renin, Col9a3, and a novel gene termed tescalcin had patterns of expression that suggested a role in testis determination. ^ Studies of the tescalcin gene revealed that it is organized into eight exons and seven introns. The gene was located at 64 cM in mouse chromosome 5, where it spans approximately 35 Kb. Three mRNA variants resulting from alternative splicing of intron 5 were identified in mouse tissues. Gel mobility shift assays demonstrated that Sp1 and Sp3 from Y-1, msc-1, and MIN-6 cells nuclear extracts bind the GC-boxes within the tescalcin proximal promoter. Bisulfite sequencing analysis of tescalcin CpG island revealed that it is differentially methylated in male and female mouse embryonic gonads, and that hypermethylation of this region represses expression of tescalcin in the β-TC3 cell line. ^ The major tescalcin mRNA encodes a protein with 214 amino acids that contains a consensus EF-hand Ca2+-binding domain and an N-myristoylation motif. The amino acid sequence of tescalcin is highly conserved among various species, and it showed the highest homology with calcineurin B homologous proteins 1 and 2, and calcineurin B. Western blot analysis using antibodies generated against the tescalcin protein confirmed its presence in specific mouse tissues and cell lines. Immunohistochemical analysis of mouse embryos confirmed the pattern of expression of tescalcin mRNA in fetal testis. Using pull-down assays, glyceraidehydes-3-phosphate dehydrogenase was identified as an interacting and potential functional partner of tescalcin. ^ The identification and characterization of tescalcin as a novel embryonic testicular marker will contribute to the elucidation of the genetic pathways involved in testis development and likely to the understanding of pathological conditions such as sex reversal and infertility. ^

The Role of TET Proteins in the Epigenetic Regulation of Neural Gene Expression and Behavior

Understanding how genes affect behavior is critical to develop precise therapies for human behavioral disorders. The ability to investigate the relationship between genes and behavior has been greatly advanced over the last few decades due to progress in gene-targeting technology. Recently, the Tet gene family was discovered and implicated in epigenetic modification of DNA methylation by converting 5-methylcytosine to 5-hydroxymethylcytosine (5hmC). 5hmC and its catalysts, the TET proteins, are highly abundant in the postnatal brain but with unclear functions. To investigate their neural functions, we generated new lines of Tet1 and Tet3 mutant mice using a gene targeting approach. We designed both mutations to cause a frameshift by deleting the largest coding exon of Tet1 (Tet1Δe4) and the catalytic domain of Tet3 (Tet3Δe7-9). As Tet1 is also highly expressed in embryonic stem cells (ESCs), we generated Tet1 homozygous deleted ESCs through sequential targeting to compare the function of Tet1 in the brain to its role in ESCs. To test our hypothesis that TET proteins epigenetically regulate transcription of key neural genes important for normal brain function, we examined transcriptional and epigenetic differences in the Tet1Δe4 mouse brain. The oxytocin receptor (OXTR), a neural gene implicated in social behaviors, is suggested to be epigenetically regulated by an unknown mechanism. Interestingly, several human studies have found associations between OXTR DNA hypermethylation and a wide spectrum of behavioral traits and neuropsychiatric disorders including autism spectrum disorders. Here we report the first evidence for an epigenetic mechanism of Oxtr transcription as expression of Oxtr is reduced in the brains of Tet1Δe4-/- mice. Likewise, the CpG island overlapping the promoter of Oxtr is hypermethylated during early embryonic development and persists into adulthood. We also discovered altered histone modifications at the hypermethylated regions, indicating the loss of TET1 has broad effects on the chromatin structure at Oxtr. Unexpectedly, we discovered an array of novel mRNA isoforms of Oxtr that are selectively reduced in Tet1Δe4-/- mice. Additionally, Tet1Δe4-/- mice display increased agonistic behaviors and impaired maternal care and short-term memory. Our findings support a novel role for TET1 in regulating Oxtr expression by preventing DNA hypermethylation and implicate TET1 in social behaviors, offering novel insight into Oxtr epigenetic regulation and its role in neuropsychiatric disorders.

De novo methylation of tumor suppressor gene p16/INK4a is a frequent finding in multiple myeloma patients at diagnosis.

The p16 gene competes with cyclin D for binding to CDK4/CDK6 and therefore inhibits CDK4/6 complex kinase activity, resulting in dephosphorylation of pRb and related G1 growth arrest. Inactivation of this gene has been involved in a variety of tumors by different mechanisms: homozygous/hemyzygous deletions, point mutations and methylation of a 5' CpG island into exon E1alpha of the p16 gene. Homozygous deletions have been rarely found in multiple myeloma (MM) and no point mutations have been reported. Two recent studies have reported a high prevalence of methylation in the exon E1alpha of the p16 gene, but included only a small number of cases. We have analyzed the methylation pattern of exon E1alpha of the p16 gene in 101 untreated MM and five primary plasma cell leukemias (PCL). A PCR assay, relying on the inability of some restriction enzymes to digest methylated sequences, was used to analyze the methylation status. Southern blot analysis was used to confirm these results. Forty-one of 101 MM patients (40.5%) as well as four of the five (80%) primary PCL patients had shown methylation of the exon E1alpha. Our study confirms that hypermethylation of the p16 gene is a frequent event in MM. Leukemia (2000) 14, 183-187.

Targeting the serrated pathway of colorectal cancer with mutation in BRAF

A recently acknowledged morphological pathway to colorectal cancer originates from precursor polyps with a serrated appearance due to branching and folding of the colon epithelium. This serrated origin accounts for up to 30% of all colorectal tumors but these are heterogeneous regarding molecular characteristics and patient outcome. Here we review the current knowledge about the classification of this tumor subtype and its association with five key features: mutation status of the BRAF or KRAS genes, the CpG island methylation phenotype, microsatellite instability, immune cell infiltration, and overexpression of GTPase RAC1b. Subsequently, available therapeutic approaches for targeting these molecular characteristics are presented and critically discussed.