Fengbiao Mao

Ras-induced epigenetic inactivation of the RRAD (Ras-related associated with diabetes) gene promotes glucose uptake in a human ovarian cancer model

Background: Increased glucose uptake is essential for carcinogenesis. Results: RasV12-induced epigenetic inactivation of RRAD promotes glucose uptake and tumor formation. Conclusion: RRAD might act as a functional tumor suppressor by inhibiting glucose uptake. Significance: Down-regulation of RRAD in tumor tissues might be associated with the Warburg effect. RRAD (Ras-related associated with diabetes) is a small Ras-related GTPase that is frequently inactivated by DNA methylation of the CpG island in its promoter region in cancer tissues. However, the role of the methylation-induced RRAD inactivation in tumorigenesis remains unclear. In this study, the Ras-regulated transcriptome and epigenome were profiled by comparing T29H (a RasV12-transformed human ovarian epithelial cell line) with T29 (an immortalized but non-transformed cell line) through reduced representation bisulfite sequencing and digital gene expression. We found that RasV12-mediated oncogenic transformation was accompanied by RRAD promoter hypermethylation and a concomitant loss of RRAD expression. In addition, we found that the RRAD promoter was hypermethylated, and its transcription was reduced in ovarian cancer versus normal ovarian tissues. Treatment with the DNA methyltransferase inhibitor 5-aza-2′-deoxycytidine resulted in demethylation in the RRAD promoter and restored RRAD expression in T29H cells. Additionally, treatment with farnesyltransferase inhibitor FTI277 resulted in restored RRAD expression and inhibited DNA methytransferase expression and activity in T29H cells. By employing knockdown and overexpression techniques in T29 and T29H, respectively, we found that RRAD inhibited glucose uptake and lactate production by repressing the expression of glucose transporters. Finally, RRAD overexpression in T29H cells inhibited tumor formation in nude mice, suggesting that RRAD is a tumor suppressor gene. Our results indicate that RasV12-mediated oncogenic transformation induces RRAD epigenetic inactivation, which in turn promotes glucose uptake and may contribute to ovarian cancer tumorigenesis.

Q-RRBS: a quantitative reduced representation bisulfite sequencing method for single-cell methylome analyses

Reduced representation bisulfite sequencing (RRBS) is a powerful method of DNA methylome profiling that can be applied to single cells. However, no previous report has described how PCR-based duplication-induced artifacts affect the accuracy of this method when measuring DNA methylation levels. For quantifying the effects of duplication-induced artifacts on methylome profiling when using ultra-trace amounts of starting material, we developed a novel method, namely quantitative RRBS (Q-RRBS), in which PCR-induced duplication is excluded through the use of unique molecular identifiers (UMIs). By performing Q-RRBS on varying amounts of starting material, we determined that duplication-induced artifacts were more severe when small quantities of the starting material were used. However, through using the UMIs, we successfully eliminated these artifacts. In addition, Q-RRBS could accurately detect allele-specific methylation in absence of allele-specific genetic variants. Our results demonstrate that Q-RRBS is an optimal strategy for DNA methylation profiling of single cells or samples containing ultra-trace amounts of cells.

Comparative RNA-seq analysis reveals potential mechanisms mediating the conversion to androgen independence in an LNCaP progression cell model

The androgen-independent phenotype is an important symptom of refractory prostate cancer. However, the molecular mechanisms underlying this phenotypic conversion remain unclear. Using RNA-seq analysis of androgen-dependent prostate cancer cells (LNCaP) vs. androgen-independent cancer cells (LNCaP-AI-F), we identified 788 differentially expressed genes, 315 alternative splicing events, and eight novel LNCaP-AI-F-specific fusion genes. The fusion genes EIF2AK1-ATR and GLYR1-SLC9A8 were predicted to be damaging and oncogenic. We also observed dramatic changes in androgen receptor (AR)-mediated pathway molecules, including prostate-specific antigen (PSA, a major biomarker of prostate cancer) and AR variants, as well as neuroendocrine-like (NE-like) and tumor stem cell-like characteristics, during androgen-independent phenotype progression. Our findings provide new insights into the regulatory complexities of refractory prostate cancers.

TET1 modulates H4K16 acetylation by controlling auto-acetylation of hMOF to affect gene regulation and DNA repair function

The Ten Eleven Translocation 1 (TET1) protein is a DNA demethylase that regulates gene expression through altering statue of DNA methylation. However, recent studies have demonstrated that TET1 could modulate transcriptional expression independent of its DNA demethylation activity; yet, the detailed mechanisms underlying TET1s role in such transcriptional regulation remain not well understood. Here, we uncovered that Tet1 formed a chromatin complex with histone acetyltransferase Mof and scaffold protein Sin3a in mouse embryonic stem cells by integrative genomic analysis using publicly available ChIP-seq data sets and a series of in vitro biochemical studies in human cell lines. Mechanistically, the TET1 facilitated chromatin affinity and enzymatic activity of hMOF against acetylation of histone H4 at lysine 16 via preventing auto-acetylation of hMOF, to regulate expression of the downstream genes, including DNA repair genes. We found that Tet1 knockout MEF cells exhibited an accumulation of DNA damage and genomic instability and Tet1 deficient mice were more sensitive to x-ray exposure. Taken together, our findings reveal that TET1 forms a complex with hMOF to modulate its function and the level of H4K16Ac ultimately affect gene expression and DNA repair.

Genetic landscape of papillary thyroid carcinoma in the Chinese population

Improvement in the clinical outcome of human cancers requires characterization of the genetic alterations underlying their pathogenesis. Large‐scale genomic and transcriptomic characterization of papillary thyroid carcinomas (PTCs) in Western populations has revealed multiple oncogenic drivers which are essential for understanding pathogenic mechanisms of this disease, while, so far, the genetic landscape in Chinese patients with PTC remains uncharacterized. Here, we conducted a large‐scale genetic analysis of PTCs from patients in China to determine the mutational landscape of this cancer. By performing targeted DNA amplicon and targeted RNA deep‐sequencing, we elucidated the landscape of somatic genetic alterations in 355 Chinese patients with PTC. A total of 88.7% of PTCs were found to harbor at least one candidate oncogenic driver genetic alteration. Among them, around 72.4% of the cases carried BRAF mutations; 2.8% of cases harbored RAS mutations; and 13.8% of cases were characterized with in‐frame gene fusions, including seven newly identified kinase gene fusions. TERT promoter mutations were likely to occur in a sub‐clonal manner in our PTC cohort. The prevalence of somatic genetic alterations in PTC was significantly different between our Chinese cohort and TCGA datasets for American patients. Additionally, combined analyses of genetic alterations and clinicopathologic features demonstrated that kinase gene fusion was associated with younger age at diagnosis, larger tumor size, and lymph node metastasis in PTC. With the analyses of DNA rearrangement sites of RET gene fusions in PTC, signatures of chromosome translocations related to RET fusion events were also depicted. Collectively, our results provide fundamental insight into the pathogenesis of PTC in the Chinese population. Copyright

Whole-Genome Sequencing Reveals Genetic Variation in the Asian House Rat

Whole-genome sequencing of wild-derived rat species can provide novel genomic resources, which may help decipher the genetics underlying complex phenotypes. As a notorious pest, reservoir of human pathogens, and colonizer, the Asian house rat, Rattus tanezumi, is successfully adapted to its habitat. However, little is known regarding genetic variation in this species. In this study, we identified over 41,000,000 single-nucleotide polymorphisms, plus insertions and deletions, through whole-genome sequencing and bioinformatics analyses. Moreover, we identified over 12,000 structural variants, including 143 chromosomal inversions. Further functional analyses revealed several fixed nonsense mutations associated with infection and immunity-related adaptations, and a number of fixed missense mutations that may be related to anticoagulant resistance. A genome-wide scan for loci under selection identified various genes related to neural activity. Our whole-genome sequencing data provide a genomic resource for future genetic studies of the Asian house rat species and have the potential to facilitate understanding of the molecular adaptations of rats to their ecological niches.

MBRidge: an accurate and cost-effective method for profiling DNA methylome at single-base resolution

Organisms and cells, in response to environmental influences or during development, undergo considerable changes in DNA methylation on a genome-wide scale, which are linked to a variety of biological processes. Using MethylC-seq to decipher DNA methylome at single-base resolution is prohibitively costly. In this study, we develop a novel approach, named MBRidge, to detect the methylation levels of repertoire CpGs, by innovatively introducing C-hydroxylmethylated adapters and bisulfate treatment into the MeDIP-seq protocol and employing ridge regression in data analysis. A systematic evaluation of DNA methylome in a human ovarian cell line T29 showed that MBRidge achieved high correlation (R > 0.90) with much less cost (∼10%) in comparison with MethylC-seq. We further applied MBRidge to profiling DNA methylome in T29H, an oncogenic counterpart of T29s. By comparing methylomes of T29H and T29, we identified 131790 differential methylation regions (DMRs), which are mainly enriched in carcinogenesis-related pathways. These are substantially different from 7567 DMRs that were obtained by RRBS and related with cell development or differentiation. The integrated analysis of DMRs in the promoter and expression of DMR-corresponding genes revealed that DNA methylation enforced reverse regulation of gene expression, depending on the distance from the proximal DMR to transcription starting sites in both mRNA and lncRNA. Taken together, our results demonstrate that MBRidge is an efficient and cost-effective method that can be widely applied to profiling DNA methylomes.

Identification of a novel missense (C7W) mutation of SOD1 in a large familial amyotrophic lateral sclerosis pedigree

Mutations of Cu-Zn superoxide dismutase (SOD1) have rarely been identified in Chinese patients with amyotrophic lateral sclerosis (ALS). We recently initiated a program to screen mutations of SOD1, TARDBP, and C9orf72 genes, the most commonly mutated genes in ALS patients in Western countries, in Chinese ALS patients. In this study, we report a novel missense SOD1 mutation with a substitution of tryptophan for cysteine at the seventh amino acid (p.C7W, traditionally named p.C6W) based on HUGO Gene Nomenclature in a familial ALS pedigree. We also found that the activities of SOD1 were significantly decreased in the C7W patient and the carriers of the family, compared with the SOD1 activities of normal family members. Compared with reported C7G and C7S patients, analysis of phenotype revealed relatively mild disease phenotypes in C7W patients, which is correlated with less deteriorated alteration in protein structure. Like those of many other familial ALS families, variable clinical phenotypes in the C7W intrafamily suggest that potential genetic modifiers may contribute to this disease.

Genome-wide identification and divergent transcriptional expression of StAR-related lipid transfer (START) genes in teleosts.

The lipid transfer reactions and the steroidogenic acute regulatory protein (StAR)-related lipid transfer (START) genes have a major role in lipid metabolism. However, START genes and their physiological functions in teleost fishes are relatively unknown. Through genome-wide screening, we identified and annotated 91 START genes in 5 teleost species. Although START domain-containing proteins are augmented in teleost genomes relative to tetrapod genomes, a similar number of genes are shared between them. Asymmetry of paralogous gene loss within the teleost START family and an extra copy of some START genes in teleosts resulting from fish-specific genome duplication have been demonstrated. A distinct transcriptional expression pattern within members of some START groups under different developmental stages suggests divergent functions within the same group in the developmental process. In addition, an asymmetric molecular evolution rate deviating from the neutral expectation has been observed in 7 of 14 teleost fish extra-duplicated pairs. The present study provides valuable information for increasing our understanding of the evolution and gene expression divergence under developmental stages of the START gene family in teleost fishes.

A Bayesian framework to identify methylcytosines from high-throughput bisulfite sequencing data.

High-throughput bisulfite sequencing technologies have provided a comprehensive and well-fitted way to investigate DNA methylation at single-base resolution. However, there are substantial bioinformatic challenges to distinguish precisely methylcytosines from unconverted cytosines based on bisulfite sequencing data. The challenges arise, at least in part, from cell heterozygosis caused by multicellular sequencing and the still limited number of statistical methods that are available for methylcytosine calling based on bisulfite sequencing data. Here, we present an algorithm, termed Bycom, a new Bayesian model that can perform methylcytosine calling with high accuracy. Bycom considers cell heterozygosis along with sequencing errors and bisulfite conversion efficiency to improve calling accuracy. Bycom performance was compared with the performance of Lister, the method most widely used to identify methylcytosines from bisulfite sequencing data. The results showed that the performance of Bycom was better than that of Lister for data with high methylation levels. Bycom also showed higher sensitivity and specificity for low methylation level samples (<1%) than Lister. A validation experiment based on reduced representation bisulfite sequencing data suggested that Bycom had a false positive rate of about 4% while maintaining an accuracy of close to 94%. This study demonstrated that Bycom had a low false calling rate at any methylation level and accurate methylcytosine calling at high methylation levels. Bycom will contribute significantly to studies aimed at recalibrating the methylation level of genomic regions based on the presence of methylcytosines.