ABSTRACT
OBJECTIVE: The ability to apply results from a study to a broader population remains a primary objective in translational science. Distinct from intrinsic elements of scientific rigor, the extrinsic concept of generalization requires there be alignment between a study cohort and population in which results are expected to be applied. Widespread efforts have been made to quantify representativeness of study cohorts. These techniques, however, often consider the study and target cohorts as monolithic collections that can be directly compared. Overlooking known impacts to health from socio-demographic and environmental factors tied to individual's geographical location, and potentially obfuscating misalignment in underrepresented population subgroups. This manuscript introduces several measures to account for geographic information in the assessment of cohort representation. METHODS: Metrics were defined across two themes: First, measures of recruitment, to assess if a study cohort is drawn at an expected rate and in an expected geographical pattern with respect to individuals in a reference cohort. Second, measures of individual characteristics, to assess if the individuals in the study cohort accurately reflect the sociodemographic, clinical, and geographic diversity observed across a reference cohort while accounting for the geospatial proximity of individuals. RESULTS: As an empirical demonstration, methods are applied to an active clinical study examining asthma in Black/African American patients at a US Midwestern pediatric hospital. Results illustrate how areas of over- and under-recruitment can be identified and contextualized in light of study recruitment patterns at an individual-level, highlighting the ability to identify a subset of features for which the study cohort closely resembled the broader population. In addition they provide an opportunity to dive deeper into misalignments, to identify study cohort members that are in some way distinct from the communities for which they are expected to represent. CONCLUSION: Together, these metrics provide a comprehensive spatial assessment of a study cohort with respect to a broader target population. Such an approach offers researchers a toolset by which to target expected generalization of results derived from a given study.
ABSTRACT
Importance: The binary classification of spina bifida lesions as myelomeningocele (with sac) or myeloschisis (without sac) belies a spectrum of morphologies, which have not been correlated to clinical characteristics and outcomes. Objective: To characterize spina bifida lesion types and correlate them with preoperative presentation and postoperative outcomes. Design: Secondary analysis of images and videos obtained during fetoscopic spina bifida repair surgery from 2020-2023. Setting: Fetal surgery was performed at a quaternary care center. Participants: A prospective cohort of patients referred for fetal spina bifida underwent fetoscopic repair under an FDA-approved protocol. Of 60 lesions repaired, 57 had available images and were included in the analysis. Interventions or Exposures: We evaluated lesion morphology on high-resolution intraoperative images and videos to categorize lesions based on placode exposure and nerve root stretching. Main Outcomes and Measures: The reproducibility of the lesion classification was assessed via Kappa interrater agreement. Preoperative characteristics analyzed include ventricle size, tonsillar herniation level, lower extremities movement, and lesion dimensions. Outcomes included surgical time, need for patch for skin closure, gestational age at delivery, preterm premature rupture of membranes (PPROM), and neonatal cerebrospinal fluid (CSF) diversion. Results: We distinguished five lesion types that differ across a range of sac sizes, nerve root stretching, and placode exposure, with 93% agreement between examiners (p<0.001). Fetal characteristics at preoperative evaluation differed significantly by lesion type, including lesion volume (p<0.001), largest ventricle size (p=0.008), tonsillar herniation (p=0.005), and head circumference (p=0.03). Lesion level, talipes, and lower extremities movement did not differ by type. Surgical and perinatal outcomes differed by lesion type, including need for patch skin closure (p<0.001), gestational age at delivery (p=0.01), and NICU length of stay (p<0.001). PPROM, CSF leakage at birth, and CSF diversion in the NICU did not differ between lesion groups. Linear regression associated severity of ventriculomegaly with lesion type, but not with tonsillar herniation level. Conclusions and Relevance: There is a distinct phenotypic spectrum in open spina bifida with differential baseline presentation and outcomes. Severity of ventriculomegaly is associated with lesion type, rather than tonsillar herniation level. Our findings expand the classification of spina bifida to reveal a spectrum that warrants further study.
ABSTRACT
BACKGROUND: Type 2 diabetes (T2D) susceptibility is influenced by genetic and environmental factors. Previous findings suggest DNA methylation as a potential mechanism in T2D pathogenesis and progression. METHODS: We profiled DNA methylation in 248 blood samples from participants of European ancestry from 7 twin cohorts using a methylation sequencing platform targeting regulatory genomic regions encompassing 2,048,698 CpG sites. FINDINGS: We find and replicate 3 previously unreported T2D differentially methylated CpG positions (T2D-DMPs) at FDR 5% in RGL3, NGB and OTX2, and 20 signals at FDR 25%, of which 14 replicated. Integrating genetic variation and T2D-discordant monozygotic twin analyses, we identify both genetic-based and genetic-independent T2D-DMPs. The signals annotate to genes with established GWAS and EWAS links to T2D and its complications, including blood pressure (RGL3) and eye disease (OTX2). INTERPRETATION: The results help to improve our understanding of T2D disease pathogenesis and progression and may provide biomarkers for its complications. FUNDING: Funding acknowledgements for each cohort can be found in the Supplementary Note.
Subject(s)
CpG Islands , DNA Methylation , Diabetes Mellitus, Type 2 , Humans , Diabetes Mellitus, Type 2/genetics , Female , Male , Genome-Wide Association Study , Genetic Predisposition to Disease , Middle Aged , Epigenesis, Genetic , Otx Transcription Factors/genetics , Otx Transcription Factors/metabolism , Diabetes Complications/genetics , Gene Expression ProfilingABSTRACT
Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.
ABSTRACT
OBJECTIVE: To examine whether participants with different levels of diabetes-related DNA methylation at ABCG1 might respond differently to dietary weight loss interventions with long-term changes in adiposity and body fat distribution. RESEARCH DESIGN AND METHODS: The current study included overweight/obese participants from the POUNDS Lost trial. Blood levels of regional DNA methylation at ABCG1 were profiled by high-resolution methylC-capture sequencing at baseline among 673 participants, of whom 598 were followed up at 6 months and 543 at 2 years. Two-year changes in adiposity and computed tomography-measured body fat distribution were calculated. RESULTS: Regional DNA methylation at ABCG1 showed significantly different associations with long-term changes in body weight and waist circumference at 6 months and 2 years in dietary interventions varying in protein intake (interaction P < 0.05 for all). Among participants assigned to an average-protein (15%) diet, lower baseline regional DNA methylation at ABCG1 was associated with greater reductions in body weight and waist circumference at 6 months and 2 years, whereas opposite associations were found among those assigned to a high-protein (25%) diet. Similar interaction patterns were also observed for body fat distribution, including visceral adipose tissue, subcutaneous adipose tissue, deep subcutaneous adipose tissue, and total adipose tissue at 6 months and 2 years (interaction P < 0.05 for all). CONCLUSIONS: Baseline DNA methylation at ABCG1 interacted with dietary protein intake on long-term decreases in adiposity and body fat distribution. Participants with lower methylation at ABCG1 benefitted more in long-term reductions in body weight, waist circumference, and body fat distribution when consuming an average-protein diet.
Subject(s)
Adiposity , DNA Methylation , Humans , Adiposity/genetics , DNA Methylation/genetics , Dietary Proteins , Diet, Reducing , Obesity/genetics , Body Weight/genetics , Waist Circumference , Body Mass Index , ATP Binding Cassette Transporter, Subfamily G, Member 1/geneticsABSTRACT
There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.
Subject(s)
Osteoarthritis , Osteoclasts , Humans , Osteoclasts/metabolism , Genome-Wide Association Study/methods , Genetic Predisposition to Disease , Gene Expression Regulation , Osteoarthritis/genetics , Osteoarthritis/metabolismABSTRACT
The extravillous trophoblast cell lineage is a key feature of placentation and successful pregnancy. Knowledge of transcriptional regulation driving extravillous trophoblast cell development is limited. Here, we map the transcriptome and epigenome landscape as well as chromatin interactions of human trophoblast stem cells and their transition into extravillous trophoblast cells. We show that integrating chromatin accessibility, long-range chromatin interactions, transcriptomic, and transcription factor binding motif enrichment enables identification of transcription factors and regulatory mechanisms critical for extravillous trophoblast cell development. We elucidate functional roles for TFAP2C, SNAI1, and EPAS1 in the regulation of extravillous trophoblast cell development. EPAS1 is identified as an upstream regulator of key extravillous trophoblast cell transcription factors, including ASCL2 and SNAI1 and together with its target genes, is linked to pregnancy loss and birth weight. Collectively, we reveal activation of a dynamic regulatory network and provide a framework for understanding extravillous trophoblast cell specification in trophoblast cell lineage development and human placentation.
Subject(s)
Chromatin , Trophoblasts , Pregnancy , Female , Humans , Trophoblasts/metabolism , Chromatin/genetics , Chromatin/metabolism , Placentation/genetics , Cell Differentiation/genetics , Transcription Factors/genetics , Transcription Factors/metabolism , Cell Lineage/genetics , Placenta/metabolism , Basic Helix-Loop-Helix Transcription Factors/genetics , Basic Helix-Loop-Helix Transcription Factors/metabolismABSTRACT
Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.
Subject(s)
DNA Methylation , Rare Diseases , Humans , Haplotypes , Rare Diseases/genetics , DNA Methylation/genetics , Sequence Analysis, DNA , Base Sequence , High-Throughput Nucleotide Sequencing , Nerve Tissue Proteins/geneticsABSTRACT
BACKGROUND: DNA methylation (DNAm) may play a critical role in bridging prenatal adverse events and cardiometabolic disorders including hypertension in later life. METHODS: We included 672 adult participants with overweight or obesity, who participated in a 2-year randomized weight-loss dietary intervention study. We defined the regional DNAm levels as the average methylation level of 5'-cytosine-phosphate-guanine-3' within 500 bp of LINC00319 (cg01820192), ATP2B1 (cg00508575), and LMNA (cg12593793), respectively. Generalized linear regression models were used to assess the association between the regional DNAm and 2-year blood pressure changes. Trajectory analysis was used to identify subgroups that shared a similar underlying trajectory of 2-year blood pressure changes. RESULTS: The regional DNAm at LINC00319, showed significantly different associations with 2-year changes in systolic blood pressure and diastolic blood pressure among participants assigned to low- or high-fat diets (P for interaction<0.05 for all). In response to the low-fat diet, per SD higher regional DNAm at LINC00319 was associated with greater reductions in both 2-year changes in systolic blood pressure (ß, -1.481; P=0.020) and diastolic blood pressure (ß, -1.096; P=0.009). Three trajectories of changes in systolic blood pressure or diastolic blood pressure were identified, and participants with higher regional DNAm at LINC00319 were more likely to experience and maintain decreased systolic blood pressure and diastolic blood pressure (odds ratio of being in decrease-stable versus stable [95% CI], 1.542 [1.146-2.076] and 1.463 [1.125-1.902]). CONCLUSIONS: Our findings suggest that DNAm could be a metabolic memory bridging early and later life, and an indicator of more benefits from eating a low-fat weight-loss diet.
Subject(s)
DNA Methylation , Obesity , Adult , Humans , Blood Pressure/genetics , Birth Weight , Obesity/genetics , Weight Loss/genetics , Diet, Fat-Restricted , Plasma Membrane Calcium-Transporting ATPasesABSTRACT
Anabolic treatment options for osteoporosis remain limited. One approach to discovering novel anabolic drug targets is to identify genetic causes of extreme high bone mass (HBM). We investigated a pedigree with unexplained HBM within the UK HBM study, a national cohort of probands with HBM and their relatives. Whole exome sequencing (WES) in a family with HBM identified a rare heterozygous missense variant (NM_004482.4:c.1657C > T, p.Arg553Trp) in GALNT3, segregating appropriately. Interrogation of data from the UK HBM study and the Anglo-Australasian Osteoporosis Genetics Consortium (AOGC) revealed an unrelated individual with HBM with another rare heterozygous variant (NM_004482.4:c.831 T > A, p.Asp277Glu) within the same gene. In silico protein modeling predicted that p.Arg553Trp would disrupt salt-bridge interactions, causing instability of GALNT3, and that p.Asp277Glu would disrupt manganese binding and consequently GALNT3 catalytic function. Bi-allelic loss-of-function GALNT3 mutations alter FGF23 metabolism, resulting in hyperphosphatemia and causing familial tumoral calcinosis (FTC). However, bone mineral density (BMD) in FTC cases, when reported, has been either normal or low. Common variants in the GALNT3 locus show genome-wide significant associations with lumbar, femoral neck, and total body BMD. However, no significant associations with BMD are observed at loci coding for FGF23, its receptor FGFR1, or coreceptor klotho. Mendelian randomization analysis, using expression quantitative trait loci (eQTL) data from primary human osteoblasts and genome-wide association studies data from UK Biobank, suggested increased expression of GALNT3 reduces total body, lumbar spine, and femoral neck BMD but has no effect on phosphate concentrations. In conclusion, rare heterozygous loss-of-function variants in GALNT3 may cause HBM without altering phosphate concentration. These findings suggest that GALNT3 may affect BMD through pathways other than FGF23 regulation, the identification of which may yield novel anabolic drug targets for osteoporosis. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).
Subject(s)
Bone Density , Osteoporosis , Humans , Bone Density/genetics , Genome-Wide Association Study , Lumbar Vertebrae/physiology , Osteoporosis/genetics , PhosphatesABSTRACT
CONTEXT: Carnitine palmitoyltransferase-1A, encoded by the CPT1A gene, plays a key role in the oxidation of long-chain fatty acids in the mitochondria and may be important in triglyceride metabolism. Previous work has shown that high fat intake was negatively associated with CPT1A methylation and positively associated with CPT1A expression. OBJECTIVE: We aim to investigate the association of DNA methylation (DNAm) at the CPT1A gene with reductions in triglycerides and triglyceride-rich lipoproteins (TRLs) in response to weight-loss diet interventions. METHODS: The current study included 538 White participants, who were randomly assigned to 1 of 4 diets varying in macronutrient components. We defined the regional DNAm at CPT1A as the average methylation level over CpGs within 500 bp of the 3 triglyceride-related DNAm sites. RESULTS: Dietary fat intake significantly modified the association between baseline DNAm at CPT1A and 2-year changes in total plasma triglycerides, independent of concurrent weight loss. Among participants assigned to a low-fat diet, a higher regional DNAm level at CPT1A was associated with a greater reduction in total plasma triglycerides at 2 years (P = .01), compared with those assigned to a high-fat diet (P = .64) (P interaction = .018). Further investigation on lipids and apolipoproteins in very low-density lipoprotein (VLDL) revealed similar interaction patterns for 2-year changes in VLDL-triglycerides, VLDL-cholesterol, and VLDL-apolipoprotein B (P interaction = .009, .002, and .016, respectively), but not for VLDL-apoC-III (P interaction = .36). CONCLUSION: Participants with a higher regional DNAm level at CPT1A benefit more in long-term improvement in triglycerides, particularly in the TRLs and related apolipoproteins when consuming a low-fat weight-loss diet.
Subject(s)
Carnitine O-Palmitoyltransferase , DNA Methylation , Humans , Apolipoproteins/genetics , Apolipoproteins/metabolism , Carnitine O-Palmitoyltransferase/genetics , Diet, Reducing , Lipoproteins , Lipoproteins, LDL , Lipoproteins, VLDL , TriglyceridesABSTRACT
PURPOSE: This study aimed to assess the amount and types of clinical genetic testing denied by insurance and the rate of diagnostic and candidate genetic findings identified through research in patients who faced insurance denials. METHODS: Analysis consisted of review of insurance denials in 801 patients enrolled in a pediatric genomic research repository with either no previous genetic testing or previous negative genetic testing result identified through cross-referencing with insurance prior-authorizations in patient medical records. Patients and denials were also categorized by type of insurance coverage. Diagnostic findings and candidate genetic findings in these groups were determined through review of our internal variant database and patient charts. RESULTS: Of the 801 patients analyzed, 147 had insurance prior-authorization denials on record (18.3%). Exome sequencing and microarray were the most frequently denied genetic tests. Private insurance was significantly more likely to deny testing than public insurance (odds ratio = 2.03 [95% CI = 1.38-2.99] P = .0003). Of the 147 patients with insurance denials, 53.7% had at least 1 diagnostic or candidate finding and 10.9% specifically had a clinically diagnostic finding. Fifty percent of patients with clinically diagnostic results had immediate medical management changes (5.4% of all patients experiencing denials). CONCLUSION: Many patients face a major barrier to genetic testing in the form of lack of insurance coverage. A number of these patients have clinically diagnostic findings with medical management implications that would not have been identified without access to research testing. These findings support re-evaluation of insurance carriers' coverage policies.
Subject(s)
Genomics , Insurance Coverage , Child , HumansABSTRACT
Hyperoxia disrupts lung development in mice and causes bronchopulmonary dysplasia (BPD) in neonates. To investigate sex-dependent molecular and cellular programming involved in hyperoxia, we surveyed the mouse lung using single cell RNA sequencing (scRNA-seq), and validated our findings in human neonatal lung cells in vitro. Hyperoxia-induced inflammation in alveolar type (AT) 2 cells gave rise to damage-associated transient progenitors (DATPs). It also induced a new subpopulation of AT1 cells with reduced expression of growth factors normally secreted by AT1 cells, but increased mitochondrial gene expression. Female alveolar epithelial cells had less EMT and pulmonary fibrosis signaling in hyperoxia. In the endothelium, expansion of Car4+ EC (Cap2) was seen in hyperoxia along with an emergent subpopulation of Cap2 with repressed VEGF signaling. This regenerative response was increased in females exposed to hyperoxia. Mesenchymal cells had inflammatory signatures in hyperoxia, with a new distal interstitial fibroblast subcluster characterized by repressed lipid biosynthesis and a transcriptomic signature resembling myofibroblasts. Hyperoxia-induced gene expression signatures in human neonatal fibroblasts and alveolar epithelial cells in vitro resembled mouse scRNA-seq data. These findings suggest that neonatal exposure to hyperoxia programs distinct sex-specific stem cell progenitor and cellular reparative responses that underpin lung remodeling in BPD.
Subject(s)
Bronchopulmonary Dysplasia , Hyperoxia , Infant, Newborn , Male , Female , Animals , Mice , Humans , Bronchopulmonary Dysplasia/metabolism , Transcriptome/genetics , Hyperoxia/metabolism , Animals, Newborn , Lung/metabolism , Disease Models, AnimalABSTRACT
Transcriptional regulation of liver transplant (LT) rejection may reveal novel predictive and therapeutic targets. The purpose of this article is to test the role of differential DNA methylation in children with biopsy-proven acute cellular rejection after LT. Methods: Paired peripheral blood DNA samples were obtained before and after LT from 17 children, including 4 rejectors (Rs) and 13 nonrejectors (NRs), and assayed with MethylC capture sequencing approach covering 5 million CpGs in immune-cell-specific regulatory elements. Differentially methylated CpGs (DMCs) were identified using generalized linear regression models adjusting for sex and age and merged into differentially methylated regions (DMRs) comprising 3 or more DMCs. Results: Contrasting Rs versus NRs, we identified 2238 DMCs in post-LT and 2620 DMCs in pre-LT samples, which clustered in 216 and 282 DMRs, respectively. DMCs associated with R were enriched in enhancers and depleted in promoters. Among DMRs, the proportion of hypomethylated DMRs increased from 61/282 (22%) in pre-LT to 103/216 (48%, P < 0.0001) in post-LT samples. The highest-ranked biological processes enriched in post-LT DMCs were antigen processing and presentation via major histocompatibility complex (MHC) class I, MHC class I complex, and peptide binding (P < 7.92 × 10-17), respectively. Top-ranked DMRs mapped to genes that mediate B-cell receptor signaling (ADAP1) or regulate several immune cells (ARRB2) (P < 3.75 × 10-08). DMRs in MHC class I genes were enriched for single nucleotide polymorphisms (SNPs), which bind transcription factors, affect gene expression and splicing, or alter peptide-binding amino acid sequences. Conclusions: Dynamic methylation in distal regulatory regions reveals known transplant-relevant MHC-dependent rejection pathways and identifies novel loci for future mechanistic evaluations in pediatric transplant subcohorts.
ABSTRACT
Studies in genetically 'identical' individuals indicate that as much as 50% of complex trait variation cannot be traced to genetics or to the environment. The mechanisms that generate this 'unexplained' phenotypic variation (UPV) remain largely unknown. Here, we identify neuronatin (NNAT) as a conserved factor that buffers against UPV. We find that Nnat deficiency in isogenic mice triggers the emergence of a bi-stable polyphenism, where littermates emerge into adulthood either 'normal' or 'overgrown'. Mechanistically, this is mediated by an insulin-dependent overgrowth that arises from histone deacetylase (HDAC)-dependent ß-cell hyperproliferation. A multi-dimensional analysis of monozygotic twin discordance reveals the existence of two patterns of human UPV, one of which (Type B) phenocopies the NNAT-buffered polyphenism identified in mice. Specifically, Type-B monozygotic co-twins exhibit coordinated increases in fat and lean mass across the body; decreased NNAT expression; increased HDAC-responsive gene signatures; and clinical outcomes linked to insulinemia. Critically, the Type-B UPV signature stratifies both childhood and adult cohorts into four metabolic states, including two phenotypically and molecularly distinct types of obesity.
Subject(s)
Membrane Proteins , Nerve Tissue Proteins , Adaptation, Physiological , Adult , Animals , Child , Histone Deacetylases , Humans , Insulin , Membrane Proteins/metabolism , Mice , Nerve Tissue Proteins/genetics , Obesity/genetics , Obesity/metabolismABSTRACT
BACKGROUND: There is considerable evidence for the importance of the DNA methylome in metabolic health, for example, a robust methylation signature has been associated with body mass index (BMI). However, visceral fat (VF) mass accumulation is a greater risk factor for metabolic disease than BMI alone. In this study, we dissect the subcutaneous adipose tissue (SAT) methylome signature relevant to metabolic health by focusing on VF as the major risk factor of metabolic disease. We integrate results with genetic, blood methylation, SAT gene expression, blood metabolomic, dietary intake and metabolic phenotype data to assess and quantify genetic and environmental drivers of the identified signals, as well as their potential functional roles. METHODS: Epigenome-wide association analyses were carried out to determine visceral fat mass-associated differentially methylated positions (VF-DMPs) in SAT samples from 538 TwinsUK participants. Validation and replication were performed in 333 individuals from 3 independent cohorts. To assess functional impacts of the VF-DMPs, the association between VF and gene expression was determined at the genes annotated to the VF-DMPs and an association analysis was carried out to determine whether methylation at the VF-DMPs is associated with gene expression. Further epigenetic analyses were carried out to compare methylation levels at the VF-DMPs as the response variables and a range of different metabolic health phenotypes including android:gynoid fat ratio (AGR), lipids, blood metabolomic profiles, insulin resistance, T2D and dietary intake variables. The results from all analyses were integrated to identify signals that exhibit altered SAT function and have strong relevance to metabolic health. RESULTS: We identified 1181 CpG positions in 788 genes to be differentially methylated with VF (VF-DMPs) with significant enrichment in the insulin signalling pathway. Follow-up cross-omic analysis of VF-DMPs integrating genetics, gene expression, metabolomics, diet, and metabolic traits highlighted VF-DMPs located in 9 genes with strong relevance to metabolic disease mechanisms, with replication of signals in FASN, SREBF1, TAGLN2, PC and CFAP410. PC methylation showed evidence for mediating effects of diet on VF. FASN DNA methylation exhibited putative causal effects on VF that were also strongly associated with insulin resistance and methylation levels in FASN better classified insulin resistance (AUC=0.91) than BMI or VF alone. CONCLUSIONS: Our findings help characterise the adiposity-associated methylation signature of SAT, with insights for metabolic disease risk.
Subject(s)
Insulin Resistance , Body Mass Index , DNA Methylation , Diet , Epigenesis, Genetic , Epigenome , Humans , Insulin Resistance/geneticsABSTRACT
Circulating levels of the amino acid glutamate are associated with central fat accumulation, yet the pathophysiology of this relationship remains unknown. We aimed to (i) refine and validate the association between circulating glutamate and abdominal obesity in a large twin cohort, and (ii) investigate whether transcriptomic profiles in adipose tissue could provide insight into the biological mechanisms underlying the association. First, in a cohort of 4665 individuals from the TwinsUK resource, we identified individuals with abdominal obesity and compared prevalence of the latter across circulating glutamate quintiles. Second, we used transcriptomic signatures generated from adipose tissue, both subcutaneous and visceral, to investigate associations with circulating glutamate levels. Individuals in the top circulating glutamate quintile had a sevenfold higher prevalence of abdominal obesity compared to those in the bottom quintile. The adipose tissue transcriptomic analyses identified GLUL, encoding Glutamate-Ammonia Ligase, as being associated with circulating glutamate and abdominal obesity, with pronounced signatures in the visceral depot. In conclusion, circulating glutamate is positively associated with the prevalence of abdominal obesity which relates to dysregulated GLUL expression specifically in visceral adipose tissue.
Subject(s)
Glutamic Acid , Obesity, Abdominal , Adipose Tissue/metabolism , Body Mass Index , Gene Expression , Humans , Obesity/metabolism , Obesity, Abdominal/geneticsABSTRACT
SARS-CoV-2 is a novel betacoronavirus that caused coronavirus disease 2019 and has resulted in millions of deaths worldwide. Novel coronavirus infections in humans have steadily become more common. Understanding antibody responses to SARS-CoV-2, and identifying conserved, cross-reactive epitopes among coronavirus strains could inform the design of vaccines and therapeutics with broad application. Here, we determined that individuals with previous SARS-CoV-2 infection or vaccinated with the Pfizer-BioNTech BNT162b2 vaccine produced antibody responses that cross-reacted with related betacoronaviruses. Moreover, we designed a peptide-conjugate vaccine with a conserved SARS-CoV-2 S2 spike epitope, immunized mice and determined cross-reactive antibody binding to SARS-CoV-2 and other related coronaviruses. This conserved spike epitope also shared sequence homology to proteins in commensal gut microbiota and could prime immune responses in humans. Thus, SARS-CoV-2 conserved epitopes elicit cross-reactive immune responses to both related coronaviruses and host bacteria that could serve as future targets for broad coronavirus therapeutics and vaccines.
Subject(s)
COVID-19 , Spike Glycoprotein, Coronavirus , Animals , Antibodies, Viral , BNT162 Vaccine , COVID-19/prevention & control , Epitopes , Humans , Mice , SARS-CoV-2 , VaccinationABSTRACT
PURPOSE: This study aimed to provide comprehensive diagnostic and candidate analyses in a pediatric rare disease cohort through the Genomic Answers for Kids program. METHODS: Extensive analyses of 960 families with suspected genetic disorders included short-read exome sequencing and short-read genome sequencing (srGS); PacBio HiFi long-read genome sequencing (HiFi-GS); variant calling for single nucleotide variants (SNV), structural variant (SV), and repeat variants; and machine-learning variant prioritization. Structured phenotypes, prioritized variants, and pedigrees were stored in PhenoTips database, with data sharing through controlled access the database of Genotypes and Phenotypes. RESULTS: Diagnostic rates ranged from 11% in patients with prior negative genetic testing to 34.5% in naive patients. Incorporating SVs from genome sequencing added up to 13% of new diagnoses in previously unsolved cases. HiFi-GS yielded increased discovery rate with >4-fold more rare coding SVs compared with srGS. Variants and genes of unknown significance remain the most common finding (58% of nondiagnostic cases). CONCLUSION: Computational prioritization is efficient for diagnostic SNVs. Thorough identification of non-SNVs remains challenging and is partly mitigated using HiFi-GS sequencing. Importantly, community research is supported by sharing real-time data to accelerate gene validation and by providing HiFi variant (SNV/SV) resources from >1000 human alleles to facilitate implementation of new sequencing platforms for rare disease diagnoses.
Subject(s)
Genomics , Rare Diseases , Child , Genome , High-Throughput Nucleotide Sequencing , Humans , Pedigree , Rare Diseases/diagnosis , Rare Diseases/genetics , Sequence Analysis, DNAABSTRACT
BACKGROUND: Thioredoxin Interacting Protein (TXNIP) functions as a master regulator for glucose homeostasis. Hypomethylation at the 5'-cytosine-phosphate-guanine-3' (CpG) site cg19693031 of TXNIP has been consistently related to islet dysfunction, hyperglycemia, and type 2 diabetes. DNA methylation (DNAm) may reveal the missing mechanistic link between obesity and type 2 diabetes. We hypothesize that baseline DNAm level at TXNIP in blood may be associated with glycemic traits and their changes in response to weight-loss diet interventions. METHODS: We included 639 adult participants with overweight or obesity, who participated in a 2-year randomized weight-loss diet intervention. Baseline blood DNAm levels were profiled by high-resolution methylC-capture sequencing. We defined the regional DNAm level of TXNIP as the average methylation level over CpGs within 500 bp of cg19693031. Generalized linear regression models were used for main analyses. RESULTS: We found that higher regional DNAm at TXNIP was significantly correlated with lower fasting glucose, HbA1c, and Homeostatic Model Assessment of Insulin Resistance (HOMA-IR) at baseline (P < 0.05 for all). Significant interactions were observed between dietary protein intake and DNAm on changes in insulin (P-interaction = 0.007) and HOMA-IR (P-interaction = 0.009) at 6 months. In participants with the highest tertile of regional DNAm at TXNIP, average protein (15%) intake was associated with a greater reduction in insulin (ß: -0.14; 95% CI: -0.24, -0.03; P = 0.011) and HOMA-IR (ß: -0.15; 95% CI: -0.26, -0.03; P = 0.014) than high protein (25%) intake, whereas no significant associations were found in those with the lower tertiles (P > 0.05). The interaction was attenuated to be non-significant at 2 years, presumably related to decreasing adherence to the diet intervention. CONCLUSIONS: Our data indicate that higher regional DNAm level at TXNIP was significantly associated with better fasting glucose, HbA1c, and HOMA-IR; and people with higher regional DNAm levels benefited more in insulin and HOMA-IR improvement by taking the average-protein weight-loss diet.