Epigenetic associations in HPA axis genes related to bronchopulmonary dysplasia and antenatal steroids.

BACKGROUND: Bronchopulmonary dysplasia (BPD), a common morbidity among very preterm infants, is associated with chronic disease and neurodevelopmental impairments. A hypothesized mechanism for these outcomes lies in altered glucocorticoid (GC) activity. We hypothesized that BPD and its treatments may result in epigenetic differences in the hypothalamic-pituitary-adrenal (HPA) axis, which is modulated by GC, and could be ascertained using an established GC risk score and DNA methylation (DNAm) of HPA axis genes. METHODS: DNAm was quantified from buccal tissue (ECHO-NOVI) and from neonatal blood spots (ELGAN ECHO) via the EPIC microarray. Prenatal maternal characteristics, pregnancy complication, and neonatal medical complication data were collected from medical record review and maternal interviews. RESULTS: The GC score was not associated with steroid exposure or BPD. However, six HPA genes involved in stress response regulation demonstrated differential methylation with antenatal steroid exposure; two CpGs within FKBP5 and POMC were differentially methylated with BPD severity. These findings were sex-specific in both cohorts; males had greater magnitude of differential methylation within these genes. CONCLUSIONS: These findings suggest that BPD severity and antenatal steroids are associated with DNAm at some HPA genes in very preterm infants and the effects appear to be sex-, tissue-, and age-specific. IMPACT: This study addresses bronchopulmonary dysplasia (BPD), an important health outcome among preterm neonates, and interrogates a commonly studied pathway, the hypothalamic-pituitary-adrenal (HPA) axis. The combination of BPD, the HPA axis, and epigenetic markers has not been previously reported. In this study, we found that BPD itself was not associated with epigenetic responses in the HPA axis in infants born very preterm; however, antenatal treatment with steroids was associated with epigenetic responses.

Seasonally variant gene expression in full-term human placenta.

Seasonal exposures influence human health and development. The placenta, as a mediator of the maternal and fetal systems and a regulator of development, is an ideal tissue to understand the biological pathways underlying relationships between season of birth and later life health outcomes. Here, we conducted a differential expression (DE) analysis of season of birth in full-term human placental tissue to evaluate whether the placenta may be influenced by seasonal cues. Of the analyzed transcripts, 583 displayed DE between summer and winter births (False Discovery Rate [FDR] q < .05); among these, BHLHE40, MIR210HG, and HILPDA had increased expression among winter births (Bonferroni P < .05). Enrichment analyses of the seasonally variant genes between summer and winter births indicated overrepresentation of transcription factors HIF1A, VDR, and CLOCK, among others, and of GO term pathways related to ribosomal activity and infection. Additionally, a cosinor analysis found rhythmic expression for approximately 11.9% of all 17 664 analyzed placental transcripts. These results suggest that the placenta responds to seasonal cues and add to the growing body of evidence that the placenta acts as a peripheral clock, which may provide a molecular explanation for the extensive associations between season of birth and health outcomes.

Epigenome-wide analysis identifies genes and pathways linked to acoustic cry variation in preterm infants.

BACKGROUND: Preterm birth places infants at higher risk of adverse long-term behavioral and cognitive outcomes. Combining biobehavioral measures and molecular biomarkers may improve tools to predict the risk of long-term developmental delays. METHODS: The Neonatal Neurobehavior and Outcomes in Very Preterm Infants study was conducted at nine neonatal intensive care units between April 2014 and June 2016. Cries were recorded and buccal swabs collected during the neurobehavioral exam. Cry episodes were extracted and analyzed using a computer system and the data were summarized using factor analysis. Genomic DNA was extracted from buccal swabs, quantified using the Qubit Fluorometer, and aliquoted into standardized concentrations. DNA methylation was measured with the Illumina MethylationEPIC BeadArray, and an epigenome-wide association study was performed using cry factors (n = 335). RESULTS: Eighteen CpGs were associated with the cry factors at genome-wide significance (α = 7.08E - 09). Two CpG sites, one intergenic and one linked to gene TCF3 (important for B and T lymphocyte development), were associated with acoustic measures of cry energy. Increased methylation of TCF3 was associated with a lower energy-related cry factor. We also found that pitch (F0) and hyperpitch (F0 > 1 kHz) were associated with DNA methylation variability at 16 CpG sites. CONCLUSIONS: Acoustic cry characteristics are related to variation in DNA methylation in preterm infants. IMPACT: Preterm birth is a major public health problem and its long-term impact on health is not well understood. Cry acoustics, related to prematurity, has been linked to a variety of medical conditions. Biobehavioral measures and molecular biomarkers can improve prediction tools for long-term developmental risks of preterm birth. Variation in epigenetic modulation in preterm infants provides a potential link between preterm birth and unfavorable developmental outcomes.

A vascular endothelial growth factor A genetic variant is associated with improved ventricular function and transplant-free survival after surgery for non-syndromic CHD.

BACKGROUND: We have previously shown that the minor alleles of vascular endothelial growth factor A (VEGFA) single-nucleotide polymorphism rs833069 and superoxide dismutase 2 (SOD2) single-nucleotide polymorphism rs2758331 are both associated with improved transplant-free survival after surgery for CHD in infants, but the underlying mechanisms are unknown. We hypothesised that one or both of these minor alleles are associated with better systemic ventricular function, resulting in improved survival. METHODS: This study is a follow-up analysis of 422 non-syndromic CHD patients who underwent neonatal cardiac surgery with cardiopulmonary bypass. Echocardiographic reports were reviewed. Systemic ventricular function was subjectively categorised as normal, or as mildly, moderately, or severely depressed. The change in function was calculated as the change from the preoperative study to the last available study. Stepwise linear regression, adjusting for covariates, was performed for the outcome of change in ventricular function. Model comparison was performed using Akaike's information criterion. Only variables that improved the model prediction of change in systemic ventricular function were retained in the final model. RESULTS: Genetic and echocardiographic data were available for 335/422 subjects (79%). Of them, 33 (9.9%) developed worse systemic ventricular function during a mean follow-up period of 13.5 years. After covariate adjustment, the presence of the VEGFA minor allele was associated with preserved ventricular function (p=0.011). CONCLUSIONS: These data support the hypothesis that the mechanism by which the VEGFA single-nucleotide polymorphism rs833069 minor allele improves survival may be the preservation of ventricular function. Further studies are needed to validate this genotype-phenotype association and to determine whether this mechanism is related to increased vascular endothelial growth factor production.

Association Between Absolute Neutrophil Count and Variation at TCIRG1: The NHLBI Exome Sequencing Project.

Neutrophils are a key component of innate immunity. Individuals with low neutrophil count are susceptible to frequent infections. Linkage and association between congenital neutropenia and a single rare missense variant in TCIRG1 have been reported in a single family. Here, we report on nine rare missense variants at evolutionarily conserved sites in TCIRG1 that are associated with lower absolute neutrophil count (ANC; p = 0.005) in 1,058 participants from three cohorts: Atherosclerosis Risk in Communities (ARIC), Coronary Artery Risk Development in Young Adults (CARDIA), and Jackson Heart Study (JHS) of the NHLBI Grand Opportunity Exome Sequencing Project (GO ESP). These results validate the effects of TCIRG1 coding variation on ANC and suggest that this gene may be associated with a spectrum of mild to severe effects on ANC.

Sequencing of sporadic Attention-Deficit Hyperactivity Disorder (ADHD) identifies novel and potentially pathogenic de novo variants and excludes overlap with genes associated with autism spectrum disorder.

Attention-Deficit Hyperactivity Disorder (ADHD) has high heritability; however, studies of common variation account for <5% of ADHD variance. Using data from affected participants without a family history of ADHD, we sought to identify de novo variants that could account for sporadic ADHD. Considering a total of 128 families, two analyses were conducted in parallel: first, in 11 unaffected parent/affected proband trios (or quads with the addition of an unaffected sibling) we completed exome sequencing. Six de novo missense variants at highly conserved bases were identified and validated from four of the 11 families: the brain-expressed genes TBC1D9, DAGLA, QARS, CSMD2, TRPM2, and WDR83. Separately, in 117 unrelated probands with sporadic ADHD, we sequenced a panel of 26 genes implicated in intellectual disability (ID) and autism spectrum disorder (ASD) to evaluate whether variation in ASD/ID-associated genes were also present in participants with ADHD. Only one putative deleterious variant (Gln600STOP) in CHD1L was identified; this was found in a single proband. Notably, no other nonsense, splice, frameshift, or highly conserved missense variants in the 26 gene panel were identified and validated. These data suggest that de novo variant analysis in families with independently adjudicated sporadic ADHD diagnosis can identify novel genes implicated in ADHD pathogenesis. Moreover, that only one of the 128 cases (0.8%, 11 exome, and 117 MIP sequenced participants) had putative deleterious variants within our data in 26 genes related to ID and ASD suggests significant independence in the genetic pathogenesis of ADHD as compared to ASD and ID phenotypes. © 2017 Wiley Periodicals, Inc.

PLTP activity inversely correlates with CAAD: effects of PON1 enzyme activity and genetic variants on PLTP activity.

Recent studies have failed to demonstrate a causal cardioprotective effect of HDL cholesterol levels, shifting focus to the functional aspects of HDL. Phospholipid transfer protein (PLTP) is an HDL-associated protein involved in reverse cholesterol transport. This study sought to determine the genetic and nongenetic predictors of plasma PLTP activity (PLTPa), and separately, to determine whether PLTPa predicted carotid artery disease (CAAD). PLTPa was measured in 1,115 European ancestry participants from a case-control study of CAAD. A multivariate logistic regression model was used to elucidate the relationship between PLTPa and CAAD. Separately, a stepwise linear regression determined the nongenetic clinical and laboratory characteristics that best predicted PLTPa. A final stepwise regression considering both nongenetic and genetic variables identified the combination of covariates that explained maximal PLTPa variance. PLTPa was significantly associated with CAAD (7.90 × 10(-9)), with a 9% decrease in odds of CAAD per 1 unit increase in PLTPa (odds ratio = 0.91). Triglyceride levels (P = 0.0042), diabetes (P = 7.28 × 10(-5)), paraoxonase 1 (PON1) activity (P = 0.019), statin use (P = 0.026), PLTP SNP rs4810479 (P = 6.38 × 10(-7)), and PCIF1 SNP rs181914932 (P = 0.041) were all significantly associated with PLTPa. PLTPa is significantly inversely correlated with CAAD. Furthermore, we report a novel association between PLTPa and PON1 activity, a known predictor of CAAD.

Rare coding variation in paraoxonase-1 is associated with ischemic stroke in the NHLBI Exome Sequencing Project.

HDL-associated paraoxonase-1 (PON1) is an enzyme whose activity is associated with cerebrovascular disease. Common PON1 genetic variants have not been consistently associated with cerebrovascular disease. Rare coding variation that likely alters PON1 enzyme function may be more strongly associated with stroke. The National Heart, Lung, and Blood Institute Exome Sequencing Project sequenced the coding regions (exomes) of the genome for heart, lung, and blood-related phenotypes (including ischemic stroke). In this sample of 4,204 unrelated participants, 496 had verified, noncardioembolic ischemic stroke. After filtering, 28 nonsynonymous PON1 variants were identified. Analysis with the sequence kernel association test, adjusted for covariates, identified significant associations between PON1 variants and ischemic stroke (P = 3.01 × 10(-3)). Stratified analyses demonstrated a stronger association of PON1 variants with ischemic stroke in African ancestry (AA) participants (P = 5.03 × 10(-3)). Ethnic differences in the association between PON1 variants with stroke could be due to the effects of PON1Val109Ile (overall P = 7.88 × 10(-3); AA P = 6.52 × 10(-4)), found at higher frequency in AA participants (1.16% vs. 0.02%) and whose protein is less stable than the common allele. In summary, rare genetic variation in PON1 was associated with ischemic stroke, with stronger associations identified in those of AA. Increased focus on PON1 enzyme function and its role in cerebrovascular disease is warranted.

Epigenetic landscape of 5-hydroxymethylcytosine and associations with gene expression in placenta.

5-hydroxymethylcystosine (5hmC), is an intermediate product in the DNA demethylation pathway, but may act as a functional epigenetic modification. We have conducted the largest study of site-specific 5hmC in placenta to date using parallel bisulphite and oxidative bisulphite modification with array-based assessment. Incorporating parallel RNA-sequencing data allowed us to assess associations between 5hmC and gene expression, using expression quantitative trait hydroxymethylation (eQTHM) analysis. We identified ~ 47,000 loci with consistently elevated (systematic) 5hmC proportions. Systematic 5hmC was significantly depleted (p < 0.0001) at CpG islands (CGI), and enriched (p < 0.0001) in 'open sea' regions (CpG >4 kb from CGI). 5hmC was most and least abundant at CpGs in enhancers and active transcription start sites (TSS), respectively (p < 0.05). We identified 499 significant (empirical-p <0.05) eQTHMs within 1 MB of the assayed gene. At most (75.4%) eQTHMs, the proportion of 5hmC was positively correlated with transcript abundance. eQTHMs were significantly enriched among enhancer CpGs and depleted among CpGs in active TSS (p < 0.05 for both). Finally, we identified 107 differentially hydroxymethylated regions (DHMRs, p < 0.05) across 100 genes. Our study provides insight into placental distribution of 5hmC, and sheds light on the functional capacity of this epigenetic modification in placenta.

Novel gene-by-environment interactions: APOB and NPC1L1 variants affect the relationship between dietary and total plasma cholesterol.

Cardiovascular disease (CVD) is the leading cause of death in developed countries. Plasma cholesterol level is a key risk factor in CVD pathogenesis. Genetic and dietary variation both influence plasma cholesterol; however, little is known about dietary interactions with genetic variants influencing the absorption and transport of dietary cholesterol. We sought to determine whether gut expressed variants predicting plasma cholesterol differentially affected the relationship between dietary and plasma cholesterol levels in 1,128 subjects (772/356 in the discovery/replication cohorts, respectively). Four single nucleotide polymorphisms (SNPs) within three genes (APOB, CETP, and NPC1L1) were significantly associated with plasma cholesterol in the discovery cohort. These were subsequently evaluated for gene-by-environment (GxE) interactions with dietary cholesterol for the prediction of plasma cholesterol, with significant findings tested for replication. Novel GxE interactions were identified and replicated for two variants: rs1042034, an APOB Ser4338Asn missense SNP and rs2072183 (in males only), a synonymous NPC1L1 SNP in linkage disequilibrium with SNPs 5' of NPC1L1. This study identifies the presence of novel GxE and gender interactions implying that differential gut absorption is the basis for the variant associations with plasma cholesterol. These GxE interactions may account for part of the "missing heritability" not accounted for by genetic associations.

Novel common and rare genetic determinants of paraoxonase activity: FTO, SERPINA12, and ITGAL.

HDL-associated paraoxonase 1 (PON1) activity is associated with cardiovascular and other human diseases. As the role of genetic variants outside of the PON gene cluster on PON1 activity is unknown, we sought to identify common and rare variants in such loci. We typed 33,057 variants on the CVD chip in 1,362 subjects to test for their effects on adjusted-PON1 activity. Three novel genes (FTO, ITGAL, and SERPINA12) and the PON gene cluster had SNPs associated with PON1 arylesterase (AREase) activity. These loci were carried forward for rare-variant analysis using Exome chip genotypes in an overlapping subset of 1,051 subjects using sequence kernel association testing. PON1 (P = 2.24 × 10(-4)), PON3 (P = 0.022), FTO (P = 0.019), and SERPINA12 (P = 0.039) had both common and rare variants associated with PON1 AREase. ITGAL variants were associated with PON1 activity when using weighted sequence kernel association testing (SKAT) analysis (P = 2.63 × 10(-3)). When adjusting for the initial common variants, SERPINA12 became marginally significant (P = 0.09), whereas all other findings remained significant (P < 0.05), suggesting independent rare-variant effects. We present novel findings that common and rare variants in FTO, SERPINA12, and ITGAL predict PON1 activity. These results further link PON1 to diabetes and inflammation and may inform the role of HDL in human disease.

Dietary fatty acid intake is associated with paraoxonase 1 activity in a cohort-based analysis of 1,548 subjects.

BACKGROUND: Paraoxonase 1 (PON1) is a cardioprotective, HDL-associated glycoprotein enzyme with broad substrate specificity. Our previous work found associations between dietary cholesterol and vitamin C with PON1 activity. The goal of this study was to determine the effect of specific dietary fatty acid (DFA) intake on PON1 activity. METHODS: 1,548 participants with paraoxonase activity measures completed the Harvard Standardized Food Frequency Questionnaire to determine their daily nutrient intake over the past year. Eight saturated, 3 monounsaturated, and 6 polyunsaturated DFAs were measured by the questionnaire. To reduce the number of observations tested, only specific fatty acids that were not highly correlated (r < 0.8) with other DFAs or that were representative of other DFAs through high correlation within each respective group (saturated, monounsaturated, or polyunsaturated) were retained for analysis. Six specific DFA intakes - myristic acid (14 carbon atoms, no double bonds - 14:0), oleic acid (18:1), gadoleic acid (20:1), α-linolenic acid (18:3), arachidonic acid (20:4), and eicosapentaenoic acid (20:5) - were carried forward to stepwise linear regression, which evaluated the effect of each specific DFA on covariate-adjusted PON1 enzyme activity. RESULTS: Four of the 6 tested DFA intakes - myristic acid (p = 0.038), gadoleic acid (p = 6.68 × 10(-7)), arachidonic acid (p = 0.0007), and eicosapentaenoic acid (p = 0.013) - were independently associated with covariate-adjusted PON1 enzyme activity. Myristic acid, a saturated fat, and gadoleic acid, a monounsaturated fat, were both positively associated with PON1 activity. Both of the tested polyunsaturated fats, arachidonic acid and eicosapentaenoic acid, were negatively associated with PON1 activity. CONCLUSIONS: This study presents the largest cohort-based analysis of the relationship between dietary lipids and PON1 enzyme activity. Further research is necessary to elucidate and understand the specific biological mechanisms, whether direct or regulatory, through which DFAs affect PON1 activity.

A multi-omic approach identifies an autism spectrum disorder (ASD) regulatory complex of functional epimutations in placentas from children born preterm.

Children born preterm are at heightened risk of neurodevelopmental impairments, including Autism Spectrum Disorder (ASD). The placenta is a key regulator of neurodevelopmental processes, though the precise underlying molecular mechanisms remain unclear. Here, we employed a multi-omic approach to identify placental transcriptomic and epigenetic modifications related to ASD diagnosis at age 10, among children born preterm. Working with the extremely low gestational age (ELGAN) cohort, we hypothesized that a pro-inflammatory placental environment would be predictive of ASD diagnosis at age 10. Placental messenger RNA (mRNA) expression, CpG methylation, and microRNA (miRNA) expression were compared among 368 ELGANs (28 children diagnosed with ASD and 340 children without ASD). A total of 111 genes displayed expression levels in the placenta that were associated with ASD. Within these ASD-associated genes is an ASD regulatory complex comprising key genes that predicted ASD case status. Genes with expression that predicted ASD case status included Ewing Sarcoma Breakpoint Region 1 (EWSR1) (OR: 6.57 (95% CI: 2.34, 23.58)) and Bromodomain Adjacent To Zinc Finger Domain 2A (BAZ2A) (OR: 0.12 (95% CI: 0.03, 0.35)). Moreover, of the 111 ASD-associated genes, nine (8.1%) displayed associations with CpG methylation levels, while 14 (12.6%) displayed associations with miRNA expression levels. Among these, LRR Binding FLII Interacting Protein 1 (LRRFIP1) was identified as being under the control of both CpG methylation and miRNAs, displaying an OR of 0.42 (95% CI: 0.17, 0.95). This gene, as well as others identified as having functional epimutations, plays a critical role in immune system regulation and inflammatory response. In summary, a multi-omic approach was used to identify functional epimutations in the placenta that are associated with the development of ASD in children born preterm, highlighting future avenues for intervention.

Sex-based differences in placental DNA methylation profiles related to gestational age: an NIH ECHO meta-analysis.

The placenta undergoes many changes throughout gestation to support the evolving needs of the foetus. There is also a growing appreciation that male and female foetuses develop differently in utero, with unique epigenetic changes in placental tissue. Here, we report meta-analysed sex-specific associations between gestational age and placental DNA methylation from four cohorts in the National Institutes of Health (NIH) Environmental influences on Child Health Outcomes (ECHO) Programme (355 females/419 males, gestational ages 23-42 weeks). We identified 407 cytosine-guanine dinucleotides (CpGs) in females and 794 in males where placental methylation levels were associated with gestational age. After cell-type adjustment, 55 CpGs in females and 826 in males were significant. These were enriched for biological processes critical to the immune system in females and transmembrane transport in males. Our findings are distinct between the sexes: in females, associations with gestational age are largely explained by differences in placental cellular composition, whereas in males, gestational age is directly associated with numerous alterations in methylation levels.

Dietary cholesterol increases paraoxonase 1 enzyme activity.

HDL-associated paraoxonase 1 (PON1) activity has been consistently associated with cardiovascular and other diseases. Vitamins C and E intake have previously been positively associated with PON1 in a subset of the Carotid Lesion Epidemiology and Risk (CLEAR) cohort. The goal of this study was to replicate these findings and determine whether other nutrient intake affected PON1 activity. To predict nutrient and mineral intake values, 1,402 subjects completed a standardized food frequency survey of their dietary habits over the past year. Stepwise regression was used to evaluate dietary and covariate effects on PON1 arylesterase activity. Five dietary components, cholesterol (P < 2.0 × 10(-16)), alcohol (P = 8.51 × 10(-8)), vitamin C (P = 7.97 × 10(-5)), iron (P = 0.0026), and folic acid (0.037) were independently predictive of PON1 activity. Dietary cholesterol was positively associated and predicted 5.5% of PON1 activity, second in variance explained. This study presents a novel finding of dietary cholesterol, iron, and folic acid predicting PON1 activity in humans and confirms prior reported associations, including that with vitamin C. Identifying and understanding environmental factors that affect PON1 activity is necessary to understand its role and that of HDL in human disease.

Developmental chronodisruption alters placental signaling in mice.

Chronodisruption has been largely overlooked as a developmental exposure. The placenta, a conduit between the maternal and fetal environments, may relay circadian cues to the fetus. We have previously shown that developmental chronodisruption causes visual impairment and increased retinal microglial and macrophage marker expression. Here, we investigated the impacts of environmental chronodisruption on fetal and placental outcomes in a C57BL/6J mouse (Mus musculus) model. Developmental chronodisruption had no effect on embryo count, placental weight, or fetal sex ratio. When measured with RNAseq, mice exposed to developmental chronodisruption (CD) had differential placental expression of several transcripts including Serpinf1, which encodes pigment epithelium-derived factor (PEDF). Immunofluorescence of microglia/macrophage markers, Iba1 and CD11b, also revealed significant upregulation of immune cell markers in CD-exposed placenta. Our results suggest that in utero chronodisruption enhances placental immune cell expression, potentially programming a pro-inflammatory tissue environment.

Maternal circadian disruption is associated with variation in placental DNA methylation.

Circadian disruption is a common environmental and occupational exposure with public health consequences, but not much is known about whether circadian disruption affects in utero development. We investigated whether maternal circadian disruption, using night shift work as a proxy, is associated with variations in DNA methylation patterns of placental tissue in an epigenome-wide association study (EWAS) of night shift work. Here, we compared cytosine-guanosine dinucleotide (CpG) specific methylation genome-wide of placental tissue (measured with the Illumina 450K array) from participants (n = 237) in the Rhode Island Child Health Study (RICHS) who did (n = 53) and did not (n = 184) report working the night shift, using robust linear modeling and adjusting for maternal age, pre-pregnancy smoking, infant sex, maternal adversity, and putative cell mixture. Statistical analyses were adjusted for multiple comparisons and results presented with Bonferroni or Benjamini and Hochberg (BH) adjustment for false discovery rate. Night shift work was associated with differential methylation in placental tissue, including CpG sites in the genes NAV1, SMPD1, TAPBP, CLEC16A, DIP2C, FAM172A, and PLEKHG6 (Bonferroni-adjusted p<0.05). CpG sites within NAV1, MXRA8, GABRG1, PRDM16, WNT5A, and FOXG1 exhibited the most hypomethylation, while CpG sites within TDO2, ADAMTSL3, DLX2, and SERPINA1 exhibited the most hypermethylation (BH q<0.10). Functional analysis indicated GO-terms associated with cell-cell adhesion and enriched GWAS results for psoriasis. Night shift work was associated with differential methylation of the placenta, which may have implications for fetal health and development. This is the first study to examine the epigenetic impacts of night shift exposure, as a proxy for circadian disruption, on placental methylation in humans, and, while results should be interpreted with caution, suggests circadian disruption may have epigenetic impacts.

Protein-coding variants implicate novel genes related to lipid homeostasis contributing to body-fat distribution.

Body-fat distribution is a risk factor for adverse cardiovascular health consequences. We analyzed the association of body-fat distribution, assessed by waist-to-hip ratio adjusted for body mass index, with 228,985 predicted coding and splice site variants available on exome arrays in up to 344,369 individuals from five major ancestries (discovery) and 132,177 European-ancestry individuals (validation). We identified 15 common (minor allele frequency, MAF ≥5%) and nine low-frequency or rare (MAF <5%) coding novel variants. Pathway/gene set enrichment analyses identified lipid particle, adiponectin, abnormal white adipose tissue physiology and bone development and morphology as important contributors to fat distribution, while cross-trait associations highlight cardiometabolic traits. In functional follow-up analyses, specifically in Drosophila RNAi-knockdowns, we observed a significant increase in the total body triglyceride levels for two genes (DNAH10 and PLXND1). We implicate novel genes in fat distribution, stressing the importance of interrogating low-frequency and protein-coding variants.

Concentration of Smaller High-Density Lipoprotein Particle (HDL-P) Is Inversely Correlated With Carotid Intima Media Thickening After Confounder Adjustment: The Multi Ethnic Study of Atherosclerosis (MESA).

BACKGROUND: Recent studies have failed to establish a causal relationship between high-density lipoprotein cholesterol levels (HDL-C) and cardiovascular disease (CVD), shifting focus to other HDL measures. We previously reported that smaller/denser HDL levels are protective against cerebrovascular disease. This study sought to determine which of small+medium HDL particle concentration (HDL-P) or large HDL-P was more strongly associated with carotid intima-media thickening (cIMT) in an ethnically diverse cohort. METHODS AND RESULTS: In cross-sectional analyses of participants from the Multi Ethnic Study of Atherosclerosis (MESA), we evaluated the associations of nuclear magnetic resonance spectroscopy-measured small+medium versus large HDL-P with cIMT measured in the common and internal carotid arteries, through linear regression. After adjustment for CVD confounders, low-density lipoprotein cholesterol (LDL-C), HDL-C, and small+medium HDL-P remained significantly and inversely associated with common (coefficient=-1.46 µm; P=0.00037; n=6512) and internal cIMT (coefficient=-3.82 µm; P=0.0051; n=6418) after Bonferroni correction for 4 independent tests (threshold for significance=0.0125; α=0.05/4). Large HDL-P was significantly and inversely associated with both cIMT outcomes before HDL-C adjustment; however, after adjustment for HDL-C, the association of large HDL-P with both common (coefficient=1.55 µm; P=0.30; n=6512) and internal cIMT (coefficient=4.84 µm; P=0.33; n=6418) was attenuated. In a separate sample of 126 men, small/medium HDL-P was more strongly correlated with paraoxonase 1 activity (rp=0.32; P=0.00023) as compared to both total HDL-P (rp=0.27; P=0.0024) and large HDL-P (rp=0.02; P=0.41) measures. CONCLUSIONS: Small+medium HDL-P is significantly and inversely correlated with cIMT measurements. Correlation of small+medium HDL-P with cardioprotective paraoxonase 1 activity may reflect a functional aspect of HDL responsible for this finding.

Burden of potentially pathologic copy number variants is higher in children with isolated congenital heart disease and significantly impairs covariate-adjusted transplant-free survival.

OBJECTIVES: Copy number variants (CNVs) are duplications or deletions of genomic regions. Large CNVs are potentially pathogenic and are overrepresented in children with congenital heart disease (CHD). We sought to determine the frequency of large CNVs in children with isolated CHD, and to evaluate the relationship of these potentially pathogenic CNVs with transplant-free survival. METHODS: These cases are derived from a prospective cohort of patients with nonsyndromic CHD (n = 422) identified before first surgery. Healthy pediatric controls (n = 500) were obtained from the electronic Medical Records and Genetic Epidemiology Network, and CNV frequency was contrasted for CHD cases and controls. CNVs were determined algorithmically; subsequently screened for >95% overlap between 2 methods, size (>300 kb), quality score, overlap with a gene, and novelty (absent from databases of known, benign CNVs); and separately validated by quantitative polymerase chain reaction. Survival likelihoods for cases were calculated using Cox proportional hazards modeling to evaluate the joint effect of CNV burden and known confounders on transplant-free survival. RESULTS: Children with nonsyndromic CHD had a higher burden of potentially pathogenic CNVs compared with pediatric controls (12.1% vs 5.0%; P = .00016). Presence of a CNV was associated with significantly decreased transplant-free survival after surgery (hazard ratio, 3.42; 95% confidence interval, 1.66-7.09; P = .00090) with confounder adjustment. CONCLUSIONS: We confirm that children with isolated CHD have a greater burden of rare/large CNVs. We report a novel finding that these CNVs are associated with an adjusted 2.55-fold increased risk of death or transplant. These data suggest that CNV burden is an important modifier of survival after surgery for CHD.