Loss of epigenetic information as a cause of mammalian aging.

All living things experience an increase in entropy, manifested as a loss of genetic and epigenetic information. In yeast, epigenetic information is lost over time due to the relocalization of chromatin-modifying proteins to DNA breaks, causing cells to lose their identity, a hallmark of yeast aging. Using a system called "ICE" (inducible changes to the epigenome), we find that the act of faithful DNA repair advances aging at physiological, cognitive, and molecular levels, including erosion of the epigenetic landscape, cellular exdifferentiation, senescence, and advancement of the DNA methylation clock, which can be reversed by OSK-mediated rejuvenation. These data are consistent with the information theory of aging, which states that a loss of epigenetic information is a reversible cause of aging.

Ebola Virus Transmission Initiated by Relapse of Systemic Ebola Virus Disease.

During the 2018-2020 Ebola virus disease (EVD) outbreak in North Kivu province in the Democratic Republic of Congo, EVD was diagnosed in a patient who had received the recombinant vesicular stomatitis virus-based vaccine expressing a ZEBOV glycoprotein (rVSV-ZEBOV) (Merck). His treatment included an Ebola virus (EBOV)-specific monoclonal antibody (mAb114), and he recovered within 14 days. However, 6 months later, he presented again with severe EVD-like illness and EBOV viremia, and he died. We initiated epidemiologic and genomic investigations that showed that the patient had had a relapse of acute EVD that led to a transmission chain resulting in 91 cases across six health zones over 4 months. (Funded by the Bill and Melinda Gates Foundation and others.).

Blood Leukocyte DNA Methylation Predicts Risk of Future Myocardial Infarction and Coronary Heart Disease.

BACKGROUND: DNA methylation is implicated in coronary heart disease (CHD), but current evidence is based on small, cross-sectional studies. We examined blood DNA methylation in relation to incident CHD across multiple prospective cohorts. METHODS: Nine population-based cohorts from the United States and Europe profiled epigenome-wide blood leukocyte DNA methylation using the Illumina Infinium 450k microarray, and prospectively ascertained CHD events including coronary insufficiency/unstable angina, recognized myocardial infarction, coronary revascularization, and coronary death. Cohorts conducted race-specific analyses adjusted for age, sex, smoking, education, body mass index, blood cell type proportions, and technical variables. We conducted fixed-effect meta-analyses across cohorts. RESULTS: Among 11 461 individuals (mean age 64 years, 67% women, 35% African American) free of CHD at baseline, 1895 developed CHD during a mean follow-up of 11.2 years. Methylation levels at 52 CpG (cytosine-phosphate-guanine) sites were associated with incident CHD or myocardial infarction (false discovery rate<0.05). These CpGs map to genes with key roles in calcium regulation (ATP2B2, CASR, GUCA1B, HPCAL1), and genes identified in genome- and epigenome-wide studies of serum calcium (CASR), serum calcium-related risk of CHD (CASR), coronary artery calcified plaque (PTPRN2), and kidney function (CDH23, HPCAL1), among others. Mendelian randomization analyses supported a causal effect of DNA methylation on incident CHD; these CpGs map to active regulatory regions proximal to long non-coding RNA transcripts. CONCLUSION: Methylation of blood-derived DNA is associated with risk of future CHD across diverse populations and may serve as an informative tool for gaining further insight on the development of CHD.

DNA Methylation Analysis Identifies Loci for Blood Pressure Regulation.

Genome-wide association studies have identified hundreds of genetic variants associated with blood pressure (BP), but sequence variation accounts for a small fraction of the phenotypic variance. Epigenetic changes may alter the expression of genes involved in BP regulation and explain part of the missing heritability. We therefore conducted a two-stage meta-analysis of the cross-sectional associations of systolic and diastolic BP with blood-derived genome-wide DNA methylation measured on the Infinium HumanMethylation450 BeadChip in 17,010 individuals of European, African American, and Hispanic ancestry. Of 31 discovery-stage cytosine-phosphate-guanine (CpG) dinucleotides, 13 replicated after Bonferroni correction (discovery: N = 9,828, p < 1.0 × 10-7; replication: N = 7,182, p < 1.6 × 10-3). The replicated methylation sites are heritable (h2 > 30%) and independent of known BP genetic variants, explaining an additional 1.4% and 2.0% of the interindividual variation in systolic and diastolic BP, respectively. Bidirectional Mendelian randomization among up to 4,513 individuals of European ancestry from 4 cohorts suggested that methylation at cg08035323 (TAF1B-YWHAQ) influences BP, while BP influences methylation at cg00533891 (ZMIZ1), cg00574958 (CPT1A), and cg02711608 (SLC1A5). Gene expression analyses further identified six genes (TSPAN2, SLC7A11, UNC93B1, CPT1A, PTMS, and LPCAT3) with evidence of triangular associations between methylation, gene expression, and BP. Additional integrative Mendelian randomization analyses of gene expression and DNA methylation suggested that the expression of TSPAN2 is a putative mediator of association between DNA methylation at cg23999170 and BP. These findings suggest that heritable DNA methylation plays a role in regulating BP independently of previously known genetic variants.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels.

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Menopause accelerates biological aging.

Although epigenetic processes have been linked to aging and disease in other systems, it is not yet known whether they relate to reproductive aging. Recently, we developed a highly accurate epigenetic biomarker of age (known as the "epigenetic clock"), which is based on DNA methylation levels. Here we carry out an epigenetic clock analysis of blood, saliva, and buccal epithelium using data from four large studies: the Women's Health Initiative (n = 1,864); Invecchiare nel Chianti (n = 200); Parkinson's disease, Environment, and Genes (n = 256); and the United Kingdom Medical Research Council National Survey of Health and Development (n = 790). We find that increased epigenetic age acceleration in blood is significantly associated with earlier menopause (P = 0.00091), bilateral oophorectomy (P = 0.0018), and a longer time since menopause (P = 0.017). Conversely, epigenetic age acceleration in buccal epithelium and saliva do not relate to age at menopause; however, a higher epigenetic age in saliva is exhibited in women who undergo bilateral oophorectomy (P = 0.0079), while a lower epigenetic age in buccal epithelium was found for women who underwent menopausal hormone therapy (P = 0.00078). Using genetic data, we find evidence of coheritability between age at menopause and epigenetic age acceleration in blood. Using Mendelian randomization analysis, we find that two SNPs that are highly associated with age at menopause exhibit a significant association with epigenetic age acceleration. Overall, our Mendelian randomization approach and other lines of evidence suggest that menopause accelerates epigenetic aging of blood, but mechanistic studies will be needed to dissect cause-and-effect relationships further.

Characterization of TCF21 Downstream Target Regions Identifies a Transcriptional Network Linking Multiple Independent Coronary Artery Disease Loci.

To functionally link coronary artery disease (CAD) causal genes identified by genome wide association studies (GWAS), and to investigate the cellular and molecular mechanisms of atherosclerosis, we have used chromatin immunoprecipitation sequencing (ChIP-Seq) with the CAD associated transcription factor TCF21 in human coronary artery smooth muscle cells (HCASMC). Analysis of identified TCF21 target genes for enrichment of molecular and cellular annotation terms identified processes relevant to CAD pathophysiology, including "growth factor binding," "matrix interaction," and "smooth muscle contraction." We characterized the canonical binding sequence for TCF21 as CAGCTG, identified AP-1 binding sites in TCF21 peaks, and by conducting ChIP-Seq for JUN and JUND in HCASMC confirmed that there is significant overlap between TCF21 and AP-1 binding loci in this cell type. Expression quantitative trait variation mapped to target genes of TCF21 was significantly enriched among variants with low P-values in the GWAS analyses, suggesting a possible functional interaction between TCF21 binding and causal variants in other CAD disease loci. Separate enrichment analyses found over-representation of TCF21 target genes among CAD associated genes, and linkage disequilibrium between TCF21 peak variation and that found in GWAS loci, consistent with the hypothesis that TCF21 may affect disease risk through interaction with other disease associated loci. Interestingly, enrichment for TCF21 target genes was also found among other genome wide association phenotypes, including height and inflammatory bowel disease, suggesting a functional profile important for basic cellular processes in non-vascular tissues. Thus, data and analyses presented here suggest that study of GWAS transcription factors may be a highly useful approach to identifying disease gene interactions and thus pathways that may be relevant to complex disease etiology.

Leveraging information from genetic risk scores of coronary atherosclerosis.

PURPOSE OF REVIEW: Genome-wide association studies (GWAS) have identified ∼60 loci for coronary artery disease (CAD). Through genetic risk scores (GRSs), investigators are leveraging this genomic information to gain insights on both the fundamental mechanisms driving these associations as well as their utility in improving risk prediction. RECENT FINDINGS: GRSs of CAD track with the earliest atherosclerosis lesions in the coronary including fatty streaks and uncomplicated raised lesions. In multiple cohort studies, they predict incident CAD events independent of all traditional and lifestyle risk factors. The incorporation of SNPs with suggestive but not genome-wide association in GWAS into GRSs often increases the strength of these associations. GRS may also predict recurrent events and identify patients most likely to respond to statins. The effect of the GRS on discrimination metrics remains modest but the minimal degree of improvement needed for clinical utility is unknown. SUMMARY: Most novel loci for CAD identified through GWAS facilitate the formation of coronary atherosclerosis and stratify individuals based on their underlying burden of coronary atherosclerosis. GRSs may one day be routinely used in clinical practice to not only assess the risk of incident events but also to predict who will respond best to established prevention strategies.

Contemporary Considerations for Constructing a Genetic Risk Score: An Empirical Approach.

Genetic risk scores are an increasingly popular tool for summarizing the cumulative risk of a set of Single Nucleotide Polymorphisms (SNPs) with disease. Typically only the set of the SNPs that have reached genome-wide significance compose these scores. However recent work suggests that including additional SNPs may aid risk assessment. In this paper, we used the Atherosclerosis Risk in Communities (ARIC) Study cohort to illustrate how one can choose the optimal set of SNPs for a genetic risk score (GRS). In addition to P-value threshold, we also examined linkage disequilibrium, imputation quality, and imputation type. We provide a variety of evaluation metrics. Results suggest that P-value threshold had the greatest impact on GRS quality for the outcome of coronary heart disease, with an optimal threshold around 0.001. However, GRSs are relatively robust to both linkage disequilibrium and imputation quality. We also show that the optimal GRS partially depends on the evaluation metric and consequently the way one intends to use the GRS. Overall the implications highlight both the robustness of GRS and a means to empirically choose the best set of GRSs.

Genetic evidence for causal relationships between age at natural menopause and the risk of ageing-associated adverse health outcomes.

BACKGROUND: A later age at natural menopause (ANM) has been linked to several ageing-associated traits including an increased risk of breast and endometrial cancer and a decreased risk of lung cancer, osteoporosis and Alzheimer disease. However, ANM is also related to several proxies for overall health that may confound these associations. METHODS: We investigated the causal association of ANM with these clinical outcomes using Mendelian randomization (MR). Participants and outcomes analysed were restricted to post-menopausal females. We conducted a one-sample MR analysis in both the Women's Health Initiative and UK Biobank. We further analysed and integrated several additional data sets of post-menopausal women using a two-sample MR design. We used ≤55 genetic variants previously discovered to be associated with ANM as our instrumental variable. RESULTS: A 5-year increase in ANM was causally associated with a decreased risk of osteoporosis [odds ratio (OR) = 0.80, 95% CI (0.70-0.92)] and fractures (OR = 0.76, 95% CI, 0.62-0.94) as well as an increased risk of lung cancer (OR = 1.35, 95% CI, 1.06-1.71). Other associations including atherosclerosis-related outcomes were null. CONCLUSIONS: Our study confirms that the decline in bone density with menopause causally translates into fractures and osteoporosis. Additionally, this is the first causal epidemiological analysis to our knowledge to find an increased risk of lung cancer with increasing ANM. This finding is consistent with molecular and epidemiological studies suggesting oestrogen-dependent growth of lung tumours.

Genetic insights into resting heart rate and its role in cardiovascular disease.

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Local monitoring of SARS-CoV-2 variants in two large California counties in 2021.

Coronavirus Disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), continues to persist due to mutations resulting in newer, more infectious variants of concern. We aimed to leverage an ongoing private SARS-CoV-2 testing laboratory's infrastructure to monitor SARS-CoV-2 variants in two large California counties. Study enrollment was offered to adults aged 18 years or older in Los Angeles County and Riverside County who recently tested positive for SARS-CoV-2 with a polymerase chain reaction (PCR) assay. A cycle threshold value less than or equal to 30 cycles was considered a positive test for sequencing purposes. Within 5 days of study enrollment, clinician-monitored, self-collected oral fluid and anterior nares swab specimens were obtained from participants. Specimens were transported and stored at 8 °C or cooler. Samples underwent ribonucleic acid extraction, library preparation, and sequencing. SARS-CoV-2 lineages were identified using sequencing data. Participant and genomic data were analyzed using statistical tools and visualized with toolkits. The study was approved by Advarra Institutional Review Board (Pro00053729). From May 27, 2021 to September 9, 2021, 503 individuals were enrolled and underwent specimen collection. Of the 503 participants, 238 (47.3%) participants were women, 329 (63.6%) participants were vaccinated, and 221 (43.9%) participants were of Hispanic or Spanish origin. Of the cohort, 496 (98.6%) participants had symptoms at the time of collection. Among the 503 samples, 443 (88.1%) nasal specimens and 353 (70.2%) oral specimens yielded positive sequencing results. Over our study period, the prevalence of the Alpha variant of SARS-CoV-2 decreased (initially 23.1% [95% confidence interval (95% CI): 0-0.49%] to 0% [95% CI 0.0-0.0%]) as the prevalence of the Delta variant of SARS-CoV-2 increased (initially 33.3% [95% CI 0.0-100.0%] to 100.0% [95% CI 100.0-100.0%]). A strain that carried mutations of both Delta and Mu was identified. We found that outpatient SARS-CoV-2 variant surveillance could be conducted in a timely and accurate manner. The prevalence of different variants changed over time. A higher proportion of nasal specimens yielded results versus oral specimens. Timely and regional outpatient SARS-CoV-2 variant surveillance could be used for public health efforts to identify changes in SARS-CoV-2 strain epidemiology.

A multi-ethnic epigenome-wide association study of leukocyte DNA methylation and blood lipids.

Here we examine the association between DNA methylation in circulating leukocytes and blood lipids in a multi-ethnic sample of 16,265 subjects. We identify 148, 35, and 4 novel associations among Europeans, African Americans, and Hispanics, respectively, and an additional 186 novel associations through a trans-ethnic meta-analysis. We observe a high concordance in the direction of effects across racial/ethnic groups, a high correlation of effect sizes between high-density lipoprotein and triglycerides, a modest overlap of associations with epigenome-wide association studies of other cardio-metabolic traits, and a largely non-overlap with lipid loci identified to date through genome-wide association studies. Thirty CpGs reached significance in at least 2 racial/ethnic groups including 7 that showed association with the expression of an annotated gene. CpGs annotated to CPT1A showed evidence of being influenced by triglycerides levels. DNA methylation levels of circulating leukocytes show robust and consistent association with blood lipid levels across multiple racial/ethnic groups.

Genotype imputation and variability in polygenic risk score estimation.

BACKGROUND: Polygenic risk scores (PRSs) are a summarization of an individual's genetic risk for a disease or trait. These scores are being generated in research and commercial settings to study how they may be used to guide healthcare decisions. PRSs should be updated as genetic knowledgebases improve; however, no guidelines exist for their generation or updating. METHODS: Here, we characterize the variability introduced in PRS calculation by a common computational process used in their generation-genotype imputation. We evaluated PRS variability when performing genotype imputation using 3 different pre-phasing tools (Beagle, Eagle, SHAPEIT) and 2 different imputation tools (Beagle, Minimac4), relative to a WGS-based gold standard. Fourteen different PRSs spanning different disease architectures and PRS generation approaches were evaluated. RESULTS: We find that genotype imputation can introduce variability in calculated PRSs at the individual level without any change to the underlying genetic model. The degree of variability introduced by genotype imputation differs across algorithms, where pre-phasing algorithms with stochastic elements introduce the greatest degree of score variability. In most cases, PRS variability due to imputation is minor (< 5 percentile rank change) and does not influence the interpretation of the score. PRS percentile fluctuations are also reduced in the more informative tails of the PRS distribution. However, in rare instances, PRS instability at the individual level can result in singular PRS calculations that differ substantially from a whole genome sequence-based gold standard score. CONCLUSIONS: Our study highlights some challenges in applying population genetics tools to individual-level genetic analysis including return of results. Rare individual-level variability events are masked by a high degree of overall score reproducibility at the population level. In order to avoid PRS result fluctuations during updates, we suggest that deterministic imputation processes or the average of multiple iterations of stochastic imputation processes be used to generate and deliver PRS results.

Re-analysis of whole-exome sequencing data uncovers novel diagnostic variants and improves molecular diagnostic yields for sudden death and idiopathic diseases.

BACKGROUND: Whole-exome sequencing (WES) has become an efficient diagnostic test for patients with likely monogenic conditions such as rare idiopathic diseases or sudden unexplained death. Yet, many cases remain undiagnosed. Here, we report the added diagnostic yield achieved for 101 WES cases re-analyzed 1 to 7 years after initial analysis. METHODS: Of the 101 WES cases, 51 were rare idiopathic disease cases and 50 were postmortem "molecular autopsy" cases of early sudden unexplained death. Variants considered for reporting were prioritized and classified into three groups: (1) diagnostic variants, pathogenic and likely pathogenic variants in genes known to cause the phenotype of interest; (2) possibly diagnostic variants, possibly pathogenic variants in genes known to cause the phenotype of interest or pathogenic variants in genes possibly causing the phenotype of interest; and (3) variants of uncertain diagnostic significance, potentially deleterious variants in genes possibly causing the phenotype of interest. RESULTS: Initial analysis revealed diagnostic variants in 13 rare disease cases (25.4%) and 5 sudden death cases (10%). Re-analysis resulted in the identification of additional diagnostic variants in 3 rare disease cases (5.9%) and 1 sudden unexplained death case (2%), which increased our molecular diagnostic yield to 31.4% and 12%, respectively. CONCLUSIONS: The basis of new findings ranged from improvement in variant classification tools, updated genetic databases, and updated clinical phenotypes. Our findings highlight the potential for re-analysis to reveal diagnostic variants in cases that remain undiagnosed after initial WES.

An integrative cross-omics analysis of DNA methylation sites of glucose and insulin homeostasis.

Despite existing reports on differential DNA methylation in type 2 diabetes (T2D) and obesity, our understanding of its functional relevance remains limited. Here we show the effect of differential methylation in the early phases of T2D pathology by a blood-based epigenome-wide association study of 4808 non-diabetic Europeans in the discovery phase and 11,750 individuals in the replication. We identify CpGs in LETM1, RBM20, IRS2, MAN2A2 and the 1q25.3 region associated with fasting insulin, and in FCRL6, SLAMF1, APOBEC3H and the 15q26.1 region with fasting glucose. In silico cross-omics analyses highlight the role of differential methylation in the crosstalk between the adaptive immune system and glucose homeostasis. The differential methylation explains at least 16.9% of the association between obesity and insulin. Our study sheds light on the biological interactions between genetic variants driving differential methylation and gene expression in the early pathogenesis of T2D.

GWAS of epigenetic aging rates in blood reveals a critical role for TERT.

DNA methylation age is an accurate biomarker of chronological age and predicts lifespan, but its underlying molecular mechanisms are unknown. In this genome-wide association study of 9907 individuals, we find gene variants mapping to five loci associated with intrinsic epigenetic age acceleration (IEAA) and gene variants in three loci associated with extrinsic epigenetic age acceleration (EEAA). Mendelian randomization analysis suggests causal influences of menarche and menopause on IEAA and lipoproteins on IEAA and EEAA. Variants associated with longer leukocyte telomere length (LTL) in the telomerase reverse transcriptase gene (TERT) paradoxically confer higher IEAA (P < 2.7 × 10-11). Causal modeling indicates TERT-specific and independent effects on LTL and IEAA. Experimental hTERT-expression in primary human fibroblasts engenders a linear increase in DNA methylation age with cell population doubling number. Together, these findings indicate a critical role for hTERT in regulating the epigenetic clock, in addition to its established role of compensating for cell replication-dependent telomere shortening.

Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms.

Coronary artery disease (CAD) is a leading cause of morbidity and mortality worldwide. Although 58 genomic regions have been associated with CAD thus far, most of the heritability is unexplained, indicating that additional susceptibility loci await identification. An efficient discovery strategy may be larger-scale evaluation of promising associations suggested by genome-wide association studies (GWAS). Hence, we genotyped 56,309 participants using a targeted gene array derived from earlier GWAS results and performed meta-analysis of results with 194,427 participants previously genotyped, totaling 88,192 CAD cases and 162,544 controls. We identified 25 new SNP-CAD associations (P < 5 × 10-8, in fixed-effects meta-analysis) from 15 genomic regions, including SNPs in or near genes involved in cellular adhesion, leukocyte migration and atherosclerosis (PECAM1, rs1867624), coagulation and inflammation (PROCR, rs867186 (p.Ser219Gly)) and vascular smooth muscle cell differentiation (LMOD1, rs2820315). Correlation of these regions with cell-type-specific gene expression and plasma protein levels sheds light on potential disease mechanisms.