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.

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.

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.

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.

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.

Mitochondrial haplogroups in patients with rheumatoid arthritis: No association with disease and disease manifestations.

OBJECTIVE: To describe the distribution of specific mitochondrial DNA (mtDNA) haplogroups (hgs) in a cohort of patients with rheumatoid arthritis (RA). METHODS: Two-hundred nineteen consecutive patients with RA had mtDNA isolated from their blood, sequenced and haplotyped. Patients were diagnosed according to the American College of Rheumatology (ACR)/European league against Rheumatism (EULAR) criteria. Demographic and clinical data were retrieved from the Danish nationwide database (DANBIO). Logistic regression analyses were performed to test for associations. RESULTS: One-hundred eighty-four patients were eligible for analysis. Haplogroup frequencies were: H (n = 88; 47.8%), U (n = 37; 20.1%), T (n = 22; 12.0%), J (n = 16; 8.7%), K (n = 11; 5.9%), HV (n = 6; 3.3%) and V (n = 4; 2.2%). The distribution of individual hgs was identical to the background population. Radiographic erosions were significantly associated with hg clusters JT (OR = 2.37, 95% confidence interval (CI): 1.07-5.53, p = 0.038). Significantly fewer patients from hg cluster JT received biological treatment (OR = 0.17, 95% CI: 0.03-0.87, p = 0.038). Albeit, none of these associations were significant when corrected for multiple tests. CONCLUSION: There was no significant association between mtDNA hgs and presence of RA or disease manifestations. There was an, albeit insignificant, overrepresentation of patients with hg JT among patients with erosive disease; however, slightly fewer patients in the JT group were treated with biological drugs.

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.

Pharmacokinetics of Isoniazid, Pyrazinamide, and Ethambutol in Newly Diagnosed Pulmonary TB Patients in Tanzania.

Exposure to lower-than-therapeutic levels of anti-tuberculosis drugs is likely to cause selection of resistant strains of Mycobacterium tuberculosis and treatment failure. The first-line anti-tuberculosis (TB) regimen consists of rifampicin, isoniazid, pyrazinamide, and ethambutol, and correct management reduces risk of TB relapse and development of drug resistance. In this study we aimed to investigate the effect of standard of care plus nutritional supplementation versus standard care on the pharmacokinetics of isoniazid, pyrazinamide and ethambutol among sputum smear positive TB patients with and without HIV. In a clinical trial in 100 Tanzanian TB patients, with or without HIV infection, drug concentrations were determined at 1 week and 2 months post initiation of anti-TB medication. Data was analysed using population pharmacokinetic modelling. The effect of body size was described using allometric scaling, and the effects of nutritional supplementation, HIV, age, sex, CD4+ count, weight-adjusted dose, NAT2 genotype, and time on TB treatment were investigated. The kinetics of all drugs was well characterised using first-order elimination and transit compartment absorption, with isoniazid and ethambutol described by two-compartment disposition models, and pyrazinamide by a one-compartment model. Patients with a slow NAT2 genotype had higher isoniazid exposure and a lower estimate of oral clearance (15.5 L/h) than rapid/intermediate NAT2 genotype (26.1 L/h). Pyrazinamide clearance had an estimated typical value of 3.32 L/h, and it was found to increase with time on treatment, with a 16.3% increase after the first 2 months of anti-TB treatment. The typical clearance of ethambutol was estimated to be 40.7 L/h, and was found to decrease with age, at a rate of 1.41% per year. Neither HIV status nor nutritional supplementations were found to affect the pharmacokinetics of these drugs in our cohort of patients.

Longitudinal changes in C-reactive protein, proform of eosinophil major basic protein, and pregnancy-associated plasma protein-A during weight changes in obese children.

BACKGROUND: Childhood obesity is associated with several complications, including cardiovascular comorbidity. Several biomarkers, such as high-sensitive C-reactive protein (hs-CRP), proform of eosinophil major basic protein (Pro-MBP) and pregnancy associated plasma protein-A (PAPP-A), have equally been linked to increased cardiovascular susceptibility. This study investigates these biomarkers during weight loss and regain in obese children. MATERIALS AND METHODS: A longitudinal study during a 12-week weight loss program with a 28 months follow-up was conducted. Anthropometrics and plasma concentrations of hs-CRP, Pro-MBP, and PAPP-A were measured at baseline; at days 14, 33 and 82 during weight loss; and at months 10, 16, and 28 during follow-up. RESULTS: Fifty-three boys and 62 girls aged 8-15 years with a median body mass index (BMI) standard deviation score (SDS) at baseline of 2.78 (boys), and 2.70 (girls) were included. Ninety children completed the weight loss program and 68 children entered the follow-up program. Pro-MBP and PAPP-A, but not hs-CRP, exhibited individual-specific levels (tracking) during weight loss and regain. The PAPP-A/Pro-MBP correlation was strong, whereas the hs-CRP/PAPP-A correlation was weak during weight fluctuations. CONCLUSION: Hs-CRP changes reflect weight changes. PAPP-A and Pro-MBP exhibited tracking during weight perturbations and may contribute as early risk markers of cardiovascular susceptibility.

Nutritional supplementation increases rifampin exposure among tuberculosis patients coinfected with HIV.

Nutritional supplementation to tuberculosis (TB) patients has been associated with increased weight and reduced mortality, but its effect on the pharmacokinetics of first-line anti-TB drugs is unknown. A cohort of 100 TB patients (58 men; median age, 35 [interquartile range {IQR}, 29 to 40] years, and median body mass index [BMI], 18.8 [17.3 to 19.9] kg/m(2)) were randomized to receive nutritional supplementation during the intensive phase of TB treatment. Rifampin plasma concentrations were determined after 1 week and 2 months of treatment. The effects of nutritional supplementation, HIV, time on treatment, body weight, and SLCO1B1 rs4149032 genotype were examined using a population pharmacokinetic model. The model adjusted for body size via allometric scaling, accounted for clearance autoinduction, and detected an increase in bioavailability (+14%) for the patients in the continuation phase. HIV coinfection in patients not receiving the supplementation was found to decrease bioavailability by 21.8%, with a median maximum concentration of drug in serum (Cmax) and area under the concentration-time curve from 0 to 24 h (AUC0-24) of 5.6 µg/ml and 28.6 µg · h/ml, respectively. HIV-coinfected patients on nutritional supplementation achieved higher Cmax and AUC0-24 values of 6.4 µg/ml and 31.6 µg · h/ml, respectively, and only 13.3% bioavailability reduction. No effect of the SLCO1B1 rs4149032 genotype was observed. In conclusion, nutritional supplementation during the first 2 months of TB treatment reduces the decrease in rifampin exposure observed in HIV-coinfected patients but does not affect exposure in HIV-uninfected patients. If confirmed in other studies, the use of defined nutritional supplementation in HIV-coinfected TB patients should be considered in TB control programs. (This study has the controlled trial registration number ISRCTN 16552219.).

Dysfunctional mitochondrial respiration in the striatum of the Huntington's disease transgenic R6/2 mouse model.

Metabolic dysfunction and mitochondrial involvement are recognised as part of the pathology in Huntington's Disease (HD). Post-mortem examinations of the striatum from end-stage HD patients have shown a decrease in the in vitro activity of complexes II, III and IV of the electron transport system (ETS). In different models of HD, evidence of enzyme defects have been reported in complex II and complex IV using enzyme assays. However, such assays are highly variable and results have been inconsistent. We investigated the integrated ETS function ex vivo using a sensitive high-resolution respirometric (HRR) method. The O2 flux in a whole-cell sample combined with the addition of mitochondrial substrates, uncouplers and inhibitors enabled us to accurately quantitate the function of individual mitochondrial complexes in intact mitochondria, while retaining mitochondrial regulation and compensatory mechanisms. We used HRR to examine the mitochondrial function in striata from 12-week old R6/2 mice expressing exon 1 of human HTT with 130 CAG repeats. A significant reduction in complex II and complex IV flux control ratios was found in the R6/2 mouse striatum at 12 weeks of age compared to controls, confirming previous findings obtained with spectrophotometric enzyme assays.