The Impact of Childhood Mental Health and Substance Use on Methylation Aging Into Adulthood.

OBJECTIVE: To test whether childhood mental health symptoms, substance use, and early adversity accelerate the rate of DNA methylation (DNAm) aging from adolescence to adulthood. METHOD: DNAm was assayed from blood samples in 381 participants in both adolescence (mean [SD] age = 13.9 [1.6] years) and adulthood (mean [SD] age = 25.9 [2.7] years). Structured diagnostic interviews were completed with participants and their parents at multiple childhood observations (1,950 total) to assess symptoms of common mental health disorders (attention-deficit/hyperactivity disorder, oppositional defiant disorder, conduct disorder, anxiety, and depression) and common types of substance use (alcohol, cannabis, nicotine) and early adversities. RESULTS: Neither childhood mental health symptoms nor substance use variables were associated with DNAm aging cross-sectionally. In contrast, the following mental health symptoms and substance variables were associated with accelerated DNAm aging from adolescence to adulthood: depressive symptoms (b = 0.314, SE = 0.127, p = .014), internalizing symptoms (b = 0.108, SE = 0.049, p = .029), weekly cannabis use (b =1.665, SE = 0.591, p = .005), and years of weekly cannabis use (b = 0.718, SE = 0.283, p = .012). In models testing all individual variables simultaneously, the combined effect of the variables was equivalent to a potential difference of 3.17 to 3.76 years in DNAm aging. A final model tested a variable assessing cumulative exposure to mental health symptoms, substance use, and early adversities. This cumulative variable was strongly associated with accelerated aging (b = 0.126, SE = 0.044, p = .005). CONCLUSION: Mental health symptoms and substance use accelerated DNAm aging into adulthood in a manner consistent with a shared risk mechanism.

Development and validation of a tumor immune cell infiltration-related gene signature for recurrence prediction by weighted gene co-expression network analysis in prostate cancer.

Introduction: Prostate cancer (PCa) is the second most common malignancy in men. Despite multidisciplinary treatments, patients with PCa continue to experience poor prognoses and high rates of tumor recurrence. Recent studies have shown that tumor-infiltrating immune cells (TIICs) are associated with PCa tumorigenesis. Methods: The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were used to derive multi-omics data for prostate adenocarcinoma (PRAD) samples. The CIBERSORT algorithm was used to calculate the landscape of TIICs. Weighted gene co-expression network analysis (WGCNA) was performed to determine the candidate module most significantly associated with TIICs. LASSO Cox regression was applied to screen a minimal set of genes and construct a TIIC-related prognostic gene signature for PCa. Then, 78 PCa samples with CIBERSORT output p-values of less than 0.05 were selected for analysis. WGCNA identified 13 modules, and the MEblue module with the most significant enrichment result was selected. A total of 1143 candidate genes were cross-examined between the MEblue module and active dendritic cell-related genes. Results: According to LASSO Cox regression analysis, a risk model was constructed with six genes (STX4, UBE2S, EMC6, EMD, NUCB1 and GCAT), which exhibited strong correlations with clinicopathological variables, tumor microenvironment context, antitumor therapies, and tumor mutation burden (TMB) in TCGA-PRAD. Further validation showed that the UBE2S had the highest expression level among the six genes in five different PCa cell lines. Discussion: In conclusion, our risk-score model contributes to better predicting PCa patient prognosis and understanding the underlying mechanisms of immune responses and antitumor therapies in PCa.

An innovative pyroptosis-related long-noncoding-RNA signature predicts the prognosis of gastric cancer via affecting immune cell infiltration landscape.

Background: Gastric cancer (GC) is a worldwide popular malignant tumor. However, the survival rate of advanced GC remains low. Pyroptosis and long non-coding RNAs (lncRNAs) are important in cancer progression. Thus, we aimed to find out a pyroptosis-related lncRNAs (PRLs) signature and use it to build a practical risk model with the purpose to predict the prognosis of patients with GC. Methods: Univariate Cox regression analysis was used to identify PRLs linked to GC patient's prognosis. Subsequently, to construct a PRLs signature, the least absolute shrinkage and selection operator regression, and multivariate Cox regression analysis were used. Kaplan-Meier analysis, principal component analysis, and receiver operating characteristic curve analysis were performed to assess our novel lncRNA signature. The correlation between risk signature and clinicopathological features was also examined. Finally, the relationship of pyroptosis and immune cells were evaluated through the CIBERSORT tool and single-sample lncRNA set enrichment analysis (ssGSEA). Results: A PRLs signature comprising eight lncRNAs was discerned as a self-determining predictor of prognosis. GC patients were sub-divided into high-risk and low-risk groups via this risk-model. Stratified analysis of different clinical factors also displayed that the PRLs signature was a good prognosis factor. According to the risk score and clinical characteristics, a nomogram was established. Moreover, the difference between the groups is significance in immune cells and immune pathways. Conclusion: This study established an effective prognostic signature consist of eight PRLs in GC, and constructed an efficient nomogram model. Further, the PRLs correlated with immune cells and immune pathways.

Genetic Alterations in Papillary Thyroid Carcinoma With Hashimoto's Thyroiditis： ANK3, an Indolent Maintainer of Papillary Thyroid Carcinoma.

Hashimoto's thyroiditis (TH) is a risk factor for the occurrence of papillary thyroid carcinoma (PTC), which is considered to be the most common type of thyroid cancer. In recent years, the prevalence of PTC with TH has been increasing, but little is known about the genetic alteration in PTC with TH. This study analyzed the mutation spectrum and mutation signature of somatic single nucleotide variants (SNV) for 10 non-tumor and tumor pair tissues of PTC with TH using whole-exome sequencing. The ANK3 protein expression was evaluated by immunohistochemistry in PTC with TH and PTC samples. Moreover, the functional role of ANK3 in PTC cells was determined by CCK-8 proliferation assay, colony formation assays, cell cycle analysis, cell invasion and migration and in vivo study through overexpression assay. Our results showed three distinct mutational signatures and the C>T/G>A substitution was the most common type of SNV. Gene-set enrichment analysis showed that most of the significantly mutated genes were enriched in the regulation of actin cytoskeleton signaling. Moreover, NCOR2, BPTF, ANK3, and PCSK5 were identified as the significantly mutated genes in PTC with TH, most of which have not been previously characterized. Unexpectedly, it was found that ANK3 was overexpressed in cytoplasm close to the membrane of PTC cells with TH and in almost all PTC cases, suggesting its role as a diagnostic marker of PTC. Ectopic expression of ANK3 suppressed invasion and migration, increased apoptosis of B-CPAP and TPC-1 cells. Moreover, our findings revealed that enhanced ANK3 expression inhibits growth of PTC cells both in vitro and in vivo. Ectopic expression of ANK3 significantly enhanced E-cadherin protein expression and inhibited PTC progression, at least in part, by suppression of epithelial-mesenchymal transition (EMT). Our study shows that ANK3 exerts an anti-oncogenic role in the development of PTC and might be an indolent maintainer of PTC.

A Prognostic Survival Model of Pancreatic Adenocarcinoma Based on Metabolism-Related Gene Expression.

Accurately predicting the survival prospects of patients suffering from pancreatic adenocarcinoma (PAAD) is challenging. In this study, we analyzed RNA matrices of 182 subjects with PAAD based on public datasets obtained from The Cancer Genome Atlas (TCGA) as training datasets and those of 63 subjects obtained from the Gene Expression Omnibus (GEO) database as the validation dataset. Genes regulating the metabolism of PAAD cells correlated with survival were identified. Furthermore, LASSO Cox regression analyses were conducted to identify six genes (XDH, MBOAT2, PTGES, AK4, PAICS, and CKB) to create a metabolic risk score. The proposed scoring framework attained the robust predictive performance, with 2-year survival areas under the curve (AUCs) of 0.61 in the training cohort and 0.66 in the validation cohort. Compared with the subjects in the low-risk cohort, subjects in the high-risk training cohort presented a worse survival outcome. The metabolic risk score increased the accuracy of survival prediction in patients suffering from PAAD.

An Early-Onset Advanced Rectal Cancer Patient With Increased KRAS Gene Copy Number Showed A Primary Resistance to Cetuximab in Combination With Chemotherapy: A Case Report.

Mutations in KRAS (codon 12/13), NRAS, BRAF V600E, and amplification of ERBB2 and MET account for 70-80% of anti-epidermal growth factor receptor (EGFR) monoclonal antibody primary resistance. However, the list of anti-EGFR monoclonal antibody primary resistance biomarkers is still incomplete. Herein, we report a case of wild-type RAS/BRAF metastatic colorectal cancer (CRC) with resistance to anti-EGFR monoclonal antibody and chemotherapy. Initially, mutation detection in postoperative tumor tissue by using amplification-refractory mutation system polymerase chain reaction indicated wild-type RAS/BRAF without point mutations, insertion deletions, or fusion mutations. Therefore, we recommended combined therapy of cetuximab and FOLFIRI after failure of platinum-based adjuvant chemotherapy, but the disease continued to progress. Next generation sequencing analysis of the postoperative tumor tissue revealed that KRAS copy number was increased and detected SMAD4, RNF43, and PREX2 mutations. This is the first case of advanced CRC with increased copy numbers of KRAS resistant to cetuximab and chemotherapy, which results in poor patient survival, and other mutated genes may be associated with the outcomes. Our findings indicate KRAS copy number alterations should also be examined, especially with anti-EGFR monoclonal antibody therapy in CRC, since it may be related with the primary resistance to these drugs.

Prognostic role and clinicopathological features of SMAD4 gene mutation in colorectal cancer: a systematic review and meta-analysis.

BACKGROUND: Approximately 5.0-24.2% of colorectal cancers (CRCs) have inactivating mutations in SMAD4, making it one of the frequently mutated genes in CRC. We thus carried out a comprehensive system review and meta-analysis investigating the prognostic significance and clinicopathological features of SMAD4 gene mutation in CRC patients. METHODS: A detailed literature search was conducted in PubMed, Web of Science and Embase databases to study the relationship between SMAD4 mutations and the demographic and clinicopathological characteristics in CRC patients. The hazard ratios (HRs) with 95% confidence intervals (CI) were used to evaluate the effect of SMAD4 mutations on overall survival (OS) and progression-free survival (PFS)/recurrence-free survival (RFS). RESULTS: Ten studies enrolling 4394 patients were eligible for inclusion. Data on OS were available from 5 studies and data on PFS/RFS were available from 3 studies. Comparing SMAD4-mutated CRC patients with SMAD4 wild-type CRC patients, the summary HR for OS was 1.46 (95% CI 1.28-1.67, P = 0.001), the summary HR for PFS/RFS was 1.59 (95% CI 1.14-2.22, P = 0.006). In terms of clinicopathology parameters, 9 studies have data that can be extracted, SMAD4 mutations were associated with tumor location (odds ratio [OR] = 1.15, colon/rectum, 95% CI 1.01-1.31, P = 0.042), TNM stage (OR = 1.28, stage IV/I-III, 95% CI 1.03-1.58, P = 0.025), lymph node metastasis (OR = 1.42, N1 + N2/N0, 95% CI 1.20-1.67, P < 0.001), mucinous differentiation (OR = 2.23, 95% CI 1.85-2.70, P < 0.001) and rat sarcoma viral oncogene homolog (RAS) mutation status (OR = 2.13, 95% CI 1.37-3.34, P = 0.001). No connection was found with age, gender, tumor grade, microsatellite instability status and b-viral oncogene homolog B1 mutation status. Besides, publication bias was not observed in any study. CONCLUSIONS: This meta-analysis suggests that SMAD4 mutation was associated with OS, PFS/RFS, and clinicopathological parameters, including tumor site, disease stage, RAS status, lymph node metastasis and mucinous differentiation. Our meta-analysis indicated that SMAD4 mutations could predict the poor prognosis and aggressive clinicopathological characteristics of CRC. More large-sample cohort studies are needed to confirm this conclusion. Since SMAD4 mutations are closely related to RAS mutations, their relationship warrants further investigation.

CD36 promotes the epithelial-mesenchymal transition and metastasis in cervical cancer by interacting with TGF-ß.

BACKGROUND: Accumulating evidence indicates that CD36 initiates metastasis and correlates with an unfavorable prognosis in cancers. However, there are few reports regarding the roles of CD36 in initiation and metastasis of cervical cancer. METHODS: Using immunohistochemistry, we analyzed 133 cervical cancer samples for CD36 protein expression levels, and then investigated the correlation between changes in its expression and clinicopathologic parameters. The effect of CD36 expression on the epithelial-mesenchymal transition (EMT) in cervical cancer cells was evaluated by Western immunoblotting analysis. In vitro invasion and in vivo metastasis assays were also used to evaluate the role of CD36 in cervical cancer metastasis. RESULTS: In the present study, we confirmed that CD36 was highly expressed in cervical cancer samples relative to normal cervical tissues. Moreover, overexpression of CD36 promoted invasiveness and metastasis of cervical cancer cells in vitro and in vivo, while CD36 knockdown suppressed proliferation, migration, and invasiveness. We demonstrated that TGF-ß treatment attenuated E-cadherin expression and enhanced the expression levels of CD36, vimentin, slug, snail, and twist in si-SiHa, si-HeLa, and C33a-CD36 cells, suggesting that TGF-ß synergized with CD36 on EMT via active CD36 expression. We also observed that the expression levels of TGF-ß in si-SiHa cells and si-HeLa cells were down-regulated, whereas the expression levels of TGF-ß were up-regulated in C33a-CD36 cells. These results imply that CD36 and TGF-ß interact with each other to promote the EMT in cervical cancer. CONCLUSIONS: Our findings suggest that CD36 is likely to be an effective target for guiding individualized clinical therapy of cervical cancer.

A Whole Methylome Study of Ethanol Exposure in Brain and Blood: An Exploration of the Utility of Peripheral Blood as Proxy Tissue for Brain in Alcohol Methylation Studies.

BACKGROUND: Recent reviews have highlighted the potential use of blood-based methylation biomarkers as diagnostic and prognostic tools of current and future alcohol use and addiction. Due to the substantial overlap that often exists between methylation patterns across different tissues, including blood and brain, blood-based methylation may track methylation changes in brain; however, little work has explored the overlap in alcohol-related methylation in these tissues. METHODS: To study the effects of alcohol on the brain methylome and identify possible biomarkers of these changes in blood, we performed a methylome-wide association study in brain and blood from 40 male DBA/2J mice that received either an acute ethanol (EtOH) or saline intraperitoneal injection. To investigate all 22 million CpGs in the mouse genome, we enriched for the methylated genomic fraction using methyl-CpG binding domain (MBD) protein capture followed by next-generation sequencing (MBD-seq). We performed association tests in blood and brain separately followed by enrichment testing to determine whether there was overlapping alcohol-related methylation in the 2 tissues. RESULTS: The top result for brain was a CpG located in an intron of Ttc39b (p = 5.65 × 10-08 ), and for blood, the top result was located in Espnl (p = 5.11 × 10-08 ). Analyses implicated pathways involved in inflammation and neuronal differentiation, such as CXCR4, IL-7, and Wnt signaling. Enrichment tests indicated significant overlap among the top results in brain and blood. Pathway analyses of the overlapping genes converge on MAPKinase signaling (p = 5.6 × 10-05 ) which plays a central role in acute and chronic responses to alcohol and glutamate receptor pathways, which can regulate neuroplastic changes underlying addictive behavior. CONCLUSIONS: Overall, we have shown some methylation changes in brain and blood after acute EtOH administration and that the changes in blood partly mirror the changes in brain suggesting the potential for DNA methylation in blood to be biomarkers of alcohol use.

Building a schizophrenia genetic network: transcription factor 4 regulates genes involved in neuronal development and schizophrenia risk.

The transcription factor 4 (TCF4) locus is a robust association finding with schizophrenia (SCZ), but little is known about the genes regulated by the encoded transcription factor. Therefore, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) of TCF4 in neural-derived (SH-SY5Y) cells to identify genome-wide TCF4 binding sites, followed by data integration with SCZ association findings. We identified 11 322 TCF4 binding sites overlapping in two ChIP-seq experiments. These sites are significantly enriched for the TCF4 Ebox binding motif (>85% having ≥1 Ebox) and implicate a gene set enriched for genes downregulated in TCF4 small-interfering RNA (siRNA) knockdown experiments, indicating the validity of our findings. The TCF4 gene set was also enriched among (1) gene ontology categories such as axon/neuronal development, (2) genes preferentially expressed in brain, in particular pyramidal neurons of the somatosensory cortex and (3) genes downregulated in postmortem brain tissue from SCZ patients (odds ratio, OR = 2.8, permutation P < 4x10-5). Considering genomic alignments, TCF4 binding sites significantly overlapped those for neural DNA-binding proteins such as FOXP2 and the SCZ-associated EP300. TCF4 binding sites were modestly enriched among SCZ risk loci from the Psychiatric Genomic Consortium (OR = 1.56, P = 0.03). In total, 130 TCF4 binding sites occurred in 39 of the 108 regions published in 2014. Thirteen genes within the 108 loci had both a TCF4 binding site ±10kb and were differentially expressed in siRNA knockdown experiments of TCF4, suggesting direct TCF4 regulation. These findings confirm TCF4 as an important regulator of neural genes and point toward functional interactions with potential relevance for SCZ.

Methyl-CpG-Binding Domain Sequencing: MBD-seq.

Detailed biological knowledge about the potential importance of the methylome is typically lacking for common diseases. Therefore, methylome-wide association studies (MWAS) are critical to detect disease relevant methylation sites. Methyl-CpG-binding domain sequencing (MBD-seq) offers potential advantages compared to antibody-based enrichment, but performance depends critically on using an optimal protocol. Using an optimized protocol, MBD-seq can approximate the sensitivity/specificity obtained with whole-genome bisulfite sequencing, but at a fraction of the costs and time to complete the project. Thus, MBD-seq offers a comprehensive first pass at the CpG methylome and is economically feasible with the samples sizes required for MWAS.

Deep Sequencing of 71 Candidate Genes to Characterize Variation Associated with Alcohol Dependence.

BACKGROUND: Previous genomewide association studies (GWASs) have identified a number of putative risk loci for alcohol dependence (AD). However, only a few loci have replicated and these replicated variants only explain a small proportion of AD risk. Using an innovative approach, the goal of this study was to generate hypotheses about potentially causal variants for AD that can be explored further through functional studies. METHODS: We employed targeted capture of 71 candidate loci and flanking regions followed by next-generation deep sequencing (mean coverage 78X) in 806 European Americans. Regions included in our targeted capture library were genes identified through published GWAS of alcohol, all human alcohol and aldehyde dehydrogenases, reward system genes including dopaminergic and opioid receptors, prioritized candidate genes based on previous associations, and genes involved in the absorption, distribution, metabolism, and excretion of drugs. We performed single-locus tests to determine if any single variant was associated with AD symptom count. Sets of variants that overlapped with biologically meaningful annotations were tested for association in aggregate. RESULTS: No single, common variant was significantly associated with AD in our study. We did, however, find evidence for association with several variant sets. Two variant sets were significant at the q-value <0.10 level: a genic enhancer for ADHFE1 (p = 1.47 × 10-5 ; q = 0.019), an alcohol dehydrogenase, and ADORA1 (p = 5.29 × 10-5 ; q = 0.035), an adenosine receptor that belongs to a G-protein-coupled receptor gene family. CONCLUSIONS: To our knowledge, this is the first sequencing study of AD to examine variants in entire genes, including flanking and regulatory regions. We found that in addition to protein coding variant sets, regulatory variant sets may play a role in AD. From these findings, we have generated initial functional hypotheses about how these sets may influence AD.

AKAP-9 promotes colorectal cancer development by regulating Cdc42 interacting protein 4.

Our previous studies have shown that PRKA kinase anchor protein 9 (AKAP-9) is involved in colorectal cancer (CRC) cell proliferation and migration in vitro. However, whether or not AKAP-9 is important for CRC development or metastasis in vivo remains unknown. In the present study, we found that AKAP-9 expression was significantly higher in human colorectal cancer tissues than the paired normal tissues. In fact, AKAP-9 level correlated with the CRC infiltrating depth and metastasis. Moreover, the higher AKAP-9 expression was associated with the lower survival rate in patients. In cultured CRC cells, knockdown of AKAP-9 inhibited cell proliferation, invasion, and migration. AKAP-9 deficiency also attenuated CRC tumor growth and metastasis in vivo. Mechanistically, AKAP-9 interacted with cdc42 interacting protein 4 (CIP4) and regulated its expression. CIP4 levels were interrelated to the AKAP-9 level in CRC cells. Functionally, AKAP-9 was essential for TGF-ß1-induced epithelial-mesenchymal transition of CRC cells, and CIP4 played a critical role in mediating the function of AKAP-9. Importantly, CIP4 expression was significantly up-regulated in human CRC tissues. Taken together, our results demonstrated that AKAP-9 facilitates CRC development and metastasis via regulating CIP4-mediated epithelial-mesenchymal transition of CRC cells.

Long non-coding RNA MALAT1 increases AKAP-9 expression by promoting SRPK1-catalyzed SRSF1 phosphorylation in colorectal cancer cells.

Our earlier findings indicate that the long non-coding RNA MALAT1 promotes colorectal cancer (CRC) cell proliferation, invasion and metastasis in vitro and in vivo by increasing expression of AKAP-9. In the present study, we investigated the molecular mechanism by which MALAT1 enhances AKAP9 expression in CRC SW480 cells. We found that MALAT1 interacts with both SRPK1 and SRSF1. MALAT1 increases AKAP-9 expression by promoting SRPK1-catalyzed SRSF1 phosphorylation. Following MALAT1 knockdown, overexpression of SRPK1 was sufficient to restore SRSF1 phosphorylation and AKAP-9 expression to a level that promoted cell proliferation, invasion and migration in vitro. Conversely, SRPK1 knockdown after overexpression of MALAT1 in SW480 cells diminished SRSF1 phosphorylation and AKAP-9 expression and suppressed cell proliferation, invasion and migration in vitro. These findings suggest MALAT1 increases AKAP-9 expression by promoting SRPK1-catalyzed SRSF1 phosphorylation in CRC cells. These results reveal a novel molecular mechanism by which MALAT1 regulates AKAP-9 expression in CRC cells.

Deep Sequencing of Three Loci Implicated in Large-Scale Genome-Wide Association Study Smoking Meta-Analyses.

INTRODUCTION: Genome-wide association study meta-analyses have robustly implicated three loci that affect susceptibility for smoking: CHRNA5\CHRNA3\CHRNB4, CHRNB3\CHRNA6 and EGLN2\CYP2A6. Functional follow-up studies of these loci are needed to provide insight into biological mechanisms. However, these efforts have been hampered by a lack of knowledge about the specific causal variant(s) involved. In this study, we prioritized variants in terms of the likelihood they account for the reported associations. METHODS: We employed targeted capture of the CHRNA5\CHRNA3\CHRNB4, CHRNB3\CHRNA6, and EGLN2\CYP2A6 loci and flanking regions followed by next-generation deep sequencing (mean coverage 78×) to capture genomic variation in 363 individuals. We performed single locus tests to determine if any single variant accounts for the association, and examined if sets of (rare) variants that overlapped with biologically meaningful annotations account for the associations. RESULTS: In total, we investigated 963 variants, of which 71.1% were rare (minor allele frequency < 0.01), 6.02% were insertion/deletions, and 51.7% were catalogued in dbSNP141. The single variant results showed that no variant fully accounts for the association in any region. In the variant set results, CHRNB4 accounts for most of the signal with significant sets consisting of directly damaging variants. CHRNA6 explains most of the signal in the CHRNB3\CHRNA6 locus with significant sets indicating a regulatory role for CHRNA6. Significant sets in CYP2A6 involved directly damaging variants while the significant variant sets suggested a regulatory role for EGLN2. CONCLUSIONS: We found that multiple variants implicating multiple processes explain the signal. Some variants can be prioritized for functional follow-up.

Combined Whole Methylome and Genomewide Association Study Implicates CNTN4 in Alcohol Use.

BACKGROUND: Methylome-wide association (MWAS) studies present a new way to advance the search for biological correlates for alcohol use. A challenge with methylation studies of alcohol involves the causal direction of significant methylation-alcohol associations. One way to address this issue is to combine MWAS data with genomewide association study (GWAS) data. METHODS: Here, we combined MWAS and GWAS results for alcohol use from 619 individuals. Our MWAS data were generated by next-generation sequencing of the methylated genomic DNA fraction, producing over 60 million reads per subject to interrogate methylation levels at ~27 million autosomal CpG sites in the human genome. Our GWAS included 5,571,786 single nucleotide polymorphisms (SNPs) imputed with 1000 Genomes. RESULTS: When combining the MWAS and GWAS data, our top finding was a region in an intron of CNTN4 (p = 2.55 × 10(-8) ), located between chr3: 2,555,403 and 2,555,524, encompassing SNPs rs1382874 and rs1382875. This finding was then replicated in an independent sample of 730 individuals. We used bisulfite pyrosequencing to measure methylation and found significant association with regular alcohol use in the same direction as the MWAS (p = 0.021). Rs1382874 and rs1382875 were genotyped and found to be associated in the same direction as the GWAS (p = 0.008 and p = 0.009). After integrating the MWAS and GWAS findings from the replication sample, we replicated our combined analysis finding (p = 0.0017) in CNTN4. CONCLUSIONS: Through combining methylation and SNP data, we have identified CNTN4 as a risk factor for regular alcohol use.

Evaluation of Methyl-Binding Domain Based Enrichment Approaches Revisited.

Methyl-binding domain (MBD) enrichment followed by deep sequencing (MBD-seq), is a robust and cost efficient approach for methylome-wide association studies (MWAS). MBD-seq has been demonstrated to be capable of identifying differentially methylated regions, detecting previously reported robust associations and producing findings that replicate with other technologies such as targeted pyrosequencing of bisulfite converted DNA. There are several kits commercially available that can be used for MBD enrichment. Our previous work has involved MethylMiner (Life Technologies, Foster City, CA, USA) that we chose after careful investigation of its properties. However, in a recent evaluation of five commercially available MBD-enrichment kits the performance of the MethylMiner was deemed poor. Given our positive experience with MethylMiner, we were surprised by this report. In an attempt to reproduce these findings we here have performed a direct comparison of MethylMiner with MethylCap (Diagenode Inc, Denville, NJ, USA), the best performing kit in that study. We find that both MethylMiner and MethylCap are two well performing MBD-enrichment kits. However, MethylMiner shows somewhat better enrichment efficiency and lower levels of background "noise". In addition, for the purpose of MWAS where we want to investigate the majority of CpGs, we find MethylMiner to be superior as it allows tailoring the enrichment to the regions where most CpGs are located. Using targeted bisulfite sequencing we confirmed that sites where methylation was detected by either MethylMiner or by MethylCap indeed were methylated.

Refinement of schizophrenia GWAS loci using methylome-wide association data.

Recent genome-wide association studies (GWAS) have made substantial progress in identifying disease loci. The next logical step is to design functional experiments to identify disease mechanisms. This step, however, is often hampered by the large size of loci identified in GWAS that is caused by linkage disequilibrium between SNPs. In this study, we demonstrate how integrating methylome-wide association study (MWAS) results with GWAS findings can narrow down the location for a subset of the putative casual sites. We use the disease schizophrenia as an example. To handle "data analytic" variation, we first combined our MWAS results with two GWAS meta-analyses (N = 32,143 and 21,953), that had largely overlapping samples but different data analysis pipelines, separately. Permutation tests showed significant overlapping association signals between GWAS and MWAS findings. This significant overlap justified prioritizing loci based on the concordance principle. To further ensure that the methylation signal was not driven by chance, we successfully replicated the top three methylation findings near genes SDCCAG8, CREB1 and ATXN7 in an independent sample using targeted pyrosequencing. In contrast to the SNPs in the selected region, the methylation sites were largely uncorrelated explaining why the methylation signals implicated much smaller regions (median size 78 bp). The refined loci showed considerable enrichment of genomic elements of possible functional importance and suggested specific hypotheses about schizophrenia etiology. Several hypotheses involved possible variation in transcription factor-binding efficiencies.

MALAT1 promotes colorectal cancer cell proliferation/migration/invasion via PRKA kinase anchor protein 9.

Our previous studies have shown that the 3' end of metastasis associated lung adenocarcinoma transcript 1 (MALAT1) is involved in colorectal cancer (CRC) cell proliferation and migration/invasion in vitro. The role and mechanism of MALAT1 in CRC metastasis in vivo, however, remain largely unknown. In the present study, we found that MALAT1 was up-regulated in human primary CRC tissues with lymph node metastasis. Overexpression of MALAT1 via RNA activation promoted CRC cell proliferation, invasion and migration in vitro, and stimulated tumor growth and metastasis in mice in vivo. Conversely, knockdown of MALAT1 inhibited CRC tumor growth and metastasis. MALAT1 regulated at least 243 genes in CRC cells in a genome-wide expression profiling. Among these genes, PRKA kinase anchor protein 9 (AKAP-9) was significantly up-regulated at both mRNA and protein levels. AKAP-9 was highly expressed in CRC cells with metastatic potential and human primary CRC tissues with lymph node metastasis, but not in normal cells or tissues. Importantly, knockdown of AKAP-9 blocked MALAT1-mediated CRC cell proliferation, migration and invasion. These data indicate that MALAT1 may promote CRC tumor development via its target protein AKAP-9.

Humanizing mouse folate metabolism: conversion of the dual-promoter mouse folylpolyglutamate synthetase gene to the human single-promoter structure.

The mouse is extensively used to model human folate metabolism and therapeutic outcomes with antifolates. However, the folylpoly-γ-glutamate synthetase (fpgs) gene, whose product determines folate/antifolate intracellular retention and antifolate antitumor activity, displays a pronounced species difference. The human gene uses only a single promoter, whereas the mouse uses two: P2, akin to the human promoter, at low levels in most tissues; and P1, an upstream promoter used extensively in liver and kidney. We deleted the mouse P1 promoter through homologous recombination to study the dual-promoter mouse system and to create a mouse with a humanized fpgs gene structure. Despite the loss of the predominant fpgs mRNA species in liver and kidney (representing 95 and 75% of fpgs transcripts in these tissues, respectively), P1-knockout mice developed and reproduced normally. The survival of these mice was explained by increased P2 transcription due to relief of transcriptional interference, by a 3-fold more efficient translation of P2-derived than P1-derived transcripts, and by 2-fold higher stability of P2-derived FPGS. In combination, all 3 effects reinstated FPGS function, even in liver. By eliminating mouse P1, we created a mouse model that mimicked the human housekeeping pattern of fpgs gene expression.