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.

Rare Variants in RPPH1 Real-Time Quantitative PCR Control Assay Binding Sites Result in Incorrect Copy Number Calls.

Real-time quantitative PCR (qPCR) using RPPH1 as a reference gene is a standard method for assessment and validation of genomic copy number variations. However, variants in the reference amplicon may cause errors, which was investigated herein. While conducting copy number variation validations for birth defects research studies, 13 of 1634 specimens with multiple loci that appeared to be present as three copies were unexpectedly detected. This apparent trisomy was hypothesized to be an amplification artifact caused by a variant in the RPPH1 amplicon. Sequencing revealed all 13 individuals carried one of the four different variants within the RPPH1 amplicon. These variants could produce allelic dropout or altered reaction efficiency, causing an inaccurate measurement of copy number. Additional genotyping predicted a low frequency of the most common variant (rs3093876; 14/3562 alleles; minor allele frequency, 0.39%). Laboratories should recognize the potential for inaccurate results when using a single qPCR control assay. Overestimated CFTR and SMN2 copy numbers identified during newborn screening that otherwise would have been incorrectly called were also detected. Variants in reference loci may produce false-negative normal results for test loci when real deletions are present. For clinical laboratories screening for heterozygous deletions for diagnostic testing or prenatal/carrier screening via qPCR, the most cost-effective solution to maximize sensitivity is to run triplex reactions targeting the region of interest with two control genes.

Validation of a Custom Next-Generation Sequencing Assay for Cystic Fibrosis Newborn Screening.

Newborn screening (NBS) for Cystic Fibrosis (CF) is associated with improved outcomes. All US states screen for CF; however, CF NBS algorithms have high false positive (FP) rates. In New York State (NYS), the positive predictive value of CF NBS improved from 3.7% to 25.2% following the implementation of a three-tier IRT-DNA-SEQ approach using commercially available tests. Here we describe a modification of the NYS CF NBS algorithm via transition to a new custom next-generation sequencing (NGS) platform for more comprehensive cystic fibrosis transmembrane conductance regulator (CFTR) gene analysis. After full gene sequencing, a tiered strategy is used to first analyze only a specific panel of 338 clinically relevant CFTR variants (second-tier), followed by unblinding of all sequence variants and bioinformatic assessment of deletions/duplications in a subset of samples requiring third-tier analysis. We demonstrate the analytical and clinical validity of the assay and the feasibility of use in the NBS setting. The custom assay has streamlined our molecular workflow, increased throughput, and allows for bioinformatic customization of second-tier variant panel content. NBS aims to identify those infants with the highest disease risk. Technological molecular improvements can be applied to NBS algorithms to reduce the burden of FP referrals without loss of sensitivity.

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.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

Copy number variants in hypoplastic right heart syndrome.

Hypoplastic right heart syndrome (HRHS) is a rare congenital defect characterized by underdeveloped and malformed structures of the right heart. Familial recurrence of HRHS indicates genetic factors contribute to its etiology. Our study investigates the presence of copy number variants (CNVs) in HRHS cases. We genotyped 42 HRHS cases identified from live births throughout California (2003-2010) using the Illumina HumanOmni2.5-8 array. We identified 14 candidate CNVs in 14 HRHS cases (33%) based on the genes included in the CNVs and their functions. Duplications overlapping part of ERBB4 were identified in two unrelated cases. ERBB4 is a neuregulin receptor with a pivotal role in cardiomyocyte differentiation and heart development. We also described a 7.5 Mb duplication at 16q11-12. Multiple genes in the duplicated region have previously been linked to heart defects and cardiac development, including RPGRIP1L, RBL2, SALL1, and MYLK3. Of the 14 validated CNVs, we identified four CNVs in close proximity to genes linked to the Wnt signaling pathway. This study expands on our previous work supporting the role of genetics in HRHS. We identified CNVs affecting crucial genes and signaling pathways involved in right heart development. ERBB4 and duplication of the 16q11-12 region are important areas for future investigation.

Rare copy number variants identified in prune belly syndrome.

Prune belly syndrome (PBS), also known as Eagle-Barrett syndrome, is a rare congenital disorder characterized by absence or hypoplasia of the abdominal wall musculature, urinary tract anomalies, and cryptorchidism in males. The etiology of PBS is largely unresolved, but genetic factors are implicated given its recurrence in families. We examined cases of PBS to identify novel pathogenic copy number variants (CNVs). A total of 34 cases (30 males and 4 females) with PBS identified from all live births in New York State (1998-2005) were genotyped using Illumina HumanOmni2.5 microarrays. CNVs were prioritized if they were absent from in-house controls, encompassed ≥10 consecutive probes, were ≥20 Kb in size, had ≤20% overlap with common variants in population reference controls, and had ≤20% overlap with any variant previously detected in other birth defect phenotypes screened in our laboratory. We identified 17 candidate autosomal CNVs; 10 cases each had one CNV and four cases each had two CNVs. The CNVs included a 158 Kb duplication at 4q22 that overlaps the BMPR1B gene; duplications of different sizes carried by two cases in the intron of STIM1 gene; a 67 Kb duplication 202 Kb downstream of the NOG gene, and a 1.34 Mb deletion including the MYOCD gene. The identified rare CNVs spanned genes involved in mesodermal, muscle, and urinary tract development and differentiation, which might help in elucidating the genetic contribution to PBS. We did not have parental DNA and cannot identify whether these CNVs were de novo or inherited. Further research on these CNVs, particularly BMP signaling is warranted to elucidate the pathogenesis of PBS.

Copy number variants in Ebstein anomaly.

BACKGROUND: Ebstein anomaly (EA) is a rare congenital defect characterized by apical displacement of the septal tricuspid leaflets and atrialization of the right ventricle. The etiology of EA is unclear; however, recurrence in families and the association of EA with genetic syndromes and copy number variants (CNVs) suggest a genetic component. OBJECTIVE: We performed a population-based study to search for recurrent and novel CNVs in a previously unreported set of EA cases. METHODS: We genotyped 60 EA cases identified from all live births (2,891,076) from selected California counties (1991-2010) using the Illumina HumanOmni2.5-8 array. We identified 38 candidate CNVs in 28 (46%) cases and prioritized and validated 11 CNVs based on the genes included. RESULTS: Five CNVs (41%) overlapped or were close to genes involved in early myocardial development, including NODAL, PDLIM5, SIX1, ASF1A and FGF12. We also replicated a previous association of EA with CNVs at 1p34.1 and AKAP12. Finally, we identified four CNVs overlapping or in close proximity to the transcription factors HES3, TRIM71, CUX1 and EIF4EBP2. CONCLUSIONS: This study supports the relationship of genetic factors to EA and demonstrates that defects in cardiomyocytes and myocardium differentiation may play a role. Abnormal differentiation of cardiomyocytes and how genetic factors contribute should be examined for their association with EA.

Copy-number variants and candidate gene mutations in isolated split hand/foot malformation.

Split hand/foot malformation (SHFM) is a congenital limb deficiency with missing or shortened central digits. Some SHFM genes have been identified but the cause of many SHFM cases is unknown. We used single-nucleotide polymorphism (SNP) microarray analysis to detect copy-number variants (CNVs) in 25 SHFM cases without other birth defects from New York State (NYS), prioritized CNVs absent from population CNV databases, and validated these CNVs using quantitative real-time polymerase chain reaction (qPCR). We tested for the validated CNVs in seven cases from Iowa using qPCR, and also sequenced 36 SHFM candidate genes in all the subjects. Seven NYS cases had a potentially deleterious variant: two had a p.R225H or p.R225L mutation in TP63, one had a 17q25 microdeletion, one had a 10q24 microduplication and three had a 17p13.3 microduplication. In addition, one Iowa case had a de novo 10q24 microduplication. The 17q25 microdeletion has not been reported previously in SHFM and included two SHFM candidate genes (SUMO2 and GRB2), while the 10q24 and 17p13.3 CNVs had breakpoints within genomic regions that contained putative regulatory elements and a limb development gene. In SHFM pathogenesis, the microdeletion may cause haploinsufficiency of SHFM genes and/or deletion of their regulatory regions, and the microduplications could disrupt regulatory elements that control transcription of limb development genes.

Rare copy number variants in a population-based investigation of hypoplastic right heart syndrome.

BACKGROUND: Hypoplastic right heart syndrome (HRHS) is a rare congenital defect characterized by underdevelopment of the right heart structures commonly accompanied by an atrial septal defect. Familial HRHS reports suggest genetic factor involvement. We examined the role of copy number variants (CNVs) in HRHS. METHODS: We genotyped 32 HRHS cases identified from all New York State live births (1998-2005) using Illumina HumanOmni2.5 microarrays. CNVs were called with PennCNV and prioritized if they were ≥20 Kb, contained ≥10 SNPs and had minimal overlap with CNVs from in-house controls, the Database of Genomic Variants, HapMap3, and Childrens Hospital of Philadelphia database. RESULTS: We identified 28 CNVs in 17 cases; several encompassed genes important for right heart development. One case had a 2p16-2p23 duplication spanning LBH, a limb and heart development transcription factor. Lbh mis-expression results in right ventricular hypoplasia and pulmonary valve defects. This duplication also encompassed SOS1, a factor associated with pulmonary valve stenosis in Noonan syndrome. Sos1-/- mice display thin and poorly trabeculated ventricles. In another case, we identified a 1.5 Mb deletion associated with Williams-Beuren syndrome, a disorder that includes valvular malformations. A third case had a 24 Kb deletion upstream of the TGFß ligand ITGB8. Embryos genetically null for Itgb8, and its intracellular interactant Band 4.1B, display lethal cardiac phenotypes. CONCLUSION: To our knowledge, this is the first study of CNVs in HRHS. We identified several rare CNVs that overlap genes related to right ventricular wall and valve development, suggesting that genetics plays a role in HRHS and providing clues for further investigation. Birth Defects Research 109:16-26, 2017. © 2016 Wiley Periodicals, Inc.

Copy number variants in a population-based investigation of Klippel-Trenaunay syndrome.

Klippel-Trenaunay syndrome (KTS) is a rare congenital vascular disorder that is thought to occur sporadically; however, reports of familial occurrence suggest a genetic component. We examined KTS cases to identify novel, potentially causal copy number variants (CNVs). We identified 17 KTS cases from all live-births occurring in New York (1998-2010). Extracted DNA was genotyped using Illumina microarrays and CNVs were called using PennCNV software. CNVs selected for follow-up had ≥10 single nucleotide polymorphisms (SNPs) and minimal overlap with in-house controls or controls from the Database of Genomic Variants. We identified 15 candidate CNVs in seven cases; among them a deletion in two cases within transcripts of HDAC9, a histone deacetylase essential for angiogenic sprouting of endothelial cells. One of them also had a duplication upstream of SALL3, a transcription factor essential for embryonic development that inhibits DNMT3A, a DNA methyltransferase responsible for embryonic de novo DNA methylation. Another case had a duplication spanning ING5, a histone acetylation regulator active during embryogenesis. We identified rare genetic variants related to chromatin modification which may have a key role in regulating vascular development during embryogenesis. Further investigation of their implications in the pathogenesis of KTS is warranted. © 2016 Wiley Periodicals, Inc.

Genetic Variants in Isolated Ebstein Anomaly Implicated in Myocardial Development Pathways.

Ebstein anomaly (EA) is a rare heart defect in which the tricuspid valve is malformed and displaced. The tricuspid valve abnormalities can lead to backflow of blood from the right ventricle to the right atrium, preventing proper circulation of blood to the lungs. Although the etiology of EA is largely unresolved, increased prevalence of EA in those with a family history of congenital heart disease suggests EA has a genetic component. Copy number variants (CNVs) are a major source of genetic variation and have been implicated in a range of congenital heart defect phenotypes. We performed a systematic, genome-wide search for CNVs in 47 isolated EA cases using genotyping microarrays. In addition, we used a custom HaloPlex panel to sequence three known EA genes and 47 candidate EA genes. We identified 35 candidate CNVs in 24 (51%) EA cases. Rare sequence variants in genes associated with cardiomyopathy were identified in 11 (23%) EA cases. Two CNVs near the transcriptional repressor HEY1, a member of the NOTCH signaling pathway, were identified in three unrelated cases. All other candidate CNVs were each identified in a single case. At least 11 of 35 candidate CNVs include genes involved in myocardial development or function, including multiple genes in the BMP signaling pathway. We identified enrichment of gene sets involved in histone modification and cardiomyocyte differentiation, supporting the involvement of the developing myocardium in the etiology of EA. Gene set enrichment analysis also identified ribosomal RNA processing, a potentially novel pathway of altered cardiac development in EA. Our results suggest an altered myocardial program may contribute to abnormal tricuspid valve development in EA. Future studies should investigate abnormal differentiation of cardiomyocytes as a potential etiological factor in EA.

Copy-number variant analysis of classic heterotaxy highlights the importance of body patterning pathways.

Classic heterotaxy consists of congenital heart defects with abnormally positioned thoracic and abdominal organs. We aimed to uncover novel, genomic copy-number variants (CNVs) in classic heterotaxy cases. A microarray containing 2.5 million single-nucleotide polymorphisms (SNPs) was used to genotype 69 infants (cases) with classic heterotaxy identified from California live births from 1998 to 2009. CNVs were identified using the PennCNV software. We identified 56 rare CNVs encompassing genes in the NODAL (NIPBL, TBX6), BMP (PPP4C), and WNT (FZD3) signaling pathways, not previously linked to classic heterotaxy. We also identified a CNV involving FGF12, a gene previously noted in a classic heterotaxy case. CNVs involving RBFOX1 and near MIR302F were detected in multiple cases. Our findings illustrate the importance of body patterning pathways for cardiac development and left/right axes determination. FGF12, RBFOX1, and MIR302F could be important in human heterotaxy, because they were noted in multiple cases. Further investigation into genes involved in the NODAL, BMP, and WNT body patterning pathways and into the dosage effects of FGF12, RBFOX1, and MIR302F is warranted.

Rare copy number variants implicated in posterior urethral valves.

The cause of posterior urethral valves (PUV) is unknown, but genetic factors are suspected given their familial occurrence. We examined cases of isolated PUV to identify novel copy number variants (CNVs). We identified 56 cases of isolated PUV from all live-births in New York State (1998-2005). Samples were genotyped using Illumina HumanOmni2.5 microarrays. Autosomal and sex-linked CNVs were identified using PennCNV and cnvPartition software. CNVs were prioritized for follow-up if they were absent from in-house controls, contained ≥ 10 consecutive probes, were ≥ 20 Kb in size, had ≤ 20% overlap with variants detected in other birth defect phenotypes screened in our lab, and were rare in population reference controls. We identified 47 rare candidate PUV-associated CNVs in 32 cases; one case had a 3.9 Mb deletion encompassing BMP7. Mutations in BMP7 have been associated with severe anomalies in the mouse urethra. Other interesting CNVs, each detected in a single PUV case included: a deletion of PIK3R3 and TSPAN1, duplication/triplication in FGF12, duplication of FAT1--a gene essential for normal growth and development, a large deletion (>2 Mb) on chromosome 17q that involves TBX2 and TBX4, and large duplications (>1 Mb) on chromosomes 3q and 6q. Our finding of previously unreported novel CNVs in PUV suggests that genetic factors may play a larger role than previously understood. Our data show a potential role of CNVs in up to 57% of cases examined. Investigation of genes in these CNVs may provide further insights into genetic variants that contribute to PUV.

Novel copy-number variants in a population-based investigation of classic heterotaxy.

PURPOSE: Heterotaxy is a clinically and genetically heterogeneous disorder. We investigated whether screening cases restricted to a classic phenotype would result in the discovery of novel, potentially causal copy-number variants. METHODS: We identified 77 cases of classic heterotaxy from all live births in New York State during 1998-2005. DNA extracted from each infant's newborn dried blood spot was genotyped with a microarray containing 2.5 million single-nucleotide polymorphisms. Copy-number variants were identified with PennCNV and cnvPartition software. Candidates were selected for follow-up if they were absent in unaffected controls, contained 10 or more consecutive probes, and had minimal overlap with variants published in the Database of Genomic Variants. RESULTS: We identified 20 rare copy-number variants including a deletion of BMP2, which has been linked to laterality disorders in mice but not previously reported in humans. We also identified a large, terminal deletion of 10q and a microdeletion at 1q23.1 involving the MNDA gene; both are rare variants suspected to be associated with heterotaxy. CONCLUSION: Our findings implicate rare copy-number variants in classic heterotaxy and highlight several candidate gene regions for further investigation. We also demonstrate the efficacy of copy-number variant genotyping in blood spots using microarrays.