Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Investigating Genetic Determinants of Plasma Inositol Status in Adult Humans.

BACKGROUND: Myo-inositol (MI) is incorporated into numerous biomolecules, including phosphoinositides and inositol phosphates. Disturbance of inositol availability or metabolism is associated with various disorders, including neurological conditions and cancers, whereas supplemental MI has therapeutic potential in conditions such as depression, polycystic ovary syndrome, and congenital anomalies. Inositol status can be influenced by diet, synthesis, transport, utilization, and catabolism. OBJECTIVES: We aimed to investigate potential genetic regulation of circulating MI status and to evaluate correlation of MI concentration with other metabolites. METHODS: GC-MS was used to determine plasma MI concentration of >2000 healthy, young adults (aged 18-28 y) from the Trinity Student Study. Genotyping data were used to test association of plasma MI with single nucleotide polymorphisms (SNPs) in candidate genes, encoding inositol transporters and synthesizing enzymes, and test for genome-wide association. We evaluated potential correlation of plasma MI with d-chiro-inositol (DCI), glucose, and other metabolites by Spearman rank correlation. RESULTS: Mean plasma MI showed a small but significant difference between males and females (28.5 and 26.9 µM, respectively). Candidate gene analysis revealed several nominally significant associations with plasma MI, most notably for SLC5A11 (solute carrier family 5 member 11), encoding a sodium-coupled inositol transporter, also known as SMIT2 (sodium-dependent myo-inositol transporter 2). However, these did not survive correction for multiple testing. Subsequent testing for genome-wide association with plasma MI did not identify associations of genome-wide significance (P < 5 × 10-8). However, 8 SNPs exceeded the threshold for suggestive significant association with plasma MI concentration (P < 1 × 10-5), 3 of which were located within or close to genes: MTDH (metadherin), LAPTM4B (lysosomal protein transmembrane 4 ß), and ZP2 (zona pellucida 2). We found significant positive correlation of plasma MI concentration with concentration of dci and several other biochemicals including glucose, methionine, betaine, sarcosine, and tryptophan. CONCLUSIONS: Our findings suggest potential for modulation of plasma MI in young adults by variation in SLC5A11, which is worthy of further investigation.

Probing the functional consequence and clinical relevance of CD320 p.E88del, a variant in the transcobalamin receptor gene.

The biological and clinical significance of the p.E88del variant in the transcobalamin receptor, CD320, is unknown. This allele is annotated in ClinVar as likely benign, pathogenic, and of uncertain significance. To determine functional consequence and clinical relevance of this allele, we employed cell culture and genetic association studies. Fibroblasts from 16 CD320 p.E88del homozygotes exhibited reduced binding and uptake of cobalamin. Complete ascertainment of newborns with transiently elevated C3 (propionylcarnitine) in New York State demonstrated that homozygosity for CD320 p.E88del was over-represented (7/348, p < 6 × 10-5 ). Using population data, we estimate that ~85% of the p.E88del homozygotes born in the same period did not have elevated C3, suggesting that cobalamin metabolism in the majority of these infants with this genotype is unaffected. Clinical follow-up of 4/9 homozygous individuals uncovered neuropsychological findings, mostly in speech and language development. None of these nine individuals exhibited perturbation of cobalamin metabolism beyond the newborn stage even during periods of acute illness. Newborns homozygous for this allele in the absence of other factors are at low risk of requiring clinical intervention, although more studies are required to clarify the natural history of various CD320 variants across patient populations.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia.

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Exome sequencing of child-parent trios with bladder exstrophy: Findings in 26 children.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

Mice lacking the transcobalamin-vitamin B12 receptor, CD320, suffer from anemia and reproductive deficits when fed vitamin B12-deficient diet.

In humans, poor nutrition, malabsorption and variation in cobalamin (vitamin B12) metabolic genes are associated with hematological, neurological and developmental pathologies. Cobalamin is transported from blood into tissues via the transcobalamin (TC) receptor encoded by the CD320 gene. We created mice carrying a targeted deletion of the mouse ortholog, Cd320. Knockout (KO) mice lacking this TC receptor have elevated levels of plasma methylmalonic acid and homocysteine but are otherwise healthy, viable, fertile and not anemic. To challenge the Cd320 KO mice we maintained them on a vitamin B12-deficient diet. After 5 weeks on this diet, reproductive failure develops in Cd320 KO females but not males. In vitro, homozygous Cd320 KO embryos from cobalamin-deficient Cd320 KO dams develop normally to embryonic day (E) 3.5, while in vivo, few uterine decidual implantation sites are observed at E7.5, suggesting that embryos perish around the time of implantation. Dietary restriction of vitamin B12 induces a severe macrocytic anemia in Cd320 KO mice after 10-12 months while control mice on this diet are anemia-free up to 2 years. Despite the severe anemia, cobalamin-deficient KO mice do not exhibit obvious neurological symptoms. Our results with Cd320 KO mice suggest that an alternative mechanism exists for mice to transport cobalamin independent of the Cd320 encoded receptor. Our findings with deficient diet are consistent with historical and epidemiological data suggesting that low vitamin B12 levels in humans are associated with infertility and developmental abnormalities. Our Cd320 KO mouse model is an ideal model system for studying vitamin B12 deficiency.

The FUT2 secretor variant p.Trp154Ter influences serum vitamin B12 concentration via holo-haptocorrin, but not holo-transcobalamin, and is associated with haptocorrin glycosylation.

Vitamin B12 deficiency is common in older individuals. Circulating vitamin B12 concentration can be used to diagnose deficiency, but this test has substantial false positive and false negative rates. We conducted genome-wide association studies (GWAS) in which we resolved total serum vitamin B12 into the fractions bound to transcobalamin and haptocorrin: two carrier proteins with very different biological properties. We replicated reported associations between total circulating vitamin B12 concentrations and a common null variant in FUT2. This allele determines the secretor phenotype in which blood group antigens are found in non-blood body fluids. Vitamin B12 bound to haptocorrin (holoHC) remained highly associated with FUT2 rs601338 (p.Trp154Ter). Transcobalamin bound vitamin B12 (holoTC) was not influenced by this variant. HoloTC is the bioactive the form of the vitamin and is taken up by all tissues. In contrast, holoHC is only taken up by the liver. Using holoHC from individuals with known FUT2 genotypes, we demonstrated that FUT2 rs601338 genotype influences the glycosylation of haptocorrin. We then developed an experimental model demonstrating that holoHC is transported into cultured hepatic cells (HepG2) via the asialoglycoprotein receptor (ASGR). Our data challenge current published hypotheses on the influence of genetic variation on this clinically important measure and are consistent with a model in which FUT2 rs601338 influences holoHC by altering haptocorrin glycosylation, whereas B12 bound to non-glycosylated transcobalamin (i.e. holoTC) is not affected. Our findings explain some of the observed disparity between use of total B12 or holoTC as first-line clinical tests of vitamin B12 status.

A Common Polymorphism in HIBCH Influences Methylmalonic Acid Concentrations in Blood Independently of Cobalamin.

Methylmalonic acid (MMA) is a by-product of propionic acid metabolism through the vitamin B12 (cobalamin)-dependent enzyme methylmalonyl CoA mutase. Elevated MMA concentrations are a hallmark of several inborn errors of metabolism and indicators of cobalamin deficiency in older persons. In a genome-wide analysis of 2,210 healthy young Irish adults (median age 22 years) we identified a strong association of plasma MMA with SNPs in 3-hydroxyisobutyryl-CoA hydrolase (HIBCH, p = 8.42 × 10(-89)) and acyl-CoA synthetase family member 3 (ACSF3, p = 3.48 × 10(-19)). These loci accounted for 12% of the variance in MMA concentration. The most strongly associated SNP (HIBCH rs291466; c:2T>C) causes a missense change of the initiator methionine codon (minor-allele frequency = 0.43) to threonine. Surprisingly, the resulting variant, p.Met1?, is associated with increased expression of HIBCH mRNA and encoded protein. These homozygotes had, on average, 46% higher MMA concentrations than methionine-encoding homozygotes in young adults with generally low MMA concentrations (0.17 [0.14-0.21] µmol/L; median [25(th)-75(th) quartile]). The association between MMA levels and HIBCH rs291466 was highly significant in a replication cohort of 1,481 older individuals (median age 79 years) with elevated plasma MMA concentrations (0.34 [0.24-0.51] µmol/L; p = 4.0 × 10(-26)). In a longitudinal study of 185 pregnant women and their newborns, the association of this SNP remained significant across the gestational trimesters and in newborns. HIBCH is unique to valine catabolism. Studies evaluating flux through the valine catabolic pathway in humans should account for these variants. Furthermore, this SNP could help resolve equivocal clinical tests where plasma MMA values have been used to diagnose cobalamin deficiency.

Genetic Risk Factors for Folate-Responsive Neural Tube Defects.

Neural tube defects (NTDs) are the most severe congenital malformations of the central nervous system. The etiology is complex, with both genetic and environmental factors having important contributions. Researchers have known for the past two decades that maternal periconceptional use of the B vitamin folic acid can prevent many NTDs. Though this finding is arguably one of the most important recent discoveries in birth defect research, the mechanism by which folic acid exerts this benefit remains unknown. Research to date has focused on the hypothesis that an underlying genetic susceptibility interacts with folate-sensitive metabolic processes at the time of neural tube closure. Little progress has been made searching for risk-causative variants in candidate genes; therefore, more complex genetic and epigenetic methodologies are now being considered. This article reviews the research to date that has been targeted on this important gene-nutrient locus.

Association of a transcobalamin II genetic variant with falsely low results for the holotranscobalamin immunoassay.

BACKGROUND: The clinical use of holotranscobalamin (holoTC) testing to evaluate vitamin B12 status has increased in recent years. We present two patients (African Caribbean and Indian heritage), in which the holoTC assay indicated severe B12 deficiency (< 5 pmol/L). Additional clinical tests revealed that these patients had normal levels of total vitamin B12 in blood and unremarkable levels of two other markers of vitamin B12 status, homocysteine and methylmalonic acid. We hypothesized that these patients carry a variant in the transcobalamin gene (TCN2) that influences the most widely commercially available holoTC test - Active-B12 (Axis-Shield Diagnostics Ltd). DESIGN: Exon sequencing of the TCN2 gene was carried out. Protein characterization included total transcobalamin (TCN2) detection by Western blot, and holoTC by (57) Co-labelled B12 binding followed by size fractionation. RESULTS: Exon sequencing of TCN2 revealed both patients were homozygous for the minor allele of rs35838082 (p.R215W). Western blot and chromatographic analyses revealed that the serum of these patients contains intact transcobalamin and that this variant-containing protein binds vitamin B12 . The variant is rare in Caucasians (minor allele frequency (MAF) < 0·01) but more common in South Asians (MAF ~ 0·02) and those of African origin (MAF ~ 0·25). CONCLUSIONS: The impeded ability to detect normal levels of holoTC in these two patients may be due to this variant interfering with the detection of holoTC by one or both of the monoclonal antibodies currently employed in the Active-B12 test. Laboratories should be aware of this variant and use confirmatory tests when applicable.

Evaluation of proton-coupled folate transporter (SLC46A1) polymorphisms as risk factors for neural tube defects and oral clefts.

Many folate-related genes have been investigated for possible causal roles in neural tube defects (NTDs) and oral clefts. However, no previous reports have examined the major gene responsible for folate uptake, the proton-coupled folate transporter (SLC46A1). We tested for association between these birth defects and single nucleotide polymorphisms in the SLC46A1 gene. The NTD study population included 549 complete and incomplete case-family triads, and 999 controls from Ireland. The oral clefts study population comprised a sample from Utah (495 complete and incomplete case-family triads and 551 controls) and 221 Filipino multiplex cleft families. There was suggestive evidence of increased NTD case risk with the rs17719944 minor allele (odds ratio (OR): 1.29; 95% confidence intervals (CI): [1.00-1.67]), and decreased maternal risk of an NTD pregnancy with the rs4795436 minor allele (OR: 0.62; [0.39-0.99]). In the Utah sample, the rs739439 minor allele was associated with decreased case risk for cleft lip with cleft palate (genotype relative risk (GRR): 0.56 [0.32-0.98]). Additionally, the rs2239907 minor allele was associated with decreased case risk for cleft lip with cleft palate in several models, and with cleft palate only in a recessive model (OR: 0.41; [0.20-0.85]). These associations did not remain statistically significant after correcting for multiple hypothesis testing. Nominal associations between SLC46A1 polymorphisms and both Irish NTDs and oral clefts in the Utah population suggest some role in the etiology of these birth defects, but further investigation in other populations is needed.

Common Variants at Putative Regulatory Sites of the Tissue Nonspecific Alkaline Phosphatase Gene Influence Circulating Pyridoxal 5'-Phosphate Concentration in Healthy Adults.

BACKGROUND: Vitamin B-6 interconversion enzymes are important for supplying pyridoxal 5'-phosphate (PLP), the co-enzyme form, to tissues. Variants in the genes for these enzymes [tissue nonspecific alkaline phosphatase (ALPL), pyridoxamine 5'-phosphate oxidase, pyridoxal kinase, and pyridoxal phosphatase] could affect enzyme function and vitamin B-6 status. OBJECTIVES: We tested whether single-nucleotide polymorphisms (SNPs) in these genes influence vitamin B-6 status markers [plasma PLP, pyridoxal (PL), and 4-pyridoxic acid (PA)], and explored potential functional effects of the SNPs. METHODS: Study subjects were young, healthy adults from Ireland (n = 2345). We measured plasma PLP, PL, and PA with liquid chromatography-tandem mass spectrometry and genotyped 66 tag SNPs in the 4 genes. We tested for associations with single SNPs in candidate genes and also performed genome-wide association study (GWAS) and gene-based analyses. RESULTS: Seventeen SNPs in ALPL were associated with altered plasma PLP in candidate gene analyses (P < 1.89 × 10(-4)). In the GWAS, 5 additional ALPL SNPs were associated with altered plasma PLP (P < 5.0 × 10(-8)). Gene-based analyses that used the functional linear model ß-spline (P = 4.04 × 10(-15)) and Fourier spline (P = 5.87 × 10(-15)) methods also showed associations between ALPL and altered plasma PLP. No SNPs in other genes were associated with plasma PLP. The association of the minor CC genotype of 1 ALPL SNP, rs1256341, with reduced ALPL expression in the HapMap Northern European ancestry population is consistent with the positive association between the CC genotype and plasma PLP in our study (P = 0.008). No SNP was associated with altered plasma PL or PA. CONCLUSIONS: In healthy adults, common variants in ALPL influence plasma PLP concentration, the most frequently used biomarker for vitamin B-6 status. Whether these associations are indicative of functional changes in vitamin B-6 status requires more investigation.

Replication and exploratory analysis of 24 candidate risk polymorphisms for neural tube defects.

BACKGROUND: Neural tube defects (NTDs), which are among the most common congenital malformations, are influenced by environmental and genetic factors. Low maternal folate is the strongest known contributing factor, making variants in genes in the folate metabolic pathway attractive candidates for NTD risk. Multiple studies have identified nominally significant allelic associations with NTDs. We tested whether associations detected in a large Irish cohort could be replicated in an independent population. METHODS: Replication tests of 24 nominally significant NTD associations were performed in racially/ethnically matched populations. Family-based tests of fifteen nominally significant single nucleotide polymorphisms (SNPs) were repeated in a cohort of NTD trios (530 cases and their parents) from the United Kingdom, and case-control tests of nine nominally significant SNPs were repeated in a cohort (190 cases, 941 controls) from New York State (NYS). Secondary hypotheses involved evaluating the latter set of nine SNPs for NTD association using alternate case-control models and NTD groupings in white, African American and Hispanic cohorts from NYS. RESULTS: Of the 24 SNPs tested for replication, ADA rs452159 and MTR rs10925260 were significantly associated with isolated NTDs. Of the secondary tests performed, ARID1A rs11247593 was associated with NTDs in whites, and ALDH1A2 rs7169289 was associated with isolated NTDs in African Americans. CONCLUSIONS: We report a number of associations between SNP genotypes and neural tube defects. These associations were nominally significant before correction for multiple hypothesis testing. These corrections are highly conservative for association studies of untested hypotheses, and may be too conservative for replication studies. We therefore believe the true effect of these four nominally significant SNPs on NTD risk will be more definitively determined by further study in other populations, and eventual meta-analysis.

Vitamin B12 status and folic acid supplementation influence mitochondrial heteroplasmy levels in mice.

One-carbon metabolism is a complex network of metabolic reactions that are essential for cellular function including DNA synthesis. Vitamin B12 and folate are micronutrients that are utilized in this pathway and their deficiency can result in the perturbation of one-carbon metabolism and subsequent perturbations in DNA replication and repair. This effect has been well characterized in nuclear DNA but to date, mitochondrial DNA (mtDNA) has not been investigated extensively. Mitochondrial variants have been associated with several inherited and age-related disease states; therefore, the study of factors that impact heteroplasmy are important for advancing our understanding of the mitochondrial genome's impact on human health. Heteroplasmy studies require robust and efficient mitochondrial DNA enrichment to carry out in-depth mtDNA sequencing. Many of the current methods for mtDNA enrichment can introduce biases and false-positive results. Here, we use a method that overcomes these limitations and have applied it to assess mitochondrial heteroplasmy in mouse models of altered one-carbon metabolism. Vitamin B12 deficiency was found to cause increased levels of mitochondrial DNA heteroplasmy across all tissues that were investigated. Folic acid supplementation also contributed to elevated mitochondrial DNA heteroplasmy across all mouse tissues investigated. Heteroplasmy analysis of human data from the Framingham Heart Study suggested a potential sex-specific effect of folate and vitamin B12 status on mitochondrial heteroplasmy. This is a novel relationship that may have broader consequences for our understanding of one-carbon metabolism, mitochondrial-related disease and the influence of nutrients on DNA mutation rates.

Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

Genotyping of a tri-allelic polymorphism by a novel melting curve assay in MTHFD1L: an association study of nonsyndromic Cleft in Ireland.

BACKGROUND: Polymorphisms within the MTHFD1L gene were previously associated with risk of neural tube defects in Ireland. We sought to test the most significant MTHFD1L polymorphisms for an association with risk of cleft in an Irish cohort. This required the development of a new melting curve assay to genotype the technically challenging MTHFD1L triallelic deletion/insertion polymorphism (rs3832406). METHODS: Melting curve analysis was used to genotype the MTHFD1L triallelic deletion/insertion polymorphism (rs3832406) and a Single Nucleotide Polymorphism rs17080476 in an Irish cohort consisting of 981 Irish case-parent trios and 1,008 controls. Tests for association with nonsyndromic cleft lip with or without cleft palate and cleft palate included case/control analysis, mother/control analysis and Transmission Disequilibrium Tests of case-parent trios. RESULTS: A successful melting curve genotyping assay was developed for the deletion/insertion polymorphism (rs3832406). The TDT analysis initially showed that the rs3832406 polymorphism was associated with isolated cleft lip with or without cleft palate. However, corrected p-values indicated that this association was not significant. CONCLUSIONS: Melting Curve Analysis can be employed to successfully genotype challenging polymorphisms such as the MTHFD1L triallelic deletion/insertion polymorphism (DIP) reported here (rs3832406) and is a viable alternative to capillary electrophoresis. Corrected p-values indicate no association between MTHFD1L and risk of cleft in an Irish cohort.

Evaluation of common genetic variants in 82 candidate genes as risk factors for neural tube defects.

BACKGROUND: Neural tube defects (NTDs) are common birth defects (~1 in 1000 pregnancies in the US and Europe) that have complex origins, including environmental and genetic factors. A low level of maternal folate is one well-established risk factor, with maternal periconceptional folic acid supplementation reducing the occurrence of NTD pregnancies by 50-70%. Gene variants in the folate metabolic pathway (e.g., MTHFR rs1801133 (677 C > T) and MTHFD1 rs2236225 (R653Q)) have been found to increase NTD risk. We hypothesized that variants in additional folate/B12 pathway genes contribute to NTD risk. METHODS: A tagSNP approach was used to screen common variation in 82 candidate genes selected from the folate/B12 pathway and NTD mouse models. We initially genotyped polymorphisms in 320 Irish triads (NTD cases and their parents), including 301 cases and 341 Irish controls to perform case-control and family based association tests. Significantly associated polymorphisms were genotyped in a secondary set of 250 families that included 229 cases and 658 controls. The combined results for 1441 SNPs were used in a joint analysis to test for case and maternal effects. RESULTS: Nearly 70 SNPs in 30 genes were found to be associated with NTDs at the p < 0.01 level. The ten strongest association signals (p-value range: 0.0003-0.0023) were found in nine genes (MFTC, CDKN2A, ADA, PEMT, CUBN, GART, DNMT3A, MTHFD1 and T (Brachyury)) and included the known NTD risk factor MTHFD1 R653Q (rs2236225). The single strongest signal was observed in a new candidate, MFTC rs17803441 (OR = 1.61 [1.23-2.08], p = 0.0003 for the minor allele). Though nominally significant, these associations did not remain significant after correction for multiple hypothesis testing. CONCLUSIONS: To our knowledge, with respect to sample size and scope of evaluation of candidate polymorphisms, this is the largest NTD genetic association study reported to date. The scale of the study and the stringency of correction are likely to have contributed to real associations failing to survive correction. We have produced a ranked list of variants with the strongest association signals. Variants in the highest rank of associations are likely to include true associations and should be high priority candidates for further study of NTD risk.

Mito-SiPE is a sequence-independent and PCR-free mtDNA enrichment method for accurate ultra-deep mitochondrial sequencing.

The analysis of somatic variation in the mitochondrial genome requires deep sequencing of mitochondrial DNA. This is ordinarily achieved by selective enrichment methods, such as PCR amplification or probe hybridization. These methods can introduce bias and are prone to contamination by nuclear-mitochondrial sequences (NUMTs), elements that can introduce artefacts into heteroplasmy analysis. We isolated intact mitochondria using differential centrifugation and alkaline lysis and subjected purified mitochondrial DNA to a sequence-independent and PCR-free method to obtain ultra-deep (>80,000X) sequencing coverage of the mitochondrial genome. This methodology avoids false-heteroplasmy calls that occur when long-range PCR amplification is used for mitochondrial DNA enrichment. Previously published methods employing mitochondrial DNA purification did not measure mitochondrial DNA enrichment or utilise high coverage short-read sequencing. Here, we describe a protocol that yields mitochondrial DNA and have quantified the increased level of mitochondrial DNA post-enrichment in 7 different mouse tissues. This method will enable researchers to identify changes in low frequency heteroplasmy without introducing PCR biases or NUMT contamination that are incorrectly identified as heteroplasmy when long-range PCR is used.

Exome sequencing identifies variants in infants with sacral agenesis.

BACKGROUND: Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA. METHODS: Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced. RESULTS: Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children. CONCLUSIONS: To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies.