DNA Methylation Age Mediates Effect of Metabolic Profile on Cardiovascular and General Aging

BACKGROUND

Alterations in lipid metabolism and DNA methylation are 2 hallmarks of aging. Connecting metabolomic, epigenomic, and aging outcomes help unravel the complex mechanisms underlying aging. We aimed to assess whether DNA methylation clocks mediate the association of circulating metabolites with incident atherosclerotic cardiovascular disease (ASCVD) and frailty.

METHODS

The China Kadoorie Biobank is a prospective cohort study with a baseline survey from 2004 to 2008 and a follow-up period until December 31, 2018. We used the Infinium Methylation EPIC BeadChip to measure the methylation levels of 988 participants' baseline blood leukocyte DNA. Metabolite profiles, including lipoprotein particles, lipid constituents, and various circulating metabolites, were measured using quantitative nuclear magnetic resonance. The pace of DNA methylation age acceleration (AA) was calculated using 5 widely used epigenetic clocks (the first generation: Horvath, Hannum, and Li; the second generation: Grim and Pheno). Incident ASCVD was ascertained through linkage with local death and disease registries and national health insurance databases, supplemented by active follow-up. The frailty index was constructed using medical conditions, symptoms, signs, and physical measurements collected at baseline.

RESULTS

A total of 508 incident cases of ASCVD were documented during a median follow-up of 9.5 years. The first generation of epigenetic clocks was associated with the risk of ASCVD (P<0.05). For each SD increment in LiAA, HorvathAA, and HannumAA, the corresponding hazard ratios for ASCVD risk were 1.16 (1.05-1.28), 1.10 (1.00-1.22), and 1.17 (1.04-1.31), respectively. Only LiAA mediated the association of various metabolites (lipids, fatty acids, histidine, and inflammatory biomarkers) with ASCVD, with the mediating proportion reaching up to 15% for the diameter of low-density lipoprotein (P=1.2×10-2). Regarding general aging, a 1-SD increase in GrimAA was associated with an average increase of 0.10 in the frailty index (P=2.0×10-3), and a 33% and 63% increased risk of prefrailty and frailty at baseline (P=1.5×10-2 and 5.8×10-2), respectively; this association was not observed with other clocks. GrimAA mediated the effect of various lipids, fatty acids, glucose, lactate, and inflammatory biomarkers on the frailty index, with the mediating proportion reaching up to 22% for triglycerides in very small-sized very low-density lipoprotein (P=6.0×10-3).

CONCLUSIONS

These findings suggest that epigenomic mechanisms may play a role in the associations between circulating metabolites and the aging process. Different mechanisms underlie the first and second generations of DNA methylation age in cardiovascular and general aging.

Genome-wide DNA methylation profiling in blood reveals epigenetic signature of incident acute coronary syndrome

DNA methylation (DNAm) has been implicated in acute coronary syndrome (ACS), but the causality remains unclear in cross-sectional studies. Here, we conduct a prospective epigenome-wide association study of incident ACS in two Chinese cohorts (discovery: 751 nested case-control pairs; replication: 476 nested case-control pairs). We identified and validated 26 differentially methylated positions (DMPs, false discovery rate [FDR] <0.05), including three mapped to known cardiovascular disease genes (PRKCZ, PRDM16, EHBP1L1) and four with causal evidence from Mendelian randomization (PRKCZ, TRIM27, EMC2, EHBP1L1). Two hypomethylated DMPs were negatively correlated with the expression in blood of their mapped genes (PIGG and EHBP1L1), which were further found to overexpress in leukocytes and/or atheroma plaques. Finally, our DMPs could substantially improve the prediction of ACS over traditional risk factors and polygenic scores. These findings demonstrate the importance of DNAm in the pathogenesis of ACS and highlight DNAm as potential predictive biomarkers and treatment targets.

Association of Computed Tomography Radiomics Signature with Progression-free Survival in Neuroblastoma Patients

AIMS

To investigate the association of computed tomography radiomics signature with progression-free survival (PFS) in neuroblastoma patients.

MATERIALS AND METHODS

We retrospectively included 167 neuroblastoma patients who were divided into a training set and a test set through stratified sampling at a ratio of 7:3. Regions of interest of the primary tumours were delineated on pretreatment contrast-enhanced computed tomography images and radiomics features were extracted from them. The intraclass correlation coefficient, Pearson correlation coefficient, and least absolute shrinkage and selection operator Cox regression algorithm were applied to select radiomics features and construct the radiomics signature. The effectiveness of the signature in predicting PFS was evaluated using the concordance index (C-index) and 95% confidence interval in both the training and the test sets. The time-dependent receiver operator characteristic curve of the radiomics signature was plotted and the area under the curve (AUC) was calculated. A calibration curve was used to assess the difference between the predicted probability of the radiomics signature and the observed probability at different time points.

RESULTS

The radiomics signature was composed of six features, which achieved a C-index of 0.733 (95% confidence interval 0.664-0.803) in the training set and 0.734 (95% confidence interval 0.608-0.861) in the test set. In the training set, the radiomics signature yielded an AUC of 0.707, 0.737, 0.788, 0.859 and 0.829 for 1-, 2-, 3-, 4- and 5-year PFS, respectively. Similarly, the radiomics signature exhibited an AUC of 0.738, 0.807, 0.761, 0.787 and 0.818 for 1-, 2-, 3-, 4- and 5-year PFS, respectively, in the test set. The calibration curves showed no significant difference between the predicted probability of the radiomics signature and the observed probability for up to 5 years.

CONCLUSIONS

Computed tomography radiomics features exhibit a significant correlation with the PFS of neuroblastoma patients, particularly in terms of long-term outcomes.

Minimal improvement in coronary artery disease risk prediction in Chinese population using polygenic risk scores: evidence from the China Kadoorie Biobank

BACKGROUND

Several studies have reported that polygenic risk scores (PRSs) can enhance risk prediction of coronary artery disease (CAD) in European populations. However, research on this topic is far from sufficient in non-European countries, including China. We aimed to evaluate the potential of PRS for predicting CAD for primary prevention in the Chinese population.

METHODS

Participants with genome-wide genotypic data from the China Kadoorie Biobank were divided into training ( n = 28,490) and testing sets ( n = 72,150). Ten previously developed PRSs were evaluated, and new ones were developed using clumping and thresholding or LDpred method. The PRS showing the strongest association with CAD in the training set was selected to further evaluate its effects on improving the traditional CAD risk-prediction model in the testing set. Genetic risk was computed by summing the product of the weights and allele dosages across genome-wide single-nucleotide polymorphisms. Prediction of the 10-year first CAD events was assessed using hazard ratios (HRs) and measures of model discrimination, calibration, and net reclassification improvement (NRI). Hard CAD (nonfatal I21-I23 and fatal I20-I25) and soft CAD (all fatal or nonfatal I20-I25) were analyzed separately.

RESULTS

In the testing set, 1214 hard and 7201 soft CAD cases were documented during a mean follow-up of 11.2 years. The HR per standard deviation of the optimal PRS was 1.26 (95% CI:1.19-1.33) for hard CAD. Based on a traditional CAD risk prediction model containing only non-laboratory-based information, the addition of PRS for hard CAD increased Harrell's C index by 0.001 (-0.001 to 0.003) in women and 0.003 (0.001 to 0.005) in men. Among the different high-risk thresholds ranging from 1% to 10%, the highest categorical NRI was 3.2% (95% CI: 0.4-6.0%) at a high-risk threshold of 10.0% in women. The association of the PRS with soft CAD was much weaker than with hard CAD, leading to minimal or no improvement in the soft CAD model.

CONCLUSIONS

In this Chinese population sample, the current PRSs minimally changed risk discrimination and offered little improvement in risk stratification for soft CAD. Therefore, this may not be suitable for promoting genetic screening in the general Chinese population to improve CAD risk prediction.

Maternal pre-pregnancy BMI, offspring epigenome-wide DNA methylation, and childhood obesity: findings from the Boston Birth Cohort

BACKGROUND

Maternal pre-pregnancy obesity is an established risk factor for childhood obesity. Investigating epigenetic alterations induced by maternal obesity during fetal development could gain mechanistic insight into the developmental origins of childhood obesity. While obesity disproportionately affects underrepresented racial and ethnic mothers and children in the USA, few studies investigated the role of prenatal epigenetic programming in intergenerational obesity of these high-risk populations.

METHODS

This study included 903 mother-child pairs from the Boston Birth Cohort, a predominantly urban, low-income minority birth cohort. Mother-infant dyads were enrolled at birth and the children were followed prospectively to age 18 years. Infinium Methylation EPIC BeadChip was used to measure epigenome-wide methylation level of cord blood. We performed an epigenome-wide association study of maternal pre-pregnancy body mass index (BMI) and cord blood DNA methylation (DNAm). To quantify the degree to which cord blood DNAm mediates the maternal BMI-childhood obesity, we further investigated whether maternal BMI-associated DNAm sites impact birthweight or childhood overweight or obesity (OWO) from age 1 to age 18 and performed corresponding mediation analyses.

RESULTS

The study sample contained 52.8% maternal pre-pregnancy OWO and 63.2% offspring OWO at age 1-18 years. Maternal BMI was associated with cord blood DNAm at 8 CpG sites (genome-wide false discovery rate [FDR] < 0.05). After accounting for the possible interplay of maternal BMI and smoking, 481 CpG sites were discovered for association with maternal BMI. Among them 123 CpGs were associated with childhood OWO, ranging from 42% decrease to 87% increase in OWO risk for each SD increase in DNAm. A total of 14 identified CpG sites showed a significant mediation effect on the maternal BMI-child OWO association (FDR < 0.05), with mediating proportion ranging from 3.99% to 25.21%. Several of these 14 CpGs were mapped to genes in association with energy balance and metabolism (AKAP7) and adulthood metabolic syndrome (CAMK2B).

CONCLUSIONS

This prospective birth cohort study in a high-risk yet understudied US population found that maternal pre-pregnancy OWO significantly altered DNAm in newborn cord blood and provided suggestive evidence of epigenetic involvement in the intergenerational risk of obesity.

Healthy lifestyle, DNA methylation age acceleration, and incident risk of coronary heart disease

BACKGROUND

DNA methylation clocks emerged as a tool to determine biological aging and have been related to mortality and age-related diseases. Little is known about the association of DNA methylation age (DNAm age) with coronary heart disease (CHD), especially in the Asian population.

RESULTS

Methylation level of baseline blood leukocyte DNA was measured by Infinium Methylation EPIC BeadChip for 491 incident CHD cases and 489 controls in the prospective China Kadoorie Biobank. We calculated the methylation age using a prediction model developed among Chinese. The correlation between chronological age and DNAm age was 0.90. DNA methylation age acceleration (Δage) was defined as the residual of regressing DNA methylation age on the chronological age. After adjustment for multiple risk factors of CHD and cell type proportion, compared with participants in the bottom quartile of Δage, the OR (95% CI) for CHD was 1.84 (1.17, 2.89) for participants in the top quartile. One SD increment in Δage was associated with 30% increased risk of CHD (OR = 1.30; 95% CI 1.09, 1.56; Ptrend = 0.003). The average number of cigarette equivalents consumed per day and waist-to-hip ratio were positively associated with Δage; red meat consumption was negatively associated with Δage, characterized by accelerated aging in those who never or rarely consumed red meat (all P < 0.05). Further mediation analysis revealed that 10%, 5% and 18% of the CHD risk related to smoking, waist-to-hip ratio and never or rarely red meat consumption was mediated through methylation aging, respectively (all P for mediation effect < 0.05).

CONCLUSIONS

We first identified the association between DNAm age acceleration and incident CHD in the Asian population, and provided evidence that unfavorable lifestyle-induced epigenetic aging may play an important part in the underlying pathway to CHD.

Temporal associations between leukocytes DNA methylation and blood lipids: a longitudinal study

BACKGROUND

The associations between blood lipids and DNA methylation have been investigated in epigenome-wide association studies mainly among European ancestry populations. Several studies have explored the direction of the association using cross-sectional data, while evidence of longitudinal data is still lacking.

RESULTS

We tested the associations between peripheral blood leukocytes DNA methylation and four lipid measures from Illumina 450 K or EPIC arrays in 1084 participants from the Chinese National Twin Registry and replicated the result in 988 participants from the China Kadoorie Biobank. A total of 23 associations of 19 CpG sites were identified, with 4 CpG sites located in or adjacent to 3 genes (TMEM49, SNX5/SNORD17 and CCDC7) being novel. Among the validated associations, we conducted a cross-lagged analysis to explore the temporal sequence and found temporal associations of methylation levels of 2 CpG sites with triglyceride and 2 CpG sites with high-density lipoprotein-cholesterol (HDL-C) in all twins. In addition, methylation levels of cg11024682 located in SREBF1 at baseline were temporally associated with triglyceride at follow-up in only monozygotic twins. We then performed a mediation analysis with the longitudinal data and the result showed that the association between body mass index and HDL-C was partially mediated by the methylation level of cg06500161 (ABCG1), with a mediation proportion of 10.1%.

CONCLUSIONS

Our study indicated that the DNA methylation levels of ABCG1, AKAP1 and SREBF1 may be involved in lipid metabolism and provided evidence for elucidating the regulatory mechanism of lipid homeostasis.

Hsa_circ_0007312 Promotes Third-Generation Epidermal Growth Factor Receptor-Tyrosine Kinase Inhibitor Resistance through Pyroptosis and Apoptosis via the MiR-764/MAPK1 Axis in Lung Adenocarcinoma Cells

Purposes: Osimertinib is a third-generation epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI) used for patients with gefitinib (first-generation EGFR-TKI) resistance, but osimertinib resistance inevitably occurs. Therefore, it is necessary to explore the mechanisms of osimertinib resistance. Materials and Methods: We performed quantitative real-time polymerase chain reaction to detect hsa_circ_0007312 (circ7312), miR-764, and MAPK1 expressions in tissues and cells. Western blotting was used to detect protein levels in cells. Cell Counting Kit-8, apoptotic, and Transwell assays were used to explore biological functions. Luciferase assays were used to identify the interactions between circ7312 and miR-764, MAPK1 and miR-764. A xenograft experiment was performed to clarify the role of circ7312 in vivo. Public datasets were used to identify the relation between circ7312 expression and the cell half maximal inhibitory concentration value of osimertinib in 41 lung adenocarcinoma cell lines. The Student t-test, Kaplan-Meier analysis, and Pearson correlation analysis were used in data analysis. Results: We found that circ7312 knockdown increased miR-764 expression and decreased MAPK1 expression, and circ7312 regulated MAPK1 by sponging miR-764. In addition, high circ7312 expression has significant positive correlation with osimertinib IC50 values, circ7312 knockdown decreased the cell half maximal inhibitory concentration value of osimertinib and increased pyroptosis and apoptosis by sponging the miR-764/MAPK1 axis. We also found that circ7312 and MAPK1 were highly expressed in tumor tissues and related to poor prognosis. Xenograft experiments revealed that circ7312 knockdown decreased osimertinib resistance in vivo. Conclusion: We demonstrated that the inhibition of circ7312 decreased osimertinib resistance by promoting pyroptosis and apoptosis via the miR-764/MAPK1 axis, providing a novel target for osimertinib resistance therapy.

circ_0000567/miR-421/TMEM100 Axis Promotes the Migration and Invasion of Lung Adenocarcinoma and Is Associated with Prognosis

Purpose: Due to the high metastatic ability and poor prognosis of lung adenocarcinoma (LUAD), we identified novel non-coding RNAs, which constitute approximately 60% of human transcripts, as prognostic biomarkers and potential therapeutic targets for LUAD.

Methods: In this study, we downloaded and analyzed microRNA (miRNA) datasets from The Cancer Genome Atlas (TCGA) to identify dysregulated miRNAs correlating with the overall survival (OS) of LUAD patients. miR-421, circ_0000567, and TMEM100 expression levels were examined by quantitative real-time polymerase chain reaction (qRT-PCR) in NSCLC tissues from 73 patients and adjacent normal tissues. Cell migration and invasion were assayed using wound healing and transwell assays. miR-421 target predictions were conducted using starBase, CircInteractome, circBank, TargetScan, miRanda, MirDB, miRpath, and Gene Expression Omnibus (GEO) databases. The circular structure and stability of circ_0000567 were verified by RNase R digestion and qRT-PCR using oligo(dT) and random primers. A luciferase reporter assay was used to evaluate the relationship between miR-421, circ_0000567, and TMEM100.

Results: The miRNA panel associated with OS in patients with LUAD was screened according to the hazard ratio (HR) of miRNAs from high to low. Based on the correlation between these miRNAs and OS, as well as miRNA expression levels, miR-421 was selected for further outcome analysis. High miR-421 expression was an independent risk factor for shorter OS in 73 patients collected from our department. Bioinformatic analyses, luciferase reporter assays, and functional assays showed that circ_0000567 could act as a sponge for miR-421 and prevent it from directly targeting the 3'-untranslated region of TMEM100 mRNA and further degrading it in LUAD. miR-421 expression was significantly upregulated, while circ_0000567 and TMEM100 were downregulated in tumor tissues of LUAD, compared to their counterparts in normal tissues. Gain- and loss-of-function assays showed that miR-421 promoted LUAD cell migration and invasion. Overexpression of circ_0000567 inhibited migration and invasion, whereas co-transfection of circ_0000567 and miR-421 mimics partly counteracted this effect. TMEM100 was upregulated by enhanced circ_0000567 in LUAD cells, and the expression of TMEM100 was inversely proportional to miR-421, whereas it was directly proportional to circ_0000567 in 73 LUAD specimens, which confirmed the competitive endogenous RNA (ceRNA) network.

Conclusion: Our findings suggest that miR-421 promotes the migration and invasion of lung adenocarcinoma via circ_0000567/miR-421/TMEM100 signaling and could be a prognostic biomarker for LUAD.

Integrated analysis of whole genome and transcriptome sequencing reveals a frameshift mutation associated with recessive embryonic lethality in Holstein cattle

Embryo loss is an important factor affecting fertility in dairy production. HH2 was identified as a haplotype on chromosome 1 associated with embryonic lethality in Holstein cattle. In the current study, both short- and long-read WGS was performed on four carriers and four non-carriers of HH2 to screen for variants in concordance with HH2 haplotype status. Sequence variation analysis revealed five putative functional variants of protein-coding genes, including a frameshift mutation (g.107172616delT) in intraflagellar transport protein 80 (IFT80) gene. Transcriptome analysis of whole blood indicated that no gene exhibited significantly differential expression or allele-specific expression between carriers and non-carriers in the candidate region. This evidence points to g.107172616delT as the highest priority causative mutation for HH2. Protein prediction reveals that the frameshift mutation results in a premature stop codon to reduce the peptide chain from 760 to 383 amino acids and greatly alters the structure and function of IFT80 protein. Our results demonstrate that the use of a combination of multiple high-throughput sequencing technologies is an efficient strategy to screen for the candidate causative mutations responsible for Mendelian traits, including genetic disorders.

HIF1A-AS2 induces osimertinib resistance in lung adenocarcinoma patients by regulating the miR-146b-5p/IL-6/STAT3 axis

Although epidermal growth factor receptor tyrosine kinase inhibitors (TKIs) show efficacy in lung adenocarcinoma (LUAD) patients, TKI resistance inevitably develops, limiting long-term results. Thus, there is an urgent need to address drug resistance in LUAD. Long non-coding RNA (lncRNA) HIF1A-AS2 could be a critical mediator in the progression of various tumor types. We examined the function of HIF1A-AS2 in modifying tumor aggravation and osimertinib resistance in lung adenocarcinoma. Using clinical samples, we showed that HIF1A-AS2 was upregulated in LUAD specimens, predicting poorer overall survival and disease-free survival. HIF1A-AS2 silencing inhibited the proliferation, migration, and tumorigenesis of LUAD cells and therapeutic efficacy of osimertinib against tumor cells in vitro and in vivo. RNA precipitation assays, western blotting, luciferase assays, and rescue experiments demonstrated that HIF1A-AS2 sponged microRNA-146b-5p (miR-146b-5p), promoting interleukin-6 (IL-6) expression, activating the IL-6/STAT3 pathway, and leading to LUAD progression. miR-146b-5p and IL-6 levels were correlated with the prognosis of LUAD patients. Our results indicated that HIF1A-AS2 functions as an oncogenic factor in adenocarcinoma cells by targeting the miR-146b-5p/IL-6/STAT3 axis and may be a prognostic indicator of survival. Moreover, it can be a potential therapeutic target to enhance the efficacy of osimertinib in LUAD patients.

Construction of a circRNA-miRNA-mRNA Regulated Pathway Involved in EGFR-TKI Lung Adenocarcinoma Resistance

Objectives: Epidermal growth factor receptor-tyrosine kinase inhibitors are widely used for lung epidermal growth factor receptor-positive lung adenocarcinomas, but acquired resistance is inevitable. Although non-coding RNAs, such as circular RNA and microRNA, are known to play vital roles in epidermal growth factor receptor-tyrosine kinase inhibitor resistance, comprehensive analysis is lacking. Thus, this study aimed to explore the circular RNA-microRNA-messenger RNA regulatory network involved in epidermal growth factor receptor-tyrosine kinase inhibitor resistance. Methods: To identify differentially expressed genes between the epidermal growth factor receptor-tyrosine kinase inhibitor sensitive cell line PC9 and resistant cell line PC9/ epidermal growth factor receptor-tyrosine kinase inhibitor resistance(PC9/ER), circular RNA, microRNA and messenger RNA microarrays were performed. Candidates were then identified to construct a circular RNA-microRNA-messenger RNA network using bioinformatics. Additionally, Gene Oncology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted to evaluate the network messenger RNA, setting up a protein-protein interaction network for hub-gene identification. Afterwards, RNA immunoprecipitation was performed to enrich microRNA, and quantitative real-time PCR was used to estimated gene expression levels. Results: In total, 603, 377, and 1863 differentially expressed circular RNA, microRNA, messenger RNAs, respectively, were identified using microarray analysis, constructing a circular RNA-microRNA-messenger RNA network containing 18 circular RNAs, 17 microRNAs and 175 messenger RNAs. Moreover, Gene Oncology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the most enriched biological process terms and pathways were related to epidermal growth factor receptor-tyrosine kinase inhibitor resistance, including Wnt and Hippo signaling pathways. Based on the competing endogenous RNA and protein-protein interaction network, circ-0007312 was showed to interact with miR-764 and both circ-0003748 and circ-0001398 were shown to interact with miR-628; both these microRNAs targeted MAPK1. Furthermore, circ-0007312, circ-0003748, circ-0001398, and MAPK1 were up-regulated, whereas miR-764 and miR-628 were downregulated in PC9/ER cells as compared to parental PC9 cells. We also found that circ-0007312 and miR-764 were positively expressed in plasma. Conclusions: Our original study associated with mechanism of target therapy in lung cancer provided a systematic and comprehensive regulation of circular RNA, microRNA and messenger RNA in epidermal growth factor receptor-tyrosine kinase inhibitor resistance. It was found that circ-0007312- miR-764-MAPK1, circ-0003748-miR-628-MAPK1, and circ-0001398-miR-628-MAPK1 axis may play key roles in epidermal growth factor receptor-tyrosine kinase inhibitor resistance.

A transdisciplinary approach to understand the epigenetic basis of race/ethnicity health disparities

Health disparities correspond to differences in disease burden and mortality among socially defined population groups. Such disparities may emerge according to race/ethnicity, socioeconomic status and a variety of other social contexts, and are documented for a wide range of diseases. Here, we provide a transdisciplinary perspective on the contribution of epigenetics to the understanding of health disparities, with a special emphasis on disparities across socially defined racial/ethnic groups. Scientists in the fields of biological anthropology, bioinformatics and molecular epidemiology provide a summary of theoretical, statistical and practical considerations for conducting epigenetic health disparities research, and provide examples of successful applications from cancer research using this approach.

A large-scale genome-wide association analysis of lung function in the Chinese population identifies novel loci and highlights shared genetic aetiology with obesity

BACKGROUND

Lung function is a heritable complex phenotype with obesity being one of its important risk factors. However, knowledge of their shared genetic basis is limited. Most genome-wide association studies (GWASs) for lung function have been based on European populations, limiting the generalisability across populations. Large-scale lung function GWASs in other populations are lacking.

METHODS

We included 100 285 subjects from the China Kadoorie Biobank (CKB). To identify novel loci for lung function, single-trait GWAS analyses were performed on forced expiratory volume in 1 s (FEV1), forced vital capacity (FVC) and FEV1/FVC in the CKB. We then performed genome-wide cross-trait analysis between lung function and obesity traits (body mass index (BMI), BMI-adjusted waist-to-hip ratio and BMI-adjusted waist circumference) to investigate the shared genetic effects in the CKB. Finally, polygenic risk scores (PRSs) of lung function were developed in the CKB and their interaction with BMI's association on lung function were examined. We also conducted cross-trait analysis in parallel with the CKB using up to 457 756 subjects from the UK Biobank (UKB) for replication and investigation of ancestry-specific effects.

RESULTS

We identified nine genome-wide significant novel loci for FEV1, six for FVC and three for FEV1/FVC in the CKB. FEV1 and FVC showed significant negative genetic correlation with obesity traits in both the CKB and UKB. Genetic loci shared between lung function and obesity traits highlighted important biological pathways, including cell proliferation, embryo, skeletal and tissue development, and regulation of gene expression. Mendelian randomisation analysis suggested significant negative causal effects of BMI on FEV1 and on FVC in both the CKB and UKB. Lung function PRSs significantly modified the effect of change in BMI on change in lung function during an average follow-up of 8 years.

CONCLUSION

This large-scale GWAS of lung function identified novel loci and shared genetic aetiology between lung function and obesity. Change in BMI might affect change in lung function differently according to a subject's polygenic background. These findings may open new avenues for the development of molecular-targeted therapies for obesity and lung function improvement.

Epigenome-wide analysis of DNA methylation and coronary heart disease: a nested case-control study

Background

Identifying environmentally responsive genetic loci where DNA methylation is associated with coronary heart disease (CHD) may reveal novel pathways or therapeutic targets for CHD. We conducted the first prospective epigenome-wide analysis of DNA methylation in relation to incident CHD in the Asian population.

Methods

We did a nested case-control study comprising incident CHD cases and 1:1 matched controls who were identified from the 10 year follow-up of the China Kadoorie Biobank. Methylation level of baseline blood leukocyte DNA was measured by Infinium Methylation EPIC BeadChip. We performed the single cytosine-phosphate-guanine (CpG) site association analysis and network approach to identify CHD-associated CpG sites and co-methylation gene module.

Results

After quality control, 982 participants (mean age 50.1 years) were retained. Methylation level at 25 CpG sites across the genome was associated with incident CHD (genome-wide false discovery rate [FDR] < 0.05 or module-specific FDR < 0.01). One SD increase in methylation level of identified CpGs was associated with differences in CHD risk, ranging from a 47 % decrease to a 118 % increase. Mediation analyses revealed 28.5 % of the excessed CHD risk associated with smoking was mediated by methylation level at the promoter region of ANKS1A gene (P for mediation effect = 0.036). Methylation level at the promoter region of SNX30 was associated with blood pressure and subsequent risk of CHD, with the mediating proportion to be 7.7 % (P = 0.003) via systolic blood pressure and 6.4 % (P = 0.006) via diastolic blood pressure. Network analysis revealed a co-methylation module associated with CHD.

Conclusions

We identified novel blood methylation alterations associated with incident CHD in the Asian population and provided evidence of the possible role of epigenetic regulations in the smoking- and blood pressure-related pathways to CHD risk.

Funding

This work was supported by National Natural Science Foundation of China (81390544 and 91846303). The CKB baseline survey and the first re-survey were supported by a grant from the Kadoorie Charitable Foundation in Hong Kong. The long-term follow-up is supported by grants from the UK Wellcome Trust (202922/Z/16/Z, 088158/Z/09/Z, 104085/Z/14/Z), grant (2016YFC0900500, 2016YFC0900501, 2016YFC0900504, 2016YFC1303904) from the National Key R&D Program of China, and Chinese Ministry of Science and Technology (2011BAI09B01).

Small Cell Size Circulating Aneuploid Cells as a Biomarker of Prognosis in Resectable Non-Small Cell Lung Cancer

Objective

The size distribution of circulating aneuploid cells (CACs) and its clinical significance were investigated in resectable non-small cell lung cancer (NSCLC).

Patients and Methods

A total of 50 patients with resectable NSCLC were enrolled in this study. Blood samples (50 pre-surgery and 35 post-surgery) were collected and used for the detection of CAC chromosome 8 heteroploidy through the subtraction enrichment and immunostaining fluorescence in situ hybridization (SE-iFISH) method.

Results

Less than 20% small cell size and more than 80% large cell size CACs were detected. Karyotypes, including triploid, tetraploid, and multiploid, had varying distributions. The triploid subtype accounted for the majority of small cell size CACs, whereas the multiploid subtype accounted for the majority of large cell size CACs. We found that total small cell size and triploid small cell size CACs, but not large cell size CACs, derived from pre-surgery samples, were associated with shorter disease-free survival. Moreover, total small cell size and triploid small cell size CACs were associated with higher TNM stage and recurrence. Nevertheless, the variation between pre- and post-surgery CACs was not related to survival among patients with resectable NSCLC.

Conclusions

Pre-surgery small cell size CACs, especially the triploid subtype, could be regarded as a potential prognostic biomarker for patients with resectable NSCLC.

Improved lipidomic profile mediates the effects of adherence to healthy lifestyles on coronary heart disease

Adherence to healthy lifestyles is associated with reduced risk of coronary heart disease (CHD), but uncertainty persists about the underlying lipid pathway. In a case-control study of 4681 participants nested in the prospective China Kadoorie Biobank, 61 lipidomic markers in baseline plasma were measured by targeted nuclear magnetic resonance spectroscopy. Baseline lifestyles included smoking, alcohol consumption, dietary habit, physical activity, and adiposity levels. Genetic instrument was used to mimic the lipid-lowering effect of statins. We found that 35 lipid metabolites showed statistically significant mediation effects in the pathway from healthy lifestyles to CHD reduction, including very low-density lipoprotein (VLDL) particles and their cholesterol, large-sized high-density lipoprotein (HDL) particle and its cholesterol, and triglyceride in almost all lipoprotein subfractions. The statins genetic score was associated with reduced intermediate- and low-density lipoprotein, but weak or no association with VLDL and HDL. Lifestyle interventions and statins may improve different components of the lipid profile.

Caspofungin combined with TMP/SMZ as a first-line therapy for moderate-to-severe PCP in patients with human immunodeficiency virus infection

OBJECTIVES

The effectiveness of trimethoprim/sulfamethoxazole (TMP/SMZ) for pneumocystis pneumonia (PCP) is limited with adverse events. Caspofungin, by inhibiting the cyst form of Pneumocystis jirovecii, may be an alternative therapy for PCP. However, the availability of clinical data about caspofungin combined with TMP/SMZ in the treatment of PCP in HIV-infected patients is limited. Thus, we aimed to examine the clinical effectiveness and safety of caspofungin combined with TMP/SMZ as a first-line therapy for moderate-to-severe PCP in HIV-infected patients.

METHODS

From January 2017 to December 2019, data of HIV-infected patients with moderate-to-severe PCP who received either TMP/SMZ alone or caspofungin combined with TMP/SMZ as first-line therapy were retrospectively reviewed to assess the effectiveness and safety of each regimen. The Kaplan-Meier curve and log-rank test were used for survival analysis.

RESULTS

A total of 278 patients met the criteria. The overall positive response rate of PCP treatment was 48.92%, and the overall all-cause in-hospital mortality rate was 33.09%. Patients who received combination therapy consisting of caspofungin and TMP/SMZ had a better positive response rate (59.44% vs. 37.78%, P < 0.001) and lower all-cause in-hospital mortality rate (24.48% vs. 42.22%, P = 0.003). Also, patients who received combination therapy had higher survival rate during a hospital stay (75.52% vs. 57.78%, P = 0.004), and those who received longer combination therapy were more likely to have higher survival rate (P = 0.042). We found that age (P = 0.019), CD4 cell count (P = 0.001) and therapeutic regimen (P = 0.002) were significant risk factors for all-cause in-hospital mortality rate in univariate analysis. In multivariate analysis, only CD4 cell count and therapeutic regimen were statistically significant factors associated with all-cause in-hospital mortality rate. Patients with a CD4 count of > 30 cells/µL and patients who received combination therapy consisting of caspofungin and TMP/SMZ were more likely to survive from PCP (P = 0.011 and P = 0.002, respectively). There were no additional severe adverse events caused by adding caspofungin.

CONCLUSIONS

For HIV-infected patients with moderate-to-severe PCP, combination therapy with caspofungin and TMP/SMZ is an effective and promising first-line therapy with no greater number of adverse events compared with TMP/SMZ monotherapy. Patients who received caspofungin had better positive response rates and lower all-cause in-hospital mortality rates. Also, we recommend early initiation of caspofungin.

Early famine exposure and adult disease risk based on a 10-year prospective study of Chinese adults

OBJECTIVE

To comprehensively examine the potential impacts of prenatal experience of the Chinese Great Famine on chronic disease risks in the middle age.

METHODS

This study included 92 284 participants aged 39-51 years from China Kadoorie Biobank born around the famine period and without major chronic diseases at baseline. We categorised participants into non-famine births (born between 1 October 1956 and 30 September 1958, and 1 October 1962 and 30 September 1964) and famine births (born between 1 October 1959 and 30 September 1961). The outcomes were incident cardiovascular disease, cancer and respiratory system disease. Cox regression was used to estimate adjusted HR and 95% CI for famine exposure. Subgroup analyses were performed according to baseline characteristics.

RESULTS

During a median 10.1 years of follow-up, we identified 4626 incident ischaemic heart disease (IHD) cases, 7332 cerebrovascular disease cases, 3111 cancer cases and 16 081 respiratory system disease cases. In the whole population, prenatal famine exposure was not statistically associated with the risks of developing any chronic diseases in adulthood. However, for urban participants, compared with non-famine births, famine births had a higher risk of cerebrovascular disease (HR 1.18; 95% CI 1.09 to 1.28); such association was not shown for rural participants (p for interaction <0.001). Also, we observed the associations of prenatal famine exposure with IHD (HR 1.15; 95% CI 1.05 to 1.26) and cerebrovascular disease (HR 1.13; 95% CI 1.05 to 1.21) in participants with lower physical activity level, but not in those with higher ones (all p for interaction=0.003).

CONCLUSION

Our findings indicate that prenatal exposure to the Chinese famine might be associated with an increased cardiovascular risk and such risk may be modified by adult lifestyle.

Shisa3 brakes resistance to EGFR-TKIs in lung adenocarcinoma by suppressing cancer stem cell properties

Background

Although EGFR tyrosine kinase inhibitors (EGFR-TKIs) are beneficial to lung adenocarcinoma patients with sensitive EGFR mutations, resistance to these inhibitors induces a cancer stem cell (CSC) phenotype. Here, we clarify the function and molecular mechanism of shisa3 as a suppressor that can reverse EGFR-TKI resistance and inhibit CSC properties.

Methods

The suppresser genes involved in EGFR-TKI resistance were identified and validated by transcriptome sequencing, quantitative real-time PCR (qRT-PCR) and immunohistochemistry. Biological function analyses, cell half maximal inhibitory concentration (IC50), self-renewal, and migration and invasion capacities, were detected by CCK8, sphere formation and Transwell assays. Tumorigenesis and therapeutic effects were investigated in nonobese diabetic/severe combined immunodeficiency (nod-scid) mice. The underlying mechanisms were explored by Western blot and immunoprecipitation analyses.

Results

We found that low expression of shisa3 was related to EGFR-TKI resistance in lung adenocarcinoma patients. Ectopic overexpression of shisa3 inhibited CSC properties and the cell cycle in the lung adenocarcinoma cells resistant to gefitinib/osimertinib. In contrast, suppression of shisa3 promoted CSC phenotypes and the cell cycle in the cells sensitive to EGFR-TKIs. For TKI-resistant PC9/ER tumors in nod-scid mice, overexpressed shisa3 had a significant inhibitory effect. In addition, we verified that shisa3 inhibited EGFR-TKI resistance by interacting with FGFR1/3 to regulate AKT/mTOR signaling. Furthermore, combinational administration of inhibitors of FGFR/AKT/mTOR and cell cycle signaling could overcome EGFR-TKI resistance associated with shisa3-mediated CSC capacities in vivo.

Conclusion

Taken together, shisa3 was identified as a brake to EGFR-TKI resistance and CSC characteristics, probably through the FGFR/AKT/mTOR and cell cycle pathways, indicating that shisa3 and concomitant inhibition of its regulated signaling may be a promising therapeutic strategy for reversing EGFR-TKI resistance.

[A single-center retrospective study on axillary evaluation in 1 557 breast ductal carcinoma in situ patients between 2006 and 2016]

Objectives: To examine the influence factors on axillary evaluation in ductal carcinoma in situ (DCIS) patients, and the prognosis of different choices of axillary evaluation in a single-center retrospective study. Methods: Totally 1 557 DCIS patients admitted in Department of Breast Surgery, Fudan University Shanghai Cancer Center from January 2006 to November 2016 were retrospectively enrolled. All patients were female. The median age was 49 years (range: 21 to 85 years). Surgical methods included modified radical mastectomy, simple mastectomy (with or without axillary evaluation) and breast conservation surgery (with or without axillary evaluation). Axillary evaluation included axillary lymph node dissection (ALND) and sentinel lymph node biopsy (SLNB). T tests, χ(2) test and Logistic regression analysis was used to analyze influence factors on axillary evaluation, respectively. Kaplan-Meier curve and Log-rank analysis were used to evaluate recurrence-free survival (RFS) and loco-regional recurrence-free survival (LRRFS) in patients with different surgical methods. Results: Among the 1 557 DCIS patients, there were 1 226 cases received axillary evaluation, while 331 cases not received axillary evaluation. Patients were separated into 3 groups by different axillary evaluation choices: SLNB group (957 cases, 61.46%), ALND group (197 cases, 12.65%) and no evaluation group (403 cases, 25.88%). The patients in SLNB group increased significantly (P=0.000), from 3.85% (60/1 557) in 2006 to 75.19% (1 170/1 557) in 2016. The independent influence factors of receiving axillary evaluation were high nuclear grade (OR=3.191, 95%CI: 1.722 to 5.912, P=0.001) and tumor size>15 mm (OR=1.698, 95%CI: 1.120 to 2.573, P=0.012). Also, patients received breast conservation surgery were more likely to refuse axillary evaluation (OR=0.155, 95%CI: 0.103 to 0.233, P=0.000). There were no significant differences in RFS and LRRFS in patients with different axillary evaluation choices. Conclusions: The investigation in trends and influence factors of different axillary evaluation choices provided basis on surgical precision medicine in DCIS patients. Patients received SLNB increased significantly. The independent influence factors of axillary evaluation were nuclear grade, tumor size and surgical methods. There was no significant differences in prognosis among the groups receiving different axillary evaluations.