A cross-disorder dosage sensitivity map of the human genome.

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies.

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

Genotype-Phenotype Correlation of Distal 2q37 Deletions.

Brachydactyly mental retardation syndrome (BDMR) typically results from large deletions (>2-9 Mb) in distal 2q37. Haploinsufficiency of HDAC4 with incomplete penetrance has been proposed as the primary genetic cause of BDMR. To date, pure 2q37 deletions distal to HDAC4 were reported only in a limited number of individuals who share a subset of the clinical manifestations seen in cases with 2q37 deletions encompassing HDAC4. Here, we present a 4-year-old African American male who carries the smallest established 2q37.3 deletion distal to HDAC4 (827.1 kb; 16 OMIM genes). His clinical features that overlap with BDMR phenotypes include expressive-receptive language delay, behavioral issues, mild facial dysmorphism such as frontal bossing, and bilateral 5th finger brachydactyly and clinodactyly. The deletion was inherited from his mother with a history of learning difficulties and similar facial dysmorphism. This case provides important genotype-phenotype correlation information and suggests a 2q37 region distal to HDAC4 encompassing the HDLBP gene may contribute to a subset of clinical features overlapping with those seen in individuals with BDMR.

Two Patients with Complex Rearrangements Suggestive of Germline Chromoanagenesis.

Chromoanagenesis, a phenomenon characterized by complex chromosomal rearrangement and reorganization events localized to a limited number of genomic regions, includes the subcategories chromothripsis, chromoanasynthesis, and chromoplexy. Although definitions of these terms are evolving, constitutional chromoanagenesis events have been reported in a limited number of patients with variable phenotypes. We report on 2 cases with complex genomic events characterized by multiple copy number gains and losses confined to a single chromosome region, which are suggestive of constitutional chromoanagenesis. Case 1 is a 43-year-old male with intellectual disability and recently developed generalized tonic-clonic seizures. Chromosomal microarray analysis identified a complex rearrangement involving chromosome region 14q31.1q32.2, consisting of 16 breakpoints ranging in size from 0.2 to 6.2 Mb, with 5 segments of normal copy number present between these alterations. Interestingly, this case represents the oldest known patient with a complex rearrangement indicative of constitutional chromoanagenesis. Case 2 is a 2-year-old female with developmental delay, speech delay, low muscle tone, and seizures. Chromosomal microarray analysis identified a complex rearrangement consisting of 28 breakpoints localized to 18q21.32q23. The size of the copy number alterations ranged from 0.042 to 5.1 Mb, flanked by 12 small segments of normal copy number. These cases add to a growing body of literature demonstrating complex chromosomal rearrangements as a disease mechanism for congenital anomalies.

Correction: A new microdeletion syndrome involving TBC1D24, ATP6V0C, and PDPK1 causes epilepsy, microcephaly, and developmental delay.

The original version of this Article contained an error in the spelling of the author Siddharth Banka, which was incorrectly given as Siddhart Banka. This has now been corrected in both the PDF and HTML versions of the Article.

A new microdeletion syndrome involving TBC1D24, ATP6V0C, and PDPK1 causes epilepsy, microcephaly, and developmental delay.

PURPOSE: Contiguous gene deletions are known to cause several neurodevelopmental syndromes, many of which are caused by recurrent events on chromosome 16. However, chromosomal microarray studies (CMA) still yield copy-number variants (CNVs) of unknown clinical significance. We sought to characterize eight individuals with overlapping 205-kb to 504-kb 16p13.3 microdeletions that are distinct from previously published deletion syndromes. METHODS: Clinical information on the patients and bioinformatic scores for the deleted genes were analyzed. RESULTS: All individuals in our cohort displayed developmental delay, intellectual disability, and various forms of seizures. Six individuals were microcephalic and two had strabismus. The deletion was absent in all 13 parents who were available for testing. The area of overlap encompasses seven genes including TBC1D24, ATP6V0C, and PDPK1 (also known as PDK1). Bi-allelic TBC1D24 pathogenic variants are known to cause nonsyndromic deafness, epileptic disorders, or DOORS syndrome (deafness, onychodystrophy, osteodystrophy, mental retardation, seizures). Sanger sequencing of the nondeleted TBC1D24 allele did not yield any additional pathogenic variants. CONCLUSIONS: We propose that 16p13.3 microdeletions resulting in simultaneous haploinsufficiencies of TBC1D24, ATP6V0C, and PDPK1 cause a novel rare contiguous gene deletion syndrome of microcephaly, developmental delay, intellectual disability, and epilepsy.

Paired-Duplication Signatures Mark Cryptic Inversions and Other Complex Structural Variation.

Copy-number variants (CNVs) have been the predominant focus of genetic studies of structural variation, and chromosomal microarray (CMA) for genome-wide CNV detection is the recommended first-tier genetic diagnostic screen in neurodevelopmental disorders. We compared CNVs observed by CMA to the structural variation detected by whole-genome large-insert sequencing in 259 individuals diagnosed with autism spectrum disorder (ASD) from the Simons Simplex Collection. These analyses revealed a diverse landscape of complex duplications in the human genome. One remarkably common class of complex rearrangement, which we term dupINVdup, involves two closely located duplications ("paired duplications") that flank the breakpoints of an inversion. This complex variant class is cryptic to CMA, but we observed it in 8.1% of all subjects. We also detected other paired-duplication signatures and duplication-mediated complex rearrangements in 15.8% of all ASD subjects. Breakpoint analysis showed that the predominant mechanism of formation of these complex duplication-associated variants was microhomology-mediated repair. On the basis of the striking prevalence of dupINVdups in this cohort, we explored the landscape of all inversion variation among the 235 highest-quality libraries and found abundant complexity among these variants: only 39.3% of inversions were canonical, or simple, inversions without additional rearrangement. Collectively, these findings indicate that dupINVdups, as well as other complex duplication-associated rearrangements, represent relatively common sources of genomic variation that is cryptic to population-based microarray and low-depth whole-genome sequencing. They also suggest that paired-duplication signatures detected by CMA warrant further scrutiny in genetic diagnostic testing given that they might mark complex rearrangements of potential clinical relevance.

Rare exonic deletions implicate the synaptic organizer Gephyrin (GPHN) in risk for autism, schizophrenia and seizures.

The GPHN gene codes for gephyrin, a key scaffolding protein in the neuronal postsynaptic membrane, responsible for the clustering and localization of glycine and GABA receptors at inhibitory synapses. Gephyrin has well-established functional links with several synaptic proteins that have been implicated in genetic risk for neurodevelopmental disorders such as autism spectrum disorder (ASD), schizophrenia and epilepsy including the neuroligins (NLGN2, NLGN4), the neurexins (NRXN1, NRXN2, NRXN3) and collybistin (ARHGEF9). Moreover, temporal lobe epilepsy has been linked to abnormally spliced GPHN mRNA lacking exons encoding the G-domain of the gephyrin protein, potentially arising due to cellular stress associated with epileptogenesis such as temperature and alkalosis. Here, we present clinical and genomic characterization of six unrelated subjects, with a range of neurodevelopmental diagnoses including ASD, schizophrenia or seizures, who possess rare de novo or inherited hemizygous microdeletions overlapping exons of GPHN at chromosome 14q23.3. The region of common overlap across the deletions encompasses exons 3-5, corresponding to the G-domain of the gephyrin protein. These findings, together with previous reports of homozygous GPHN mutations in connection with autosomal recessive molybdenum cofactor deficiency, will aid in clinical genetic interpretation of the GPHN mutation spectrum. Our data also add to the accumulating evidence implicating neuronal synaptic gene products as key molecular factors underlying the etiologies of a diverse range of neurodevelopmental conditions.

Rare deletions at the neurexin 3 locus in autism spectrum disorder.

The three members of the human neurexin gene family, neurexin 1 (NRXN1), neurexin 2 (NRXN2), and neurexin 3 (NRXN3), encode neuronal adhesion proteins that have important roles in synapse development and function. In autism spectrum disorder (ASD), as well as in other neurodevelopmental conditions, rare exonic copy-number variants and/or point mutations have been identified in the NRXN1 and NRXN2 loci. We present clinical characterization of four index cases who have been diagnosed with ASD and who possess rare inherited or de novo microdeletions at 14q24.3-31.1, a region that overlaps exons of the alpha and/or beta isoforms of NRXN3. NRXN3 deletions were found in one father with subclinical autism and in a carrier mother and father without formal ASD diagnoses, indicating issues of penetrance and expressivity at this locus. Notwithstanding these clinical complexities, this report on ASD-affected individuals who harbor NRXN3 exonic deletions advances the understanding of the genetic etiology of autism, further enabling molecular diagnoses.

Recurrent duplications of 17q12 associated with variable phenotypes.

The ability to identify the clinical nature of the recurrent duplication of chromosome 17q12 has been limited by its rarity and the diverse range of phenotypes associated with this genomic change. In order to further define the clinical features of affected patients, detailed clinical information was collected in the largest series to date (30 patients and 2 of their siblings) through a multi-institutional collaborative effort. The majority of patients presented with developmental delays varying from mild to severe. Though dysmorphic features were commonly reported, patients do not have consistent and recognizable features. Cardiac, ophthalmologic, growth, behavioral, and other abnormalities were each present in a subset of patients. The newly associated features potentially resulting from 17q12 duplication include height and weight above the 95th percentile, cataracts, microphthalmia, coloboma, astigmatism, tracheomalacia, cutaneous mosaicism, pectus excavatum, scoliosis, hypermobility, hypospadias, diverticulum of Kommerell, pyloric stenosis, and pseudohypoparathryoidism. The majority of duplications were inherited with some carrier parents reporting learning disabilities or microcephaly. We identified additional, potentially contributory copy number changes in a subset of patients, including one patient each with 16p11.2 deletion and 15q13.3 deletion. Our data further define and expand the clinical spectrum associated with duplications of 17q12 and provide support for the role of genomic modifiers contributing to phenotypic variability.

Analysis of MDM2 Amplification in 43 Endometrial Stromal Tumors: A Potential Diagnostic Pitfall.

MDM2 amplification is known to occur in a variety of neoplasms and its detection by fluorescence in situ hybridization is helpful in distinguishing well-differentiated and dedifferentated liposarcoma from classic lipoma. We recently evaluated a mesenteric mass initially diagnosed as dedifferentiated liposarcoma, largely due to the neoplasm's myxoid morphology and MDM2 expression by immunohistochemistry, from a 46-yr-old woman with a history of uterine low-grade endometrial stromal sarcoma (LG-ESS) with a JAZF1 rearrangement. Our workup of the mesenteric mass revealed a JAZF1 rearrangement and a revised diagnosis of metastatic LG-ESS with myxoid change was rendered. Retrospective testing of the mesenteric mass was negative for MDM2 amplification, an uncommon, but known diagnostic pitfall in MDM2 expression by immunohistochemistry. As MDM2 amplification is not specific for the diagnosis of liposarcoma, we investigated its occurrence in 43 cases of endometrial stromal tumors: 14 uterine LG-ESS, 11 metastatic or recurrent uterine LG-ESS, 8 undifferentiated uterine sarcomas, 5 endometrial stromal nodules, and 4 high-grade ESS with YHWAE rearrangement. In addition, 40 of the 43 cases had previously undergone fluorescence in situ hybridization analysis of JAZF1, PHF1, and YHWAE. Two of the 43 cases (5%) had MDM2 amplification: one was a uterine LG-ESS (JAZF1 rearrangement) and the other was a undifferentiated uterine sarcoma (polysomy intact JAZF1, PHF1, and YHWAE), both metastatic to the lung. Both cases positive for MDM2 amplification showed MDM2 expression by immunohistochemistry. At last follow-up, both patients had died of disease (19 and 60 mo). Our study is the first to demonstrate MDM2 amplification in endometrial stromal tumor. Awareness of MDM2 amplification in endometrial stromal tumor is critical; particularly in locations more common to liposarcoma, to avoid diagnostic errors.

The phenotype spectrum of X-linked ichthyosis identified by chromosomal microarray.

BACKGROUND: Steroid sulfatase (STS) gene disruption causes X-linked ichthyosis (XLI). Interrogating the entire genome through chromosomal microarray (CMA), a test primarily used to screen patients with noncutaneous congenital anomalies, may detect STS deletions incidentally. OBJECTIVE: We sought to determine the variability of skin features associated with STS deletions diagnosed through CMA and to compare these findings with XLI cases reported in the literature and recognized in a dermatology clinic. METHODS: Male patients with an STS deletion were identified from 23,172 consecutive postnatal blood samples tested with CMA at Mayo Clinic. A comparison group of male patients with biochemically confirmed XLI was ascertained in the dermatology clinic. The available patient medical records, skin histopathology, and photographs were evaluated and a literature search of patients with XLI was conducted. RESULTS: Children whose diagnosis was made incidentally through CMA had milder skin phenotypes, including dryness or eczema, or both, and did not manifest the polygonal or "dirty" scale described as typical of XLI in the literature. LIMITATIONS: The small sample size, limited clinical information, and assessment by nondermatologists in a subset of cases may have influenced the results. CONCLUSION: STS deletions may cause a milder skin phenotype than the typical presentation of XLI.

Expression profiling of uterine leiomyomata cytogenetic subgroups reveals distinct signatures in matched myometrium: transcriptional profilingof the t(12;14) and evidence in support of predisposing genetic heterogeneity.

Uterine leiomyomata (UL), the most common neoplasm in reproductive-age women, are classified into distinct genetic subgroups based on recurrent chromosome abnormalities. To develop a molecular signature of UL with t(12;14)(q14-q15;q23-q24), we took advantage of the multiple UL arising as independent clonal lesions within a single uterus. We compared genome-wide expression levels of t(12;14) UL to non-t(12;14) UL from each of nine women in a paired analysis, with each sample weighted for the percentage of t(12;14) cells to adjust for mosaicism with normal cells. This resulted in a transcriptional profile that confirmed HMGA2, known to be overexpressed in t(12;14) UL, as the most significantly altered gene. Pathway analysis of the differentially expressed genes showed significant association with cell proliferation, particularly G1/S checkpoint regulation. This is consistent with the known larger size of t(12;14) UL relative to karyotypically normal UL or to UL in the deletion 7q22 subgroup. Unsupervised hierarchical clustering demonstrated that patient variability is relatively dominant to the distinction of t(12;14) UL compared with non-t(12;14) UL or of t(12;14) UL compared with del(7q) UL. The paired design we employed is therefore important to produce an accurate t(12;14) UL-specific gene list by removing the confounding effects of genotype and environment. Interestingly, myometrium not only clustered away from the tumors, but generally separated based on associated t(12;14) versus del(7q) status. Nine genes were identified whose expression can distinguish the myometrium origin. This suggests an underlying constitutional genetic predisposition to these somatic changes which could potentially lead to improved personalized management and treatment.

Molecular cytogenetic analysis for TFE3 rearrangement in Xp11.2 renal cell carcinoma and alveolar soft part sarcoma: validation and clinical experience with 75 cases.

Renal cell carcinoma with TFE3 rearrangement at Xp11.2 is a distinct subtype manifesting an indolent clinical course in children, with recent reports suggesting a more aggressive entity in adults. This subtype is morphologically heterogeneous and can be misclassified as clear cell or papillary renal cell carcinoma. TFE3 is also rearranged in alveolar soft part sarcoma. To aid in diagnosis, a break-apart strategy fluorescence in situ hybridization (FISH) probe set specific for TFE3 rearrangement and a reflex dual-color, single-fusion strategy probe set involving the most common TFE3 partner gene, ASPSCR1, were validated on formalin-fixed, paraffin-embedded tissues from nine alveolar soft part sarcoma, two suspected Xp11.2 renal cell carcinoma, and nine tumors in the differential diagnosis. The impact of tissue cut artifact was reduced through inclusion of a chromosome X centromere control probe. Analysis of the UOK-109 renal carcinoma cell line confirmed the break-apart TFE3 probe set can distinguish the subtle TFE3/NONO fusion-associated inversion of chromosome X. Subsequent extensive clinical experience was gained through analysis of 75 cases with an indication of Xp11.2 renal cell carcinoma (n=54), alveolar soft part sarcoma (n=13), perivascular epithelioid cell neoplasms (n=2), chordoma (n=1), or unspecified (n=5). We observed balanced and unbalanced chromosome X;17 translocations in both Xp11.2 renal cell carcinoma and alveolar soft part sarcoma, supporting a preference but not a necessity for the translocation to be balanced in the carcinoma and unbalanced in the sarcoma. We further demonstrate the unbalanced separation is atypical, with TFE3/ASPSCR1 fusion and loss of the derivative X chromosome but also an unanticipated normal X chromosome gain in both males and females. Other diverse sex chromosome copy number combinations were observed. Our TFE3 FISH assay is a useful adjunct to morphologic analysis of such challenging cases and will be applicable to assess the growing spectrum of TFE3-rearranged tumors.

Recurrent 8q13.2-13.3 microdeletions associated with branchio-oto-renal syndrome are mediated by human endogenous retroviral (HERV) sequence blocks.

BACKGROUND: Human endogenous retroviral (HERV) sequences are the remnants of ancient retroviral infection and comprise approximately 8% of the human genome. The high abundance and interspersed nature of homologous HERV sequences make them ideal substrates for genomic rearrangements. A role for HERV sequences in mediating human disease-associated rearrangement has been reported but is likely currently underappreciated. METHODS AND RESULTS: In the present study, two independent de novo 8q13.2-13.3 microdeletion events were identified in patients with clinical features of Branchio-Oto-Renal (BOR) syndrome. Nucleotide-level mapping demonstrated the identical breakpoints, suggesting a recurrent microdeletion including multiple genes such as EYA1, SULF1, and SLCO5A1, which is mediated by HERV1 homologous sequences. CONCLUSIONS: These findings raise the potential that HERV sequences may more commonly underlie recombination of dosage sensitive regions associated with recurrent syndromes.

Uterine cellular leiomyomata with chromosome 1p deletions represent a distinct entity.

OBJECTIVE: This study aimed to determine whether 1p deletion defines a subset of cellular leiomyomata (CL), which is a hypercellular variant of uterine leiomyomata that may have delayed malignant potential, and to correlate this genetic change with clinical and pathologic characteristics including those present in uterine sarcomas. STUDY DESIGN: Available CL cases at the Mayo Clinic (n = 101) and variant cases reported in another article (n = 16) were identified. Each case with sufficient tissue that met histologic criteria for CL when reviewed by a single pathologist underwent interphase fluorescence in situ hybridization to determine the presence of 1p deletion. Clinical characteristics of women with confirmed CL were compared on the basis of 1p deletion status using univariate analysis. RESULTS: Of the Mayo Clinic cohort of histologically confirmed CL, 23% had deletion of 1p. Women with this subset of CL, when compared to those without 1p deletion, were more likely to be postmenopausal (P = .049) and their uteri tended to be heavier (P = .039) with a larger dominant leiomyoma (P = .030). The pathologic features associated with 1p deletion were high cellularity (P = .036) and hyaline necrosis (P = .047), which remained significant after inclusion of the CL cases from a previously published series. CONCLUSION: Deletion of 1p occurs in approximately one-quarter of CL cases. This genetic alteration is potentially associated with clinicopathologic features that are present in uterine sarcomas, which suggests a distinct clinical entity that may have malignant potential. Our findings are particularly pertinent considering the increased preference for uterine-sparing options in leiomyoma treatment, suggesting assessment of 1p deletion status in CL may influence clinical surveillance decisions.

Deletions in 16q24.2 are associated with autism spectrum disorder, intellectual disability and congenital renal malformation.

BACKGROUND: The contribution of copy-number variation (CNV) to disease has been highlighted with the widespread adoption of array-based comparative genomic hybridisation (aCGH) and microarray technology. Contiguous gene deletions involving ANKRD11 in 16q24.3 are associated with autism spectrum disorder (ASD) and intellectual disability (ID), while 16q24.1 deletions affecting FOXF1 are associated with congenital renal malformations, alveolar capillary dysplasia, and various other abnormalities. The disease associations of deletions in the intervening region, 16q24.2, have only been defined to a limited extent. AIM: To determine whether deletions affecting 16q24.2 are correlated with congenital anomalies. METHODS: 35 individuals, each having a deletion in 16q24.2, were characterised clinically and by aCGH and/or SNP-genotyping microarray. RESULTS: Several of the 35 16q24.2 deletions identified here closely abut or overlap the coding regions of FOXF1 and ANKRD11, two genes that have been previously associated with the disease. 25 patients were reported to have ASD/ID, and three were found to have bilateral hydronephrosis. 14 of the deletions associated with ASD/ID overlap the coding regions of FBXO31 and MAP1LC3B. These same genes and two others, C16orf95 and ZCCHC14, are also included in the area of minimal overlap of the three deletions associated with hydronephrosis. CONCLUSIONS: Our data highlight 16q24.2 as a region of interest for ASD, ID and congenital renal malformations. These conditions are associated, albeit without complete penetrance, with deletions affecting C16orf95, ZCCHC14, MAP1LC3B and FBXO31. The function of each gene in development and disease warrants further investigation.

Development of a Multifaceted Program for Pharmacogenetics Adoption at an Academic Medical Center: Practical Considerations and Lessons Learned.

In 2019, Indiana University launched the Precision Health Initiative to enhance the institutional adoption of precision medicine, including pharmacogenetics (PGx) implementation, at university-affiliated practice sites across Indiana. The overarching goal of this PGx implementation program was to facilitate the sustainable adoption of genotype-guided prescribing into routine clinical care. To accomplish this goal, we pursued the following specific objectives: (i) to integrate PGx testing into existing healthcare system processes; (ii) to implement drug-gene pairs with high-level evidence and educate providers and pharmacists on established clinical management recommendations; (iii) to engage key stakeholders, including patients to optimize the return of results for PGx testing; (iv) to reduce health disparities through the targeted inclusion of underrepresented populations; (v) and to track third-party reimbursement. This tutorial details our multifaceted PGx implementation program, including descriptions of our interventions, the critical challenges faced, and the major program successes. By describing our experience, we aim to assist other clinical teams in achieving sustainable PGx implementation in their health systems.

A Multicenter Analysis of Abnormal Chromosomal Microarray Findings in Congenital Heart Disease.

Background Chromosomal microarray analysis (CMA) provides an opportunity to understand genetic causes of congenital heart disease (CHD). The methods for describing cardiac phenotypes in patients with CMA abnormalities have been inconsistent, which may complicate clinical interpretation of abnormal testing results and hinder a more complete understanding of genotype-phenotype relationships. Methods and Results Patients with CHD and abnormal clinical CMA were accrued from 9 pediatric cardiac centers. Highly detailed cardiac phenotypes were systematically classified and analyzed for their association with CMA abnormality. Hierarchical classification of each patient into 1 CHD category facilitated broad analyses. Inclusive classification allowing multiple CHD types per patient provided sensitive descriptions. In 1363 registry patients, 28% had genomic disorders with well-recognized CHD association, 67% had clinically reported copy number variants (CNVs) with rare or no prior CHD association, and 5% had regions of homozygosity without CNV. Hierarchical classification identified expected CHD categories in genomic disorders, as well as uncharacteristic CHDs. Inclusive phenotyping provided sensitive descriptions of patients with multiple CHD types, which occurred commonly. Among CNVs with rare or no prior CHD association, submicroscopic CNVs were enriched for more complex types of CHD compared with large CNVs. The submicroscopic CNVs that contained a curated CHD gene were enriched for left ventricular obstruction or septal defects, whereas CNVs containing a single gene were enriched for conotruncal defects. Neuronal-related pathways were over-represented in single-gene CNVs, including top candidate causative genes NRXN3, ADCY2, and HCN1. Conclusions Intensive cardiac phenotyping in multisite registry data identifies genotype-phenotype associations in CHD patients with abnormal CMA.

Computational pharmacogenotype extraction from clinical next-generation sequencing.

Background: Next-generation sequencing (NGS), including whole genome sequencing (WGS) and whole exome sequencing (WES), is increasingly being used for clinic care. While NGS data have the potential to be repurposed to support clinical pharmacogenomics (PGx), current computational approaches have not been widely validated using clinical data. In this study, we assessed the accuracy of the Aldy computational method to extract PGx genotypes from WGS and WES data for 14 and 13 major pharmacogenes, respectively. Methods: Germline DNA was isolated from whole blood samples collected for 264 patients seen at our institutional molecular solid tumor board. DNA was used for panel-based genotyping within our institutional Clinical Laboratory Improvement Amendments- (CLIA-) certified PGx laboratory. DNA was also sent to other CLIA-certified commercial laboratories for clinical WGS or WES. Aldy v3.3 and v4.4 were used to extract PGx genotypes from these NGS data, and results were compared to the panel-based genotyping reference standard that contained 45 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, TPMT, and VKORC1. Results: Mean WGS read depth was >30x for all variant regions except for G6PD (average read depth was 29 reads), and mean WES read depth was >30x for all variant regions. For 94 patients with WGS, Aldy v3.3 diplotype calls were concordant with those from the genotyping reference standard in 99.5% of cases when excluding diplotypes with additional major star alleles not tested by targeted genotyping, ambiguous phasing, and CYP2D6 hybrid alleles. Aldy v3.3 identified 15 additional clinically actionable star alleles not covered by genotyping within CYP2B6, CYP2C19, DPYD, SLCO1B1, and NUDT15. Within the WGS cohort, Aldy v4.4 diplotype calls were concordant with those from genotyping in 99.7% of cases. When excluding patients with CYP2D6 copy number variation, all Aldy v4.4 diplotype calls except for one CYP3A4 diplotype call were concordant with genotyping for 161 patients in the WES cohort. Conclusion: Aldy v3.3 and v4.4 called diplotypes for major pharmacogenes from clinical WES and WGS data with >99% accuracy. These findings support the use of Aldy to repurpose clinical NGS data to inform clinical PGx.