👤 David Roberts

Mild traumatic brain injury can disrupt brain function and is associated with high morbidity and healthcare utilization. While many individuals recover from mild traumatic brain injury, a significant proportion experience long-term sequelae, collectively known as post-concussion syndrome. Symptoms of post-concussion syndrome include headache, dizziness, insomnia, cognitive processing difficulties and mental health disturbances. The disease burden is augmented by the current lack of objective measures to accurately predict long-term symptoms and deficits, providing an opportunity to utilize biomarkers in biofluids. A large proportion of available diagnostic clinical tools are subjective symptom scores. This review aims to explore current fluid biomarkers, grouped by clinical symptoms. With the available literature, we have discovered a wide range of fluid biomarkers that have been investigated for predicting post-traumatic headache, including neuropeptides; sleep disturbances, such as cortisol and melatonin; vestibular disturbances, including interleukin-6 and neurone-specific enolase; and vomiting, such as S100B. Along with physical symptoms, biomarkers investigated for predicting cognitive disturbances include inflammatory markers, S100B, neurofilament light chain, tau, microRNA and hormones. Biomarkers to predict mental health disturbances may include brain-derived neurotrophic factor, tau and cortisol. By utilizing such biomarkers, there is capacity to adopt a personalized medicine approach to facilitate early interventions for those most in need while also identifying individuals with a favourable prognosis who can safely return to their normal activities. Show less

To understand the relative contributions of 5' UTR elements to translation output, we performed a comprehensive analysis of upstream open reading frames (uORFs) in the 5' UTRs of the BDNF (brain-derived neurotrophic factor) transcripts. Predicted uORFs were identified in 14 out of 17 BDNF RefSeq transcript isoforms, and we experimentally validated five of these transcripts as being uORF-repressed, suggesting that uORF elements play an important role in shaping the protein output from this locus. We explored several approaches to disrupt BDNF uORF function. Deletion of a 5' UTR exon in BDNF v11 (containing eight predicted uORFs), in order to simulate an exon skipping outcome, resulted in pronounced upregulation in a reporter construct system. This effect was found to be partially uORF-dependent but was also dependent on the disruption of an RNA secondary structure element. However, this transcript variant was found to not be expressed in human brain. Conversely, direct disruption of a single uORF start codon in the widely expressed BDNF v4 transcript variant using an adenine base editing approach resulted in a ∼1.8-fold upregulation of endogenous BDNF protein expression in cell culture. This study characterizes uORF-mediated regulation of the BDNF locus and demonstrates the potential for BDNF protein upregulation via base editing-mediated uORF disruption. Show less

Single-cell RNA-sequencing has identified that Alzheimer's disease (AD) pathology in humans is associated with activation of disease-associated microglia (DAM). Microglial signatures of human AD have not been consistently identified in AD mouse models. Since the inflammatory response of rats is more like humans, we profiled microglial transcriptomes in aging TgF344-AD rats, which overexpress two human AD risk genes. Classic DAM gene activation ( Show less

Fatty acids are important as structural components, energy sources, and signaling mediators. While studies have extensively explored genetic regulation of fatty acids in serum and other bodily fluids, their regulation within adipose tissue, a crucial regulator of cardiovascular and metabolic health, remains unclear. Here, we investigated the genetic regulation of 18 fatty acids in subcutaneous adipose tissue from 569 female twins from TwinsUK. Using twin models, the heritability of fatty acids ranged from 5% to 59%, indicating a substantial genetic regulation of fatty acid levels within adipose tissue, which was also tissue specific in many cases. Genome-wide association studies identified 10 significant loci, in SCD, ADAMTSL1, ZBTB41, SNTB1, EXOC6B, ACSL3, LINC02055, MKRN2/TSEN2, FADS1, and HAPLN across 13 fatty acids or fatty acid product-to-precursor ratios. Using adipose gene expression and methylation, which were concurrently measured in these samples, we detected five fatty acid-associated signals that colocalized with expression quantitative trait locus (eQTL) and methylation quantitative trait locus (meQTL) signals, highlighting fatty acids that are regulated by molecular processes within adipose tissue. We explored links between polygenic scores of common metabolic traits and adipose fatty acid levels and identified associations between polygenic scores of BMI, body-fat distribution, and triglycerides and several fatty acids, indicating these risk scores impact local adipose tissue content. Overall, our results identified local genetic regulation of fatty acids within adipose tissue and highlighted their links with renal and cardio-metabolic health. Show less

Despite a vast literature on the role of macrophages in wound healing, the role of dermal monocyte (Mo)-derived antigen presenting cells (APC) has received scant attention. Using scRNAseq and flow cytometry, we identify a population of APC that is prominent in wounds of non-diabetic mice but is reduced in wounds of diabetic mice. Using adoptive transfer experiments and Show less

Dual agonists targeting glucagon-like peptide-1 receptor (GLP1R) and glucose-dependent insulinotropic polypeptide receptor (GIPR) are breakthrough treatments for patients with type 2 diabetes and obesity. Compared to GLP1R agonists, dual agonists show superior efficacy for glucose lowering and weight reduction. However, delineation of dual agonist cell targets remains challenging. Here, we develop and test daLUXendin and daLUXendin+, non-lipidated and lipidated fluorescent GLP1R/GIPR dual agonist probes, and use them to visualize cellular targets. daLUXendins are potent GLP1R/GIPR dual agonists that advantageously show less functional selectivity for mouse GLP1R over mouse GIPR. daLUXendins label rodent and human pancreatic islet cells, with a signal intensity of β cells > α cells = δ cells. Systemic administration of daLUXendin strongly labels GLP1R Show less

The next generation of obesity medicines harness the activity of the glucose-dependent insulinotropic polypeptide and glucagon-like peptide 1 receptors (GIPR and GLP-1R), but their mechanism of action remains unclear. Here, we report that the GIPR is enriched in oligodendrocytes and GIPR signaling bidirectionally regulates oligodendrogenesis. In mice with adult-onset deletion of GIPR in oligodendrocytes, GIPR agonism fails to enhance the weight-loss effects of GLP-1R agonism. Mechanistically, GIPR agonism increases brain access of GLP-1R agonists, and GIPR signaling in oligodendrocytes is required for this effect. In addition, we show that vasopressin neurons of the paraventricular hypothalamus are necessary for the weight-loss response to GLP-1R activation, targeted by peripherally administered GLP-1R agonists via their axonal compartment, and this access is increased by activation of the GIPR in oligodendrocytes. Collectively, our findings identify a novel mechanism by which incretin therapies may function to promote synergistic weight loss in the management of excess adiposity. Show less

Glucose-dependent insulinotropic polypeptide (GIP) was the first incretin identified and plays an essential role in the maintenance of glucose tolerance in healthy humans. Until recently GIP had not been developed as a therapeutic and thus has been overshadowed by the other incretin, glucagon-like peptide 1 (GLP-1), which is the basis for several successful drugs to treat diabetes and obesity. However, there has been a rekindling of interest in GIP biology in recent years, in great part due to pharmacology demonstrating that both GIPR agonism and antagonism may be beneficial in treating obesity and diabetes. This apparent paradox has reinvigorated the field, led to new lines of investigation, and deeper understanding of GIP. In this review, we provide a detailed overview on the multifaceted nature of GIP biology and discuss the therapeutic implications of GIPR signal modification on various diseases. Following its classification as an incretin hormone, GIP has emerged as a pleiotropic hormone with a variety of metabolic effects outside the endocrine pancreas. The numerous beneficial effects of GIPR signal modification render the peptide an interesting candidate for the development of pharmacotherapies to treat obesity, diabetes, drug-induced nausea and both bone and neurodegenerative disorders. Show less

The apolipoprotein E (APOE) gene is the best established genetic risk factor for Alzheimer's disease in later life, with the ε4 allele conferring higher risk. APOE disclosure is becoming increasingly common in the clinical care of people with Alzheimer's disease and in cognitively unimpaired adults. In this study, we aimed to describe changes in measures of genetic disease knowledge and psychiatric symptoms following APOE disclosure to cognitively unimpaired adults. Data were collected as part of the screening phase of the global, multicentre, Alzheimer's Prevention Initiative Generation Study 1 (NCT02565511). Eligible individuals were cognitively unimpaired (Mini-Mental State Exam total score ≥24), aged 60-75 years, and psychologically pre-screened for readiness (by measures of depressive symptoms and anxiety) to receive their APOE genotype from a health-care provider. Participants were assessed before disclosure, and 2-7 days, 6 weeks, 6 months, and 12 months after disclosure. Multivariable linear and ordinal logistic regressions were used to compare changes in genetic disease knowledge, anxiety, depression, and distress by APOE4 genotype status, adjusting for key covariates, with a focus on 2-7 days after disclosure. Multiple imputation by chained equations methods was used to account for missing outcome data. The trial took place between Nov 30, 2015, and Sept 23, 2019. In total, 9496 participants (including 790 APOE4 homozygotes, 4869 heterozygotes, and 3837 non-carriers) learned their APOE genotype from a health-care provider as part of Generation Study 1 screening. 4038 (42·5%) participants were in the 65-69-year age group, 5790 (61·0%) were female, 3706 (39·0%) were male, and 8862 (93·3%) self-identified as White. Increase in genetic disease knowledge 2-7 days after disclosure was greater in APOE4 homozygotes (mean 1·19 [SD 3·95]) than in heterozygotes (0·78 [3·95], p=0·042) and non-carriers (0·29 [3·96], p=0·0002). Disease-specific distress 2-7 days after disclosure increased more in homozygotes (2·25 [6·42]) than in heterozygotes (0·53 [5·08], p<0·0001) and non-carriers (0·79 [4·95], p<0·0001). Levels of anxiety 2-7 days after disclosure increased in homozygotes (0·17 [2·95]) but decreased in heterozygotes (-0·67 [2·68], p<0·0001) and non-carriers (-0·66 [2·67], p<0·0001). There were no significant changes in depressive symptoms following disclosure for any APOE4 group. Notably, for all APOE4 groups, increases in distress and anxiety were small and did not reach predefined levels of clinical concern. In cognitively unimpaired, psychologically pre-screened adults, APOE disclosure by a trained health-care provider was generally safe and well tolerated, consistent with results from previous studies. To our knowledge, this is the largest study experience of APOE disclosure to date, especially for homozygotes, and is notable for the older age of participants compared with previous research. These results are timely and important given anticipated increases in APOE disclosure to guide clinical decision making once an Alzheimer's disease prevention treatment is approved for cognitively unimpaired adults or if patients' family members are interested in genetic testing. Scalable approaches for returning Alzheimer's disease risk information are critical to meeting anticipated demand. Results from this study may be useful to bolster clinical translatability of disclosure programmes. The National Institute on Aging, Alzheimer's Association, Banner Alzheimer's Foundation, GHR Foundation, F-Prime Biomedical Research Initiative (FBRI), and Novartis Pharma. Show less

The increased prevalence of opioid use disorder (OUD) makes it imperative to disentangle the biological mechanisms contributing to individual differences in OUD vulnerability. OUD shows strong heritability, however genetic variants contributing to vulnerability remain poorly defined. We performed a genome-wide association study using over 850 male and female heterogeneous stock (HS) rats to identify genes underlying behaviors associated with OUD such as nociception, as well as heroin-taking, extinction and seeking behaviors. By using an animal model of OUD, we were able to identify genetic variants associated with distinct OUD behaviors while maintaining a uniform environment, an experimental design not easily achieved in humans. Furthermore, we used a novel non-linear network-based clustering approach to characterize rats based on OUD vulnerability to assess genetic variants associated with OUD susceptibility. Our findings confirm the heritability of several OUD-like behaviors, including OUD susceptibility. Additionally, several genetic variants associated with nociceptive threshold prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin were identified. Tom1, a microglial component, was implicated for nociception. Several genes involved in dopaminergic signaling, neuroplasticity and substance use disorders, including Brwd1, Pcp4, Phb1l2 and Mmp15 were implicated for the heroin traits. Additionally, an OUD vulnerable phenotype was associated with genetic variants for consumption and break point, suggesting a specific genetic contribution for OUD-like traits contributing to vulnerability. Together, these findings identify novel genetic markers related to the susceptibility to OUD-relevant behaviors in HS rats. Show less

Hypertrophic cardiomyopathy (HCM) is an important cause of morbidity and mortality with both monogenic and polygenic components. Here, we report results from a large genome-wide association study and multitrait analysis including 5,900 HCM cases, 68,359 controls and 36,083 UK Biobank participants with cardiac magnetic resonance imaging. We identified 70 loci (50 novel) associated with HCM and 62 loci (20 novel) associated with relevant left ventricular traits. Among the prioritized genes in the HCM loci, we identify a novel HCM disease gene, SVIL, which encodes the actin-binding protein supervillin, showing that rare truncating SVIL variants confer a roughly tenfold increased risk of HCM. Mendelian randomization analyses support a causal role of increased left ventricular contractility in both obstructive and nonobstructive forms of HCM, suggesting common disease mechanisms and anticipating shared response to therapy. Taken together, these findings increase our understanding of the genetic basis of HCM, with potential implications for disease management. Show less

Atrial fibrillation (AF) is a prevalent and morbid abnormality of the heart rhythm with a strong genetic component. Here, we meta-analyzed genome and exome sequencing data from 36 studies that included 52,416 AF cases and 277,762 controls. In burden tests of rare coding variation, we identified novel associations between AF and the genes MYBPC3, LMNA, PKP2, FAM189A2 and KDM5B. We further identified associations between AF and rare structural variants owing to deletions in CTNNA3 and duplications of GATA4. We broadly replicated our findings in independent samples from MyCode, deCODE and UK Biobank. Finally, we found that CRISPR knockout of KDM5B in stem-cell-derived atrial cardiomyocytes led to a shortening of the action potential duration and widespread transcriptomic dysregulation of genes relevant to atrial homeostasis and conduction. Our results highlight the contribution of rare coding and structural variants to AF, including genetic links between AF and cardiomyopathies, and expand our understanding of the rare variant architecture for this common arrhythmia. Show less

Variants with large effect contribute to congenital heart disease (CHD). To date, recessive genotypes (RGs) have commonly been implicated through anecdotal ascertainment of consanguineous families and candidate gene-based analysis; the recessive contribution to the broad range of CHD phenotypes has been limited. We analyzed whole exome sequences of 5,424 CHD probands. Rare damaging RGs were estimated to contribute to at least 2.2% of CHD, with greater enrichment among laterality phenotypes (5.4%) versus other subsets (1.4%). Among 108 curated human recessive CHD genes, there were 66 RGs, with 54 in 11 genes with >1 RG, 12 genes with 1 RG, and 85 genes with zero. RGs were more prevalent among offspring of consanguineous union (4.7%, 32/675) than among nonconsanguineous probands (0.7%, 34/4749). Founder variants in Show less

To broaden our understanding of bradyarrhythmias and conduction disease, we performed common variant genome-wide association analyses in up to 1.3 million individuals and rare variant burden testing in 460,000 individuals for sinus node dysfunction (SND), distal conduction disease (DCD) and pacemaker (PM) implantation. We identified 13, 31 and 21 common variant loci for SND, DCD and PM, respectively. Four well-known loci (SCN5A/SCN10A, CCDC141, TBX20 and CAMK2D) were shared for SND and DCD, while others were more specific for SND or DCD. SND and DCD showed a moderate genetic correlation (r Show less

The increased prevalence of opioid use disorder (OUD) makes it imperative to disentangle the biological mechanisms contributing to individual differences in OUD vulnerability. OUD shows strong heritability, however genetic variants contributing toward vulnerability remain poorly defined. We performed a genome-wide association study using over 850 male and female heterogeneous stock (HS) rats to identify genes underlying behaviors associated with OUD such as nociception, as well as heroin-taking, extinction and seeking behaviors. By using an animal model of OUD, we were able to identify genetic variants associated with distinct OUD behaviors while maintaining a uniform environment, an experimental design not easily achieved in humans. Furthermore, we used a novel non-linear network-based clustering approach to characterize rats based on OUD vulnerability to assess genetic variants associated with OUD susceptibility. Our findings confirm the heritability of several OUD-like behaviors, including OUD susceptibility. Additionally, several genetic variants associated with nociceptive threshold prior to heroin experience, heroin consumption, escalation of intake, and motivation to obtain heroin were identified. Tom1, a microglial component, was implicated for nociception. Several genes involved in dopaminergic signaling, neuroplasticity and substance use disorders, including Brwd1, Pcp4, Phb1l2 and Mmp15 were implicated for the heroin traits. Additionally, an OUD vulnerable phenotype was associated with genetic variants for consumption and break point, suggesting a specific genetic contribution for OUD-like traits contributing to vulnerability. Together, these findings identify novel genetic markers related to the susceptibility to OUD-relevant behaviors in HS rats. Show less

Diabetes mellitus (DM) increases heart failure incidence and worsens prognosis, but its molecular basis is poorly defined in humans. We aimed to define the diabetic myocardial transcriptome and validate hits in their circulating protein form to define disease mechanisms and biomarkers. RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project was used to define differentially expressed genes (DEGs) in right atrial (RA) and left ventricular (LV) myocardium from people with vs. without DM (type 1 or 2). DEGs were validated as plasma proteins in the UK Biobank cohort, searching for directionally concordant differential expression. Validated plasma proteins were characterized in UK Biobank participants, irrespective of diabetes status, using cardiac magnetic resonance imaging, incident heart failure, and cardiovascular mortality. We found 32 and 32 DEGs associated with DM in the RA and LV, respectively, with no overlap between these. Plasma proteomic data were available for 12, with ERBB3, NRXN3, and HSPA2 (all LV hits) exhibiting directional concordance. Irrespective of DM status, lower circulating ERBB3 and higher HSPA2 were associated with impaired LV contractility and higher LV mass. Participants in the lowest quartile of circulating ERBB3 or highest quartile of circulating HSPA2 had increased incident heart failure and cardiovascular death vs. all other quartiles. DM is characterized by lower Erbb3 and higher Hspa2 expression in the myocardium, with directionally concordant differences in their plasma protein concentration. These are associated with LV dysfunction, incident heart failure, and cardiovascular mortality. Show less

Transitioning from a genetic association signal to an effector gene and a targetable molecular mechanism requires the application of functional fine-mapping tools such as reporter assays and genome editing. In this report, we undertook such studies on the osteoarthritis (OA) risk that is marked by single nucleotide polymorphism (SNP) rs34195470 (A > G). The OA risk-conferring G allele of this SNP associates with increased DNA methylation (DNAm) at two CpG dinucleotides within WWP2. This gene encodes a ubiquitin ligase and is the host gene of microRNA-140 (miR-140). WWP2 and miR-140 are both regulators of TGFβ signaling. Nucleic acids were extracted from adult OA (arthroplasty) and foetal cartilage. Samples were genotyped and DNAm quantified by pyrosequencing at the two CpGs plus 14 flanking CpGs. CpGs were tested for transcriptional regulatory effects using a chondrocyte cell line and reporter gene assay. DNAm was altered using epigenetic editing, with the impact on gene expression determined using RT-qPCR. In silico analysis complemented laboratory experiments. rs34195470 genotype associates with differential methylation at 14 of the 16 CpGs in OA cartilage, forming a methylation quantitative trait locus (mQTL). The mQTL is less pronounced in foetal cartilage (5/16 CpGs). The reporter assay revealed that the CpGs reside within a transcriptional regulator. Epigenetic editing to increase their DNAm resulted in altered expression of the full-length and N-terminal transcript isoforms of WWP2. No changes in expression were observed for the C-terminal isoform of WWP2 or for miR-140. As far as we are aware, this is the first experimental demonstration of an OA association signal targeting specific transcript isoforms of a gene. The WWP2 isoforms encode proteins with varying substrate specificities for the components of the TGFβ signaling pathway. Future analysis should focus on the substrates regulated by the two WWP2 isoforms that are the targets of this genetic risk. Show less

Newborn genome sequencing (NBSeq) can detect infants at risk for treatable disorders currently undetected by conventional newborn screening. Despite broad stakeholder support for NBSeq, the perspectives of rare disease experts regarding which diseases should be screened have not been ascertained. To query rare disease experts about their perspectives on NBSeq and which gene-disease pairs they consider appropriate to evaluate in apparently healthy newborns. This survey study, designed between November 2, 2021, and February 11, 2022, assessed experts' perspectives on 6 statements related to NBSeq. Experts were also asked to indicate whether they would recommend including each of 649 gene-disease pairs associated with potentially treatable conditions in NBSeq. The survey was administered between February 11 and September 23, 2022, to 386 experts, including all 144 directors of accredited medical and laboratory genetics training programs in the US. Expert perspectives on newborn screening using genome sequencing. The proportion of experts indicating agreement or disagreement with each survey statement and those who selected inclusion of each gene-disease pair were tabulated. Exploratory analyses of responses by gender and age were conducted using t and χ2 tests. Of 386 experts invited, 238 (61.7%) responded (mean [SD] age, 52.6 [12.8] years [range 27-93 years]; 126 [52.9%] women and 112 [47.1%] men). Among the experts who responded, 161 (87.9%) agreed that NBSeq for monogenic treatable disorders should be made available to all newborns; 107 (58.5%) agreed that NBSeq should include genes associated with treatable disorders, even if those conditions were low penetrance; 68 (37.2%) agreed that actionable adult-onset conditions should be sequenced in newborns to facilitate cascade testing in parents, and 51 (27.9%) agreed that NBSeq should include screening for conditions with no established therapies or management guidelines. The following 25 genes were recommended by 85% or more of the experts: OTC, G6PC, SLC37A4, CYP11B1, ARSB, F8, F9, SLC2A1, CYP17A1, RB1, IDS, GUSB, DMD, GLUD1, CYP11A1, GALNS, CPS1, PLPBP, ALDH7A1, SLC26A3, SLC25A15, SMPD1, GATM, SLC7A7, and NAGS. Including these, 42 gene-disease pairs were endorsed by at least 80% of experts, and 432 genes were endorsed by at least 50% of experts. In this survey study, rare disease experts broadly supported NBSeq for treatable conditions and demonstrated substantial concordance regarding the inclusion of a specific subset of genes in NBSeq. Show less

Genetic screening for youth with obesity in the absence of syndromic findings has not been part of obesity management. For children with early onset obesity, genetic screening is recommended for those having clinical features of genetic obesity syndromes (including hyperphagia). The overarching goal of this work is to report the findings and experiences from one pediatric weight management program that implemented targeted sequencing analysis for genes known to cause rare genetic disorders of obesity. This exploratory study evaluated youth tested over an 18-month period using a panel of 40-genes in the melanocortin 4 receptor pathway. Medical records were reviewed for demographic and visit information, including body mass index (BMI) percent of 95th percentile (%BMIp95) and two eating behaviors. Of 117 subjects: 51.3% were male; 53.8% Hispanic; mean age 10.2 years (SD 3.8); mean %BMIp95 157% (SD 29%). Most subjects were self- or caregiver-reported to have overeating to excess or binge eating (80.3%) and sneaking food or eating in secret (59.0%). Among analyzed genes, 72 subjects (61.5%) had at least one variant reported; 50 (42.7%) had a single variant reported; 22 (18.8%) had 2-4 variants reported; most variants were rare (<0.05% minor allele frequency [MAF]), and of uncertain significance; all variants were heterozygous. Nine subjects (7.7%) had a variant reported as PSCK1 "risk" or MC4R "likely pathogenic"; 39 (33.3%) had a Bardet-Biedl Syndrome (BBS) gene variant (4 with "pathogenic" or "likely pathogenic" variants). Therefore, 9 youth (7.7%) had gene variants previously identified as increasing risk for obesity and 4 youth (3.4%) had BBS carrier status. Panel testing identified rare variants of uncertain significance in most youth tested, and infrequently identified variants previously reported to increase the risk for obesity. Further research in larger cohorts is needed to understand how genetic variants influence the expression of non-syndromic obesity. Show less

Lipid transfer proteins mediate the exchange of lipids between closely apposed membranes at organelle contact sites and play key roles in lipid metabolism, membrane homeostasis, and cellular signaling. A recently discovered novel family of lipid transfer proteins, which includes the VPS13 proteins (VPS13A-D), adopt a rod-like bridge conformation with an extended hydrophobic groove that enables the bulk transfer of membrane lipids for membrane growth. Loss of function mutations in VPS13A and VPS13C cause chorea acanthocytosis and Parkinson's disease, respectively. VPS13A and VPS13C localize to multiple organelle contact sites, including endoplasmic reticulum (ER) - lipid droplet (LD) contact sites, but the functional roles of these proteins in LD regulation remains mostly unexplored. Here we employ CRISPR-Cas9 genome editing to generate VPS13A and VPS13C knockout cell lines in U-2 OS cells via deletion of exon 2 and introduction of an early frameshift. Analysis of LD content in these cell lines revealed that loss of either VPS13A or VPS13C results in reduced LD abundance under oleate-stimulated conditions. These data implicate two lipid transfer proteins, VPS13A and VPS13C, in LD regulation. Show less

The heart muscle diseases hypertrophic (HCM) and dilated (DCM) cardiomyopathies are leading causes of sudden death and heart failure in young, otherwise healthy, individuals. We conducted genome-wide association studies and multi-trait analyses in HCM (1,733 cases), DCM (5,521 cases) and nine left ventricular (LV) traits (19,260 UK Biobank participants with structurally normal hearts). We identified 16 loci associated with HCM, 13 with DCM and 23 with LV traits. We show strong genetic correlations between LV traits and cardiomyopathies, with opposing effects in HCM and DCM. Two-sample Mendelian randomization supports a causal association linking increased LV contractility with HCM risk. A polygenic risk score explains a significant portion of phenotypic variability in carriers of HCM-causing rare variants. Our findings thus provide evidence that polygenic risk score may account for variability in Mendelian diseases. More broadly, we provide insights into how genetic pathways may lead to distinct disorders through opposing genetic effects. Show less

Circulating metabolite levels are biomarkers for cardiovascular disease (CVD). Here we studied, association of rare variants and 226 serum lipoproteins, lipids and amino acids in 7,142 (discovery plus follow-up) healthy participants. We leveraged the information from multiple metabolite measurements on the same participants to improve discovery in rare variant association analyses for gene-based and gene-set tests by incorporating correlated metabolites as covariates in the validation stage. Gene-based analysis corrected for the effective number of tests performed, confirmed established associations at APOB, APOC3, PAH, HAL and PCSK (p<1.32x10-7) and identified novel gene-trait associations at a lower stringency threshold with ACSL1, MYCN, FBXO36 and B4GALNT3 (p<2.5x10-6). Regulation of the pyruvate dehydrogenase (PDH) complex was associated for the first time, in gene-set analyses also corrected for effective number of tests, with IDL and LDL parameters, as well as circulating cholesterol (pMETASKAT<2.41x10-6). In conclusion, using an approach that leverages metabolite measurements obtained in the same participants, we identified novel loci and pathways involved in the regulation of these important metabolic biomarkers. As large-scale biobanks continue to amass sequencing and phenotypic information, analytical approaches such as ours will be useful to fully exploit the copious amounts of biological data generated in these efforts. Show less

In humans, abnormal thickening of the left ventricle of the heart clinically defines hypertrophic cardiomyopathy (HCM), a common inherited cardiovascular disorder that can precede a sudden cardiac death event. The wide range of clinical presentations in HCM obscures genetic variants that may influence an individual's susceptibility to sudden cardiac death. Although exon sequencing of major sarcomere genes can be used to detect high-impact causal mutations, this strategy is successful in only half of patient cases. The incidence of left ventricular hypertrophy (LVH) in a managed research colony of rhesus macaques provides an excellent comparative model in which to explore the genomic etiology of severe HCM and sudden cardiac death. Because no rhesus HCM-associated mutations have been reported, we used a next-generation genotyping assay that targets 7 sarcomeric rhesus genes within 63 genomic sites that are orthologous to human genomic regions known to harbor HCM disease variants. Amplicon sequencing was performed on 52 macaques with confirmed LVH and 42 unrelated, unaffected animals representing both the Indian and Chinese rhesus macaque subspecies. Bias-reduced logistic regression uncovered a risk haplotype in the rhesus Show less

Whereas the metabolic consequences of obesity have been studied extensively in the rhesus macaque, corollary genetic studies of obesity are nonexistent. This study assessed genetic contributions to spontaneous adiposity in this species. Phenotypic variation by age class and sex for BMI, waist to height ratio, waist to thigh ratio, and waist circumference was assessed in 583 macaques. Total and sex-specific heritability for all traits was estimated, including waist to thigh ratio adjusted for BMI, as well as genotypic and phenotypic correlations. In addition, functional genetic variation at BDNF, FTO, LEP, LEPR, MC4R, PCSK1, POMC, and SIM1 was assessed in four animals with extreme spontaneous adiposity. Trait heritability in the combined sample was low to moderate (0.14-0.32), whereas sex-specific heritability was more substantial (0.20-0.67). Heritability was greater in females for all traits except BMI. All traits were robustly correlated, with genetic correlations of 0.63 to 0.93 indicating substantial pleiotropy. Likely functional variants were discovered in the four macaques at all eight human obesity genes, including six missense mutations in BDNF, FTO, LEP, LEPR, and PCSK1 and, notably, one nonsense mutation in LEPR. A moderate polygenic contribution to adiposity in rhesus macaques was found, as well as mutations with potentially larger effects in multiple genes that influence obesity in humans. Show less

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity. Show less

The omnigenic model of complex disease stipulates that the majority of the heritability will be explained by the effects of common variation on genes in the periphery of core disease pathways. Rare variant associations, expected to explain far less of the heritability, may be enriched in core disease genes and thus will be instrumental in the understanding of complex disease pathogenesis and their potential therapeutic targets. Here, using complementary whole-exome sequencing, high-density imputation, and in vitro cellular assays, we identify candidate core genes in the pathogenesis of systemic lupus erythematosus (SLE). Using extreme-phenotype sampling, we sequenced the exomes of 30 SLE parent-affected-offspring trios and identified 14 genes with missense de novo mutations (DNM), none of which are within the >80 SLE susceptibility loci implicated through genome-wide association studies. In a follow-up cohort of 10, 995 individuals of matched European ancestry, we imputed genotype data to the density of the combined UK10K-1000 genomes Phase III reference panel across the 14 candidate genes. Gene-level analyses indicate three functional candidates: DNMT3A, PRKCD, and C1QTNF4. We identify a burden of rare variants across PRKCD associated with SLE risk (P = 0.0028), and across DNMT3A associated with two severe disease prognosis sub-phenotypes (P = 0.0005 and P = 0.0033). We further characterise the TNF-dependent functions of the third candidate gene C1QTNF4 on NF-κB activation and apoptosis, which are inhibited by the p.His198Gln DNM. Our results identify three novel genes in SLE susceptibility and support extreme-phenotype sampling and DNM gene discovery to aid the search for core disease genes implicated through rare variation. Show less

Metabolites derived from dietary choline and L-carnitine, such as trimethylamine N-oxide and betaine, have recently been identified as novel risk factors for atherosclerosis in mice and humans. We sought to identify genetic factors associated with plasma betaine levels and determine their effect on risk of coronary artery disease (CAD). A two-stage genome-wide association study (GWAS) identified two significantly associated loci on chromosomes 2q34 and 5q14.1. The lead variant on 2q24 (rs715) localizes to carbamoyl-phosphate synthase 1 (CPS1), which encodes a mitochondrial enzyme that catalyses the first committed reaction and rate-limiting step in the urea cycle. Rs715 is also significantly associated with decreased levels of urea cycle metabolites and increased plasma glycine levels. Notably, rs715 yield a strikingly significant and protective association with decreased risk of CAD in only women. These results suggest that glycine metabolism and/or the urea cycle represent potentially novel sex-specific mechanisms for the development of atherosclerosis. Show less

The liver X receptors (LXRs, NR1H2 and NR1H3) and peroxisome proliferator-activated receptor gamma (PPARG, NR1C3) nuclear receptor transcription factors (TFs) are master regulators of energy homeostasis. Intriguingly, recent studies suggest that these metabolic regulators also impact tumor cell proliferation. However, a comprehensive temporal molecular characterization of the LXR and PPARG gene regulatory responses in tumor cells is still lacking. To better define the underlying molecular processes governing the genetic control of cellular growth in response to extracellular metabolic signals, we performed a comprehensive, genome-wide characterization of the temporal regulatory cascades mediated by LXR and PPARG signaling in HT29 colorectal cancer cells. For this analysis, we applied a multi-tiered approach that incorporated cellular phenotypic assays, gene expression profiles, chromatin state dynamics, and nuclear receptor binding patterns. Our results illustrate that the activation of both nuclear receptors inhibited cell proliferation and further decreased glutathione levels, consistent with increased cellular oxidative stress. Despite a common metabolic reprogramming, the gene regulatory network programs initiated by these nuclear receptors were widely distinct. PPARG generated a rapid and short-term response while maintaining a gene activator role. By contrast, LXR signaling was prolonged, with initial, predominantly activating functions that transitioned to repressive gene regulatory activities at late time points. Through the use of a multi-tiered strategy that integrated various genomic datasets, our data illustrate that distinct gene regulatory programs elicit common phenotypic effects, highlighting the complexity of the genome. These results further provide a detailed molecular map of metabolic reprogramming in cancer cells through LXR and PPARG activation. As ligand-inducible TFs, these nuclear receptors can potentially serve as attractive therapeutic targets for the treatment of various cancers. Show less