Cross-trait prediction accuracy of summary statistics in genome-wide association studies.

In the era of big data, univariate models have widely been used as a workhorse tool for quickly producing marginal estimators; and this is true even when in a high-dimensional dense setting, in which many features are "true," but weak signals. Genome-wide association studies (GWAS) epitomize this type of setting. Although the GWAS marginal estimator is popular, it has long been criticized for ignoring the correlation structure of genetic variants (i.e., the linkage disequilibrium [LD] pattern). In this paper, we study the effects of LD pattern on the GWAS marginal estimator and investigate whether or not additionally accounting for the LD can improve the prediction accuracy of complex traits. We consider a general high-dimensional dense setting for GWAS and study a class of ridge-type estimators, including the popular marginal estimator and the best linear unbiased prediction (BLUP) estimator as two special cases. We show that the performance of GWAS marginal estimator depends on the LD pattern through the first three moments of its eigenvalue distribution. Furthermore, we uncover that the relative performance of GWAS marginal and BLUP estimators highly depends on the ratio of GWAS sample size over the number of genetic variants. Particularly, our finding reveals that the marginal estimator can easily become near-optimal within this class when the sample size is relatively small, even though it ignores the LD pattern. On the other hand, BLUP estimator has substantially better performance than the marginal estimator as the sample size increases toward the number of genetic variants, which is typically in millions. Therefore, adjusting for the LD (such as in the BLUP) is most needed when GWAS sample size is large. We illustrate the importance of our results by using the simulated data and real GWAS.

Maximal Information Coefficient-Based Testing to Identify Epistasis in Case-Control Association Studies.

Interactions between genetic variants (epistasis) are ubiquitous in the model system and can significantly affect evolutionary adaptation, genetic mapping, and precision medical efforts. In this paper, we proposed a method for epistasis detection, called EpiMIC (epistasis detection through a maximal information coefficient (MIC)). MIC is a promising bivariate dependence measure explicitly designed for rapidly exploring various function types equally and for interpreting and comparing them on the same scale. Most epistasis detection approaches make assumptions about the form of the association between genetic variants, resulting in limited statistical performance. Based on the notion that if two SNPs do not interact, their joint distribution in all samples and in only cases should not be substantially different. We developed a statistic that utilizes the difference of MIC as a signal of epistasis and combined it with a permutation resampling strategy to estimate the empirical distribution of our statistic. Results of simulation and real-world data set showed that EpiMIC outperformed previous approaches for identifying epistasis at varying degrees of heredity.

Aberrant epigenetic and transcriptional events associated with breast cancer risk.

BACKGROUND: Genome-wide association studies have identified several breast cancer susceptibility loci. However, biomarkers for risk assessment are still missing. Here, we investigated cancer-related molecular changes detected in tissues from women at high risk for breast cancer prior to disease manifestation. Disease-free breast tissue cores donated by healthy women (N = 146, median age = 39 years) were processed for both methylome (MethylCap) and transcriptome (Illumina's HiSeq4000) sequencing. Analysis of tissue microarray and primary breast epithelial cells was used to confirm gene expression dysregulation. RESULTS: Transcriptomic analysis identified 69 differentially expressed genes between women at high and those at average risk of breast cancer (Tyrer-Cuzick model) at FDR < 0.05 and fold change ≥ 2. Majority of the identified genes were involved in DNA damage checkpoint, cell cycle, and cell adhesion. Two genes, FAM83A and NEK2, were overexpressed in tissue sections (FDR < 0.01) and primary epithelial cells (p < 0.05) from high-risk breasts. Moreover, 1698 DNA methylation changes were identified in high-risk breast tissues (FDR < 0.05), partially overlapped with cancer-related signatures, and correlated with transcriptional changes (p < 0.05, r ≤ 0.5). Finally, among the participants, 35 women donated breast biopsies at two time points, and age-related molecular alterations enhanced in high-risk subjects were identified. CONCLUSIONS: Normal breast tissue from women at high risk of breast cancer bears molecular aberrations that may contribute to breast cancer susceptibility. This study is the first molecular characterization of the true normal breast tissues, and provides an opportunity to investigate molecular markers of breast cancer risk, which may lead to new preventive approaches.

Prediction of New Risk Genes and Potential Drugs for Rheumatoid Arthritis from Multiomics Data.

Rheumatoid arthritis (RA) is an autoimmune and inflammatory disease for which there is a lack of therapeutic options. Genome-wide association studies (GWASs) have identified over 100 genetic loci associated with RA susceptibility; however, the most causal risk genes (RGs) associated with, and molecular mechanism underlying, RA remain unknown. In this study, we collected 95 RA-associated loci from multiple GWASs and detected 87 candidate high-confidence risk genes (HRGs) from these loci via integrated multiomics data (the genome-scale chromosome conformation capture data, enhancer-promoter linkage data, and gene expression data) using the Bayesian integrative risk gene selector (iRIGS). Analysis of these HRGs indicates that these genes were indeed, markedly associated with different aspects of RA. Among these, 36 and 46 HRGs have been reported to be related to RA and autoimmunity, respectively. Meanwhile, most novel HRGs were also involved in the significantly enriched RA-related biological functions and pathways. Furthermore, drug repositioning prediction of the HRGs revealed three potential targets (ERBB2, IL6ST, and MAPK1) and nine possible drugs for RA treatment, of which two IL-6 receptor antagonists (tocilizumab and sarilumab) have been approved for RA treatment and four drugs (trastuzumab, lapatinib, masoprocol, and arsenic trioxide) have been reported to have a high potential to ameliorate RA. In summary, we believe that this study provides new clues for understanding the pathogenesis of RA and is important for research regarding the mechanisms underlying RA and the development of therapeutics for this condition.

SNP characteristics and validation success in genome wide association studies.

Genome wide association studies (GWASs) have identified tens of thousands of single nucleotide polymorphisms (SNPs) associated with human diseases and characteristics. A significant fraction of GWAS findings can be false positives. The gold standard for true positives is an independent validation. The goal of this study was to identify SNP features associated with validation success. Summary statistics from the Catalog of Published GWASs were used in the analysis. Since our goal was an analysis of reproducibility, we focused on the diseases/phenotypes targeted by at least 10 GWASs. GWASs were arranged in discovery-validation pairs based on the time of publication, with the discovery GWAS published before validation. We used four definitions of the validation success that differ by stringency. Associations of SNP features with validation success were consistent across the definitions. The strongest predictor of SNP validation was the level of statistical significance in the discovery GWAS. The magnitude of the effect size was associated with validation success in a non-linear manner. SNPs with risk allele frequencies in the range 30-70% showed a higher validation success rate compared to rarer or more common SNPs. Missense, 5'UTR, stop gained, and SNPs located in transcription factor binding sites had a higher validation success rate compared to intergenic, intronic and synonymous SNPs. There was a positive association between validation success and the level of evolutionary conservation of the sites. In addition, validation success was higher when discovery and validation GWASs targeted the same ethnicity. All predictors of validation success remained significant in a multivariate logistic regression model indicating their independent contribution. To conclude, we identified SNP features predicting validation success of GWAS hits. These features can be used to select SNPs for validation and downstream functional studies.

A genetic risk score and diabetes predict development of alcohol-related cirrhosis in drinkers.

BACKGROUND & AIMS: Only a minority of excess alcohol drinkers develop cirrhosis. We developed and evaluated risk stratification scores to identify those at highest risk. METHODS: Three cohorts (GenomALC-1: n = 1,690, GenomALC-2: n = 3,037, UK Biobank: relevant n = 6,898) with a history of heavy alcohol consumption (≥80 g/day (men), ≥50 g/day (women), for ≥10 years) were included. Cases were participants with alcohol-related cirrhosis. Controls had a history of similar alcohol consumption but no evidence of liver disease. Risk scores were computed from up to 8 genetic loci identified previously as associated with alcohol-related cirrhosis and 3 clinical risk factors. Score performance for the stratification of alcohol-related cirrhosis risk was assessed and compared across the alcohol-related liver disease spectrum, including hepatocellular carcinoma (HCC). RESULTS: A combination of 3 single nucleotide polymorphisms (SNPs) (PNPLA3:rs738409, SUGP1-TM6SF2:rs10401969, HSD17B13:rs6834314) and diabetes status best discriminated cirrhosis risk. The odds ratios (ORs) and (95% CIs) between the lowest (Q1) and highest (Q5) score quintiles of the 3-SNP score, based on independent allelic effect size estimates, were 5.99 (4.18-8.60) (GenomALC-1), 2.81 (2.03-3.89) (GenomALC-2), and 3.10 (2.32-4.14) (UK Biobank). Patients with diabetes and high risk scores had ORs of 14.7 (7.69-28.1) (GenomALC-1) and 17.1 (11.3-25.7) (UK Biobank) compared to those without diabetes and with low risk scores. Patients with cirrhosis and HCC had significantly higher mean risk scores than patients with cirrhosis alone (0.76 ± 0.06 vs. 0.61 ± 0.02, p = 0.007). Score performance was not significantly enhanced by information on additional genetic risk variants, body mass index or coffee consumption. CONCLUSIONS: A risk score based on 3 genetic risk variants and diabetes status enables the stratification of heavy drinkers based on their risk of cirrhosis, allowing for the provision of earlier preventative interventions. LAY SUMMARY: Excessive chronic drinking leads to cirrhosis in some people, but so far there is no way to identify those at high risk of developing this debilitating disease. We developed a genetic risk score that can identify patients at high risk. The risk of cirrhosis is increased >10-fold with just two risk factors - diabetes and a high genetic risk score. Risk assessment using this test could enable the early and personalised management of this disease in high-risk patients.

Rare GATA6 variants associated with risk of congenital heart disease phenotypes in 200,000 UK Biobank exomes.

Congenital heart disease (CHD) has a complex and largely uncharacterised genetic etiology. Using 200,000 UK Biobank (UKB) exomes, we assess the burden of ultra-rare, potentially pathogenic variants in the largest case/control cohort of predominantly mild CHD to date. We find an association with GATA6, a member of the GATA family of transcription factors that play an important role during heart development and has been linked with several CHD phenotypes previously. Several identified GATA6 variants are previously unreported and their roles in conferring risk to CHD warrants further study. We demonstrate that despite limitations regarding detailed familial phenotype information in large-scale biobank projects, through careful consideration of case and control cohorts it is possible to derive important associations.

Analyses of biomarker traits in diverse UK biobank participants identify associations missed by European-centric analysis strategies.

Despite the dramatic underrepresentation of non-European populations in human genetics studies, researchers continue to exclude participants of non-European ancestry, as well as variants rare in European populations, even when these data are available. This practice perpetuates existing research disparities and can lead to important and large effect size associations being missed. Here, we conducted genome-wide association studies (GWAS) of 31 serum and urine biomarker quantitative traits in African (n = 9354), East Asian (n = 2559), and South Asian (n = 9823) ancestry UK Biobank (UKBB) participants. We adjusted for all known GWAS catalog variants for each trait, as well as novel signals identified in a recent European ancestry-focused analysis of UKBB participants. We identify 7 novel signals in African ancestry and 2 novel signals in South Asian ancestry participants (p < 1.61E-10). Many of these signals are highly plausible, including a cis pQTL for the gene encoding gamma-glutamyl transferase and PIEZO1 and G6PD variants with impacts on HbA1c through likely erythrocytic mechanisms. This work illustrates the importance of using the genetic data we already have in diverse populations, with novel discoveries possible in even modest sample sizes.

The genetics of testosterone contributes to "femaleness/maleness" of cardiometabolic traits and type 2 diabetes.

The genetic architecture of testosterone is highly distinct between sexes. Moreover, obesity is associated with higher testosterone in females but lower testosterone in males. Here, we ask whether male-specific testosterone variants are associated with a male pattern of obesity and type 2 diabetes (T2D) in females, and vice versa. In the UK Biobank, we conducted sex-specific genome-wide association studies and computed polygenic scores for total (PGSTT) and bioavailable testosterone (PGSBT). We tested sex-congruent and sex-incongruent associations between sex-specific PGSTs and metabolic traits, as well as T2D diagnosis. Female-specific PGSBT was associated with an elevated cardiometabolic risk and probability of T2D, in both sexes. Male-specific PGSTT was associated with traits conferring a lower cardiometabolic risk and probability of T2D, in both sexes. We demonstrate the value in considering polygenic testosterone as sex-related continuous traits, in each sex.

Accurate prediction of acute pancreatitis severity based on genome-wide cell free DNA methylation profiles.

BACKGROUND: Patients with severe acute pancreatitis (SAP) have a high mortality, thus early diagnosis and interventions are critical for improving survival. However, conventional tests are limited in acute pancreatitis (AP) stratification. We aimed to assess AP severity by integrating the informative clinical measurements with cell free DNA (cfDNA) methylation markers. METHODS: One hundred and seventy-five blood samples were collected from 61 AP patients at multiple time points, plus 24 samples from healthy individuals. Genome-wide cfDNA methylation profiles of all samples were characterized with reduced representative bisulfite sequencing. Clinical blood tests covering 93 biomarkers were performed on AP patients within 24 h. SAP predication models were built based on cfDNA methylation and conventional blood biomarkers separately and in combination. RESULTS: We identified 565 and 59 cfDNA methylation markers informative for acute pancreatitis and its severity. These markers were used to develop prediction models for AP and SAP with area under the receiver operating characteristic of 0.92 and 0.81, respectively. Twelve blood biomarkers were systematically screened for a predictor of SAP with a sensitivity of 87.5% for SAP, and a specificity of 100% in mild acute pancreatitis, significantly higher than existing blood tests. An expanded model integrating 12 conventional blood biomarkers with 59 cfDNA methylation markers further improved the SAP prediction sensitivity to 92.2%. CONCLUSIONS: These findings have demonstrated that accurate prediction of SAP by the integration of conventional and novel blood molecular markers, paving the way for early and effective SAP intervention through a non-invasive rapid diagnostic test.

Genome-wide methylation patterns in Marfan syndrome.

BACKGROUND: Marfan syndrome (MFS) is a connective tissue disorder caused by mutations in the Fibrillin-1 gene (FBN1). Here, we undertook the first epigenome-wide association study (EWAS) in patients with MFS aiming at identifying DNA methylation loci associated with MFS phenotypes that may shed light on the disease process. METHODS: The Illumina 450 k DNA-methylation array was used on stored peripheral whole-blood samples of 190 patients with MFS originally included in the COMPARE trial. An unbiased genome-wide approach was used, and methylation of CpG-sites across the entire genome was evaluated. Additionally, we investigated CpG-sites across the FBN1-locus (15q21.1) more closely, since this is the gene defective in MFS. Differentially Methylated Positions (DMPs) and Differentially Methylated Regions (DMRs) were identified through regression analysis. Associations between methylation levels and aortic diameters and presence or absence of 21 clinical features of MFS at baseline were analyzed. Moreover, associations between aortic diameter change, and the occurrence of clinical events (death any cause, type-A or -B dissection/rupture, or aortic surgery) and methylation levels were analyzed. RESULTS: We identified 28 DMPs that are significantly associated with aortic diameters in patients with MFS. Seven of these DMPs (25%) could be allocated to a gene that was previously associated with cardiovascular diseases (HDAC4, IGF2BP3, CASZ1, SDK1, PCDHGA1, DIO3, PTPRN2). Moreover, we identified seven DMPs that were significantly associated with aortic diameter change and five DMP's that associated with clinical events. No significant associations at p < 10-8 or p < 10-6 were found with any of the non-cardiovascular phenotypic MFS features. Investigating DMRs, clusters were seen mostly on X- and Y, and chromosome 18-22. The remaining DMRs indicated involvement of a large family of protocadherins on chromosome 5, which were not reported in MFS before. CONCLUSION: This EWAS in patients with MFS has identified a number of methylation loci significantly associated with aortic diameters, aortic dilatation rate and aortic events. Our findings add to the slowly growing literature on the regulation of gene expression in MFS patients.

Simultaneous estimation of bi-directional causal effects and heritable confounding from GWAS summary statistics.

Mendelian Randomisation (MR) is an increasingly popular approach that estimates the causal effect of risk factors on complex human traits. While it has seen several extensions that relax its basic assumptions, most suffer from two major limitations; their under-exploitation of genome-wide markers, and sensitivity to the presence of a heritable confounder of the exposure-outcome relationship. To overcome these limitations, we propose a Latent Heritable Confounder MR (LHC-MR) method applicable to association summary statistics, which estimates bi-directional causal effects, direct heritabilities, and confounder effects while accounting for sample overlap. We demonstrate that LHC-MR outperforms several existing MR methods in a wide range of simulation settings and apply it to summary statistics of 13 complex traits. Besides several concordant results with other MR methods, LHC-MR unravels new mechanisms (how disease diagnosis might lead to improved lifestyle) and reveals new causal effects (e.g. HDL cholesterol being protective against high systolic blood pressure), hidden from standard MR methods due to a heritable confounder of opposite effect direction.

CDK5 Regulatory Subunit-Associated Protein 1-Like 1 Gene Polymorphisms and Gestational Diabetes Mellitus Risk: A Trial Sequential Meta-Analysis of 13,306 Subjects.

Objectives: The CDK5 regulatory subunit-associated protein 1-like 1 (CDKAL1) contributes to islet ß-cell function and insulin secretion by inhibiting the activation of CDK5. The current studies on the relationship between CDKAL1 polymorphisms rs7756992 A>G and rs7754840 C>G and the risk of gestational diabetes mellitus (GDM) have drawn contradictory conclusions. Materials and Methods: A meta-analysis with a fixed- or random-effects model was conducted to estimate the correlation between studied CDKAL1 polymorphisms and GDM risk with the summary odds ratio (OR) and 95% confidence interval (CI). In addition, trial sequential analysis (TSA) and false-positive report probability (FPRP) analysis were performed to confirm the study findings. Results: A total of 13,306 subjects were included in the present study. Meta-analysis results showed that the variant heterozygous and homozygous genotypes of the two polymorphisms were associated with increased GDM risk in comparison with the wild-type AA genotype (AG vs. AA: OR = 1.23, 95% CI = 1.08, 1.41, p = 0.002; GG vs. AA: OR = 1.47, 95% CI = 1.05, 2.05, p = 0.024 for rs7756992; and CG vs. GG: OR = 1.36, 95% CI = 1.13, 1.65, p = 0.002; CC vs. GG: OR = 1.76, 95% CI = 1.37, 2.26, p < 0.001 for rs7754840). The TSA confirmed a significant association between rs7754840 and the susceptibility to GDM because the cumulative Z-curve crossed both the conventional cutoff value and the TSA boundaries under the heterozygote and homozygote models. Conclusions: This study supported the finding that rs7756992 and rs7754840 are associated with susceptibility to GDM. However, further functional studies are warranted to clarify the mechanism.

Combining the strengths of inverse-variance weighting and Egger regression in Mendelian randomization using a mixture of regressions model.

With the increasing availability of large-scale GWAS summary data on various traits, Mendelian randomization (MR) has become commonly used to infer causality between a pair of traits, an exposure and an outcome. It depends on using genetic variants, typically SNPs, as instrumental variables (IVs). The inverse-variance weighted (IVW) method (with a fixed-effect meta-analysis model) is most powerful when all IVs are valid; however, when horizontal pleiotropy is present, it may lead to biased inference. On the other hand, Egger regression is one of the most widely used methods robust to (uncorrelated) pleiotropy, but it suffers from loss of power. We propose a two-component mixture of regressions to combine and thus take advantage of both IVW and Egger regression; it is often both more efficient (i.e. higher powered) and more robust to pleiotropy (i.e. controlling type I error) than either IVW or Egger regression alone by accounting for both valid and invalid IVs respectively. We propose a model averaging approach and a novel data perturbation scheme to account for uncertainties in model/IV selection, leading to more robust statistical inference for finite samples. Through extensive simulations and applications to the GWAS summary data of 48 risk factor-disease pairs and 63 genetically uncorrelated trait pairs, we showcase that our proposed methods could often control type I error better while achieving much higher power than IVW and Egger regression (and sometimes than several other new/popular MR methods). We expect that our proposed methods will be a useful addition to the toolbox of Mendelian randomization for causal inference.

A generalized linear mixed model association tool for biobank-scale data.

Compared with linear mixed model-based genome-wide association (GWA) methods, generalized linear mixed model (GLMM)-based methods have better statistical properties when applied to binary traits but are computationally much slower. In the present study, leveraging efficient sparse matrix-based algorithms, we developed a GLMM-based GWA tool, fastGWA-GLMM, that is severalfold to orders of magnitude faster than the state-of-the-art tools when applied to the UK Biobank (UKB) data and scalable to cohorts with millions of individuals. We show by simulation that the fastGWA-GLMM test statistics of both common and rare variants are well calibrated under the null, even for traits with extreme case-control ratios. We applied fastGWA-GLMM to the UKB data of 456,348 individuals, 11,842,647 variants and 2,989 binary traits (full summary statistics available at http://fastgwa.info/ukbimpbin ), and identified 259 rare variants associated with 75 traits, demonstrating the use of imputed genotype data in a large cohort to discover rare variants for binary complex traits.

An open approach to systematically prioritize causal variants and genes at all published human GWAS trait-associated loci.

Genome-wide association studies (GWASs) have identified many variants associated with complex traits, but identifying the causal gene(s) is a major challenge. In the present study, we present an open resource that provides systematic fine mapping and gene prioritization across 133,441 published human GWAS loci. We integrate genetics (GWAS Catalog and UK Biobank) with transcriptomic, proteomic and epigenomic data, including systematic disease-disease and disease-molecular trait colocalization results across 92 cell types and tissues. We identify 729 loci fine mapped to a single-coding causal variant and colocalized with a single gene. We trained a machine-learning model using the fine-mapped genetics and functional genomics data and 445 gold-standard curated GWAS loci to distinguish causal genes from neighboring genes, outperforming a naive distance-based model. Our prioritized genes were enriched for known approved drug targets (odds ratio = 8.1, 95% confidence interval = 5.7, 11.5). These results are publicly available through a web portal ( http://genetics.opentargets.org ), enabling users to easily prioritize genes at disease-associated loci and assess their potential as drug targets.

Validating biomarkers and models for epigenetic inference of alcohol consumption from blood.

BACKGROUND: Information on long-term alcohol consumption is relevant for medical and public health research, disease therapy, and other areas. Recently, DNA methylation-based inference of alcohol consumption from blood was reported with high accuracy, but these results were based on employing the same dataset for model training and testing, which can lead to accuracy overestimation. Moreover, only subsets of alcohol consumption categories were used, which makes it impossible to extrapolate such models to the general population. By using data from eight population-based European cohorts (N = 4677), we internally and externally validated the previously reported biomarkers and models for epigenetic inference of alcohol consumption from blood and developed new models comprising all data from all categories. RESULTS: By employing data from six European cohorts (N = 2883), we empirically tested the reproducibility of the previously suggested biomarkers and prediction models via ten-fold internal cross-validation. In contrast to previous findings, all seven models based on 144-CpGs yielded lower mean AUCs compared to the models with less CpGs. For instance, the 144-CpG heavy versus non-drinkers model gave an AUC of 0.78 ± 0.06, while the 5 and 23 CpG models achieved 0.83 ± 0.05, respectively. The transportability of the models was empirically tested via external validation in three independent European cohorts (N = 1794), revealing high AUC variance between datasets within models. For instance, the 144-CpG heavy versus non-drinkers model yielded AUCs ranging from 0.60 to 0.84 between datasets. The newly developed models that considered data from all categories showed low AUCs but gave low AUC variation in the external validation. For instance, the 144-CpG heavy and at-risk versus light and non-drinkers model achieved AUCs of 0.67 ± 0.02 in the internal cross-validation and 0.61-0.66 in the external validation datasets. CONCLUSIONS: The outcomes of our internal and external validation demonstrate that the previously reported prediction models suffer from both overfitting and accuracy overestimation. Our results show that the previously proposed biomarkers are not yet sufficient for accurate and robust inference of alcohol consumption from blood. Overall, our findings imply that DNA methylation prediction biomarkers and models need to be improved considerably before epigenetic inference of alcohol consumption from blood can be considered for practical applications.

High-throughput and affordable genome-wide methylation profiling of circulating cell-free DNA by methylated DNA sequencing (MeD-seq) of LpnPI digested fragments.

BACKGROUND: DNA methylation detection in liquid biopsies provides a highly promising and much needed means for real-time monitoring of disease load in advanced cancer patient care. Compared to the often-used somatic mutations, tissue- and cancer-type specific epigenetic marks affect a larger part of the cancer genome and generally have a high penetrance throughout the tumour. Here, we describe the successful application of the recently described MeD-seq assay for genome-wide DNA methylation profiling on cell-free DNA (cfDNA). The compatibility of the MeD-seq assay with different types of blood collection tubes, cfDNA input amounts, cfDNA isolation methods, and vacuum concentration of samples was evaluated using plasma from both metastatic cancer patients and healthy blood donors (HBDs). To investigate the potential value of cfDNA methylation profiling for tumour load monitoring, we profiled paired samples from 8 patients with resectable colorectal liver metastases (CRLM) before and after surgery. RESULTS: The MeD-seq assay worked on plasma-derived cfDNA from both EDTA and CellSave blood collection tubes when at least 10 ng of cfDNA was used. From the 3 evaluated cfDNA isolation methods, both the manual QIAamp Circulating Nucleic Acid Kit (Qiagen) and the semi-automated Maxwell® RSC ccfDNA Plasma Kit (Promega) were compatible with MeD-seq analysis, whereas the QiaSymphony DSP Circulating DNA Kit (Qiagen) yielded significantly fewer reads when compared to the QIAamp kit (p < 0.001). Vacuum concentration of samples before MeD-seq analysis was possible with samples in AVE buffer (QIAamp) or water, but yielded inconsistent results for samples in EDTA-containing Maxwell buffer. Principal component analysis showed that pre-surgical samples from CRLM patients were very distinct from HBDs, whereas post-surgical samples were more similar. Several described methylation markers for colorectal cancer monitoring in liquid biopsies showed differential methylation between pre-surgical CRLM samples and HBDs in our data, supporting the validity of our approach. Results for MSC, ITGA4, GRIA4, and EYA4 were validated by quantitative methylation specific PCR. CONCLUSIONS: The MeD-seq assay provides a promising new method for cfDNA methylation profiling. Potential future applications of the assay include marker discovery specifically for liquid biopsy analysis as well as direct use as a disease load monitoring tool in advanced cancer patients.

Impact of BMI and waist circumference on epigenome-wide DNA methylation and identification of epigenetic biomarkers in blood: an EWAS in multi-ethnic Asian individuals.

BACKGROUND: The prevalence of obesity and its related chronic diseases have been increasing especially in Asian countries. Obesity-related genetic variants have been identified, but these explain little of the variation in BMI. Recent studies reported associations between DNA methylation and obesity, mostly in non-Asian populations. METHODS: We performed an epigenome-wide association study (EWAS) on general adiposity (body mass index, BMI) and abdominal adiposity (waist circumference, WC) in 409 multi-ethnic Asian individuals and replicated BMI and waist-associated DNA methylation CpGs identified in other populations. The cross-lagged panel model and Mendelian randomization were used to assess the temporal relationship between methylation and BMI. The temporal relationship between the identified CpGs and inflammation and metabolic markers was also examined. RESULTS: EWAS identified 116 DNA methylation CpGs independently associated with BMI and eight independently associated with WC at false discovery rate PFDR < 0.05 in 409 Asian samples. We replicated 110 BMI-associated CpGs previously reported in Europeans and identified six novel BMI-associated CpGs and two novel WC-associated CpGs. We observed high consistency in association direction of effect compared to studies in other populations. Causal relationship analyses indicated that BMI was more likely to be the cause of DNA methylation alteration, rather than the consequence. The causal analyses using BMI-associated methylation risk score also suggested that higher levels of the inflammation marker IL-6 were likely the consequence of methylation change. CONCLUSION: Our study provides evidence of an association between obesity and DNA methylation in multi-ethnic Asians and suggests that obesity can drive methylation change. The results also suggested possible causal influence that obesity-related methylation changes might have on inflammation and lipoprotein levels.

Whole-genome sequencing as a tool for studying the microevolution of drug-resistant serial Mycobacterium tuberculosis isolates.

Treatment of drug-resistant tuberculosis requires extended use of more toxic and less effective drugs and may result in retreatment cases due to failure, abandonment or disease recurrence. It is therefore important to understand the evolutionary process of drug resistance in Mycobacterium tuberculosis. We here in describe the microevolution of drug resistance in serial isolates from six previously treated patients. Drug resistance was initially investigated through phenotypic methods, followed by genotypic approaches. The use of whole-genome sequencing allowed the identification of mutations in the katG, rpsL and rpoB genes associated with drug resistance, including the detection of rare mutations in katG and mixed populations of strains. Molecular docking simulation studies of the impact of observed mutations on isoniazid binding were also performed. Whole-genome sequencing detected 266 single nucleotide polymorphisms between two isolates obtained from one patient, suggesting a case of exogenous reinfection. In conclusion, sequencing technologies can detect rare mutations related to drug resistance, identify subpopulations of resistant strains, and identify diverse populations of strains due to exogenous reinfection, thus improving tuberculosis control by guiding early implementation of appropriate clinical and therapeutic interventions.