Combinatorial batching of DNA for ultralow-cost detection of pathogenic variants.

BACKGROUND: Next-generation sequencing (NGS) based population screening holds great promise for disease prevention and earlier diagnosis, but the costs associated with screening millions of humans remain prohibitive. New methods for population genetic testing that lower the costs of NGS without compromising diagnostic power are needed. METHODS: We developed double batched sequencing where DNA samples are batch-sequenced twice - directly pinpointing individuals with rare variants. We sequenced batches of at-birth blood spot DNA using a commercial 113-gene panel in an explorative (n = 100) and a validation (n = 100) cohort of children who went on to develop pediatric cancers. All results were benchmarked against individual whole genome sequencing data. RESULTS: We demonstrated fully replicable detection of cancer-causing germline variants, with positive and negative predictive values of 100% (95% CI, 0.91-1.00 and 95% CI, 0.98-1.00, respectively). Pathogenic and clinically actionable variants were detected in RB1, TP53, BRCA2, APC, and 19 other genes. Analyses of larger batches indicated that our approach is highly scalable, yielding more than 95% cost reduction or less than 3 cents per gene screened for rare disease-causing mutations. We also show that double batched sequencing could cost-effectively prevent childhood cancer deaths through broad genomic testing. CONCLUSIONS: Our ultracheap genetic diagnostic method, which uses existing sequencing hardware and standard newborn blood spots, should readily open up opportunities for population-wide risk stratification using genetic screening across many fields of clinical genetics and genomics.

Genetic liability to posttraumatic stress disorder and its association with postpartum depression.

BACKGROUND: Childbirth may be a traumatic experience and vulnerability to posttraumatic stress disorder (PTSD) may increase the risk of postpartum depression (PPD). We investigated whether genetic vulnerability to PTSD as measured by polygenic score (PGS) increases the risk of PPD and whether a predisposition to PTSD in PPD cases exceeds that of major depressive disorder (MDD) outside the postpartum period. METHODS: This case-control study included participants from the iPSYCH2015, a case-cohort of all singletons born in Denmark between 1981 and 2008. Restricting to women born between 1981 and 1997 and excluding women with a first diagnosis other than depression (N = 22 613), 333 were identified with PPD. For each PPD case, 999 representing the background population and 993 with MDD outside the postpartum were matched by calendar year at birth, cohort selection, and age. PTSD PGS was calculated from summary statistics from the Psychiatric Genomics Consortium with LDpred2-auto. Odds ratios (ORs) were estimated using conditional logistic regression adjusted for parental psychiatric history and country of origin, PGS for MDD and age at first birth, and the first 10 principal components. RESULTS: The PTSD PGS was significantly associated with PPD (OR 1.42, 95% CI 1.20-1.68 per standard deviation increase in PTSD PGS) compared to healthy female controls. Genetic PTSD vulnerability in PPD cases did not exceed that of matched female depression cases outside the postpartum period (OR 1.10, 95% CI 0.94-1.30 per standard deviation increase). CONCLUSIONS: Genetic vulnerability to PTSD increased the risk of PPD but did not differ between PPD cases and women with depression at other times.

Neonatal Screening for Congenital Adrenal Hyperplasia in Denmark: 10 Years of Experience.

INTRODUCTION: Early detection of salt-wasting congenital adrenal hyperplasia (SW-CAH) is important to reduce CAH-related morbidity. However, neonatal screening has shown to have a low positive predictive value (PPV), especially among preterm newborns. Here, the Danish CAH screening is evaluated by comparing incidence and morbidity of SW-CAH 10 years before and after introduction of screening. Furthermore, sensitivity, specificity, and PPV are determined. METHODS: All newborns in Denmark born during 1999-2018 and diagnosed with SW-CAH were identified in the Danish National Patient Registry and/or at the Department of Clinical Genetics, Rigshospitalet. Newborns with a positive neonatal CAH screening were identified at Statens Serum Institut. Correct diagnosis was evaluated by medical record review. RESULTS: A total of 65 newborns with SW-CAH were identified. The incidence of SW-CAH was 5:100,000 both before and after introduction of screening. Performance of sensitivity and specificity of the screening were 97% and 100%, respectively, and the PPV was 55% for the given period. Stratified according to gestational age, the PPV was 33% and 61% for pre -and fullterm newborns, respectively. Though not significant, the proportion of newborns presenting with SW-crisis decreased after introduction of screening from 29% versus 10% (p = 0.07). DISCUSSION AND CONCLUSION: Neonatal screening for SW-CAH has not led to an increase in the incidence of newborns diagnosed with SW-CAH. The screening algorithm has effectively identified newborns with SW-CAH. After 2009, there was a tendency toward a lower proportion of newborns with SW-crisis at diagnosis. Finally, the study emphasizes the benefits of using second-tier screening as well as repeated screening of premature newborns.

Implementation of SCID Screening in Denmark.

Screening for SCID was added to the Danish Neonatal Screening Program in February 2020. The screening uses a RealtimePCR kit and we here present the results and experiences with the validation of the kit and the first 10 months of screening.

Co-occurring hydrocephalus in autism spectrum disorder: a Danish population-based cohort study.

BACKGROUND: The association between autism spectrum disorder and hydrocephalus is not well understood, despite demonstrated links between autism spectrum disorder and cerebrospinal fluid abnormalities. Based on the hypothesis that autism spectrum disorder and hydrocephalus may, at least in some cases, be two manifestations of a shared congenital brain pathology, we investigated the potential association between autism spectrum disorder and hydrocephalus in a large Danish population-based cohort. METHODS: Patients and controls were obtained from the Lundbeck Foundation Initiative for Integrative Psychiatric Research iPSYCH2012 case-cohort, which includes all patients with selected psychiatric disorders born in Denmark 1981-2005 along with randomly selected population controls (end of follow-up, December 31, 2016). The associations between individual psychiatric disorders and hydrocephalus were estimated using binary logistic regression with adjustment for age and sex. RESULTS: The cohort consisted of 86,571 individuals, of which 14,654 were diagnosed with autism spectrum disorder, 28,606 were population controls, and the remaining were diagnosed with other psychiatric disorders. We identified 201 hydrocephalus cases; 68 among autism spectrum disorder patients and 40 among controls (OR 3.77, 95% CI 2.48-5.78), which corresponds to an absolute risk of 0.46 % (i.e. approximately one in 217 children with autism spectrum disorder had co-occurring hydrocephalus). The autism spectrum disorder-hydrocephalus association was significant over the entire subgroup spectrum of autism spectrum disorder. CONCLUSIONS: Given the considerable risk of hydrocephalus among patients with autism spectrum disorder, we suggest that patients with autism spectrum disorder should be evaluated for co-occurring hydrocephalus on a routine basis as timely neurosurgical intervention is important. Likewise, attention must be paid to traits of autism spectrum disorder in children with hydrocephalus. The results of this study call for future investigations on a potential shared aetiology between hydrocephalus and autism spectrum disorder, including the role abnormal CSF dynamics in the pathogenesis of autism spectrum disorder.

Identification of genetic loci associated with nocturnal enuresis: a genome-wide association study.

BACKGROUND: Nocturnal enuresis (bedwetting) is a common disorder affecting 10-16% of 7-year-old children globally. Nocturnal enuresis is highly heritable, but its genetic determinants remain unknown. We aimed to identify genetic variants associated with nocturnal enuresis and explore its genetic architecture and underlying biology. METHODS: We did a genome-wide association study (GWAS) of nocturnal enuresis. Nocturnal enuresis cases were identified in iPSYCH2012, a large Danish population-based case cohort established to investigate mental disorders, on the basis of 10th revision of the International Statistical Classification of Diseases (ICD-10) diagnoses and redeemed desmopressin prescriptions in Danish registers. The GWAS was done in a genetically homogeneous sample of unrelated individuals using logistic regression with relevant covariates. All genome-wide significant variants were analysed for their association with nocturnal enuresis in an independent Icelandic sample from deCODE genetics. Standardised polygenic risk scores for attention-deficit hyperactivity disorder (ADHD) and autism spectrum disorder were constructed from summary statistics of large GWASs and analysed for association with nocturnal enuresis. FINDINGS: The GWAS included 3882 nocturnal enuresis cases and 31 073 controls. We found two loci at chromosome 6 and chromosome 13 significantly associated with nocturnal enuresis. Six genetic variants at the two loci (five variants at chromosome 6q16.2 and one variant at chromosome 13q22.3) surpassed the threshold for genome-wide significance (p<5 × 10-8). There were two lead variants: rs9376454 (chromosome 6q16.2), with an odds ratio (OR) of 1·199 (95% CI 1·135-1·267; p=9·91 × 10-11), and rs60721117 (chromosome 13q22.3), with an OR of 1·149 (1·095-1·205; p=1·21 × 10-8). All associated variants in the chromosome 6 locus were replicated (p<8 × 10-3) in the independent Icelandic cohort of 5475 nocturnal enuresis cases and 303 996 controls, whereas the associated variant in the chromosome 13 locus showed nominal significant association (p=0·031). The percentage of nocturnal enuresis phenotypic variance explained by the common genetic variants was 23·9-30·4%. Polygenic risk for ADHD was associated with nocturnal enuresis (OR 1·06, 95% CI, 1·01-1·10; p=0·011). Among the potential nocturnal enuresis risk genes mapped, PRDM13 and EDNRB have biological functions associated with known pathophysiological mechanisms in nocturnal enuresis, and SIM1 regulates the formation of the hypothalamic neuroendocrine lineage that produces arginine vasopressin, a well known nocturnal enuresis drug target. INTERPRETATION: This study shows that common genetic variants contribute considerably to nocturnal enuresis, and it identifies potential nocturnal enuresis risk genes with roles in sleep, urine production, and bladder function. Given that available treatments target these mechanisms, any of the identified genes and their functional gene networks are potential drug targets. FUNDING: The Lundbeck Foundation Initiative for Integrative Psychiatric Research (iPSYCH), Stanley Foundation.

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.

Genetic predisposition to hypertension is associated with preeclampsia in European and Central Asian women.

Preeclampsia is a serious complication of pregnancy, affecting both maternal and fetal health. In genome-wide association meta-analysis of European and Central Asian mothers, we identify sequence variants that associate with preeclampsia in the maternal genome at ZNF831/20q13 and FTO/16q12. These are previously established variants for blood pressure (BP) and the FTO variant has also been associated with body mass index (BMI). Further analysis of BP variants establishes that variants at MECOM/3q26, FGF5/4q21 and SH2B3/12q24 also associate with preeclampsia through the maternal genome. We further show that a polygenic risk score for hypertension associates with preeclampsia. However, comparison with gestational hypertension indicates that additional factors modify the risk of preeclampsia.

Association Between Childhood Green Space, Genetic Liability, and the Incidence of Schizophrenia.

Childhood exposure to green space has previously been associated with lower risk of developing schizophrenia later in life. It is unclear whether this association is mediated by genetic liability or whether the 2 risk factors work additively. Here, we investigate possible gene-environment associations with the hazard ratio (HR) of schizophrenia by combining (1) an estimate of childhood exposure to residential-level green space based on the normalized difference vegetation index (NDVI) from Landsat satellite images, with (2) genetic liability estimates based on polygenic risk scores for 19 746 genotyped individuals from the Danish iPSYCH sample. We used information from the Danish registers of health, residential address, and socioeconomic status to adjust HR estimates for established confounders, ie, parents' socioeconomic status, and family history of mental illness. The adjusted HRs show that growing up surrounded by the highest compared to the lowest decile of NDVI was associated with a 0.52-fold (95% confidence interval [CI]: 0.40 to 0.66) lower schizophrenia risk, and children with the highest polygenic risk score had a 1.24-fold (95% CI: 1.18 to 1.30) higher schizophrenia risk. We found that NDVI explained 1.45% (95% CI: 1.07 to 1.90) of the variance on the liability scale, while polygenic risk score for schizophrenia explained 1.01% (95% CI: 0.77 to 1.46). Together they explained 2.40% (95% CI: 1.99 to 3.07) with no indication of a gene-environment interaction (P = .29). Our results suggest that risk of schizophrenia is associated additively with green space exposure and genetic liability, and provide no support for an environment-gene interaction between NDVI and schizophrenia.

Association of Childhood Exposure to Nitrogen Dioxide and Polygenic Risk Score for Schizophrenia With the Risk of Developing Schizophrenia.

Importance: Schizophrenia is a highly heritable psychiatric disorder, and recent studies have suggested that exposure to nitrogen dioxide (NO2) during childhood is associated with an elevated risk of subsequently developing schizophrenia. However, it is not known whether the increased risk associated with NO2 exposure is owing to a greater genetic liability among those exposed to highest NO2 levels. Objective: To examine the associations between childhood NO2 exposure and genetic liability for schizophrenia (as measured by a polygenic risk score), and risk of developing schizophrenia. Design, Setting, and Participants: Population-based cohort study including individuals with schizophrenia (International Statistical Classification of Diseases and Related Health Problems, Tenth Revision code F20) and a randomly selected subcohort. Using national registry data, all individuals born in Denmark between May 1, 1981, and December 31, 2002, were followed up from their 10th birthday until the first occurrence of schizophrenia, emigration, death, or December 31, 2012, whichever came first. Statistical analyses were conducted between October 24, 2018, and June 17, 2019. Exposures: Individual exposure to NO2 during childhood estimated as mean daily exposure to NO2 at residential addresses from birth to the 10th birthday. Polygenic risk scores were calculated as the weighted sum of risk alleles at selected single-nucleotide polymorphisms based on genetic material obtained from dried blood spot samples from the Danish Newborn Screening Biobank and on the Psychiatric Genomics Consortium genome-wide association study summary statistics file. Main Outcomes and Measures: The main outcome was schizophrenia. Weighted Cox proportional hazards regression models were fitted to estimate adjusted hazard ratios (AHRs) for schizophrenia with 95% CIs according to the exposures. Results: Of a total of 23 355 individuals, 11 976 (51.3%) were male and all had Danish-born parents. During the period of the study, 3531 were diagnosed with schizophrenia. Higher polygenic risk scores were correlated with higher childhood NO2 exposure (ρ = 0.0782; 95% CI, 0.065-0.091; P < .001). A 10-µg/m3 increase in childhood daily NO2 exposure (AHR, 1.23; 95% CI, 1.15-1.32) and a 1-SD increase in polygenic risk score (AHR, 1.29; 95% CI, 1.23-1.35) were independently associated with increased schizophrenia risk. Conclusions and Relevance: These findings suggest that the apparent association between NO2 exposure and schizophrenia is only slightly confounded by a higher polygenic risk score for schizophrenia among individuals living in areas with greater NO2. The findings demonstrate the utility of including polygenic risk scores in epidemiologic studies.

Reduced neonatal brain-derived neurotrophic factor is associated with autism spectrum disorders.

Mental disorders have for the majority of cases an unknown etiology, but several studies indicate that neurodevelopmental changes happen in utero or early after birth. We performed a nested case-control study of the relation between blood levels of neuro-developmental (S100B, BDNF, and VEGF-A) and inflammatory (MCP-1, TARC, IL-8, IL-18, CRP, and IgA) biomarkers in newborns, and later development of autism spectrum disorders (ASD, N = 751), attention deficit hyperactivity disorders (ADHD, N = 801), schizophrenia (N = 1969), affective (N = 641) or bipolar disorders (N = 641). Samples and controls were obtained as part of the iPSYCH Danish Case-Cohort Study using dried blood spot samples collected between 1981 and 2004, and stored frozen at the Danish National Biobank. In newborns lower blood level of BDNF was significantly associated with increased odds (OR 1.15) of developing ASD (p = 0.001). This difference could not be explained by genetic variation in the BDNF coding gene region. A tendency of decreased levels of all the neurotrophic markers and increased levels of all inflammatory markers was noted. The low newborn blood levels of BDNF in children developing ASD is an important finding, suggesting that lower BDNF levels in newborns contributes to the etiology of ASD and indicates new directions for further research. It may also help identifying a long-sought marker for high-ASD risk in, e.g., younger siblings of ASD children.

Polygenic risk for psychiatric disorder and singleness in patients with severe mental illness and controls.

We aimed to investigate whether the polygenic risk score (PRS) for schizophrenia influences time in couple relationships for patients with severe mental illness and controls. We combined the nationwide Danish registers with genetic information from dried neonatal blood spots. We included 2,599 individuals with schizophrenia, 1,446 with bipolar disorder, 20,315 with depression, and 6,963 controls. PRS for schizophrenia, depression, and bipolar disorder were estimated using data from the Psychiatric Genetics Consortium and analyzed both as a scale-predictor and as highest versus other deciles. The main outcome was number of days in couple relationships. Patients with schizophrenia had markedly fewer days/year in couple relationships: 64 (95% CI; 61-69) than patients with depression: 119 (95% CI; 117-121), bipolar disorder: 103 (95% CI 97-110), and controls: 136 (95% CI 133-139). PRS for schizophrenia was associated with fewer days in couple relationships in patients with schizophrenia (scale-PRS: IRR = 0.95 (0.93-0.97)) or depression (highest decile: IRR = 0.93 (0.87-0.98)). PRS for bipolar disorder (as scale) was also associated with fewer days in couple relationships in patients with depression (IRR = 0.99 (0.99-1.00)) or bipolar disorder (IRR = 0.96 (0.94-0.99)) and controls (IRR = 0.99 (0.97-1.00), and IRR = 0.89 (0.81-0.98) for the highest decile). Due to the number of statistical tests, however, it cannot be concluded definitely that some of these may not be spurious findings. In conclusion, our findings implicate high genetic loading for schizophrenia as a predisposing factor to singleness in patients with schizophrenia or depression, and genetic loading for bipolar disorder a similar predisposing factor in patients with depression, bipolar disorder or controls.

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

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

Exploring Cuba's population structure and demographic history using genome-wide data.

Cuba is the most populated country in the Caribbean and has a rich and heterogeneous genetic heritage. Here, we take advantage of dense genomic data from 860 Cuban individuals to reconstruct the genetic structure and ancestral origins of this population. We found distinct admixture patterns between and within the Cuban provinces. Eastern provinces have higher African and Native American ancestry contributions (average 26% and 10%, respectively) than the rest of the Cuban provinces (average 17% and 5%, respectively). Furthermore, in the Eastern Cuban region, we identified more intense sex-specific admixture patterns, strongly biased towards European male and African/Native American female ancestries. Our subcontinental ancestry analyses in Cuba highlight the Iberian population as the best proxy European source population, South American and Mesoamerican populations as the closest Native American ancestral component, and populations from West Central and Central Africa as the best proxy sources of the African ancestral component. Finally, we found complex admixture processes involving two migration pulses from both Native American and African sources. Most of the inferred Native American admixture events happened early during the Cuban colonial period, whereas the African admixture took place during the slave trade and more recently as a probable result of large-scale migrations from Haiti.

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

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

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.

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.

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.

Alcohol-related breast cancer in postmenopausal women - effect of CYP19A1, PPARG and PPARGC1A polymorphisms on female sex-hormone levels and interaction with alcohol consumption and NSAID usage in a nested case-control study and a randomised controlled trial.

BACKGROUND: Alcohol consumption is associated with increased risk of breast cancer (BC), and the underlying mechanism is thought to be sex-hormone driven. In vitro and observational studies suggest a mechanism involving peroxisome proliferator-activated receptor gamma (PPARγ) in a complex with peroxisome proliferator-activated receptor gamma coactivator 1-α (PGC-1α) and interaction with aromatase (encoded by CYP19A1). Use of non-steroidal anti-inflammatory drugs (NSAID) may also affect circulating sex-hormone levels by modifying PPARγ activity. METHODS: In the present study we assessed whether genetic variation in CYP19A1 is associated with risk of BC in a case-control study group nested within the Danish "Diet, Cancer and Health" cohort (ncases = 687 and ncontrols = 687) and searched for gene-gene interaction between CYP19A1 and PPARGC1A, and CYP19A1 and PPARG, and gene-alcohol and gene-NSAID interactions. Association between the CYP19A1 polymorphisms and hormone levels was also examined among 339 non-HRT users. Incidence rate ratios were calculated based on Cox' proportional hazards model. Furthermore, we performed a pilot randomised controlled trial to determine the effect of the PPARG Pro(12)Ala polymorphism and the PPARγ stimulator Ibuprofen on sex-hormone levels following alcohol intake in postmenopausal women (n = 25) using linear regression. RESULTS: Genetic variations in CYP19A1 were associated with hormone levels (estrone: P rs11070844 = 0.009, estrone sulphate: P rs11070844 = 0.01, P rs749292 = 0.004, P rs1062033 = 0.007 and P rs10519297 = 0.03, and sex hormone-binding globulin (SHBG): P rs3751591 = 0.03) and interacted with alcohol intake in relation to hormone levels (estrone sulphate: P interaction/rs2008691 = 0.02 and P interaction/rs1062033= 0.03, and SHBG: P interaction/rs11070844 = 0.03). CYP19A1/rs3751591 was both associated with SHBG levels (P = 0.03) and with risk of BC (Incidence Rate Ratio = 2.12; 95 % Confidence Interval: 1.02-4.43) such that homozygous variant allele carriers had increased levels of serum SHBG and were at increased risk of BC. Acute intake of alcohol decreased blood estrone (P = <0.0001), estrone sulphate (P = <0.0001), and SHBG (P = 0.009) levels, whereas Ibuprofen intake and PPARG Pro(12)Ala genotype had no effect on hormone levels. CONCLUSIONS: Our results suggest that genetically determined variation in CYP19A1 is associated with differences in sex hormone levels. However, the genetically determined differences in sex hormone levels were not convincingly associated with BC risk. The results therefore indicate that the genetically determined variation in CYP19A1 contributes little to BC risk and to alcohol-mediated BC risk. TRIAL REGISTRATION: NCT02463383, June 3, 2015.