Genome-wide scans using archived neonatal dried blood spot samples

Genome-wide Association Study of Susceptibility to Respiratory Syncytial Virus Hospitalization in Young Children <5 Years of age

BACKGROUND

Worldwide, respiratory syncytial virus (RSV) infections are among the most common causes of infant hospitalization. Host genetic factors influencing the risk and severity of RSV infection are not well known.

METHODS

We conducted a genome-wide association study (GWAS) to investigate single-nucleotide polymorphisms (SNPs) associated with severe RSV infections using a nested case-control design based on 2 Danish cohorts. We compared SNPs from 1786 children hospitalized with RSV to 45 060 controls without an RSV-coded hospitalization. We performed gene-based testing, tissue enrichment, gene-set enrichment, and a meta-analysis of the 2 cohorts. Finally, an analysis of potential associations between the severity of RSV infection and genetic markers was performed.

RESULTS

We did not detect any significant genome-wide associations between SNPs and RSV infection or the severity of RSV. We did find potential loci associated with RSV infections on chromosome 5 in 1 cohort but failed to replicate any signals in both cohorts.

CONCLUSIONS

Despite being the largest GWAS of severe RSV infection, we did not detect any genome-wide significant loci. This may be an indication of a lack of power or an absence of signal. Future studies might include mild illness and need to be larger to detect any significant associations.

Methods applied to neonatal dried blood spot samples for secondary research purposes: a scoping review

This scoping review aimed to synthesize the analytical techniques used and methodological limitations encountered when undertaking secondary research using residual neonatal dried blood spot (DBS) samples. Studies that used residual neonatal DBS samples for secondary research (i.e. research not related to newborn screening for inherited genetic and metabolic disorders) were identified from six electronic databases: Cochrane Library, Cumulative Index to Nursing and Allied Health Literature (CINAHL), Embase, Medline, PubMed and Scopus. Inclusion was restricted to studies published from 1973 and written in or translated into English that reported the storage, extraction and testing of neonatal DBS samples. Sixty-seven studies were eligible for inclusion. Included studies were predominantly methodological in nature and measured various analytes, including nucleic acids, proteins, metabolites, environmental pollutants, markers of prenatal substance use and medications. Neonatal DBS samples were stored over a range of temperatures (ambient temperature, cold storage or frozen) and durations (two weeks to 40.5 years), both of which impacted the recovery of some analytes, particularly amino acids, antibodies and environmental pollutants. The size of DBS sample used and potential contamination were also cited as methodological limitations. Residual neonatal DBS samples retained by newborn screening programs are a promising resource for secondary research purposes, with many studies reporting the successful measurement of analytes even from neonatal DBS samples stored for long periods of time in suboptimal temperatures and conditions.

Identification of women at high risk of postpartum psychiatric episodes: A population-based study quantifying relative and absolute risks following exposure to selected risk factors and genetic liability

BACKGROUND

We quantified relative and absolute risks of postpartum psychiatric episodes (PPE) following risk factors: Young age, past personal or family history of psychiatric disorders, and genetic liability.

METHODS

We conducted a register-based study using the iPSYCH2012 case-cohort sample. Exposures were personal history of psychiatric episodes prior to childbirth, being a young mother (giving birth before the age of 21.5 years), having a family history of psychiatric disorders, and a high (highest quartile) polygenic score (PGS) for major depression. PPE was defined within 12 months postpartum by prescription of psychotropic medication or in- and outpatient contact to a psychiatric facility. We included primiparous women born 1981-1999, giving birth before January 1st, 2016. We conducted Cox regression to calculate hazard ratios (HRs) of PPE, absolute risks were calculated using cumulative incidence functions.

RESULTS

We included 8174 primiparous women, and the estimated baseline PPE risk was 6.9% (95% CI 6.0%-7.8%, number of PPE cases: 2169). For young mothers with a personal and family history of psychiatric disorders, the absolute risk of PPE was 21.6% (95% CI 15.9%-27.8%). Adding information on high genetic liability to depression, the risk increased to 29.2% (95% CI 21.3%-38.4%) for PPE.

CONCLUSIONS

Information on prior personal and family psychiatric episodes as well as age may assist in estimating a personalized risk of PPE. Furthermore, additional information on genetic liability could add even further to this risk assessment.

Genetic predictors of neurocognitive outcomes in survivors of pediatric brain tumors

BACKGROUND

Neurocognitive deficits are common in pediatric brain tumor survivors. The use of single nucleotide polymorphism (SNP) analysis in DNA repair genes may identify children treated with radiation therapy for brain tumors at increased risk for treatment toxicity and adverse neurocognitive outcomes.

MATERIALS

The Human 660W-Quad v1.0 DNA BeadChip analysis (Illumina) was used to evaluate 1048 SNPs from 59 DNA repair genes in 46 subjects. IQ testing was measured by the Wechsler Intelligence Scale for Children. Linear regression was used to identify the 10 SNPs with the strongest association with IQ scores while adjusting for radiation type.

RESULTS

The low vs high IQ patient cohorts were well matched for time from first treatment to most recent IQ, first treatment age, sex, and treatments received. 5 SNPs on 3 different genes (CYP29, XRCC1, and BRCA1) and on 3 different chromosomes (10, 19, and 17) had the strongest association with most recent IQ score that was not modified by radiation type. Furthermore, 5 SNPs on 4 different genes (WRN, NR3C1, ERCC4, RAD51L1) on 4 different chromosomes (8, 5, 16, 14) had the strongest association with change in IQ independent of radiation type, first IQ, and years between IQ measures.

CONCLUSIONS

SNPs offer the potential to predict adverse neurocognitive outcomes in pediatric brain tumor survivors. Our results require validation in a larger patient cohort. Improving the ability to identify children at risk of treatment related neurocognitive deficits could allow for better treatment stratification and early cognitive interventions.

Preparation and processing of dried blood spots for microRNA sequencing

Dried blood spots (DBS) are biological samples commonly collected from newborns and in geographic areas distanced from laboratory settings for the purposes of disease testing and identification. MicroRNAs (miRNAs)-small non-coding RNAs that regulate gene activity at the post-transcriptional level-are emerging as critical markers and mediators of disease, including cancer, infectious diseases, and mental disorders. This protocol describes optimized procedural steps for utilizing DBS as a reliable source of biological material for obtaining peripheral miRNA expression profiles. We outline key practices, such as the method of DBS rehydration that maximizes RNA extraction yield, and the use of degenerate oligonucleotide adapters to mitigate ligase-dependent biases that are associated with small RNA sequencing. The standardization of miRNA readout from DBS offers numerous benefits: cost-effectiveness in sample collection and processing, enhanced reliability and consistency of miRNA profiling, and minimal invasiveness that facilitates repeated testing and retention of participants. The use of DBS-based miRNA sequencing is a promising method to investigate disease mechanisms and to advance personalized medicine.

Multi-ancestry genome-wide association study of 4069 children with glioma identifies 9p21.3 risk locus

BACKGROUND

Although recent sequencing studies have revealed that 10% of childhood gliomas are caused by rare germline mutations, the role of common variants is undetermined and no genome-wide significant risk loci for pediatric central nervous system tumors have been identified to date.

METHODS

Meta-analysis of 3 population-based genome-wide association studies comprising 4069 children with glioma and 8778 controls of multiple genetic ancestries. Replication was performed in a separate case-control cohort. Quantitative trait loci analyses and a transcriptome-wide association study were conducted to assess possible links with brain tissue expression across 18 628 genes.

RESULTS

Common variants in CDKN2B-AS1 at 9p21.3 were significantly associated with astrocytoma, the most common subtype of glioma in children (rs573687, P-value of 6.974e-10, OR 1.273, 95% CI 1.179-1.374). The association was driven by low-grade astrocytoma (P-value of 3.815e-9) and exhibited unidirectional effects across all 6 genetic ancestries. For glioma overall, the association approached genome-wide significance (rs3731239, P-value of 5.411e-8), while no significant association was observed for high-grade tumors. Predicted decreased brain tissue expression of CDKN2B was significantly associated with astrocytoma (P-value of 8.090e-8).

CONCLUSIONS

In this population-based genome-wide association study meta-analysis, we identify and replicate 9p21.3 (CDKN2B-AS1) as a risk locus for childhood astrocytoma, thereby establishing the first genome-wide significant evidence of common variant predisposition in pediatric neuro-oncology. We furthermore provide a functional basis for the association by showing a possible link to decreased brain tissue CDKN2B expression and substantiate that genetic susceptibility differs between low- and high-grade astrocytoma.

Polygenic liability, stressful life events and risk for secondary-treated depression in early life: a nationwide register-based case-cohort study

BACKGROUND

In this study, we examined the relationship between polygenic liability for depression and number of stressful life events (SLEs) as risk factors for early-onset depression treated in inpatient, outpatient or emergency room settings at psychiatric hospitals in Denmark.

METHODS

Data were drawn from the iPSYCH2012 case-cohort sample, a population-based sample of individuals born in Denmark between 1981 and 2005. The sample included 18 532 individuals who were diagnosed with depression by a psychiatrist by age 31 years, and a comparison group of 20 184 individuals. Information on SLEs was obtained from nationwide registers and operationalized as a time-varying count variable. Hazard ratios and cumulative incidence rates were estimated using Cox regressions.

RESULTS

Risk for depression increased by 35% with each standard deviation increase in polygenic liability (p < 0.0001), and 36% (p < 0.0001) with each additional SLE. There was a small interaction between polygenic liability and SLEs (β = -0.04, p = 0.0009). The probability of being diagnosed with depression in a hospital-based setting between ages 15 and 31 years ranged from 1.5% among males in the lowest quartile of polygenic liability with 0 events by age 15, to 18.8% among females in the highest quartile of polygenic liability with 4+ events by age 15.

CONCLUSIONS

These findings suggest that although there is minimal interaction between polygenic liability and SLEs as risk factors for hospital-treated depression, combining information on these two important risk factors could potentially be useful for identifying high-risk individuals.

A method to correct for the influence of bovine serum albumin-associated vitamin D metabolites in protein extracts from neonatal dried blood spots

Objective

We developed an assay to measure the concentration of 25 hydroxyvitamin D₂ and D₃ in protein extracts derived from stored neonatal dried blood spots. During this study, we postulated that these samples had been contaminated with exogenous vitamin D metabolites because of the addition of bovine serum albumin (BSA) as part of an extraction step undertaken 7 years earlier. The aim of the current study was to develop methods in order to adjust for this contamination.

Results

We identified between-plate variations in 25 hydroxyvitamin D₂ and D₃ concentrations which suggested the presence of three different BSA batches. Based on repeat extraction (without the addition of BSA) and testing of 395 samples, we developed models to correct for the exogenous 25 hydroxyvitamin D₂ and D_3. The regression models were Diff_25OHD3 = − 8.2 + 1.8* Diff_25OHD2 for low contamination, Diff_25OHD3 = 23.8 + 1.7* Diff_25OHD2 for middle contamination, and Diff_25OHD3 = 14.3 + 3.0* Diff_25OHD2 for high contamination. After these corrections, the three subsamples had comparable distributions within the expected range for both 25 hydroxyvitamin D₂ and D₃.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-06077-1.

Sensitive and Robust LC-MS/MS Assay to Quantify 25-Hydroxyvitamin D in Leftover Protein Extract from Dried Blood Spots

Neonatal dried blood spots (DBS) provide a remarkable resource for biobanks. These microsamples can provide information related to the genetic correlates of disease and can be used to quantify a range of analytes, such as proteins and small molecules. However, after routine neonatal screening, the amount of DBS sample available is limited. To optimize the use of these samples, there is a need for sensitive assays which are integrated across different analytic platforms. For example, after DNA extraction, protein extracts are available for additional analyses. We describe a sensitive and robust LC-MS/MS method for 25-hydroxyvitamin D₂ and 25-hydroxyvitamin D₃ optimized for leftover protein extracts from DBS, which has excellent recovery, precision, and accuracy.

The Eating Disorders Genetics Initiative (EDGI): study protocol

BACKGROUND

The Eating Disorders Genetics Initiative (EDGI) is an international investigation exploring the role of genes and environment in anorexia nervosa, bulimia nervosa, and binge-eating disorder.

METHODS

A total of 14,500 individuals with eating disorders and 1500 controls will be included from the United States (US), Australia (AU), New Zealand (NZ), and Denmark (DK). In the US, AU, and NZ, participants will complete comprehensive online phenotyping and will submit a saliva sample for genotyping. In DK, individuals with eating disorders will be identified by the National Patient Register, and genotyping will occur using bloodspots archived from birth. A genome-wide association study will be conducted within EDGI and via meta-analysis with other data from the Eating Disorders Working Group of the Psychiatric Genomics Consortium (PGC-ED).

DISCUSSION

EDGI represents the largest genetic study of eating disorders ever to be conducted and is designed to rapidly advance the study of the genetics of the three major eating disorders (anorexia nervosa, bulimia nervosa, and binge-eating disorder). We will explicate the genetic architecture of eating disorders relative to each other and to other psychiatric and metabolic disorders and traits. Our goal is for EDGI to deliver "actionable" findings that can be transformed into clinically meaningful insights.

TRIAL REGISTRATION

EDGI is a registered clinical trial: clinicaltrials.gov NCT04378101 .

Genome-wide association study across pediatric central nervous system tumors implicates shared predisposition and points to 1q25.2 (PAPPA2) and 11p12 (LRRC4C) as novel candidate susceptibility loci

INTRODUCTION

Central nervous system (CNS) tumors constitute the most common form of solid neoplasms in children, but knowledge on genetic predisposition is sparse. In particular, whether susceptibility attributable to common variants is shared across CNS tumor types in children has not been investigated. The purpose of this study was to explore potential common genetic risk variants exhibiting pleiotropic effects across pediatric CNS tumors. We also investigated whether such susceptibility differs between early and late onset of disease.

METHOD

A Danish nationwide genome-wide association study (GWAS) of 1,097 consecutive patients (< 15 years of age) with CNS tumors and a cohort of 4,745 population-based controls.

RESULTS

For both the overall cohort and patients diagnosed after the age of four, the strongest association was rs12064625 which maps to PAPPA2 at 1q25.2 (p = 3.400 × 10^-7 and 9.668 × 10^-8, respectively). PAPPA2 regulates local bioavailability of insulin-like growth factor I (IGF-I). IGF-I is fundamental to CNS development and is involved in tumorigenesis across a wide range of different cancers. For the younger children, the strongest association was provided by rs11036373 mapping to LRRC4C at 11p12 (p = 7.620 × 10^-7), which encoded protein acts as an axon guidance molecule during CNS development and has not formerly been associated with brain tumors.

DISCUSSION

This GWAS indicates shared susceptibility attributable to common variants across pediatric CNS tumor types. Variations in genetic loci with roles in CNS development appear to be involved, possibly via altered IGF-I related pathways.

A genome-wide association study on medulloblastoma

INTRODUCTION

Medulloblastoma is a malignant embryonal tumor of the cerebellum that occurs predominantly in children. To find germline genetic variants associated with medulloblastoma risk, we conducted a genome-wide association study (GWAS) including 244 medulloblastoma cases and 247 control subjects from Sweden and Denmark.

METHODS

Genotyping was performed using Illumina BeadChips, and untyped variants were imputed using IMPUTE2.

RESULTS

Fifty-nine variants in 11 loci were associated with increased medulloblastoma risk (p < 1 × 10-5), but none were statistically significant after adjusting for multiple testing (p < 5 × 10-8). Thirteen of these variants were genotyped, whereas 46 were imputed. Genotyped variants were further investigated in a validation study comprising 249 medulloblastoma cases and 629 control subjects. In the validation study, rs78021424 (18p11.23, PTPRM) was associated with medulloblastoma risk with OR in the same direction as in the discovery cohort (ORT = 1.59, pvalidation = 0.02). We also selected seven medulloblastoma predisposition genes for investigation using a candidate gene approach: APC, BRCA2, PALB2, PTCH1, SUFU, TP53, and GPR161. The strongest evidence for association was found for rs201458864 (PALB2, ORT = 3.76, p = 3.2 × 10-4) and rs79036813 (PTCH1, ORA = 0.42, p = 2.6 × 10-3).

CONCLUSION

The results of this study, including a novel potential medulloblastoma risk loci at 18p11.23, are suggestive but need further validation in independent cohorts.

Utilization of archived neonatal dried blood spots for genome-wide genotyping

Introduction

Heel pricks are performed on newborns for diagnostic screenings of various pre-symptomatic metabolic and genetic diseases. Excess blood is spotted on Guthrie cards and archived by many states in biobanks for follow-up diagnoses and public health research. However, storage environment may vary across biobanks and across time within biobanks. With increased applications of DNA extracted from spots for genetic studies, identifying factors associated with genotyping success is critical to maximize DNA quality for future studies.

Method

We evaluated 399 blood spots, which were part of a genome-wide association study of childhood leukemia risk in children with Down syndrome, archived at the Michigan Neonatal Biobank between 1992 and 2008. High quality DNA was defined as having post-quality control call rate ≥ 99.0% based on the Illumina GenomeStudio 2.0 GenCall algorithm after processing the samples on the Illumina Infinium Global Screening Array. Bivariate analyses and multivariable logistic regression models were applied to evaluate effects of storage environment and storage duration on DNA genotyping quality.

Results

Both storage environment and duration were associated with sample genotyping call rates (p-values < 0.001). Sample call rates were associated with storage duration independent of storage environment (p-trend = 0.006 for DBS archived in an uncontrolled environment and p-trend = 0.002 in a controlled environment). However, 95% of the total sample had high genotyping quality with a call rate ≥ 95.0%, a standard threshold for acceptable sample quality in many genetic studies.

Conclusion

Blood spot DNA quality was lower in samples archived in uncontrolled storage environments and for samples archived for longer durations. Still, regardless of storage environment or duration, neonatal biobanks including the Michigan Neonatal Biobanks can provide access to large collections of spots with DNA quality acceptable for most genotyping studies.

Genetics of suicide attempts in individuals with and without mental disorders: a population-based genome-wide association study

Family studies have shown an aggregation of suicidal behavior in families. Yet, molecular studies are needed to identify loci accounting for genetic heritability. We conducted a genome-wide association study and estimated single nucleotide polymorphisms (SNP) heritability for a suicide attempt. In a case-cohort study, national data on all individuals born in Denmark after 1981 and diagnosed with severe mental disorders prior to 2013 (n = 57,377) and individuals from the general population (n = 30,000) were obtained. After quality control, the sample consisted of 6024 cases with an incidence of suicide attempt and 44,240 controls with no record of a suicide attempt. Suggestive associations between SNPs, rs6880062 (p-value: 5.4 × 10^-8) and rs6880461 (p-value: 9.5 × 10^-8), and suicide attempt were identified when adjusting for socio-demographics. Adjusting for mental disorders, three significant associations, all on chromosome 20, were identified: rs4809706 (p-value: 2.8 × 10^-8), rs4810824 (p-value: 3.5 × 10^-8), and rs6019297 (p-value: 4.7 × 10⁸). Sub-group analysis of cases with affective disorders revealed SNPs associated with suicide attempts when compared to the general population for gene PDE4B. All SNPs explained 4.6% [CI-95: 2.9-6.3%] of the variation in suicide attempt. Controlling for mental disorders reduced the heritability to 1.9% [CI-95: 0.3-3.5%]. Affective and autism spectrum disorders exhibited a SNP heritability of 5.6% [CI-95: 1.9-9.3%] and 9.6% [CI-95: 1.1-18.1%], respectively. Using the largest sample to date, we identified significant SNP associations with suicide attempts and support for a genetic transmission of suicide attempt, which might not solely be explained by mental disorders.

Optimization of extraction of genomic DNA from archived dried blood spot (DBS): potential application in epidemiological research & bio banking

Background: Limited infrastructure is available to collect, store and transport venous blood in field epidemiological studies. Dried blood spot (DBS) is a robust potential alternative sample source for epidemiological studies & bio banking. A stable source of genomic DNA (gDNA) is required for long term storage in bio bank for its downstream applications. Our objective is to optimize the methods of gDNA extraction from stored DBS and with the aim of revealing its utility in large scale epidemiological studies.  Methods: The purpose of this study was to extract the maximum amount of gDNA from DBS on Whatman 903 protein saver card. gDNA was extracted through column  (Qiagen) & magnetic bead based (Invitrogen) methods. Quantification of extracted gDNA was performed with a spectrophotometer, fluorometer, and integrity analyzed by agarose gel electrophoresis.  Result: Large variation was observed in quantity & purity (260/280 ratio, 1.8-2.9) of the extracted gDNA. The intact gDNA bands on the electrophoresis gel reflect the robustness of DBS for gDNA even after prolonged storage time. The extracted gDNA amount 2.16 - 24 ng/µl is sufficient for its PCR based downstream application, but unfortunately it can't be used for whole genome sequencing or genotyping from extracted gDNA. Sequencing or genotyping can be achieved by after increasing template copy number through whole genome amplification of extracted gDNA. The obtained results create a base for future research to develop high-throughput research and extraction methods from blood samples. Conclusion: The above results reveal, DBS can be utilized as a potential and robust sample source for bio banking in field epidemiological studies.

Optimization of extraction of genomic DNA from archived dried blood spot (DBS): potential application in epidemiological research & bio banking

Background: Limited infrastructure is available to collect, store and transport venous blood in field epidemiological studies. Dried blood spot (DBS) is a robust potential alternative sample source for epidemiological studies & bio banking. A stable source of genomic DNA (gDNA) is required for long term storage in bio bank for its downstream applications. Our objective is to optimize the methods of gDNA extraction from stored DBS and with the aim of revealing its utility in large scale epidemiological studies. 

Methods: The purpose of this study was to extract the maximum amount of gDNA from DBS on Whatman 903 protein saver card. gDNA was extracted through column  (Qiagen) & magnetic bead based (Invitrogen) methods. Quantification of extracted gDNA was performed with a spectrophotometer, fluorometer, and integrity analyzed by agarose gel electrophoresis. 

Result: Large variation was observed in quantity & purity (260/280 ratio, 1.8-2.9) of the extracted gDNA. The intact gDNA bands on the electrophoresis gel reflect the robustness of DBS for gDNA even after prolonged storage time. The extracted gDNA amount 2.16 – 24 ng/µl is sufficient for its PCR based downstream application, but unfortunately it can’t be used for whole genome sequencing or genotyping from extracted gDNA. Sequencing or genotyping can be achieved by after increasing template copy number through whole genome amplification of extracted gDNA. The obtained results create a base for future research to develop high-throughput research and extraction methods from blood samples.

Conclusion: The above results reveal, DBS can be utilized as a potential and robust sample source for bio banking in field epidemiological studies.

Genetic risk scores for major psychiatric disorders and the risk of postpartum psychiatric disorders

Postpartum psychiatric disorders are heritable, but how genetic liability varies by other significant risk factors is unknown. We aimed to (1) estimate associations of genetic risk scores (GRS) for major depression (MD), bipolar disorder (BD), and schizophrenia (SCZ) with postpartum psychiatric disorders, (2) examine differences by prior psychiatric history, and (3) compare genetic and familial risk of postpartum psychiatric disorders. We conducted a nested case-control study based on Danish population-based registers of all women in the iPSYCH2012 cohort who had given birth before December 31, 2015 (n = 8850). Cases were women with a diagnosed psychiatric disorder or a filled psychotropic prescription within one year after delivery (n = 5829 cases, 3021 controls). Association analyses were conducted between GRS calculated from Psychiatric Genomics Consortium discovery meta-analyses for MD, BD, and SCZ and case-control status of a postpartum psychiatric disorder. Parental psychiatric history was associated with postpartum psychiatric disorders among women with previous psychiatric history (OR, 1.14; 95% CI 1.02-1.28) but not without psychiatric history (OR, 1.08; 95% CI: 0.81-1.43). GRS for MD was associated with an increased risk of postpartum psychiatric disorders in both women with (OR, 1.44; 95% CI: 1.19-1.74) and without (OR, 1.88; 95% CI: 1.26-2.81) personal psychiatric history. SCZ GRS was only minimally associated with postpartum disorders and BD GRS was not. Results suggest GRS of lifetime psychiatric illness can be applied to the postpartum period, which may provide clues about distinct environmental or genetic elements of postpartum psychiatric disorders and ultimately help identify vulnerable groups.

Schizophrenia polygenic risk scores, urbanicity and treatment-resistant schizophrenia

INTRODUCTION

To investigate the impact of a polygenic risk score for schizophrenia (PRS-SZ) and urbanicity on the risk of treatment-resistant schizophrenia (TRS) in people diagnosed with schizophrenia and to evaluate the association between PRS-SZ and TRS across levels of urbanicity.

METHODS

Cohort study of people born after 1981 with a first registered diagnosis of schizophrenia between 1996 and 2012 using Danish population registry data. Through linkage to genome-wide data, we calculated PRS-SZ based on a Psychiatric Genomics Consortium meta-analysis. We assessed urbanicity at birth (capital, provincial and rural areas). TRS was defined using prescription and hospital data. Performing Cox regression analysis, we calculated hazard rate ratios (HRs) and 95% confidence intervals (CI).

RESULTS

Among 4475 people with schizophrenia, we identified 593 (13.3%) with TRS during 17,558 person years of follow-up. The adjusted HR for TRS associated with one standard deviation (SD) increase in the PRS-SZ was 1.11 (95% CI: 1.00-1.24). The adjusted HRs for TRS across levels of urbanicity were 1.20 (95% CI: 0.98-1.47) for provincial areas and 1.19 (95% CI 0.96-1.47) for rural areas compared with the capital area. Within strata of urbanicity, the adjusted HR for TRS was 1.39 (95% CI: 1.14-1.70) in the capital area with 1 SD increase in the PRS-SZ, 0.99 (95% CI 0.84-1.17) in provincial areas, and 1.03 (95% CI: 0.86-1.25) in rural areas.

CONCLUSION

The effect of genetic liability (i.e. PRS) on risk of TRS varied across urbanicity levels and was highest for people with schizophrenia born in the capital areas.

Optimization of extraction of genomic DNA from archived dried blood spot (DBS): potential application in epidemiological research & bio banking

Background:Limited infrastructure is available to collect, store and transport venous blood in field epidemiological studies. Dried blood spot (DBS) is a robust potential alternative sample source for epidemiological studies & bio banking. A stable source of genomic DNA (gDNA) is required for long term storage in bio bank for its downstream applications. Our objective is to optimize the methods of gDNA extraction from stored DBS and with the aim of revealing its utility in large scale epidemiological studies. Methods:The purpose of this study was to extract the maximum amount of gDNA from DBS on Whatman 903 protein saver card. gDNA was extracted through column  (Qiagen) & magnetic bead based (Invitrogen) methods. Quantification of extracted gDNA was performed with a spectrophotometer, fluorometer, and integrity analyzed by agarose gel electrophoresis. Result:Large variation was observed in quantity & purity (260/280 ratio, 1.8-2.9) of the extracted gDNA. The intact gDNA bands on the electrophoresis gel reflect the robustness of DBS for gDNA even after prolonged storage time. The extracted gDNA amount 2.16 – 24 ng/µl is sufficient for its PCR based downstream application, but unfortunately it can’t be used for whole genome sequencing or genotyping from extracted gDNA. Sequencing or genotyping can be achieved by after increasing template copy number through whole genome amplification of extracted gDNA. The obtained results create a base for future research to develop high-throughput research and extraction methods from blood samples.Conclusion:The above results reveal, DBS can be utilized as a potential and robust sample source for bio banking in field epidemiological studies.

Genetic Variants in the 9p21.3 Locus Associated with Glioma Risk in Children, Adolescents, and Young Adults: A Case-Control Study

BACKGROUND

Genome-wide association studies have identified germline genetic variants in 25 genetic loci that increase the risk of developing glioma in adulthood. It is not known if these variants increase the risk of developing glioma in children and adolescents and young adults (AYA). To date, no studies have performed genome-wide analyses to find novel genetic variants associated with glioma risk in children and AYA.

METHODS

We investigated the association between 8,831,628 genetic variants and risk of glioma in 854 patients diagnosed up to the age of 29 years and 3,689 controls from Sweden and Denmark. Recruitment of patients and controls was population based. Genotyping was performed using Illumina BeadChips, and untyped variants were imputed with IMPUTE2. We selected 41 established adult glioma risk variants for detailed investigation.

RESULTS

Three adult glioma risk variants, rs634537, rs2157719, and rs145929329, all mapping to the 9p21.3 (CDKN2B-AS1) locus, were associated with glioma risk in children and AYA. The strongest association was seen for rs634537 (odds ratioG = 1.21; 95% confidence interval = 1.09-1.35; P = 5.8 × 10-4). In genome-wide analysis, an association with risk was suggested for 129 genetic variants (P <1 × 10-5).

CONCLUSIONS

Carriers of risk alleles in the 9p21.3 locus have an increased risk of glioma throughout life. The results from genome-wide association analyses require validation in independent cohorts.

IMPACT

Our findings line up with existing evidence that some, although not all, established adult glioma risk variants are associated with risk of glioma in children and AYA. Validation of results from genome-wide analyses may reveal novel susceptibility loci for glioma in children and AYA.

Post-traumatic stress following military deployment: Genetic associations and cross-disorder genetic correlations

BACKGROUND

Post-traumatic stress disorder (PTSD) is a complex psychiatric disorder that occurs with relatively high frequency after deployment to warzones (∼10%). While twin studies have estimated the heritability to be up to 40%, thus indicating a considerable genetic component in the etiology, the biological mechanisms underlying risk and development of PTSD remain unknown.

METHODS

Here, we conduct a genome-wide association study (GWAS; N = 2,481) to identify genome regions that associate with PTSD in a highly homogenous, trauma-exposed sample of Danish soldiers deployed to war and conflict zones. We perform integrated analyses of our results with gene-expression and chromatin-contact datasets to prioritized genes. We also leverage on other large GWAS (N>300,000) to investigate genetic correlations between PTSD and other psychiatric disorders and traits.

RESULTS

We discover, but do not replicate, one region, 4q31, close to the IL15 gene, which is genome-wide significantly associated with PTSD. We demonstrate that gene-set enrichment, polygenic risk score and genetic correlation analyses show consistent and significant genetic correlations between PTSD and depression, insomnia and schizophrenia.

LIMITATIONS

The limited sample size, the lack of replication, and the PTSD case definition by questionnaire are limitations to the study.

CONCLUSIONS

Our results suggest that genetic perturbations of inflammatory response may contribute to the risk of PTSD. In addition, shared genetic components contribute to observed correlations between PTSD and depression, insomnia and schizophrenia.

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.

The Anorexia Nervosa Genetics Initiative (ANGI): Overview and methods

BACKGROUND

Genetic factors contribute to anorexia nervosa (AN); and the first genome-wide significant locus has been identified. We describe methods and procedures for the Anorexia Nervosa Genetics Initiative (ANGI), an international collaboration designed to rapidly recruit 13,000 individuals with AN and ancestrally matched controls. We present sample characteristics and the utility of an online eating disorder diagnostic questionnaire suitable for large-scale genetic and population research.

METHODS

ANGI recruited from the United States (US), Australia/New Zealand (ANZ), Sweden (SE), and Denmark (DK). Recruitment was via national registers (SE, DK); treatment centers (US, ANZ, SE, DK); and social and traditional media (US, ANZ, SE). All cases had a lifetime AN diagnosis based on DSM-IV or ICD-10 criteria (excluding amenorrhea). Recruited controls had no lifetime history of disordered eating behaviors. To assess the positive and negative predictive validity of the online eating disorder questionnaire (ED100K-v1), 109 women also completed the Structured Clinical Interview for DSM-IV (SCID), Module H.

RESULTS

Blood samples and clinical information were collected from 13,363 individuals with lifetime AN and from controls. Online diagnostic phenotyping was effective and efficient; the validity of the questionnaire was acceptable.

CONCLUSIONS

Our multi-pronged recruitment approach was highly effective for rapid recruitment and can be used as a model for efforts by other groups. High online presence of individuals with AN rendered the Internet/social media a remarkably effective recruitment tool in some countries. ANGI has substantially augmented Psychiatric Genomics Consortium AN sample collection. ANGI is a registered clinical trial: clinicaltrials.govNCT01916538; https://clinicaltrials.gov/ct2/show/NCT01916538?cond=Anorexia+Nervosa&draw=1&rank=3.

Metabolic profiles derived from residual blood spot samples: A longitudinal analysis

Background: Secondary use of newborn screening dried blood spot samples include use for biomedical or epidemiological research. However, the effects of storage conditions on archival samples requires further examination. The objective of this study was to determine the utility of residual newborn samples for deriving reliable metabolic gestational age estimates. Methods: Residual newborn dried blood spot samples that had been stored for 2-, 4-, 6-, or 12-months in temperature controlled (21°C) conditions were re-analyzed for the full panel of newborn screening analytes offered by a provincial newborn screening lab in Ottawa, Canada. Data from re-analyzed samples were compared to corresponding baseline newborn screening values for absolute agreement, and Pearson and intraclass correlation. Performance of a gestational age estimation algorithm originally developed from baseline newborn screening values was then validated on data derived from stored samples. Results: A total of 307 samples were used for this study. 17-hydroxyprogesterone and newborn hemoglobin profiles measured by immunoassay and high-performance liquid chromatography, respectively, were among the most stable markers across all time points of analysis. Acylcarnitines exhibited the greatest degree of variation in stability upon repeat measurement. The largest shifts in newborn analyte profiles and the poorest performance of metabolic gestational age algorithms were observed when samples were analyzed 12-months after sample collection. Conclusions: Duration of sample storage, independent of temperature and humidity, affects newborn screening profiles and gestational age estimates derived from metabolic gestational dating algorithms. When considering use of dried blood spot samples either for clinical or research purposes, care should be taken when interpreting data stemming from secondary use.

Genome-wide association study of Hirschsprung disease detects a novel low-frequency variant at the RET locus

Hirschsprung disease (HSCR) is a congenital disorder with a population incidence of ~1/5000 live births, defined by an absence of enteric ganglia along variable lengths of the colon. HSCR genome-wide association studies (GWAS) have found common associated variants at RET, SEMA3, and NRG1, but they still fail to explain all of its heritability. To enhance gene discovery, we performed a GWAS of 170 cases identified from the Danish nationwide pathology registry with 4717 controls, based on 6.2 million variants imputed from the haplotype reference consortium panel. We found a novel low-frequency variant (rs144432435), which, when conditioning on the lead RET single-nucleotide polymorphism (SNP), was of genome-wide significance in the discovery analysis. This conditional association signal was replicated in a Swedish HSCR cohort with discovery plus replication meta-analysis conditional odds ratio of 6.6 (P = 7.7 × 10^-10; 322 cases and 4893 controls). The conditional signal was, however, not replicated in two HSCR cohorts from USA and Finland, leading to the hypothesis that rs144432435 tags a rare haplotype present in Denmark and Sweden. Using the genome-wide complex trait analysis method, we estimated the SNP heritability of HSCR to be 88%, close to estimates based on classical family studies. Moreover, by using Lasso (least absolute shrinkage and selection operator) regression we were able to construct a genetic HSCR predictor with a area under the receiver operator characteristics curve of 76% in an independent validation set. In conclusion, we combined the largest collection of sporadic Hirschsprung cases to date (586 cases) to further elucidate HSCR's genetic architecture.

The iPSYCH2012 case-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.

COMT Val158Met and MTHFR C677T moderate risk of schizophrenia in response to childhood adversity

OBJECTIVE

Mesolimbic dopamine sensitization has been hypothesized to be a mediating factor of childhood adversity (CA) on schizophrenia risk. Activity of catechol-O-methyltransferase (COMT) Val158Met increases mesolimbic dopamine signaling and may be further regulated by methylenetetrahydrofolate reductase (MTHFR) C677T. This study investigates the three-way interaction between CA, COMT, and MTHFR.

METHODS

We conducted a nested case-control study on individuals born after 1981, linking population-based registers to study the three-way interaction. We included 1699 schizophrenia cases and 1681 controls, and used conditional logistic regression to report incidence rate ratios (IRRs).

RESULTS

Childhood adversity was robustly associated with schizophrenia. No main genetic effects were observed. MTHFR C677T increased schizophrenia risk in a dose-dependent manner per MTHFR T allele (P = 0.005) consequent upon CA exposure. After inclusion of the significant (P = 0.03) COMT × MTHFR × CA interaction, the risk was further increased per high-activity COMT Val allele. Hence, exposed COMT Val/Val and MTHFR T/T carriers had an IRR of 2.76 (95% CI, 1.66-4.61). Additional adjustments for ancestry and parental history of mental illness attenuated the results with the interaction being only marginally significant.

CONCLUSION

MTHFR C677T and COMT Val158Met interact with CA to increase risk of schizophrenia.

Evaluation of whole genome amplified DNA to decrease material expenditure and increase quality

Aim

The overall aim of this study is to evaluate whole genome amplification of DNA extracted from dried blood spot samples. We wish to explore ways of optimizing the amplification process, while decreasing the amount of input material and inherently the cost. Our primary focus of optimization is on the amount of input material, the amplification reaction volume, the number of replicates and amplification time and temperature. Increasing the quality of the amplified DNA and the subsequent results of array genotyping is a secondary aim of this project.

Methods

This study is based on DNA extracted from dried blood spot samples. The extracted DNA was subsequently whole genome amplified using the REPLIg kit and genotyped on the PsychArray BeadChip (assessing > 570,000 SNPs genome wide). We used Genome Studio to evaluate the quality of the genotype data by call rates and log R ratios.

Results

The whole genome amplification process is robust and does not vary between replicates. Altering amplification time, temperature or number of replicates did not affect our results. We found that spot size i.e. amount of input material could be reduced without compromising the quality of the array genotyping data. We also showed that whole genome amplification reaction volumes can be reduced by a factor of 4, without compromising the DNA quality.

Discussion

Whole genome amplified DNA samples from dried blood spots is well suited for array genotyping and produces robust and reliable genotype data. However, the amplification process introduces additional noise to the data, making detection of structural variants such as copy number variants difficult. With this study, we explore ways of optimizing the amplification protocol in order to reduce noise and increase data quality. We found, that the amplification process was very robust, and that changes in amplification time or temperature did not alter the genotyping calls or quality of the array data. Adding additional replicates of each sample also lead to insignificant changes in the array data. Thus, the amount of noise introduced by the amplification process was consistent regardless of changes made to the amplification protocol. We also explored ways of decreasing material expenditure by reducing the spot size or the amplification reaction volume. The reduction did not affect the quality of the genotyping data.

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.

Advantages and Challenges of Dried Blood Spot Analysis by Mass Spectrometry Across the Total Testing Process

INTRODUCTION

Through the introduction of advanced analytical techniques and improved throughput, the scope of dried blood spot testing utilising mass spectrometric methods, has broadly expanded. Clinicians and researchers have become very enthusiastic about the potential applications of dried blood spot based mass spectrometric applications. Analysts on the other hand face challenges of sensitivity, reproducibility and overall accuracy of dried blood spot quantification. In this review, we aim to bring together these two facets to discuss the advantages and current challenges of non-newborn screening applications of dried blood spot quantification by mass spectrometry.

METHODS

To address these aims we performed a key word search of the PubMed and MEDLINE online databases in conjunction with individual manual searches to gather information. Keywords for the initial search included; "blood spot" and "mass spectrometry"; while excluding "newborn"; and "neonate". In addition, databases were restricted to English language and human specific. There was no time period limit applied.

RESULTS

As a result of these selection criteria, 194 references were identified for review. For presentation, this information is divided into: 1) clinical applications; and 2) analytical considerations across the total testing process; being pre-analytical, analytical and post-analytical considerations.

CONCLUSIONS

DBS analysis using MS applications is now broadly applied, with drug monitoring for both therapeutic and toxicological analysis being the most extensively reported. Several parameters can affect the accuracy of DBS measurement and further bridge experiments are required to develop adjustment rules for comparability between dried blood spot measures and the equivalent serum/plasma values. Likewise, the establishment of independent reference intervals for dried blood spot sample matrix is required.

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.

Single Lysis-Salting Out Method of Genomic DNA Extraction From Dried Blood Spots

BACKGROUND

Dried blood spots (DBS) are an important form of bio-sampling and valuable approach for storing blood samples for genetic studies. This has necessitated in developing an effective protocol to isolate genomic DNA (gDNA) from DBS samples.In this study, we have elucidated a dependable and non-hazardous "single lysis-salting out" (SLSO) protocol of gDNA extraction from DBS and compared against the available commercial kits.

METHODS

For the purpose of this study, blood spots were collected on S&S 903 filter cards from 10 healthy volunteers and 30 patients with glutaric aciduria type I (GA-I). The gDNA was extracted from theseDBS samples by SLSO, QIAamp® gDNA Micro kit and innuPREP forensic kit methods. The quantity and quality of gDNA obtained from these methods were determined by measuring the absorbance using a Nanodrop spectrophotometer.

RESULTS

The SLSO method showed four-fold and eight-fold increased yield of gDNA in healthy volunteers and patient samples, respectively, compared to commercial kits (p<0.0001). The protocol was also found to be cost efficient, reducing the per sample cost to almost half. The suitability of this method for genetic studies was confirmed by performing R402W genotyping by RFLP in GA-I patients. The genotyping results showed the presence of R402W mutation in 20% (6/30) of patients.

CONCLUSION

The SLSO method was found to be inexpensive, non-hazardous and a suitable technique for isolating gDNA from DBS samples for genetic studies.

DNA Methylation at the Neonatal State and at the Time of Diagnosis: Preliminary Support for an Association with the Estrogen Receptor 1, Gamma-Aminobutyric Acid B Receptor 1, and Myelin Oligodendrocyte Glycoprotein in Female Adolescent Patients with OCD

Obsessive-compulsive disorder (OCD) is a neuropsychiatric disorder. Non-genetic factors and their interaction with genes have attracted increasing attention. Epigenetics is regarded an important interface between environmental signals and activation/repression of genomic responses. Epigenetic mechanisms have not previously been examined in OCD in children and adolescents. The aim of the present study was to examine the DNA methylation profile of selected genes in blood spots from neonates later diagnosed with OCD and in the same children/adolescents at the time of diagnosis compared with age- and sex-matched controls. Furthermore, we wanted to characterize the association of the differential methylation profiles with the severity of OCD and treatment outcome. Dried and new blood spot samples were obtained from 21 female children/adolescents with verified OCD and 12 female controls. The differential methylation was analyzed using a linear model and the correlation with the severity of OCD and treatment outcome was analyzed using the Pearson correlation. We evaluated selected Illumina Infinium HumanMethylation450 BeadChip probes within and up to 100,000 bp up- and downstream of 14 genes previously associated with OCD (SLC1A1, SLC25A12, GABBR1, GAD1, DLGAP1, MOG, BDNF, OLIG2, NTRK2 and 3, ESR1, SL6A4, TPH2, and COMT). The study found no significantly differential methylation. However, preliminary support for a difference was found for the gamma-aminobutyric acid (GABA) B receptor 1 (cg10234998, cg17099072) in blood samples at birth and for the estrogen receptor 1 (ESR1) (cg10939667), the myelin oligodendrocyte glycoprotein (MOG) (cg16650906), and the brain-derived neurotrophic factor (BDNF) (cg14080521) in blood samples at the time of diagnosis. Preliminary support for an association was observed between the methylation profiles of GABBR1 and MOG and baseline severity, treatment effect, and responder status; and between the methylation profile of ESR1 and baseline severity. To our knowledge, this is the first study to examine the DNA methylation profiles in OCD. The study points towards possible differences in the methylation profiles and suggests a correlation with the severity of OCD. However, the results warrant further studies in larger sample sets.

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.

Systematic Integration of Brain eQTL and GWAS Identifies ZNF323 as a Novel Schizophrenia Risk Gene and Suggests Recent Positive Selection Based on Compensatory Advantage on Pulmonary Function

Genome-wide association studies have identified multiple risk variants and loci that show robust association with schizophrenia. Nevertheless, it remains unclear how these variants confer risk to schizophrenia. In addition, the driving force that maintains the schizophrenia risk variants in human gene pool is poorly understood. To investigate whether expression-associated genetic variants contribute to schizophrenia susceptibility, we systematically integrated brain expression quantitative trait loci and genome-wide association data of schizophrenia using Sherlock, a Bayesian statistical framework. Our analyses identified ZNF323 as a schizophrenia risk gene (P = 2.22×10(-6)). Subsequent analyses confirmed the association of the ZNF323 and its expression-associated single nucleotide polymorphism rs1150711 in independent samples (gene-expression: P = 1.40×10(-6); single-marker meta-analysis in the combined discovery and replication sample comprising 44123 individuals: P = 6.85×10(-10)). We found that the ZNF323 was significantly downregulated in hippocampus and frontal cortex of schizophrenia patients (P = .0038 and P = .0233, respectively). Evidence for pleiotropic effects was detected (association of rs1150711 with lung function and gene expression of ZNF323 in lung: P = 6.62×10(-5) and P = 9.00×10(-5), respectively) with the risk allele (T allele) for schizophrenia acting as protective allele for lung function. Subsequent population genetics analyses suggest that the risk allele (T) of rs1150711 might have undergone recent positive selection in human population. Our findings suggest that the ZNF323 is a schizophrenia susceptibility gene whose expression may influence schizophrenia risk. Our study also illustrates a possible mechanism for maintaining schizophrenia risk variants in the human gene pool.

Investigating interactions between early life stress and two single nucleotide polymorphisms in HSD11B2 on the risk of schizophrenia

BACKGROUND

To examine the risk of schizophrenia in a Danish population after exposure to early life stress, and whether this risk is modified by DNA sequence variation, specifically two single nucleotide polymorphisms (SNPs) (rs5479 and rs56303414) from the gene HSD11B2. This gene encodes the enzyme 11-β hydroxysteroid dehydrogenase type 2 which converts active cortisol into inactive cortisone.

METHODS

A two-stage analysis involving (1) a population-based cohort study, and (2) a nested case-control study using genotype information. Stage 1 included 1,141,447 people; here, we calculated incidence rate ratios (IRR) for the risk of schizophrenia among children of mothers who experienced loss or serious illness of close relatives before, during, and after pregnancy. In stage 2, we genotyped rs5479 and rs56303414 among 1275 schizophrenia cases and 1367 controls, and investigated interactions between genotypes and early life stress on the risk of schizophrenia.

RESULTS

In stage 1, no increased risk of schizophrenia was found in offspring after exposure during pregnancy, but offspring exposed to early life stress at age 0-2 years had a significantly increased risk of schizophrenia (adjusted IRR 1.18, 95% confidence interval 1.07-1.31). For rs5479, the minor allele was nucleotide A, and the major allele was nucleotide C. No interaction was found between rs5479 and exposure during pregnancy. Individuals with the minor A allele of rs5479, however, had a significantly increased risk of schizophrenia after exposure to early life stress at age 3-9 years (adjusted IRR 2.06, 1.04-4.06). No interaction was found between rs56303414 and exposure in any of the time periods.

CONCLUSION

No association was found between exposure to early life stress during pregnancy and schizophrenia in the offspring investigated, whereas individuals exposed to early life stress within the first two years of life had an increased risk. No interaction was found between HSD11B2 and exposure during pregnancy, but individuals with the A allele of rs5479 had an increased risk of schizophrenia after exposure at age 3-9 years.

CCND2, CTNNB1, DDX3X, GLI2, SMARCA4, MYC, MYCN, PTCH1, TP53, and MLL2 gene variants and risk of childhood medulloblastoma

Recent studies have described a number of genes that are frequently altered in medulloblastoma tumors and that have putative key roles in the development of the disease. We hypothesized that common germline genetic variations in these genes may be associated with medulloblastoma development. Based on recent publications, we selected 10 genes that were frequently altered in medulloblastoma: CCND2, CTNNB1, DDX3X, GLI2, SMARCA4, MYC, MYCN, PTCH1, TP53, and MLL2 (now renamed as KMT2D). Common genetic variants (single nucleotide polymorphisms) annotating these genes (n = 221) were genotyped in germline DNA (neonatal dried blood spot samples) from 243 childhood medulloblastoma cases and 247 control subjects from Sweden and Denmark. Eight genetic variants annotating three genes in the sonic hedgehog signaling pathway; CCND2, PTCH1, and GLI2, were found to be associated with the risk of medulloblastoma (P(combined) < 0.05). The findings were however not statistically significant following correction for multiple testing by the very stringent Bonferroni method. The results do not support our hypothesis that common germline genetic variants in the ten studied genes are associated with the risk of developing medulloblastoma.

Whole Genome Amplification from Blood Spot Samples

Whole genome amplification is an invaluable technique when working with DNA extracted from blood spots, as the DNA obtained from this source often is too limited for extensive genetic analysis. Two techniques that amplify the entire genome are common. Here, both are described with focus on the benefits and drawbacks of each system. However, in order to obtain the best possible WGA result the quality of input DNA extracted from the blood spot is essential, but also time consumption, flexibility in format and elution volume and price of the technology are factors influencing system choice. Here, three DNA extraction techniques are described and the above aspects are compared between the systems.

UGT1A1*28 polymorphism and acute lymphoblastic leukemia in children: a Danish case-control study

BACKGROUND

Oxidative stress is a possible risk factor in the development of acute lymphoblastic leukemia (ALL) in children. Bilirubin is a potent endogenous antioxidant, and the UGT1A1*28 polymorphism is the main genetic cause of variation in plasma bilirubin in Western Europe.

METHODS

In a case-control study of 665 incident cases of ALL in childhood in Denmark 1982-2010 and 1,379 controls, associations between UGT1A1*28 genotypes and ALL in childhood were estimated as odds ratios by logistic regression with adjustment for sex and birth decade. Subgroup analyses were carried out by age at onset in three groups, and on the ALL subtypes precursor B-cell, T-cell, and t(12;21) positive status. Cases were identified in The Danish Registry of Childhood Cancer, and genotypes were estimated from dried blood spots stored in The Danish Neonatal Screening Biobank. Controls were newborns with blood spots taken right before and after a case.

RESULTS

We found no association between ALL in childhood and UGT1A1*28 genotypes. The odds ratio was 1.01 (0.88-1.17) for heterozygotes and 1.03 (0.78-1.36) for homozygotes. Also, no associations were found in the subgroup analyses.

CONCLUSION

We found no association between the UGT1A1*28 genotypes and ALL in children.

Common variants associated with general and MMR vaccine-related febrile seizures

Febrile seizures represent a recognized serious adverse event following measles, mumps, and rubella (MMR) vaccination. We conducted a series of genome-wide association scans comparing children with MMR-related febrile seizures, children with febrile seizures unrelated to vaccination, and controls with no history of febrile seizures. Two loci were distinctly associated with MMR-related febrile seizures, harboring the interferon-stimulated gene IFI44L (rs273259; P = 5.9×10⁻¹² vs. controls; P =1.2×10⁻⁹ vs. MMR-unrelated febrile seizures) and the measles virus receptor CD46 (rs1318653; P = 9.6×10⁻¹¹ vs. controls; P = 1.6×10⁻⁹ vs. MMR-unrelated febrile seizures). Furthermore, four loci were associated with febrile seizures in general implicating the sodium channel genes SCN1A (rs6432860; P = 2.2×10⁻¹⁶) and SCN2A (rs3769955; P = 3.1×10⁻¹⁰), a TMEM16 family gene (TMEM16C; rs114444506; P = 3.7×10⁻²⁰), and a region associated with magnesium levels (12q21.33; rs11105468; P = 3.4×10⁻¹¹). Finally, functional relevance of TMEM16C was demonstrated with electrophysiological experiments in wild-type and knockout rats.

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.

The Nordic medical birth registers--a potential goldmine for clinical research

The Nordic medical birth registers have long been used for valuable clinical research. Their collection of data for more than four decades offers unusual possibilities for research across generations. At the same time, serum and blotting paper blood samples have been stored from most neonates. Two large cohorts (approximately 100 000 births) in Denmark and Norway have been described by questionnaires, interviews and collection of biological samples (blood, urine and milk teeth), as well as a systematic prospective follow-up of the offspring. National patient registers provide information on preceding, underlying and present health problems of the parents and their offspring. Researchers may, with permission from the national authorities, obtain access to individualized or anonymized data from the registers and tissue-banks. These data allow for multivariate analyses but their usefulness depends on knowledge of the specific registers and biological sample banks and on proper validation of the registers.

Current and future use of "dried blood spot" analyses in clinical chemistry

The analysis of blood spotted and dried on a matrix (i.e., "dried blood spot" or DBS) has been used since the 1960s in clinical chemistry; mostly for neonatal screening. Since then, many clinical analytes, including nucleic acids, small molecules and lipids, have been successfully measured using DBS. Although this pre-analytical approach represents an interesting alternative to classical venous blood sampling, its routine use is limited. Here, we review the application of DBS technology in clinical chemistry, and evaluate its future role supported by new analytical methods such as mass spectrometry.

Comparison of manual methods of extracting genomic DNA from dried blood spots collected on different cards: implications for clinical practice

Isolating genomic DNA from blood samples is essential when studying the associations between genetic variants and susceptibility to a given clinical condition, or its severity. This study of three extraction techniques and two types of commercially available cards involved 219 children attending our outpatient pediatric clinic for follow-up laboratory tests after they had been hospitalised. An aliquot of venous blood was drawn into plastic tubes without additives and, after several inversions, 80 microL were put on circles of common paper cards and Whatman FTA-treated cards. Three extraction methods were compared: the Qiagen Investigator, Gensolve, and Masterpure. The best method in terms of final DNA yield was Masterpure, which led to a significantly higher yield regardless of the type of card (p less than 0.001), followed by Qiagen Investigator and Gensolve. Masterpure was also the best in terms of price, seemed to be simple and reliable, and required less hands-on time than other techniques. These conclusions support the use of Masterpure in studies that evaluate the associations between genetic variants and the severity or prevalence of infectious diseases.

Genome-wide study of association and interaction with maternal cytomegalovirus infection suggests new schizophrenia loci

Genetic and environmental components as well as their interaction contribute to the risk of schizophrenia, making it highly relevant to include environmental factors in genetic studies of schizophrenia. This study comprises genome-wide association (GWA) and follow-up analyses of all individuals born in Denmark since 1981 and diagnosed with schizophrenia as well as controls from the same birth cohort. Furthermore, we present the first genome-wide interaction survey of single nucleotide polymorphisms (SNPs) and maternal cytomegalovirus (CMV) infection. The GWA analysis included 888 cases and 882 controls, and the follow-up investigation of the top GWA results was performed in independent Danish (1396 cases and 1803 controls) and German-Dutch (1169 cases, 3714 controls) samples. The SNPs most strongly associated in the single-marker analysis of the combined Danish samples were rs4757144 in ARNTL (P=3.78 × 10(-6)) and rs8057927 in CDH13 (P=1.39 × 10(-5)). Both genes have previously been linked to schizophrenia or other psychiatric disorders. The strongest associated SNP in the combined analysis, including Danish and German-Dutch samples, was rs12922317 in RUNDC2A (P=9.04 × 10(-7)). A region-based analysis summarizing independent signals in segments of 100 kb identified a new region-based genome-wide significant locus overlapping the gene ZEB1 (P=7.0 × 10(-7)). This signal was replicated in the follow-up analysis (P=2.3 × 10(-2)). Significant interaction with maternal CMV infection was found for rs7902091 (P(SNP × CMV)=7.3 × 10(-7)) in CTNNA3, a gene not previously implicated in schizophrenia, stressing the importance of including environmental factors in genetic studies.

A genome-wide association study identifies CDHR3 as a susceptibility locus for early childhood asthma with severe exacerbations

Asthma exacerbations are among the most frequent causes of hospitalization during childhood, but the underlying mechanisms are poorly understood. We performed a genome-wide association study of a specific asthma phenotype characterized by recurrent, severe exacerbations occurring between 2 and 6 years of age in a total of 1,173 cases and 2,522 controls. Cases were identified from national health registries of hospitalization, and DNA was obtained from the Danish Neonatal Screening Biobank. We identified five loci with genome-wide significant association. Four of these, GSDMB, IL33, RAD50 and IL1RL1, were previously reported as asthma susceptibility loci, but the effect sizes for these loci in our cohort were considerably larger than in the previous genome-wide association studies of asthma. We also obtained strong evidence for a new susceptibility gene, CDHR3 (encoding cadherin-related family member 3), which is highly expressed in airway epithelium. These results demonstrate the strength of applying specific phenotyping in the search for asthma susceptibility genes.

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.

A comparative analysis of the governance and use of residual dried blood spots from state newborn screening programs and neonatal biobanks

In contrast to prior research, which has focused mainly on legislative content, this study aimed to characterize the current governance structures and procedures used by state newborn screening programs in the United States regarding the research use of residual dried blood spots. Content analysis was performed on newborn screening laws, program policies, survey responses, and online material, and was compared to information from two neonatal biobanks. Important differences between newborn screening programs and neonatal biobanks included the types of permissible research with blood spots, the evaluation criteria used for research requests, and characteristics of the research proposal reviewers. These findings can inform ongoing policy conversations with respect to the governance and use of residual dried blood spots.

High quality genome-wide genotyping from archived dried blood spots without DNA amplification

Spots of blood are routinely collected from newborn babies onto filter paper called Guthrie cards and used to screen for metabolic and genetic disorders. The archived dried blood spots are an important and precious resource for genomic research. Whole genome amplification of dried blood spot DNA has been used to provide DNA for genome-wide SNP genotyping. Here we describe a 96 well format procedure to extract DNA from a portion of a dried blood spot that provides sufficient unamplified genomic DNA for genome-wide single nucleotide polymorphism (SNP) genotyping. We show that SNP genotyping of the unamplified DNA is more robust than genotyping amplified dried blood spot DNA, is comparable in cost, and can be done with thousands of samples. This procedure can be used for genome-wide association studies and other large-scale genomic analyses that require robust, high-accuracy genotyping of dried blood spot DNA.