Marie Bækvad-Hansen

An epigenetic clock for gestational age at birth based on blood methylation data

BackgroundGestational age is often used as a proxy for developmental maturity by clinicians and researchers alike. DNA methylation has previously been shown to be associated with age and has been used to accurately estimate chronological age in children and adults. In the current study, we examine whether DNA methylation in cord blood can be used to estimate gestational age at birth.ResultsWe find that gestational age can be accurately estimated from DNA methylation of neonatal cord blood and blood spot samples. We calculate a DNA methylation gestational age using 148 CpG sites selected through elastic net regression in six training datasets. We evaluate predictive accuracy in nine testing datasets and find that the accuracy of the DNA methylation gestational age is consistent with that of gestational age estimates based on established methods, such as ultrasound. We also find that an increased DNA methylation gestational age relative to clinical gestational age is associated with birthweight independent of gestational age, sex, and ancestry.ConclusionsDNA methylation can be used to accurately estimate gestational age at or near birth and may provide additional information relevant to developmental stage. Further studies of this predictor are warranted to determine its utility in clinical settings and for research purposes. When clinical estimates are available this measure may increase accuracy in the testing of hypotheses related to developmental age and other early life circumstances.

Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection

Schizophrenia is a debilitating psychiatric condition often associated with poor quality of life and decreased life expectancy. Lack of progress in improving treatment outcomes has been attributed to limited knowledge of the underlying biology, although large-scale genomic studies have begun to provide insights. We report a new genome-wide association study of schizophrenia (11,260 cases and 24,542 controls), and through meta-analysis with existing data we identify 50 novel associated loci and 145 loci in total. Through integrating genomic fine-mapping with brain expression and chromosome conformation data, we identify candidate causal genes within 33 loci. We also show for the first time that the common variant association signal is highly enriched among genes that are under strong selective pressures. These findings provide new insights into the biology and genetic architecture of schizophrenia, highlight the importance of mutation-intolerant genes and suggest a mechanism by which common risk variants persist in the population.A new GWAS of schizophrenia (11,260 cases and 24,542 controls) and meta-analysis identifies 50 new associated loci and 145 loci in total. The common variant association signal is highly enriched in mutation-intolerant genes and in regions under strong background selection.

The iPSYCH2012 case|[ndash]|cohort sample: new directions for unravelling genetic and environmental architectures of severe mental disorders

The Integrative Psychiatric Research (iPSYCH) consortium has established a large Danish population-based Case–Cohort sample (iPSYCH2012) aimed at unravelling the genetic and environmental architecture of severe mental disorders. The iPSYCH2012 sample is nested within the entire Danish population born between 1981 and 2005, including 1 472 762 persons. This paper introduces the iPSYCH2012 sample and outlines key future research directions. Cases were identified as persons with schizophrenia (N=3540), autism (N=16 146), attention-deficit/hyperactivity disorder (N=18 726) and affective disorder (N=26 380), of which 1928 had bipolar affective disorder. Controls were randomly sampled individuals (N=30 000). Within the sample of 86 189 individuals, a total of 57 377 individuals had at least one major mental disorder. DNA was extracted from the neonatal dried blood spot samples obtained from the Danish Neonatal Screening Biobank and genotyped using the Illumina PsychChip. Genotyping was successful for 90% of the sample. The assessments of exome sequencing, methylation profiling, metabolome profiling, vitamin-D, inflammatory and neurotrophic factors are in progress. For each individual, the iPSYCH2012 sample also includes longitudinal information on health, prescribed medicine, social and socioeconomic information, and analogous information among relatives. To the best of our knowledge, the iPSYCH2012 sample is the largest and most comprehensive data source for the combined study of genetic and environmental aetiologies of severe mental disorders.

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.

Differential DNA methylation at birth associated with mental disorder in individuals with 22q11.2 deletion syndrome

Individuals with 22q11.2 deletion syndrome (DS) have an increased risk of comorbid mental disorders including schizophrenia, attention deficit hyperactivity disorder, depression, as well as intellectual disability. Although most 22q11.2 deletion carriers have the long 3-Mb form of the hemizygous deletion, there remains a large variation in the development and progression of psychiatric disorders, which suggests that alternative factors contribute to the pathogenesis. In this study we investigated whether neonatal DNA methylation signatures in individuals with the 22q11.2 deletion associate with mental disorder later in life. DNA methylation was measured genome-wide from neonatal dried blood spots in a cohort of 164 individuals with 22q11.2DS, including 48 individuals diagnosed with a psychiatric disorder. Among several CpG sites with P-value<10−6, we identified cg23546855 (P-value=2.15 × 10−7) mapping to STK32C to be associated with a later psychiatric diagnosis. Pathway analysis of the top findings resulted in the identification of several Gene Ontology pathways to be significantly enriched (P-value<0.05 after Benjamini–Hochberg correction); among them are the following: neurogenesis, neuron development, neuron projection development, astrocyte development, axonogenesis and axon guidance. In addition, we identified differentially methylated CpG sites in LRP2BP (P-value=5.37 × 10−8) to be associated with intellectual disability (F70–79), in TOP1 (P-value=1.86 × 10−7) with behavioral disorders (F90–98), in NOSIP (P-value=5.12 × 10−8) with disorders of psychological development (F80–89) and in SEMA4B (P-value=4.02 × 10−7) with schizophrenia spectrum disorders (F20–29). In conclusion, our study suggests an association of DNA methylation differences at birth with development of mental disorder later in life in 22q11.2DS individuals.

Gene expression profiling of archived dried blood spot samples from the Danish Neonatal Screening Biobank.

A large part of the human genome is transcribed into various forms of RNA, and the global gene expression profile (GEP) has been studied for several years using technology such as RNA-microarrays. In this study, we evaluate whether neonatal dried blood spot (DBS) samples stored in the Danish Neonatal Screening Biobank (DNSB) can be used for GEP. This paper is divided into sub-studies examining the effects of: 1) different whole transcriptome amplification kits (WTA); 2) years of storage and storage in room temperature (RT) versus freezers (-20°C) on DNSB DBS samples; 3) effects of RT storage vs freezer storage on DBS samples from the USA and DNSB, and 4) using smaller disc sizes, thereby decreasing DBS use. We present evidence that reliable and reproducible GEPs can be obtained using neonatal DBS samples. The main source of variation is the storage condition. When samples are stored at -20°C, the dynamic range is increased, and Pearson correlations are higher. Differential analysis reveals no statistically significant differences between samples collected a decade apart and stored at -20°C. However, samples stored at RT show differential expression for a third of the gene-specific probes. Our data also suggests that using alternate WTA kits significantly changes the GEP. Finally, the amount of input material, i.e., the size and number of DBS discs used, can be reduced to preserve this valuable and limited material. We conclude that DNSB DBS samples provide a reproducible resource for GEP. Results are improved if the cards are stored at -20°C. Furthermore, it is important to use a single type of kit for analysis because using alternate kits introduces differential expression.

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

BackgroundAutism 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.MethodsWe 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.ResultsWe 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.ConclusionsThis 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.

RNA sequencing of archived neonatal dried blood spots.

Neonatal dried blood spots (DBS) are routinely collected on standard Guthrie cards for all-comprising national newborn screening programs for inborn errors of metabolism, hypothyroidism and other diseases. In Denmark, the Guthrie cards are stored at − 20 °C in the Danish Neonatal Screening Biobank and each sample is linked to elaborate social and medical registries. This provides a unique biospecimen repository to enable large population research at a perinatal level. Here, we demonstrate the feasibility to obtain gene expression data from DBS using next-generation RNA sequencing (RNA-seq). RNA-seq was performed on five males and five females. Sequencing results have an average of > 30 million reads per sample. 26,799 annotated features can be identified with 64% features detectable without fragments per kilobase of transcript per million mapped reads (FPKM) cutoff; number of detectable features dropped to 18% when FPKM ≥ 1. Sex can be discriminated using blood-based sex-specific gene set identified by the Genotype-Tissue Expression consortium. Here, we demonstrate the feasibility to acquire biologically-relevant gene expression from DBS using RNA-seq which provide a new avenue to investigate perinatal diseases in a high throughput manner.

ASD and ADHD have a similar burden of rare protein-truncating variants

Autism spectrum disorder (ASD) and attention-deficit/hyperactivity disorder (ADHD) are substantially heritable, but individuals with psychiatric diagnoses often do not have blood drawn as part of routine medical procedure, making it difficult to collect large cohorts for genetic study. To overcome this challenge, we drew upon two Danish national resources: the Danish Neonatal Screening Biobank (DNSB) and the Danish national psychiatric registry. We have previously validated the use of archived bloodspots from the DNSB for genotyping and sequencing, and we recently performed common variant analysis on dried bloodspot material in both ASD and ADHD. Here, we present exome sequences from over 13,000 DNSB samples, finding that ASD and ADHD show a strikingly similar burden of rare protein-truncating variants, both significantly higher than controls. Additionally, the distributions of genes hit by these variants are not distinguishable between the two disorders, suggesting that many risk genes may be shared between them. These results motivate a combined analysis across ASD and ADHD, which--in conjunction with incorporation of the gnomAD reference database as additional population controls--leads to the identification of genes conferring general risk for childhood psychiatric disorders, including the novel gene MAP1A.

Prevalence of rearrangements in the 22q11.2 region and population-based risk of neuropsychiatric and developmental disorders in a Danish population: a case-cohort study

BACKGROUND Although the pathogenic nature of copy number variants (CNVs) on chromosome 22q11.2 has been recognised for decades, unbiased estimates of their population prevalence, mortality, disease risks, and diagnostic trajectories are absent. We aimed to provide the true population prevalence of 22q11.2 CNVs and associated trajectory of disease risk and mortality by use of the unbiased, representative Danish iPSYCH population case cohort. METHODS This case-cohort study was done on a population of 86 189 individuals selected from the iPSYCH case cohort of 1 472 762 singletons born in Denmark between May 1, 1981, and Dec 31, 2005, who have a known mother from the Danish Civil Registration System, were residents in Denmark at 1 year of age, and enrolled in the iPSYCH Initiative. We used epidemiological methods in conjunction with nationwide hospital registers to analyse the iPSYCH case cohort of individuals with attention-deficit hyperactivity disorder (ADHD), major depressive disorder, schizophrenia, autism, or bipolar disorder and a random population-based sample. The main outcomes assessed were the population prevalence of 22q11.2 rearrangements, and associated unbiased, population-adjusted estimates and 31-year disease risk trajectories for major neuropsychiatric disorders. FINDINGS Population prevalence in the Danish population was one in 3672 (seven of 25 704 [0·027%; 95% CI 0·012-0·057]) for deletions and one in 1606 (17 of 25 704 [0·066%; 0·040-0·107]) for duplications. Mortality after the age of 1 year among carriers was zero, and hazard ratios for neuropsychiatric disorders ranged from 2·60 to 82·44 for both rearrangements. By the age of 32 years, about 10% of individuals with deletions or duplications had developed ADHD, autism, or intellectual disability, and deletion carriers had higher probability than duplication carriers of co-occurring intellectual disability or epilepsy. INTERPRETATION The significantly different prevalence of 22q11.2 duplications and deletions indicates distinct selective pressures on these rearrangements. Although risk of congenital abnormalities, developmental delay, and intellectual disability is elevated in deletion carriers, the overall prevalence of neuropsychiatric disorders is higher in duplication carriers, which implies that identification and clinical monitoring should extend beyond congenital traits and into child and adolescent psychiatry. FUNDING Capital Regions Research Foundation for Mental Health Research, The Lundbeck Foundation, and US National Institutes of Health.