Variable DNA methylation in neonates mediates the association between prenatal smoking and birth weight.

There is great interest in the role epigenetic variation induced by non-genetic exposures may play in the context of health and disease. In particular, DNA methylation has previously been shown to be highly dynamic during the earliest stages of development and is influenced by in utero exposures such as maternal smoking and medication. In this study we sought to identify the specific DNA methylation differences in blood associated with prenatal and birth factors, including birth weight, gestational age and maternal smoking. We quantified neonatal methylomic variation in 1263 infants using DNA isolated from a unique collection of archived blood spots taken shortly after birth (mean = 6.08 days; s.d. = 3.24 days). An epigenome-wide association study (EWAS) of gestational age and birth weight identified 4299 and 18 differentially methylated positions (DMPs) respectively, at an experiment-wide significance threshold of p < 1 × 10-7. Our EWAS of maternal smoking during pregnancy identified 110 DMPs in neonatal blood, replicating previously reported genomic loci, including AHRR. Finally, we tested the hypothesis that DNA methylation mediates the relationship between maternal smoking and lower birth weight, finding evidence that methylomic variation at three DMPs may link exposure to outcome. These findings complement an expanding literature on the epigenomic consequences of prenatal exposures and obstetric factors, confirming a link between the maternal environment and gene regulation in neonates. This article is part of the theme issue 'Developing differences: early-life effects and evolutionary medicine'.

Elevated polygenic burden for autism is associated with differential DNA methylation at birth.

BACKGROUND: Autism spectrum disorder (ASD) is a severe neurodevelopmental disorder characterized by deficits in social communication and restricted, repetitive behaviors, interests, or activities. The etiology of ASD involves both inherited and environmental risk factors, with epigenetic processes hypothesized as one mechanism by which both genetic and non-genetic variation influence gene regulation and pathogenesis. The aim of this study was to identify DNA methylation biomarkers of ASD detectable at birth. METHODS: We quantified neonatal methylomic variation in 1263 infants-of whom ~ 50% went on to subsequently develop ASD-using DNA isolated from archived blood spots taken shortly after birth. We used matched genotype data from the same individuals to examine the molecular consequences of ASD-associated genetic risk variants, identifying methylomic variation associated with elevated polygenic burden for ASD. In addition, we performed DNA methylation quantitative trait loci (mQTL) mapping to prioritize target genes from ASD GWAS findings. RESULTS: We identified robust epigenetic signatures of gestational age and prenatal tobacco exposure, confirming the utility of DNA methylation data generated from neonatal blood spots. Although we did not identify specific loci showing robust differences in neonatal DNA methylation associated with later ASD, there was a significant association between increased polygenic burden for autism and methylomic variation at specific loci. Each unit of elevated ASD polygenic risk score was associated with a mean increase in DNA methylation of - 0.14% at two CpG sites located proximal to a robust GWAS signal for ASD on chromosome 8. CONCLUSIONS: This study is the largest analysis of DNA methylation in ASD undertaken and the first to integrate genetic and epigenetic variation at birth. We demonstrate the utility of using a polygenic risk score to identify molecular variation associated with disease, and of using mQTL to refine the functional and regulatory variation associated with ASD risk variants.

High-Quality Exome Sequencing of Whole-Genome Amplified Neonatal Dried Blood Spot DNA.

Stored neonatal dried blood spot (DBS) samples from neonatal screening programmes are a valuable diagnostic and research resource. Combined with information from national health registries they can be used in population-based studies of genetic diseases. DNA extracted from neonatal DBSs can be amplified to obtain micrograms of an otherwise limited resource, referred to as whole-genome amplified DNA (wgaDNA). Here we investigate the robustness of exome sequencing of wgaDNA of neonatal DBS samples. We conducted three pilot studies of seven, eight and seven subjects, respectively. For each subject we analysed a neonatal DBS sample and corresponding adult whole-blood (WB) reference sample. Different DNA sample types were prepared for each of the subjects. Pilot 1: wgaDNA of 2x3.2mm neonatal DBSs (DBS_2x3.2) and raw DNA extract of the WB reference sample (WB_ref). Pilot 2: DBS_2x3.2, WB_ref and a WB_ref replica sharing DNA extract with the WB_ref sample. Pilot 3: DBS_2x3.2, WB_ref, wgaDNA of 2x1.6 mm neonatal DBSs and wgaDNA of the WB reference sample. Following sequencing and data analysis, we compared pairwise variant calls to obtain a measure of similarity--the concordance rate. Concordance rates were slightly lower when comparing DBS vs WB sample types than for any two WB sample types of the same subject before filtering of the variant calls. The overall concordance rates were dependent on the variant type, with SNPs performing best. Post-filtering, the comparisons of DBS vs WB and WB vs WB sample types yielded similar concordance rates, with values close to 100%. WgaDNA of neonatal DBS samples performs with great accuracy and efficiency in exome sequencing. The wgaDNA performed similarly to matched high-quality reference--whole-blood DNA--based on concordance rates calculated from variant calls. No differences were observed substituting 2x3.2 with 2x1.6 mm discs, allowing for additional reduction of sample material in future projects.

UGT1A1*28 Genotypes and Respiratory Disease in Very Preterm Infants: A Cohort Study.

BACKGROUND: Respiratory disease in the very preterm infant is frequent and often severe. Bilirubin is both a potent neurotoxin and antioxidant, and may have a clinical impact on preterm respiratory disease. The Gilbert genotype, the UGT1A1*28 allele, is the major known genetic cause of variation in bilirubin. OBJECTIVES: To study the association between respiratory disease in the very preterm infant and the UGT1A1*28 allele. METHODS: This is a cohort study of 1,354 very preterm infants (gestational age <32 weeks) born in Jutland, Denmark in 1997-2011. Genotypes were obtained from the Danish Neonatal Screening Biobank, and clinical information was obtained from the databases of two tertiary neonatal intensive care units. Outcomes were the need for surfactant therapy, any need for and duration of supplementary oxygen and bronchopulmonary dysplasia (BPD). RESULTS: Per UGT1A1*28 allele, odds were increased for any need of supplementary oxygen (odds ratio 1.26; 1.05-1.50) and for BPD (odds ratio 1.71; 1.23-2.39), the need of supplementary oxygen increased by 6.38 days (1.87-10.89), and chance per day of no longer needing supplementary oxygen was reduced (hazard rate 0.84; 0.76-0.93). No effect was observed for need of surfactant treatment (odds ratio 1.08; 0.91-1.28). Hardy-Weinberg equilibrium was unlikely for the cohort (p < 0.012). This could be explained by death prior to genotype sampling. In tests of robustness this failed to explain the primary results. CONCLUSIONS: Compared to the common genotype, UGT1A1*28 genotypes were associated with an increased need of oxygen supplementation and risk of BPD in very preterm newborns.

The role of the sodium current complex in a nonreferred nationwide cohort of sudden infant death syndrome.

BACKGROUND: Sudden infant death syndrome (SIDS) is the most common cause of death in infants between the age of 1 month and 1 year. Rare variants in Nav1.5 encoded by SCN5A are known to play a role in SIDS; however, the combined role of the sodium current complex is unknown. OBJECTIVE: The purpose of this study was to investigate the role of the sodium current complex in a nonreferred nationwide cohort of SIDS cases. METHODS: DNA was extracted from dried blood spot samples from the Danish Neonatal Screening Biobank. In total, 66 non-referred SIDS cases born in Denmark in the period of 2000-2006 were screened for genetic variants in the 8 major genes involved in the regulation of the Nav1.5 channel complex: SCN5A, SCN1B, SCN2B, SCN3B, SCN4B, GPD1L, SNTA1, and CAV3. Patch-clamp analyses were performed on variants not previously characterized. RESULTS: In total, 8 patients (12%) had nonsynonymous rare variants in the sodium current genes. SCN5A harbored 6 rare variants (R458C, R535*, S1103Y, R1193Q, S1609L, and Q1909R); CAV3, 1 rare variant (T78M); GPD1L, 1 rare variant (R220H); and SCN3B, 1 rare variant (L10P). Four variants were considered likely pathogenic and 5 variants of unknown significance. SCN5A R1193Q and GPD1L R220H (both considered variants of unknown significance) were present in the same infant. Functional analysis of variants not previously characterized (R458C, S1609L, and Q1909R in SCN5A) predominantly revealed increased transient and sustained sodium current. CONCLUSION: In a nonreferred nationwide Danish cohort of SIDS cases, up to 5/66 (7.5%) of SIDS cases can be explained by genetic variants in the sodium channel complex genes.

Extreme neonatal hyperbilirubinemia and a specific genotype: a population-based case-control study.

OBJECTIVES: Extreme hyperbilirubinemia (plasma bilirubin ≥ 24.5 mg/dL) is an important risk factor for severe bilirubin encephalopathy. Several risk factors for hyperbilirubinemia are known, but in a large number of patients, a causal factor is never established. UGT1A1 is the rate-limiting enzyme in bilirubin's metabolism. The genotype of Gilbert syndrome, the UGT1A1*28 allele, causes markedly reduced activity of this enzyme, but its association with neonatal hyperbilirubinemia is uncertain and its relationship with extreme hyperbilirubinemia has not been studied. We examined whether the UGT1A1*28 allele is associated with extreme hyperbilirubinemia. METHODS: The UGT1A1*28 allele was assessed in a case-control study of 231 white infants who had extreme hyperbilirubinemia in Denmark from 2000 to 2007 and 432 white controls. Cases were identified in the Danish Extreme Hyperbilirubinemia Database that covers the entire population. Genotypes were obtained through the Danish Neonatal Screening Biobank. Subgroup analysis was done for AB0 incompatible cases. RESULTS: No association was found between the UGT1A1*28 allele and extreme hyperbilirubinemia. With the common genotype as reference, the odds ratio of extreme hyperbilirubinemia was 0.87 (range, 0.68-1.13) for UGT1A1*28 heterozygotes and 0.77 (range, 0.46-1.27) for homozygotes. Also, no association was found for AB0 incompatible cases. CONCLUSIONS: The UGT1A1*28 allele was not associated with risk for extreme hyperbilirubinemia in this study.

Mannose-binding lectin gene, MBL2, polymorphisms are not associated with susceptibility to invasive pneumococcal disease in children.

BACKGROUND: Most children are transiently colonized with Streptococcus pneumoniae, but very few develop invasive pneumococcal disease (IPD). Host genetic variation of innate immunity may predispose to IPD. We investigated the effect of genetic variation in the mannose-binding lectin gene, MBL2, on susceptibility and disease severity of IPD in previously healthy children aged <5 years. METHODS: IPD cases were identified through national registries. DNA was obtained from the Danish Neonatal Screening Biobank. Pneumococcal serotypes were determined by Quellung reaction. The associations between MBL2 diplotypes and IPD susceptibility, serotypes, and outcome were investigated using logistic regression analysis. RESULTS: We included 372 cases with meningitis, 907 with bacteremia, and 1263 age- and sex-matched controls; 2372 individuals were successfully genotyped and assigned MBL2 diplotypes. The median age in our combined case series was 13 months. Children with defective diplotypes were not at higher risk for meningitis than children with other diplotypes (odds ratio [OR], 0.85; 95% confidence interval [CI], .56-1.28). Similar results were found for bacteremia (OR, 0.89; 95% CI, .68-1.15) as well as for all cases (OR, 0.87; 95% CI, .70-1.09). There was no association with susceptibility to recurrent IPD (n = 12) for children with defective diplotypes compared with cases with a single episode (OR, 0.53; 95% CI, .07-4.13) and with all controls (OR, 0.46; 95% CI, .06-3.56). There was no association between diplotypes and mortality or between diplotypes and pneumococcal serotypes. CONCLUSIONS: Defective MBL2 polymorphisms did not predict increased IPD susceptibility in children born in Northern Europe.

Archived neonatal dried blood spot samples can be used for accurate whole genome and exome-targeted next-generation sequencing.

Dried blood spot samples (DBSS) have been collected and stored for decades as part of newborn screening programmes worldwide. Representing almost an entire population under a certain age and collected with virtually no bias, the Newborn Screening Biobanks are of immense value in medical studies, for example, to examine the genetics of various disorders. We have previously demonstrated that DNA extracted from a fraction (2×3.2mm discs) of an archived DBSS can be whole genome amplified (wgaDNA) and used for accurate array genotyping. However, until now, it has been uncertain whether wgaDNA from DBSS can be used for accurate whole genome sequencing (WGS) and exome sequencing (WES). This study examined two individuals represented by three different types of samples each: whole-blood (reference samples), 3-year-old DBSS spotted with reference material (refDBSS), and 27- to 29-year-old archived neonatal DBSS (neoDBSS) stored at -20°C in the Danish Newborn Screening Biobank. The reference samples were genotyped using an Illumina Omni2.5M array, and all samples were sequenced on a HighSeq2000 Paired-End flow cell. First, we compared the array single nucleotide polymorphism (SNP) genotype data to the single nucleotide variation (SNV) calls from the WGS and WES SNV calls. We also compared the WGS and WES reference sample SNV calls to the DBSS SNV calls. The overall performance of the archived DBSS was similar to the whole blood reference sample. Plotting the error rates relative to coverage revealed that the error rates of DBSS were similar to that of their reference samples. SNVs called with a coverage<×8 had error rates between 1.5 and 35%, whereas the error rates of SNVs called with a coverage≥8 were <1.5%. In conclusion, the wgaDNA amplified from both new and old neonatal DBSS perform as well as their whole-blood reference samples with regards to error rates, strongly indicating that neonatal DBSS collected shortly after birth and stored for decades comprise an excellent resource for NGS studies of disease.

DNA methylome profiling using neonatal dried blood spot samples: a proof-of-principle study.

DNA methylation is the most common DNA modification and perhaps the best described epigenetic modification. It is believed to be important for genomic imprinting and gene regulation and has been associated with the development of diseases such as schizophrenia and some types of cancer. Neonatal dried blood spot samples, commonly known as Guthrie cards, are routinely collected worldwide to screen newborns for diseases. Some countries, including Denmark, have been storing the excess neonatal dried blood spot samples in biobanks for decades. Representing a high percentage of the population under a certain age, the neonatal dried blood spot samples are a potential alternative to collecting new samples to study diseases. As such, neonatal dried blood spot samples have previously been used for DNA genotyping studies with excellent results. However, the amount of material available for research is often limited, challenging researchers to generate the most data from a limited quantity of material. In this proof-of-principle study, we address whether two 3.2mm disks punched from a neonatal dried blood spot sample contain enough DNA for genome-wide methylome profiling, measuring 27,578 loci at the same time. We selected two subjects and carried out the following with each: 1) collected an adult whole-blood sample as reference, 2) spotted a fraction of the whole-blood sample onto a similar type of filter paper as used in the newborn screening and stored it for 3years to serve as a dried blood spot reference, and 3) identified the archived neonatal dried blood spot samples, stored for 26-28years, in the Danish Newborn Screening Biobank as a representative of the archived samples. For comparison, we used two different kits for DNA extraction. The DNA, extracted using the Extract-N-Amp Blood PCR kit, was analyzed, and no statistically significant differences were observed (P<0.001) when we compared the methylation profile of the reference whole-blood samples to the dried blood spot references. This indicates that two 3.2mm disks contain enough material for reliable methylome profiling and that storing the whole-blood sample on neonatal dried blood spot filter paper for 3years does not interfere with the outcome of the analysis. Furthermore, we compared the adult DNA methylation profile to the neonatal dried blood spot sample profile. Approximately 50 sites in the subjects were significantly (P<0.001) different in the newborn sample compared with the adult sample. Both being healthy adults and the high quality of the DNA methylation array led to the conclusion that the archived neonatal dried blood spot samples can be used for methylome profiling, despite decades of storage and DNA degradation. In conclusion, we show that reliable methylome data can be obtained from old neonatal dried blood spot samples, by using a reasonable amount of the limited resource. This further adds to the use of neonatal dried blood spot samples in genetic research and screening and paves the way for unique population-based studies of epigenetic modifications after birth.

Joint analysis of SNPs and proteins identifies regulatory IL18 gene variations decreasing the chance of spastic cerebral palsy.

Cerebral palsy (CP) is a permanent disorder, affecting 2-3 per 1,000 live born children, disturbing movement and posture. Spastic limbs affects about 70-80% of the CP children, and this group is the target of our study. CP is considered a multifactorial condition believed to be provoked by, for example, preterm birth, infection during pregnancy, neural disorders, and genetics, to mention some. Interestingly, the cytokine network is believed to be involved in many of these disorders. In this study, including 203 spastic CP cases and 167 controls, we measured the levels of 25 cytokine proteins, and genotyped 159 SNPs in their gene loci. Using logistic regression, we estimated the genetic association of SNP genotypes to spastic CP. In addition, fitting a Tobit regression model for each protein and each SNP in the respective gene loci, we estimated three regression coefficients corresponding three different effects of the genetic variation on the protein level. Intriguingly, two IL18 loci SNPs (rs549908:A>C and rs1290349:C>A) showed a protective effect against spastic CP, and interestingly both were associated to a decreased epidemiological expression of IL-18 protein. By joining protein data to genetic information, we have provided new data suggesting IL18's involvement in the pathogenesis of spastic CP.

The prevalence of mutations in KCNQ1, KCNH2, and SCN5A in an unselected national cohort of young sudden unexplained death cases.

INTRODUCTION: Sudden unexplained death account for one-third of all sudden natural deaths in the young (1-35 years). Hitherto, the prevalence of genopositive cases has primarily been based on deceased persons referred for postmortem genetic testing. These deaths potentially may represent the worst of cases, thus possibly overestimating the prevalence of potentially disease causing mutations in the 3 major long-QT syndrome (LQTS) genes in the general population. We therefore wanted to investigate the prevalence of mutations in an unselected population of sudden unexplained deaths in a nationwide setting. METHODS: DNA for genetic testing was available for 44 cases of sudden unexplained death in Denmark in the period 2000-2006 (equaling 33% of all cases of sudden unexplained death in the age group). KCNQ1, KCNH2, and SCN5A were sequenced and in vitro electrophysiological studies were performed on novel mutations. RESULTS: In total, 5 of 44 cases (11%) carried a mutation in 1 of the 3 genes corresponding to 11% of all investigated cases (R190W KCNQ1, F29L KCNH2 (2 cases), P297S KCNH2 and P1177L SCN5A). P1177L SCN5A has not been reported before. In vitro electrophysiological studies of P1177L SCN5A revealed an increased sustained current suggesting a LQTS phenotype. CONCLUSION: In a nationwide setting, the genetic investigation of an unselected population of sudden unexplained death cases aged 1-35 years finds a lower than expected number of mutations compared to referred populations previously reported. We therefore conclude that the prevalence of mutations in the 3 major LQTS associated genes may not be as abundant as previously estimated.

Polymorphisms in the promoter region of relaxin-2 and preterm birth: involvement of relaxin in the etiology of preterm birth.

BACKGROUND: Circulating relaxin levels have repeatedly been associated with preterm birth (PTB). Our aim was to investigate if mothers carrying promoter single nucleotide polymorphisms (SNPs) in one of the three relaxin genes (RLN1, RLN2, RLN3) are predisposed for PTB. PATIENTS AND METHODS: Maternal DNA from 80 preterm cases (40 very preterm births (24-34 weeks) and 40 moderate preterm (34-36 weeks) and 40 controls (term delivery)) nested in the Danish National Birth Cohort were examined for nine SNPs. RESULTS: Maternal homozygosity of the rarer allele in the relaxin 2 gene (RLN2, rs10115467 and rs4742076) had an increased risk of moderate PTB (odds ratio 4.1 [95% CI 1.4-12] and 8.8 [95% CI 1.03-75] respectively). Only rs10115467 remained significant after correction for multiple testing. CONCLUSION: Women homozygous for prevalent SNPs in the RLN2 gene may have a genetic susceptibility for PTB.

Genotyping whole-genome-amplified DNA from 3- to 25-year-old neonatal dried blood spot samples with reference to fresh genomic DNA.

Stored surplus of dried blood spot (DBS) samples from neonatal screening programs constitute a vast potential for large genetic epidemiological studies. However, age of the samples and the small amounts of DNA available may limit their usage. In this study we validate genotyping accuracy and efficiency of whole-genome-amplified DNA (wgaDNA) obtained from stored DBS samples, with reference to fresh genomic DNA from the same individuals. DBS samples from 29 volunteers, stored for up to 25 years, in the Danish Neonatal Screening Biobank were included and three DNA extraction methods, each using one 3.2 mm disk, were evaluated. Four whole-genome amplification kits, and one re-amplification kit, were used. Thirty-one SNPs were genotyped using the Sequenom platform and the wgaDNA samples calls were compared with their references for accuracy and efficiency evaluation. The genotype calls done blinded by the user had in many setups a 100% call- and concordance rate. Our results showed that genotyping performance is dependent on the combination of extraction procedure and amplification method, whereas years of storage did not seem to influence in this study. Based on these results we conclude that DBS samples should be considered a reliable and potential resource for future genotyping studies.

Polymorphisms in the tumor necrosis factor alpha and interleukin 1-beta promoters with possible gene regulatory functions increase the risk of preterm birth.

OBJECTIVE: To investigate the relation between 19 selected single nucleotide polymorphisms in three cytokine genes, tumor necrosis factor alpha (TNFA), interleukin 1-beta (IL1B) and interleukin 6 (IL6) and preterm birth (<37 weeks' gestation). DESIGN: Case-control association study. SAMPLE: A total of 117 singleton pregnant Danish Caucasian women, including 62 preterm birth cases and 55 controls (birth>or=37 weeks). METHODS: Genotyping was performed using TaqMan probes and traditional sequencing. Descriptive statistics were carried out with Fisher's exact test and Wilcoxon rank-sum test. All genetic data were tested for Hardy-Weinberg equilibrium and analyzed using logistic regression, 2x2 proportions or chi(2). Haplotypes were estimated for each gene and permutation used for association testing. RESULTS: Women carrying the TNFA -857 C>T rare allele (T) and those homozygous for the IL1B -31 T>C and IL1B -511 C>T rare alleles (C and T) have an increased risk of preterm birth with OR 3.1 (95% CI: 1.0-10.3) and OR 6.4 (95% CI: 1.3-60.5), respectively. Two estimated TNFA haplotypes were associated with preterm birth with OR 3.1 (p=0.037) and OR 2.7 (p=0.045). CONCLUSION: Polymorphisms in the cytokine genes TNFA and IL1B may increase the risk of preterm birth, possibly by a dysregulation of the immune system in pregnancy.