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.

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.

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.

Investigation of the involvement of MIR185 and its target genes in the development of schizophrenia.

BACKGROUND: Schizophrenia is a complex neuropsychiatric disorder of unclear etiology. The strongest known genetic risk factor is the 22q11.2 microdeletion. Research has yet to confirm which genes within the deletion region are implicated in schizophrenia. The minimal 1.5 megabase deletion contains MIR185, which encodes microRNA 185. METHODS: We determined miR-185 expression in embryonic and adult mouse brains. Common and rare variants at this locus were then investigated using a human genetics approach. First, we performed gene-based analyses for MIR185 common variants and target genes using Psychiatric Genomics Consortium genome-wide association data. Second, MIR185 was resequenced in German patients (n = 1000) and controls (n = 500). We followed up promising variants by genotyping an additional European sample (patients, n = 3598; controls, n = 4082). RESULTS: In situ hybridization in mice revealed miR-185 expression in brain regions implicated in schizophrenia. Gene-based tests revealed association between common variants in 3 MIR185 target genes (ATAT1, SH3PXD2A, NTRK3) and schizophrenia. Further analyses in mice revealed overlapping expression patterns for these target genes and miR-185. Resequencing identified 2 rare patient-specific novel variants flanking MIR185. However, follow-up genotyping provided no further evidence of their involvement in schizophrenia. LIMITATIONS: Power to detect rare variant associations was limited. CONCLUSION: Human genetic analyses generated no evidence of the involvement of MIR185 in schizophrenia. However, the expression patterns of miR-185 and its target genes in mice, and the genetic association results for the 3 target genes, suggest that further research into the involvement of miR-185 and its downstream pathways in schizophrenia is warranted.

Common variants at VRK2 and TCF4 conferring risk of schizophrenia.

Common sequence variants have recently joined rare structural polymorphisms as genetic factors with strong evidence for association with schizophrenia. Here we extend our previous genome-wide association study and meta-analysis (totalling 7 946 cases and 19 036 controls) by examining an expanded set of variants using an enlarged follow-up sample (up to 10 260 cases and 23 500 controls). In addition to previously reported alleles in the major histocompatibility complex region, near neurogranin (NRGN) and in an intron of transcription factor 4 (TCF4), we find two novel variants showing genome-wide significant association: rs2312147[C], upstream of vaccinia-related kinase 2 (VRK2) [odds ratio (OR) = 1.09, P = 1.9 × 10(-9)] and rs4309482[A], between coiled-coiled domain containing 68 (CCDC68) and TCF4, about 400 kb from the previously described risk allele, but not accounted for by its association (OR = 1.09, P = 7.8 × 10(-9)).

Whole-genome amplified DNA from stored dried blood spots is reliable in high resolution melting curve and sequencing analysis.

BACKGROUND: The use of dried blood spots (DBS) samples in genomic workup has been limited by the relative low amounts of genomic DNA (gDNA) they contain. It remains to be proven that whole genome amplified DNA (wgaDNA) from stored DBS samples, constitutes a reliable alternative to gDNA.We wanted to compare melting curves and sequencing results from wgaDNA derived from DBS samples with gDNA derived from whole blood. METHODS: gDNA was extracted from whole blood obtained from 10 patients with lone atrial fibrillation (mean age 22.3 years). From their newborn DBS samples, stored at -24°C, genomic DNA was extracted and whole-genome amplified in triplicates. Using high resolution melting curve analysis and direct sequencing in both wgaDNA and gDNA samples, all coding regions and adjacent intron regions of the genes SCN5A and KCNA5 were investigated. RESULTS: Altered melting curves was present in 85 of wgaDNA samples and 81 of gDNA samples. Sequence analysis identified a total of 31 variants in the 10 wgaDNA samples. The same 31 variants were found in the exact same pattern of samples in the gDNA group. There was no false positive or negative sequence variation in the wgaDNA group. CONCLUSIONS: The use of DNA amplified in triplicates from DBS samples is reliable and can be used both for high resolution curve melting analysis as well as direct sequence analysis. DBS samples therefore can serve as an alternative to whole blood in sequence analysis.

Robustness of genome-wide scanning using archived dried blood spot samples as a DNA source.

BACKGROUND: The search to identify disease-susceptible genes requires access to biological material from numerous well-characterized subjects. Archived residual dried blood spot (DBS) samples, also known as Guthrie cards, from national newborn screening programs may provide a DNA source for entire populations. Combined with clinical information from medical registries, DBS samples could provide a rich source for productive research. However, the amounts of DNA which can be extracted from these precious samples are minute and may be prohibitive for numerous genotypings. Previously, we demonstrated that DBS DNA can be whole-genome amplified and used for reliable genetic analysis on different platforms, including genome-wide scanning arrays. However, it remains unclear whether this approach is workable on a large sample scale. We examined the robustness of using DBS samples for whole-genome amplification following genome-wide scanning, using arrays from Illumina and Affymetrix. RESULTS: This study is based on 4,641 DBS samples from the Danish Newborn Screening Biobank, extracted for three separate genome-wide association studies. The amount of amplified DNA was significantly (P < 0.05) affected by the year of storage and storage conditions. Nine (0.2%) DBS samples failed whole-genome amplification. A total of 4,586 (98.8%) samples met our criterion of success of a genetic call-rate above 97%. The three studies used different arrays, with mean genotyping call-rates of 99.385% (Illumina Infinium Human610-Quad), 99.722% (Illumina Infinium HD HumanOmni1-Quad), and 99.206% (Affymetrix Axiom Genome-Wide CEU). We observed a concordance rate of 99.997% in the 38 methodological replications, and 99.999% in the 27 technical replications. Handling variables such as time of storage, storage conditions and type of filter paper were shown too significantly (P < 0.05) affect the genotype call-rates in some of the arrays, although the effect was minimal. CONCLUSION: Our study indicates that archived DBS samples from the Danish Newborn Screening Biobank represent a reliable resource of DNA for whole-genome amplification and subsequent genome-wide association studies. With call-rates equivalent to high quality DNA samples, our results point to new opportunities for using the neonatal biobanks available worldwide in the hunt for genetic components of disease.

Genome-wide scans using archived neonatal dried blood spot samples.

BACKGROUND: Identification of disease susceptible genes requires access to DNA from numerous well-characterised subjects. Archived residual dried blood spot samples from national newborn screening programs may provide DNA from entire populations and medical registries the corresponding clinical information. The amount of DNA available in these samples is however rarely sufficient for reliable genome-wide scans, and whole-genome amplification may thus be necessary. This study assess the quality of DNA obtained from different amplification protocols by evaluating fidelity and robustness of the genotyping of 610,000 single nucleotide polymorphisms, using the Illumina Infinium HD Human610-Quad BeadChip. Whole-genome amplified DNA from 24 neonatal dried blood spot samples stored between 15 to 25 years was tested, and high-quality genomic DNA from 8 of the same individuals was used as reference. RESULTS: Using 3.2 mm disks from dried blood spot samples the optimal DNA-extraction and amplification protocol resulted in call-rates between 99.15% - 99.73% (mean 99.56%, N = 16), and conflicts with reference DNA in only three per 10,000 genotype calls. CONCLUSION: Whole-genome amplified DNA from archived neonatal dried blood spot samples can be used for reliable genome-wide scans and is a cost-efficient alternative to collecting new samples.

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.

Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder.

Attention deficit/hyperactivity disorder (ADHD) is a highly heritable childhood behavioral disorder affecting 5% of children and 2.5% of adults. Common genetic variants contribute substantially to ADHD susceptibility, but no variants have been robustly associated with ADHD. We report a genome-wide association meta-analysis of 20,183 individuals diagnosed with ADHD and 35,191 controls that identifies variants surpassing genome-wide significance in 12 independent loci, finding important new information about the underlying biology of ADHD. Associations are enriched in evolutionarily constrained genomic regions and loss-of-function intolerant genes and around brain-expressed regulatory marks. Analyses of three replication studies: a cohort of individuals diagnosed with ADHD, a self-reported ADHD sample and a meta-analysis of quantitative measures of ADHD symptoms in the population, support these findings while highlighting study-specific differences on genetic overlap with educational attainment. Strong concordance with GWAS of quantitative population measures of ADHD symptoms supports that clinical diagnosis of ADHD is an extreme expression of continuous heritable traits.

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.

Genetic Variation in NFKBIE Is Associated With Increased Risk of Pneumococcal Meningitis in Children.

BACKGROUND: Streptococcus pneumoniae and Neisseria meningitidis are frequent pathogens in life-threatening infections. Genetic variation in the immune system may predispose to these infections. Nuclear factor-κB is a key component of the TLR-pathway, controlled by inhibitors, encoded by the genes NFKBIA, NFKBIE and NFKBIZ. We aimed to replicate previous findings of genetic variation associated with invasive pneumococcal disease (IPD), and to assess whether similar associations could be found in invasive meningococcal disease (IMD). METHODS: Cases with IPD and IMD and controls were identified by linking Danish national registries. DNA was obtained from the Danish Neonatal Screening Biobank. The association between SNPs and susceptibility to IPD and IMD, mortality and pneumococcal serotypes was investigated. RESULTS: 372 children with pneumococcal meningitis, 907 with pneumococcal bacteremia and 1273 controls were included. We included 406 cases with meningococcal meningitis, 272 with meningococcal bacteremia, and 672 controls. The NFKBIE SNP was associated with increased risk of pneumococcal meningitis (aOR 1.68; 95% CI: 1.20-2.36), but not bacteremia (aOR 1.08; 95% CI: 0.86-1.35). The remaining SNPs were not associated with susceptibility to invasive disease. None of the SNPs were associated with risk of IMD or mortality. CONCLUSIONS: A NFKBIE polymorphism was associated with increased risk of pneumococcal meningitis.

Genetic risk for autism spectrum disorders and neuropsychiatric variation in the general population.

Almost all genetic risk factors for autism spectrum disorders (ASDs) can be found in the general population, but the effects of this risk are unclear in people not ascertained for neuropsychiatric symptoms. Using several large ASD consortium and population-based resources (total n > 38,000), we find genome-wide genetic links between ASDs and typical variation in social behavior and adaptive functioning. This finding is evidenced through both LD score correlation and de novo variant analysis, indicating that multiple types of genetic risk for ASDs influence a continuum of behavioral and developmental traits, the severe tail of which can result in diagnosis with an ASD or other neuropsychiatric disorder. A continuum model should inform the design and interpretation of studies of neuropsychiatric disease biology.

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

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

Mannose-Binding Lectin Gene, MBL2, Polymorphisms Do Not Increase Susceptibility to Invasive Meningococcal Disease in a Population of Danish Children.

Background. Neisseria meningitidis is the cause of meningococcal bacteremia and meningitis, and nasopharyngeal colonization with this pathogen is common. The incidence of invasive disease is highest in infants, whereas adolescents more often are carriers. Altered regulation or dysfunction of the innate immune system may predispose to invasive meningococcal disease (IMD). In this study, we investigated the effect of genetic variation in the mannose-binding lectin gene, MBL2, and its promoter on susceptibility to IMD and IMD-associated mortality among children. Methods. Children (<5 years) diagnosed during 1982-2007 with IMD and controls were identified through Danish national registries. DNA was obtained from the Danish Neonatal Screening Biobank. The associations between MBL2 diplotypes and IMD susceptibility and 30- and 90-day mortality were investigated using logistic regression analysis. Results. We included 1351 children: 406 with meningitis, 272 with bacteremia, and 673 age- and sex-matched controls. Of the children studied, 1292 (96%) were successfully genotyped and assigned MBL2 diplotypes. The median age in IMD cases was 19.1 months (interquartile range [IQR], 8.8-32.2 months). Children with defective MBL2 diplotypes were not at higher risk for meningococcal meningitis than children with intermediate and normal diplotypes (odds ratio [OR] = 0.69; 95% confidence interval [CI], .47-1.02). Similar results were found for children with bacteremia and defective diplotypes (OR = 0.84; 95% CI, .53-1.32) as well as for all cases (OR = 0.75; 95% CI, .56-1.01). There was no association between MBL2 diplotypes and mortality. Conclusions. Defective MBL2 diplotypes did not predict either an increased IMD susceptibility or mortality in a Danish population of children.

Polygenic Risk Score, Parental Socioeconomic Status, Family History of Psychiatric Disorders, and the Risk for Schizophrenia: A Danish Population-Based Study and Meta-analysis.

IMPORTANCE: Schizophrenia has a complex etiology influenced both by genetic and nongenetic factors but disentangling these factors is difficult. OBJECTIVE: To estimate (1) how strongly the risk for schizophrenia relates to the mutual effect of the polygenic risk score, parental socioeconomic status, and family history of psychiatric disorders; (2) the fraction of cases that could be prevented if no one was exposed to these factors; (3) whether family background interacts with an individual's genetic liability so that specific subgroups are particularly risk prone; and (4) to what extent a proband's genetic makeup mediates the risk associated with familial background. DESIGN, SETTINGS, AND PARTICIPANTS: We conducted a nested case-control study based on Danish population-based registers. The study consisted of 866 patients diagnosed as having schizophrenia between January 1, 1994, and December 31, 2006, and 871 matched control individuals. Genome-wide data and family psychiatric and socioeconomic background information were obtained from neonatal biobanks and national registers. Results from a separate meta-analysis (34,600 cases and 45,968 control individuals) were applied to calculate polygenic risk scores. EXPOSURES: Polygenic risk scores, parental socioeconomic status, and family psychiatric history. MAIN OUTCOMES AND MEASURES: Odds ratios (ORs), attributable risks, liability R2 values, and proportions mediated. RESULTS: Schizophrenia was associated with the polygenic risk score (OR, 8.01; 95% CI, 4.53-14.16 for highest vs lowest decile), socioeconomic status (OR, 8.10; 95% CI, 3.24-20.3 for 6 vs no exposures), and a history of schizophrenia/psychoses (OR, 4.18; 95% CI, 2.57-6.79). The R2 values were 3.4% (95% CI, 2.1-4.6) for the polygenic risk score, 3.1% (95% CI, 1.9-4.3) for parental socioeconomic status, and 3.4% (95% CI, 2.1-4.6) for family history. Socioeconomic status and psychiatric history accounted for 45.8% (95% CI, 36.1-55.5) and 25.8% (95% CI, 21.2-30.5) of cases, respectively. There was an interaction between the polygenic risk score and family history (P = .03). A total of 17.4% (95% CI, 9.1-26.6) of the effect associated with family history of schizophrenia/psychoses was mediated through the polygenic risk score. CONCLUSIONS AND RELEVANCE: Schizophrenia was associated with the polygenic risk score, family psychiatric history, and socioeconomic status. Our study demonstrated that family history of schizophrenia/psychoses is partly mediated through the individual's genetic liability.

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

Febrile seizures represent a 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(-12) versus controls, P = 1.2 × 10(-9) versus MMR-unrelated febrile seizures) and the measles virus receptor CD46 (rs1318653: P = 9.6 × 10(-11) versus controls, P = 1.6 × 10(-9) versus 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(-16)) and SCN2A (rs3769955: P = 3.1 × 10(-10)), a TMEM16 family gene (ANO3; rs114444506: P = 3.7 × 10(-20)) and a region associated with magnesium levels (12q21.33; rs11105468: P = 3.4 × 10(-11)). Finally, we show the functional relevance of ANO3 (TMEM16C) with electrophysiological experiments in wild-type and knockout rats.

Genome-wide association analyses identify variants in developmental genes associated with hypospadias.

Hypospadias is a common congenital condition in boys in which the urethra opens on the underside of the penis. We performed a genome-wide association study on 1,006 surgery-confirmed hypospadias cases and 5,486 controls from Denmark. After replication genotyping of an additional 1,972 cases and 1,812 controls from Denmark, the Netherlands and Sweden, 18 genomic regions showed independent association with P < 5 × 10(-8). Together, these loci explain 9% of the liability to developing this condition. Several of the identified regions harbor genes with key roles in embryonic development (including HOXA4, IRX5, IRX6 and EYA1). Subsequent pathway analysis with GRAIL and DEPICT provided additional insight into possible genetic mechanisms causing hypospadias.

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.