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.

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.

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.

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.

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.

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 (F₀) and hyperpitch (F₀ > 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.

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.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

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.

[The Medical Social Efficiency of Joint Confinement of Mother and Child in Place of Imprisonment]

The article, on the basis of the results of study of 2016, analyzes a comparative medical social effect of conjoint and separate residing of mother and child in the Children's Homes of the Penal Enforcement System. The materials of study consisted of original observations while working in the given system, data of reports of the Federal Penitentiary Service of Russia concerning medical care of women and children in 2011-2015 and data of address requests to the Children's homes under correctional facility in 2013 and 2015. The experience of functioning of the Children's homes of penitentiary system is summarized. The positive effect of conjoint keeping to indices of coverage with breast feeding is established. The coverage of natural and mixed feeding amounted to 65% in Children's Homes, including 83% in Children's Homes with conjoint residing of mother and child in Children's Homes and in other Children's Homes -- 43% (p<0.05). In 2015, a single refusal of child by condemned mother occurred in Children's Homes, while in 2013 there were 4 of such refusals (p<0.05). In Children homes with extended department of conjoint residence primary morbidity of children made up to 1,679‰ and in comparison, group -- 2014‰ that is reliably higher (p<0.05). According data of primary morbidity occurred decreasing of pool of diseases of respiratory organs from 1,002 to 872 per 1,000 children due to brining into operation of dormitory of conjoint residence of mother and child. This occurrence resulted in decreasing per 12% of number of hospitalizations by emergency indications. The conjoint residence of condemned mother and newborn results in harmonious growth, physical, neuropsychic development of child and also in development of motherhood feelings that is an important factor of re-socialization of the condemned woman. It is proposed to introduce normatively obligatory conjoint residence of mother and child in Children's homes of the Penal Enforcement System.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

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.

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.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

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.

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 (r_p=0.32; P=0.00023) as compared to both total HDL‐P (r_p=0.27; P=0.0024) and large HDL‐P (r_p=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.

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.

A genome-wide map of adeno-associated virus-mediated human gene targeting

To determine which genomic features promote homologous recombination, we created a genome-wide map of gene targeting sites. We used an adeno-associated virus vector to target identical loci introduced as transcriptionally active retroviral vectors. A comparison of ~2,000 targeted and untargeted sites showed that targeting occurred throughout the human genome and was not influenced by the presence of nearby CpG islands, sequence repeats or DNase I-hypersensitive sites. Targeted sites were preferentially located within transcription units, especially when the target loci were transcribed in the opposite orientation to their surrounding chromosomal genes. We determined the impact of DNA replication by mapping replication forks, which revealed a preference for recombination at target loci transcribed toward an incoming fork. Our results constitute the first genome-wide screen of gene targeting in mammalian cells and demonstrate a strong recombinogenic effect of colliding polymerases.

Effects of dietary components on high-density lipoprotein measures in a cohort of 1,566 participants

Background

Recent data suggest that an increased level of high-density lipoprotein cholesterol (HDL-C) is not causally protective against heart disease, shifting focus to other sub-phenotypes of HDL. Prior work on the effects of dietary intakes has focused largely on HDL-C. The goal of this study was to identify the dietary intakes that affect HDL-related measures: HDL-C, HDL-2, HDL-3, and apoA1 using data from a carotid artery disease case–control cohort.

Methods

A subset of 1,566 participants with extensive lipid phenotype data completed the Harvard Standardized Food Frequency Questionnaire to determine their daily micronutrient intake over the past year. Stepwise linear regression was used to separately evaluate the effects of dietary covariates on adjusted levels of HDL-C, HDL-2, HDL-3, and apoA1.

Results

Dietary folate intake was positively associated with HDL-C (p = 0.007), HDL-2 (p = 0.0011), HDL-3 (p = 0.0022), and apoA1 (p = 0.001). Alcohol intake and myristic acid (14:0), a saturated fat, were each significantly associated with increased levels of all HDL-related measures studied. Dietary carbohydrate and iron intake were significantly associated with decreased levels of all HDL-related measures. Magnesium intake was positively associated with HDL-C, HDL-2, and HDL-3 levels, but not apoA1 levels, while vitamin C was only associated with apoA1 levels. Dietary fiber and protein intake were both associated with HDL-3 levels alone.

Conclusions

This study is the first to report that dietary folate intake is associated with HDL-C, HDL-2, HDL-3, and apoA1 levels in humans. We further identify numerous dietary intake associations with apoA1, HDL-2, and HDL-3 levels. Given the shifting focus away from HDL-C, these data will prove valuable for future epidemiologic investigation of the role of diet and multiple HDL phenotypes in heart disease.

Validation of association of the apolipoprotein E ε2 allele with neurodevelopmental dysfunction after cardiac surgery in neonates and infants

OBJECTIVE

Apolipoprotein E (APOE) genotype is a determinant of neurologic recovery after brain ischemia and traumatic brain injury. The APOE ε2 allele has been associated with worse neurodevelopmental (ND) outcome after repair of congenital heart defects (CHD) in infancy. Replication of this finding in an independent cohort is essential to validate the observed genotype-phenotype association.

METHODS

The association of APOE genotype with ND outcomes was assessed in a combined cohort of patients with single-ventricle CHD enrolled in the Single Ventricle Reconstruction and Infant Single Ventricle trials. ND outcome was assessed at 14 months using the Psychomotor Development Index (PDI) and Mental Development Index (MDI) of the Bayley Scales of Infant Development-II. Stepwise multivariable regression was performed to develop predictive models for PDI and MDI scores.

RESULTS

Complete data were available for 298 of 435 patients. After adjustment for preoperative and postoperative covariates, the APOE ε2 allele was associated with a lower PDI score (P = .038). Patients with the ε2 allele had a PDI score approximately 6 points lower than those without the risk allele, explaining 1.04% of overall PDI variance, because the ε2 allele was present in only 11% of the patients. There was a marginal effect of the ε2 allele on MDI scores (P = .058).

CONCLUSIONS

These data validate the association of the APOE ε2 allele with adverse early ND outcomes after cardiac surgery in infants, independent of patient and operative factors. Genetic variants that decrease neuroresilience and impair neuronal repair after brain injury are important risk factors for ND dysfunction after surgery for CHD.

HDL-3 is a superior predictor of carotid artery disease in a case-control cohort of 1725 participants

Background

Recent data suggest that high‐density lipoprotein cholesterol (HDL‐C) levels are likely not in the causative pathway of atheroprotection, shifting focus from HDL‐C to its subfractions and associated proteins. This study's goal was to determine which HDL phenotype was the better predictor of carotid artery disease (CAAD).

Methods and Results

HDL‐2 and HDL‐3 were measured in 1725 participants of European ancestry in a prevalent case‐control cohort study of CAAD. Stratified analyses were conducted for men (n=1201) and women (n=524). Stepwise linear regression was used to determine whether HDL‐C, HDL‐2, HDL‐3, or apolipoprotein A1 was the best predictor of CAAD, while adjusting for the confounders of censored age, diabetes, and current smoking status. In both men and women, HDL‐3 was negatively associated with CAAD (P=0.0011 and 0.033 for men and women, respectively); once HDL‐3 was included in the model, no other HDL phenotype was significantly associated with CAAD. Addition of paraoxonase 1 activity to the aforementioned regression model showed a significant and independent (of HDL‐3) association with CAAD in men (P=0.001) but not in the smaller female subgroup.

Conclusions

This study is the first to contrast the associations of HDL‐2 and HDL‐3 with CAAD. We found that HDL‐3 levels were more predictive of CAAD status than HDL‐2, HDL‐C, or apolipoprotein A1. In addition, for men, paraoxonase 1 activity improved the overall model prediction for CAAD independently and additively with HDL‐3 levels. Further investigation into the molecular mechanisms through which HDL‐3 is associated with protection from CAAD is warranted.

Patient genotypes impact survival after surgery for isolated congenital heart disease

BACKGROUND

Survival after cardiac surgery in infancy requires adaptive responses from oxidative stress management and vascular regulation pathways. We tested the hypothesis that genetic variation in these pathways influences postoperative survival in nonsyndromic congenital heart disease children.

METHODS

This is an analysis of a cohort of nonsyndromic congenital heart disease patients who underwent cardiac surgery with cardiopulmonary bypass before 6 months of age (n=422). Six single nucleotide polymorphisms (SNPs) in six genes involved in oxidative stress and vascular response pathways, identified through a priori literature search, were tested for effects on transplant-free survival. Survival curves, adjusting for confounding covariates, were calculated using the Cox proportional hazard models.

RESULTS

Long-term survival was strongly associated with vascular endothelial growth factor A gene SNP rs833069 (p=7.03×10(-4)) and superoxide dismutase 2 gene SNP rs2758331 (p=0.019). To test for joint effects of the two SNPs on transplant-free survival, the genotypes were grouped to form a risk score reflecting the cumulative number of risk alleles (0 to 4 alleles per patient). A higher risk score based on the VEGFA and SOD2 SNP genotypes was associated with worse transplant-free survival (p=3.02×10(-4)) after confounder adjustment. The total burden of risk alleles was additive; subjects with the highest risk score of 4 (n=59 subjects, 14.2% of the cohort) had a total covariate-adjusted hazard ratio of 15.64 for worse transplant-free survival.

CONCLUSIONS

After cardiac surgery, infants who are homozygous for the high-risk alleles for both the VEGFA and SOD2 SNPs have an approximately 16-fold increased risk of death or heart transplant, suggesting that genetic variants are important modifiers of survival after surgery for congenital heart disease.

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.

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.

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 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.

Prognostic value of cortical magnetic stimulation in spinal cord injury

Cortical magnetic stimulation was performed in a consecutive series of 10 patients presenting within 15 days of traumatic spinal cord injury. In those patients with complete paraplegia or quadriplegia, motor evoked potentials at presentation were absent below the level of the lesion. Six months after the injury, potentials had returned in the biceps brachii and abductor pollicis brevis muscles in some quadriplegic cases, but remained absent from the tibialis anterior in all of this group. None of those with a complete lesion made a significant functional recovery. Of the three patients with incomplete quadriplegia, two showed a significant recovery after 6 months. Motor evoked potentials were recordable below the level of the lesion at presentation in these cases, although the latencies were prolonged. In the remaining patient who failed to improve, potentials were unrecordable throughout the study. This small pilot study suggests that cortical magnetic stimulation may be useful in refining the prognosis in patients with an incomplete spinal cord injury.

Prophylactic low dose heparin anticoagulant therapy in patients with spinal cord injuries: a retrospective study

A retrospective study was performed of 100 consecutive acute traumatic spinal injury patients with neurological deficit admitted to the Yorkshire regional spinal injuries centre prior to May 1990. Ninety-seven of these patients received prophylactic low dose subcutaneous heparin and 3 patients were excluded because of noted complications. Twenty-six patients developed thromboembolic complications. Apart from the fact that acute spinal cord injury patients are considered to be in an hypercoaguable state, the occurrence of deep vein thrombosis/pulmonary embolism (DVT/PE) in our study, even though these patients were on prophylactic low dose subcutaneous heparin, was related to delay in transfers, operative intervention, level of spinal cord injury and was possibly due to loss of some amount of heparin solution from the prefilled syringe during removal of air bubbles prior to the subcutaneous heparin injection.

Myotonic dystrophy and traumatic quadriplegia: case report

A thirty-year-old man with traumatic quadriplegia, was also found to have weakness above the level of the injury. He had facial weakness, difficulty in swallowing, and recurrent respiratory problems. A diagnosis of myotonic dystrophy was supported by examination of his sister. The problems of diagnosis, and the implications of the diagnosis on the management of the patient with myotonic dystrophy and a spinal injury are discussed.

Status epilepticus complicating migraine

We describe 2 patients who developed status epilepticus in association with migraine attacks. One patient died and at the autopsy had ischaemic changes confined to the occipital lobes. The second patient made an uneventful recovery.

Successful rehabilitation in conversion paralysis

A rehabilitation programme for patients with conversion paralysis has been introduced in which they are offered physical rehabilitation. During an eight month period between October 1984 and May 1985 six patients who had been diagnosed as dependent on wheelchairs owing to conversion paralysis for a mean of 3 years (range 1-6 years) were entered into the inpatient neurorehabilitation programme. All six patients were able to walk within a mean of 41 days (range 10-70 days), and then relinquished a variety of aids and allowances as a result of their regained mobility. They continued to be independent at outpatient review for a mean of 10 months (range 8-15 months). Successful rehabilitation from wheelchair dependency can be achieved by a cost effective, prolonged, inpatient neurorehabilitation programme.

Brain stem signs in acute cervical cord trauma

Four cases are reported which illustrate mechanisms by which patients with cervical cord trauma may develop signs of acute brain stem damage. Two patients suffered vertebro-basilar ischaemia as a result of trauma to the vertebral arteries. The remaining two patients developed cavitating lesions within the cervical cord which later dissected rostrally into the brain stem.

Are spinal injuries more common in rugby union than in rugby league football?

Images

Suprapubic catheterisation after spinal cord injury: a follow-up report

The spinal injury patients who were initially treated by suprapubic catheterisation and reported from this unit in 1976 have been reviewed. Fifteen had died by early 1982; only in two cases from renal causes. All but one of the 23 survivors seen has a normal blood urea, and 15 have normal intravenous pyelograms. Eight patients have abnormal IVPs; the abnormalities were insignificant in four, and have been treated in two. Two patients have unilateral nonfunctioning kidneys. These results suggest that no long term ill effects result from the technique. In view of its considerable administrative advantages suprapubic urinary drainage should become more widely used.

Some extrapyramidal motor signs in cases of undifferentiated amentia

A group of 50 mentally handicapped patients who were without gross motor deficit and had never been on anti-psychotic drugs were examined for signs of Parkinsonism and compared with a matched control group of non-handicapped people who were also examined. There was a very significant increase in masking, positive glabella tap, absence of arm swing, rigidity and hypokinesis.

Acute spinal epidural abscess

A case of acute spinal epidural abscess is reported demonstrating many of the important features in the natural history of the disorder. The pathological changes in the cord at the level of the lesion are described. A series of 12 patients is reviewed. The average interval between the initial consultation and the onset of complete paralysis is 5 days and it is essential that the diagnosis is made during this period. Two indicators of prognosis are discussed--the duration of complete paralysis and the extent of sensory disturbance.

Percutaneous cystostomy in paraplegia--a follow-up of 41 patients

Percutaneous cystostomy using a fine calibre (9 Charriere) plastic trocar and cannula has been used in 51 male patients in the early weeks following spinal cord injury. The results are satisfactory. Of the last ten patients none has been catheterised per urethram; eight of these patients are now passing urine satisfactorily using a condom urinal.

Peripheral motor nerve function in diabetic autonomic neuropathy

Motor conduction velocity was measured in the median, ulnar and common peroneal nerves of 32 diabetics with clinical features of autonomic neuropathy. The responses to the Valsalva manoeuvre and sustained handgrip, and the postural fall in blood pressure were used to assess the integrity of the autonomic nervous system. Abnormalities in the three autonomic function tests were significantly correlated with the forearm conduction velocity of the ulnar nerve, the conduction velocity and motor latency of the common peroneal nerve, and the H reflex. These results show that, in diabetics with autonomic neuropathy, abnormalities in the autonomic nervous system parallel changes in the peripheral nerves. Any diabetic with peripheral neuropathy should be examined for evidence of autonomic nervous system involvement.