An algorithm to identify patients aged 0-3 with rare genetic disorders.

BACKGROUND: With over 7000 Mendelian disorders, identifying children with a specific rare genetic disorder diagnosis through structured electronic medical record data is challenging given incompleteness of records, inaccurate medical diagnosis coding, as well as heterogeneity in clinical symptoms and procedures for specific disorders. We sought to develop a digital phenotyping algorithm (PheIndex) using electronic medical records to identify children aged 0-3 diagnosed with genetic disorders or who present with illness with an increased risk for genetic disorders. RESULTS: Through expert opinion, we established 13 criteria for the algorithm and derived a score and a classification. The performance of each criterion and the classification were validated by chart review. PheIndex identified 1,088 children out of 93,154 live births who may be at an increased risk for genetic disorders. Chart review demonstrated that the algorithm achieved 90% sensitivity, 97% specificity, and 94% accuracy. CONCLUSIONS: The PheIndex algorithm can help identify when a rare genetic disorder may be present, alerting providers to consider ordering a diagnostic genetic test and/or referring a patient to a medical geneticist.

Identification of copy number variants with genome sequencing: Clinical experiences from the NYCKidSeq program.

Copy number variations (CNVs) play a significant role in human disease. While chromosomal microarray has traditionally been the first-tier test for CNV detection, use of genome sequencing (GS) is increasing. We report the frequency of CNVs detected with GS in a diverse pediatric cohort from the NYCKidSeq program and highlight specific examples of its clinical impact. A total of 1052 children (0-21 years) with neurodevelopmental, cardiac, and/or immunodeficiency phenotypes received GS. Phenotype-driven analysis was used, resulting in 183 (17.4%) participants with a diagnostic result. CNVs accounted for 20.2% of participants with a diagnostic result (37/183) and ranged from 0.5 kb to 16 Mb. Of participants with a diagnostic result (n = 183) and phenotypes in more than one category, 5/17 (29.4%) were solved by a CNV finding, suggesting a high prevalence of diagnostic CNVs in participants with complex phenotypes. Thirteen participants with a diagnostic CNV (35.1%) had previously uninformative genetic testing, of which nine included a chromosomal microarray. This study demonstrates the benefits of GS for reliable detection of CNVs in a pediatric cohort with variable phenotypes.

Clinical Pharmacogenomic MT-RNR1 Screening for Aminoglycoside-Induced Ototoxicity and the Post-Test Counseling Conundrum.

Aminoglycoside antibiotic exposure can result in ototoxicity and irreversible hearing loss among individuals that harbor the m.1555A>G variant in the mitochondrial 12S rRNA gene, MT-RNR1. Importantly, pre-emptive m.1555A>G screening has been shown to reduce the prevalence of pediatric aminoglycoside-induced ototoxicity; however, professional guidelines to support and guide post-test pharmacogenomic counseling in this context are not currently available. This Perspective highlights key issues with delivering MT-RNR1 results, including longitudinal familial care considerations and communicating m.1555A>G heteroplasmy.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

PURPOSE: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. METHODS: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. RESULTS: A total of 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses (P < .001). Yield was greater for GS vs TGPs in Hispanic/Latino(a) (17.2% vs 9.5%, P < .001) and White/European American (19.8% vs 7.9%, P < .001) but not in Black/African American (11.5% vs 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. CONCLUSION: GS may yield up to twice as many diagnoses in pediatric patients compared with TGP testing but not yet across all population groups.

Molecular diagnostic yield of genome sequencing versus targeted gene panel testing in racially and ethnically diverse pediatric patients.

Purpose: Adoption of genome sequencing (GS) as a first-line test requires evaluation of its diagnostic yield. We evaluated the GS and targeted gene panel (TGP) testing in diverse pediatric patients (probands) with suspected genetic conditions. Methods: Probands with neurologic, cardiac, or immunologic conditions were offered GS and TGP testing. Diagnostic yield was compared using a fully paired study design. Results: 645 probands (median age 9 years) underwent genetic testing, and 113 (17.5%) received a molecular diagnosis. Among 642 probands with both GS and TGP testing, GS yielded 106 (16.5%) and TGPs yielded 52 (8.1%) diagnoses ( P < .001). Yield was greater for GS vs . TGPs in Hispanic/Latino(a) (17.2% vs . 9.5%, P < .001) and White/European American (19.8% vs . 7.9%, P < .001), but not in Black/African American (11.5% vs . 7.7%, P = .22) population groups by self-report. A higher rate of inconclusive results was seen in the Black/African American (63.8%) vs . White/European American (47.6%; P = .01) population group. Most causal copy number variants (17 of 19) and mosaic variants (6 of 8) were detected only by GS. Conclusion: GS may yield up to twice as many diagnoses in pediatric patients compared to TGP testing, but not yet across all population groups.

Detection of mosaic variants using genome sequencing in a large pediatric cohort.

The increased use of next-generation sequencing has expanded our understanding of the involvement and prevalence of mosaicism in genetic disorders. We describe a total of eleven cases: nine in which mosaic variants detected by genome sequencing (GS) and/or targeted gene panels (TGPs) were considered to be causative for the proband's phenotype, and two of apparent parental mosaicism. Variants were identified in the following genes: PHACTR1, SCN8A, KCNT1, CDKL5, NEXMIF, CUX1, TSC2, GABRB2, and SMARCB1. In addition, we identified one large duplication including three genes, UBE3A, GABRB3, and MAGEL2, and one large deletion including deletion of ARFGAP1, EEF1A2, CHRNA4, and KCNQ2. All patients were enrolled in the NYCKidSeq study, a research program studying the communication of genomic information in clinical care, as well as the clinical utility and diagnostic yield of GS for children with suspected genetic disorders in diverse populations in New York City. We observed variability in the correlation between reported variant allele fraction and the severity of the patient's phenotype, although we were not able to determine the mosaicism percentage in clinically relevant tissue(s). Although our study was not sufficiently powered to assess differences in mosaicism detection between the two testing modalities, we saw a trend toward better detection by GS as compared with TGP testing. This case series supports the importance of mosaicism in childhood-onset genetic conditions and informs guidelines for laboratory and clinical interpretation of mosaic variants detected by GS.

Long-read HiFi sequencing of NUDT15: Phased full-gene haplotyping and pharmacogenomic allele discovery.

To determine the phase of NUDT15 sequence variants for more comprehensive star (*) allele diplotyping, we developed a novel long-read single-molecule real-time HiFi amplicon sequencing method. A 10.5 kb NUDT15 amplicon assay was validated using reference material positive controls and additional samples for specimen type and blinded accuracy assessment. Triplicate NUDT15 HiFi sequencing of two reference material samples had nonreference genotype concordances of >99.9%, indicating that the assay is robust. Notably, short-read genome sequencing of a subset of samples was unable to determine the phase of star (*) allele-defining NUDT15 variants, resulting in ambiguous diplotype results. In contrast, long-read HiFi sequencing phased all variants across the NUDT15 amplicons, including a *2/*9 diplotype that previously was characterized as *1/*2 in the 1000 Genomes Project v3 data set. Assay throughput was also tested using 8.5 kb amplicons from 100 Ashkenazi Jewish individuals, which identified a novel NUDT15 *1 suballele (c.-121G>A) and a rare likely deleterious coding variant (p.Pro129Arg). Both novel alleles were Sanger confirmed and assigned as *1.007 and *20, respectively, by the PharmVar Consortium. Taken together, NUDT15 HiFi amplicon sequencing is an innovative method for phased full-gene characterization and novel allele discovery, which could improve NUDT15 pharmacogenomic testing and subsequent phenotype prediction.

Characterization of Reference Materials for TPMT and NUDT15: A GeT-RM Collaborative Project.

Pharmacogenetic testing is increasingly provided by clinical and research laboratories; however, only a limited number of quality control and reference materials are currently available for many of the TPMT and NUDT15 variants included in clinical tests. To address this need, the Division of Laboratory Systems, Centers for Disease Control and Prevention-based Genetic Testing Reference Material (GeT-RM) coordination program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Institute for Medical Research, has characterized 19 DNA samples derived from Coriell cell lines. DNA samples were distributed to four volunteer testing laboratories for genotyping using a variety of commercially available and laboratory developed tests and/or Sanger sequencing. Of the 12 samples characterized for TPMT, newly identified variants include TPMT∗2, ∗6, ∗12, ∗16, ∗21, ∗24, ∗32, ∗33, and ∗40; for the 7 NUDT15 reference material samples, newly identified variants are NUDT15∗2, ∗3, ∗4, ∗5, ∗6, and ∗9. In addition, a novel haplotype, TPMT∗46, was identified in this study. Preexisting data on an additional 11 Coriell samples, as well as some supplemental testing, were used to create comprehensive reference material panels for TPMT and NUDT15. These publicly available and well-characterized materials can be used to support the quality assurance and quality control programs of clinical laboratories performing clinical pharmacogenetic testing.

Deletion of ERF and CIC causes abnormal skull morphology and global developmental delay.

The ETS2 repressor factor (ERF) is a transcription factor in the RAS-MEK-ERK signal transduction cascade that regulates cell proliferation and differentiation, and pathogenic sequence variants in the ERF gene cause variable craniosynostosis inherited in an autosomal dominant pattern. The reported ERF variants are largely loss-of-function, implying haploinsufficiency as a primary disease mechanism; however, ERF gene deletions have not been reported previously. Here we describe three probands with macrocephaly, craniofacial dysmorphology, and global developmental delay. Clinical genetic testing for fragile X and other relevant sequencing panels were negative; however, chromosomal microarray identified heterozygous deletions (63.7-583.2 kb) on Chromosome 19q13.2 in each proband that together included five genes associated with Mendelian diseases (ATP1A3, ERF, CIC, MEGF8, and LIPE). Parental testing indicated that the aberrations were apparently de novo in two of the probands and were inherited in the one proband with the smallest deletion. Deletion of ERF is consistent with the reported loss-of-function ERF variants, prompting clinical copy-number-variant classifications of likely pathogenic. Moreover, the recent characterization of heterozygous loss-of-function CIC sequence variants as a cause of intellectual disability and neurodevelopmental disorders inherited in an autosomal dominant pattern is also consistent with the developmental delays and intellectual disabilities identified among the two probands with CIC deletions. Taken together, this case series adds to the previously reported patients with ERF and/or CIC sequence variants and supports haploinsufficiency of both genes as a mechanism for a variable syndromic cranial phenotype with developmental delays and intellectual disability inherited in an autosomal dominant pattern.

Development and Analytical Validation of a 29 Gene Clinical Pharmacogenetic Genotyping Panel: Multi-Ethnic Allele and Copy Number Variant Detection.

To develop a novel pharmacogenetic genotyping panel, a multidisciplinary team evaluated available evidence and selected 29 genes implicated in interindividual drug response variability, including 130 sequence variants and additional copy number variants (CNVs). Of the 29 genes, 11 had guidelines published by the Clinical Pharmacogenetics Implementation Consortium. Targeted genotyping and CNV interrogation were accomplished by multiplex single-base extension using the MassARRAY platform (Agena Biosciences) and multiplex ligation-dependent probe amplification (MRC Holland), respectively. Analytical validation of the panel was accomplished by a strategic combination of > 500 independent tests performed on 170 unique reference material DNA samples, which included sequence variant and CNV accuracy, reproducibility, and specimen (blood, saliva, and buccal swab) controls. Among the accuracy controls were 32 samples from the 1000 Genomes Project that were selected based on their enrichment of sequence variants included in the pharmacogenetic panel (VarCover.org). Coupled with publicly available samples from the Genetic Testing Reference Materials Coordination Program (GeT-RM), accuracy validation material was available for the majority (77%) of interrogated sequence variants (100% with average allele frequencies > 0.1%), as well as additional structural alleles with unique copy number signatures (e.g., CYP2D6*5, *13, *36, *68; CYP2B6*29; and CYP2C19*36). Accuracy and reproducibility for both genotyping and copy number were > 99.9%, indicating that the optimized panel platforms were precise and robust. Importantly, multi-ethnic allele frequencies of the interrogated variants indicate that the vast majority of the general population carries at least one of these clinically relevant pharmacogenetic variants, supporting the implementation of this panel for pharmacogenetic research and/or clinical implementation programs.

Haploinsufficiency of the basic helix-loop-helix transcription factor HAND2 causes congenital heart defects.

Congenital heart defects (CHDs) are caused by a disruption in heart morphogenesis, which is dependent, in part, on a network of transcription factors (TFs) that regulate myocardial development. Heterozygous sequence variants in the basic helix-loop-helix TF gene heart and neural crest derivatives expressed 2 (HAND2) have been reported among some patients with CHDs; however, HAND2 has not yet been established as a Mendelian disease gene. We report a 31-month-old male with unicommissural unicuspid aortic valve, moderate aortic stenosis, and mild pulmonic stenosis. Chromosome analysis revealed a normal 46,XY karyotype, and a CHD sequencing panel was negative for pathogenic variants in NKX2.5, GATA4, TBX5, and CHD7. However, chromosomal microarray (CMA) testing identified a heterozygous 546.0-kb deletion on chromosome 4q34.1 (174364195_174910239[GRCh37/hg19]) that included exons 1 and 2 of SCRG1, HAND2, and HAND2-AS1. Familial CMA testing determined that the deletion was paternally inherited, which supported a likely pathogenic classification as the proband's father had previously undergone surgery for Tetralogy of Fallot. The family history was also notable for a paternal uncle who had previously died from complications related to an unknown heart defect. Taken together, this first report of a HAND2 and HAND2-AS1 deletion in a family with CHDs strongly supports haploinsufficiency of HAND2 as an autosomal dominant cause of CHD.

Lessons learned from expanded reproductive carrier screening in self-reported Ashkenazi, Sephardi, and Mizrahi Jewish patients.

BACKGROUND: Next-generation sequencing (NGS)-based panels have gained traction as a strategy for reproductive carrier screening. Their value for screening Ashkenazi Jewish (AJ) individuals, who have benefited greatly from population-wide targeted testing, as well as Sephardi/Mizrahi Jewish (SMJ) individuals (an underserved population), has not been fully explored. METHODS: The clinical utilization by 6,805 self-reported Jewish individuals of an expanded NGS panel, along with several ancillary assays, was assessed retrospectively. Data were extracted for a subset of 96 diseases that, during the panel design phase, were classified as being AJ-, SMJ-, or pan-Jewish/pan-ethnic-relevant. RESULTS: 64.6% of individuals were identified as carriers of one or more of these 96 diseases. Over 80% of the reported variants would have been missed by following recommended AJ screening guidelines. 10.7% of variants reported for AJs were in "SMJ-relevant genes," and 31.2% reported for SMJs were in "AJ-relevant genes." Roughly 2.5% of individuals carried a novel, likely pathogenic variant. One in 16 linked cohort couples was identified as a carrier couple for at least one of these 96 diseases. CONCLUSION: For maximal carrier identification, this study supports using expanded NGS panels for individuals of all Jewish backgrounds. This approach can better empower at-risk couples for reproductive decision making.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles.

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Prenatal cytogenomic identification and molecular refinement of compound heterozygous STRC deletion breakpoints.

Here, we report the prenatal detection of a compound heterozygous deletion at chromosome 15q15.3 by clinical chromosomal microarray (CMA) testing that included the CATSPER2 male infertility gene. However, given the low resolution of CMA at this homologous locus, it was unclear if the neighboring STRC hearing loss gene was also affected. Therefore, we developed a novel allele-specific PCR strategy, which narrowed the proximal breakpoint of the maternally inherited deletion to a 310 bp interval that was 440 bp upstream from the STRC transcription start site.

Clinical Genetic Testing for Fragile X Syndrome by Polymerase Chain Reaction Amplification and Southern Blot Analyses.

Fragile X syndrome (FXS) is characterized by mental retardation and in the vast majority of cases it is caused by expansion of CGG trinucleotide repeats in the 5' untranslated region (or UTR) in the FMR1 gene on the X chromosome. The size and methylation status of CGG repeats are correlated with the clinical phenotype of FMR1-related disorders. The methods used for clinical genetic testing of FXS include polymerase chain reaction (PCR) amplification and Southern blot analyses, which effectively detect alleles with repeats in the normal, intermediate, premutation, and full mutation size ranges, as well as the methylation status of FMR1 promoter region.

Integrated CYP2D6 interrogation for multiethnic copy number and tandem allele detection.

AIM: To comprehensively interrogate CYP2D6 by integrating genotyping, copy number analysis and novel strategies to identify CYP2D6*36 and characterize CYP2D6 duplications. METHODS: Genotyping of 16 CYP2D6 alleles, multiplex ligation-dependent probe amplification (MLPA) and CYP2D6*36 and duplication allele-specific genotyping were performed on 427 African-American, Asian, Caucasian, Hispanic, and Ashkenazi Jewish individuals. RESULTS: A novel PCR strategy determined that almost half of all CYP2D6*10 (100C>T) alleles are actually *36 (isolated or in tandem with *10) and all identified duplication alleles were characterized. Integrated results from all testing platforms enabled the refinement of genotype frequencies across all studied populations. CONCLUSION: The polymorphic CYP2D6 gene requires comprehensive interrogation to characterize allelic variation across ethnicities, which was enabled in this study by integrating multiplexed genotyping, MLPA copy number analysis, novel PCR strategies and duplication allele-specific genotyping.

The New York pilot newborn screening program for lysosomal storage diseases: Report of the First 65,000 Infants.

PURPOSE: We conducted a consented pilot newborn screening (NBS) for Pompe, Gaucher, Niemann-Pick A/B, Fabry, and MPS 1 to assess the suitability of these lysosomal storage disorders (LSDs) for public health mandated screening. METHODS: At five participating high-birth rate, ethnically diverse New York City hospitals, recruiters discussed the study with postpartum parents and documented verbal consent. Screening on consented samples was performed using multiplexed tandem mass spectrometry. Screen-positive infants underwent confirmatory enzymology, DNA testing, and biomarker quantitation when available. Affected infants are being followed for clinical management and long-term outcome. RESULTS: Over 4 years, 65,605 infants participated, representing an overall consent rate of 73%. Sixty-nine infants were screen-positive. Twenty-three were confirmed true positives, all of whom were predicted to have late-onset phenotypes. Six of the 69 currently have undetermined disease status. CONCLUSION: Our results suggest that NBS for LSDs is much more likely to detect individuals at risk for late-onset disease, similar to results from other NBS programs. This work has demonstrated the feasibility of using a novel consented pilot NBS study design that can be modified to include other disorders under consideration for public health implementation as a means to gather critical evidence for evidence-based NBS practices.

Prenatal Diagnosis of Cystic Fibrosis.

Cystic fibrosis (CF) is an inherited disease characterized by the accumulation of thick, sticky mucus which damages epithelia in organs such as the lungs, pancreas, liver, intestines, and other parts of the body. The most common symptoms are sinopulmonary disease and chronic gastrointestinal tract problems resulting from decreased mucociliary clearance and inflammation. CF is the most common life-limiting autosomal recessive disorder in people of northern European ancestry and it affects other populations with different prevalence. CF can be diagnosed by many methods including testing for blood immunoreactive trypsin, sweat chloride, transepithelial nasal potential difference, and molecular genetic testing.

Identification of rare de novo epigenetic variations in congenital disorders.

Certain human traits such as neurodevelopmental disorders (NDs) and congenital anomalies (CAs) are believed to be primarily genetic in origin. However, even after whole-genome sequencing (WGS), a substantial fraction of such disorders remain unexplained. We hypothesize that some cases of ND-CA are caused by aberrant DNA methylation leading to dysregulated genome function. Comparing DNA methylation profiles from 489 individuals with ND-CAs against 1534 controls, we identify epivariations as a frequent occurrence in the human genome. De novo epivariations are significantly enriched in cases, while RNAseq analysis shows that epivariations often have an impact on gene expression comparable to loss-of-function mutations. Additionally, we detect and replicate an enrichment of rare sequence mutations overlapping CTCF binding sites close to epivariations, providing a rationale for interpreting non-coding variation. We propose that epivariations contribute to the pathogenesis of some patients with unexplained ND-CAs, and as such likely have diagnostic relevance.

Robust identification of mosaic variants in congenital heart disease.

Mosaicism due to somatic mutations can cause multiple diseases including cancer, developmental and overgrowth syndromes, neurodevelopmental disorders, autoinflammatory diseases, and atrial fibrillation. With the increased use of next generation sequencing technology, multiple tools have been developed to identify low-frequency variants, specifically from matched tumor-normal tissues in cancer studies. To investigate whether mosaic variants are implicated in congenital heart disease (CHD), we developed a pipeline using the cancer somatic variant caller MuTect to identify mosaic variants in whole-exome sequencing (WES) data from a cohort of parent/affected child trios (n = 715) and a cohort of healthy individuals (n = 416). This is a novel application of the somatic variant caller designed for cancer to WES trio data. We identified two cases with mosaic KMT2D mutations that are likely pathogenic for CHD, but conclude that, overall, mosaicism detectable in peripheral blood or saliva does not account for a significant portion of CHD etiology.