Determining Cost-Optimal Next-Generation Sequencing Panels for Rare Disease and Pharmacogenomics Testing.

BACKGROUND: Multi-gene panel sequencing using next-generation sequencing (NGS) methods is a key tool for genomic medicine. However, with an estimated 140 000 genomic tests available, current system inefficiencies result in high genetic-testing costs. Reduced testing costs are needed to expand the availability of genomic medicine. One solution to improve efficiency and lower costs is to calculate the most cost-effective set of panels for a typical pattern of test requests. METHODS: We compiled rare diseases, associated genes, point prevalence, and test-order frequencies from a representative laboratory. We then modeled the costs of the relevant steps in the NGS process in detail. Using a simulated annealing-based optimization procedure, we determined panel sets that were more cost-optimal than whole exome sequencing (WES) or clinical exome sequencing (CES). Finally, we repeated this methodology to cost-optimize pharmacogenomics (PGx) testing. RESULTS: For rare disease testing, we show that an optimal choice of 4-6 panels, uniquely covering genes that comprise 95% of the total prevalence of monogenic diseases, saves $257-304 per sample compared with WES, and $66-135 per sample compared with CES. For PGx, we show that the optimal multipanel solution saves $6-7 (27%-40%) over a single panel covering all relevant gene-drug associations. CONCLUSIONS: Laboratories can reduce costs using the proposed method to obtain and run a cost-optimal set of panels for specific test requests. In addition, payers can use this method to inform reimbursement policy.

De novo, deleterious sequence variants that alter the transcriptional activity of the homeoprotein PBX1 are associated with intellectual disability and pleiotropic developmental defects.

We present eight patients with de novo, deleterious sequence variants in the PBX1 gene. PBX1 encodes a three amino acid loop extension (TALE) homeodomain transcription factor that forms multimeric complexes with TALE and HOX proteins to regulate target gene transcription during development. As previously reported, Pbx1 homozygous mutant mice (Pbx1-/-) develop malformations and hypoplasia or aplasia of multiple organs, including the craniofacial skeleton, ear, branchial arches, heart, lungs, diaphragm, gut, kidneys, and gonads. Clinical findings similar to those in Pbx mutant mice were observed in all patients with varying expressivity and severity, including external ear anomalies, abnormal branchial arch derivatives, heart malformations, diaphragmatic hernia, renal hypoplasia and ambiguous genitalia. All patients but one had developmental delays. Previously reported patients with congenital anomalies affecting the kidney and urinary tract exhibited deletions and loss of function variants in PBX1. The sequence variants in our cases included missense substitutions adjacent to the PBX1 homeodomain (p.Arg184Pro, p.Met224Lys, and p.Arg227Pro) or within the homeodomain (p.Arg234Pro, and p.Arg235Gln), whereas p.Ser262Glnfs*2, and p.Arg288* yielded truncated PBX1 proteins. Functional studies on five PBX1 sequence variants revealed perturbation of intrinsic, PBX-dependent transactivation ability and altered nuclear translocation, suggesting abnormal interactions between mutant PBX1 proteins and wild-type TALE or HOX cofactors. It is likely that the mutations directly affect the transcription of PBX1 target genes to impact embryonic development. We conclude that deleterious sequence variants in PBX1 cause intellectual disability and pleiotropic malformations resembling those in Pbx1 mutant mice, arguing for strong conservation of gene function between these two species.

Germline De Novo Mutations in GNB1 Cause Severe Neurodevelopmental Disability, Hypotonia, and Seizures.

Whole-exome sequencing of 13 individuals with developmental delay commonly accompanied by abnormal muscle tone and seizures identified de novo missense mutations enriched within a sub-region of GNB1, a gene encoding the guanine nucleotide-binding protein subunit beta-1, Gß. These 13 individuals were identified among a base of 5,855 individuals recruited for various undiagnosed genetic disorders. The probability of observing 13 or more de novo mutations by chance among 5,855 individuals is very low (p = 7.1 × 10(-21)), implicating GNB1 as a genome-wide-significant disease-associated gene. The majority of these 13 mutations affect known Gß binding sites, which suggests that a likely disease mechanism is through the disruption of the protein interface required for Gα-Gßγ interaction (resulting in a constitutively active Gßγ) or through the disruption of residues relevant for interaction between Gßγ and certain downstream effectors (resulting in reduced interaction with the effectors). Strikingly, 8 of the 13 individuals recruited here for a neurodevelopmental disorder have a germline de novo GNB1 mutation that overlaps a set of five recurrent somatic tumor mutations for which recent functional studies demonstrated a gain-of-function effect due to constitutive activation of G protein downstream signaling cascades for some of the affected residues.

TANGO2: expanding the clinical phenotype and spectrum of pathogenic variants.

PURPOSE: TANGO2-related disorders were first described in 2016 and prior to this publication, only 15 individuals with TANGO2-related disorder were described in the literature. Primary features include metabolic crisis with rhabdomyolysis, encephalopathy, intellectual disability, seizures, and cardiac arrhythmias. We assess whether genotype and phenotype of TANGO2-related disorder has expanded since the initial discovery and determine the efficacy of exome sequencing (ES) as a diagnostic tool for detecting variants. METHODS: We present a series of 14 individuals from 11 unrelated families with complex medical and developmental histories, in whom ES or microarray identified compound heterozygous or homozygous variants in TANGO2. RESULTS: The initial presentation of patients with TANGO2-related disorders can be variable, including primarily neurological presentations. We expand the phenotype and genotype for TANGO2, highlighting the variability of the disorder. CONCLUSION: TANGO2-related disorders can have a more diverse clinical presentation than previously anticipated. We illustrate the utility of routine ES data reanalysis whereby discovery of novel disease genes can lead to a diagnosis in previously unsolved cases and the need for additional copy-number variation analysis when ES is performed.

Correction: TANGO2: expanding the clinical phenotype and spectrum of pathogenic variants.

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

Spectrin mutations cause spinocerebellar ataxia type 5.

We have discovered that beta-III spectrin (SPTBN2) mutations cause spinocerebellar ataxia type 5 (SCA5) in an 11-generation American kindred descended from President Lincoln's grandparents and two additional families. Two families have separate in-frame deletions of 39 and 15 bp, and a third family has a mutation in the actin/ARP1 binding region. Beta-III spectrin is highly expressed in Purkinje cells and has been shown to stabilize the glutamate transporter EAAT4 at the surface of the plasma membrane. We found marked differences in EAAT4 and GluRdelta2 by protein blot and cell fractionation in SCA5 autopsy tissue. Cell culture studies demonstrate that wild-type but not mutant beta-III spectrin stabilizes EAAT4 at the plasma membrane. Spectrin mutations are a previously unknown cause of ataxia and neurodegenerative disease that affect membrane proteins involved in glutamate signaling.

A review of genetic counseling for Charcot Marie Tooth disease (CMT).

Charcot Marie Tooth disease (CMT) encompasses the inherited peripheral neuropathies. While four genes have been found to cause over 90 % of genetically identifiable causes of CMT (PMP22, GJB1, MPZ, MFN2), at least 51 genes and loci have been found to cause CMT when mutated, creating difficulties for clinicians to find a genetic subtype for families. Here, the classic features of CMT as well as characteristic features of the most common subtypes of CMT are described, as well as methods for narrowing down the possible subtypes. Psychosocial concerns particular to the CMT population are identified. This is the most inclusive publication for CMT-specific genetic counseling.

Comparison of Germline Genetic Testing Before and After a Medical Policy Covering Universal Testing Among Patients With Colorectal Cancer.

Importance: In 2020, some health insurance plans updated their medical policy to cover germline genetic testing for all patients diagnosed with colorectal cancer (CRC). Guidelines for universal tumor screening via microsatellite instability and/or immunohistochemistry (MSI/IHC) for mismatch repair protein expression for patients with CRC have been in place since 2009. Objectives: To examine whether uptake of MSI/IHC screening and germline genetic testing in patients with CRC has improved under these policies and to identify actionable findings and management implications for patients referred for germline genetic testing. Design, Setting, and Participants: The multicenter, retrospective cohort study comprised 2 analyses of patients 18 years or older who were diagnosed with CRC between January 1, 2017, and December 31, 2020. The first analysis used an insurance claims data set to examine use of MSI/IHC screening and germline genetic testing for patients diagnosed with CRC between 2017 and 2020 and treated with systemic therapy. The second comprised patients with CRC who had germline genetic testing performed in 2020 that was billed under a universal testing policy. Main Outcomes and Measures: Patient demographic characteristics, clinical information, and use of MSI/IHC screening and germline genetic testing were analyzed. Results: For 9066 patients with newly diagnosed CRC (mean [SD] age, 64.2 [12.7] years; 4964 [54.8%] male), administrative claims data indicated that MSI/IHC was performed in 6645 eligible patients (73.3%) during the study period, with 2288 (25.2%) not receiving MSI/IHC despite being eligible for coverage. Analysis of a second cohort of 55 595 patients with CRC diagnosed in 2020 and covered by insurance found that only 1675 (3.0%) received germline genetic testing. In a subset of patients for whom germline genetic testing results were available, 1 in 6 patients had pathogenic or likely pathogenic variants, with most of these patients having variants with established clinical actionability. Conclusions and Relevance: This nationwide cohort study found suboptimal rates of MSI/IHC screening and germline genetic testing uptake, resulting in clinically actionable genetic data being unavailable to patients diagnosed with CRC, despite universal eligibility. Effective strategies are required to address barriers to implementation of evidence-based universal testing policies that support precision treatment and optimal care management for patients with CRC.

Cobblestone Malformation in LAMA2 Congenital Muscular Dystrophy (MDC1A).

Congenital muscular dystrophy type 1A (MDC1A) is caused by recessive variants in laminin α2 (LAMA2). Patients have been found to have white matter signal abnormalities on magnetic resonance imaging (MRI) but rarely structural brain abnormalities. We describe the autopsy neuropathology in a 17-year-old with white matter signal abnormalities on brain MRI. Dystrophic pathology was observed in skeletal muscle, and the sural nerve manifested a mild degree of segmental demyelination and remyelination. A diffuse, bilateral cobblestone appearance, and numerous points of fusion between adjacent gyri were apparent on gross examination of the cerebrum. Brain histopathology included focal disruptions of the glia limitans associated with abnormal cerebral cortical lamination or arrested cerebellar granule cell migration. Subcortical nodular heterotopia was present within the cerebellar hemispheres. Sampling of the centrum semiovale revealed no light microscopic evidence of leukoencephalopathy. Three additional MDC1A patients were diagnosed with cobblestone malformation on brain MRI. Unlike the autopsied patient whose brain had a symmetric distribution of cobblestone pathology, the latter patients had asymmetric involvement, most severe in the occipital lobes. These cases demonstrate that cobblestone malformation may be an important manifestation of the brain pathology in MDC1A and can be present even when patients have a structurally normal brain MRI.

The genetic landscape of axonal neuropathies in the middle-aged and elderly: Focus on MME.

OBJECTIVE: To test the hypothesis that monogenic neuropathies such as Charcot-Marie-Tooth disease (CMT) contribute to frequent but often unexplained neuropathies in the elderly, we performed genetic analysis of 230 patients with unexplained axonal neuropathies and disease onset ≥35 years. METHODS: We recruited patients, collected clinical data, and conducted whole-exome sequencing (WES; n = 126) and MME single-gene sequencing (n = 104). We further queried WES repositories for MME variants and measured blood levels of the MME-encoded protein neprilysin. RESULTS: In the WES cohort, the overall detection rate for assumed disease-causing variants in genes for CMT or other conditions associated with neuropathies was 18.3% (familial cases 26.4%, apparently sporadic cases 12.3%). MME was most frequently involved and accounted for 34.8% of genetically solved cases. The relevance of MME for late-onset neuropathies was further supported by detection of a comparable proportion of cases in an independent patient sample, preponderance of MME variants among patients compared to population frequencies, retrieval of additional late-onset neuropathy patients with MME variants from WES repositories, and low neprilysin levels in patients' blood samples. Transmission of MME variants was often consistent with an incompletely penetrant autosomal-dominant trait and less frequently with autosomal-recessive inheritance. CONCLUSIONS: A detectable fraction of unexplained late-onset axonal neuropathies is genetically determined, by variants in either CMT genes or genes involved in other conditions that affect the peripheral nerves and can mimic a CMT phenotype. MME variants can act as completely penetrant recessive alleles but also confer dominantly inherited susceptibility to axonal neuropathies in an aging population.

Mutational spectrum of DMD mutations in dystrophinopathy patients: application of modern diagnostic techniques to a large cohort.

Mutations in the DMD gene, encoding the dystrophin protein, are responsible for the dystrophinopathies Duchenne Muscular Dystrophy (DMD), Becker Muscular Dystrophy (BMD), and X-linked Dilated Cardiomyopathy (XLDC). Mutation analysis has traditionally been challenging, due to the large gene size (79 exons over 2.2 Mb of genomic DNA). We report a very large aggregate data set comprised of DMD mutations detected in samples from patients enrolled in the United Dystrophinopathy Project, a multicenter research consortium, and in referral samples submitted for mutation analysis with a diagnosis of dystrophinopathy. We report 1,111 mutations in the DMD gene, including 891 mutations with associated phenotypes. These results encompass 506 point mutations (including 294 nonsense mutations) and significantly expand the number of mutations associated with the dystrophinopathies, highlighting the utility of modern diagnostic techniques. Our data supports the uniform hypermutability of CGA>TGA mutations, establishes the frequency of polymorphic muscle (Dp427m) protein isoforms and reveals unique genomic haplotypes associated with "private" mutations. We note that 60% of these patients would be predicted to benefit from skipping of a single DMD exon using antisense oligonucleotide therapy, and 62% would be predicted to benefit from an inclusive multiexonskipping approach directed toward exons 45 through 55.

MBNL Sequestration by Toxic RNAs and RNA Misprocessing in the Myotonic Dystrophy Brain.

For some neurological disorders, disease is primarily RNA mediated due to expression of non-coding microsatellite expansion RNAs (RNA(exp)). Toxicity is thought to result from enhanced binding of proteins to these expansions and depletion from their normal cellular targets. However, experimental evidence for this sequestration model is lacking. Here, we use HITS-CLIP and pre-mRNA processing analysis of human control versus myotonic dystrophy (DM) brains to provide compelling evidence for this RNA toxicity model. MBNL2 binds directly to DM repeat expansions in the brain, resulting in depletion from its normal RNA targets with downstream effects on alternative splicing and polyadenylation. Similar RNA processing defects were detected in Mbnl compound-knockout mice, highlighted by dysregulation of Mapt splicing and fetal tau isoform expression in adults. These results demonstrate that MBNL proteins are directly sequestered by RNA(exp) in the DM brain and introduce a powerful experimental tool to evaluate RNA-mediated toxicity in other expansion diseases.

A focal domain of extreme demethylation within D4Z4 in FSHD2.

OBJECTIVE: Facioscapulohumeral muscular dystrophy (FSHD) is a neuromuscular disease with an unclear genetic mechanism. Most patients have a contraction of the D4Z4 macrosatellite repeat array at 4qter, which is thought to cause partial demethylation (FSHD1) of the contracted allele. Demethylation has been surveyed at 3 restriction enzyme sites in the first repeat and only a single site across the entire array, and current models postulate that a generalized D4Z4 chromatin alteration causes FSHD. The background of normal alleles has confounded the study of epigenetic alterations; however, rare patients (FSHD2) have a form of the disease in which demethylation is global, i.e., on all D4Z4 elements throughout the genome. Our objective was to take advantage of the global nature of FSHD2 to identify where disease-relevant methylation changes occur within D4Z4. METHODS: Using bisulfite sequencing of DNA from blood and myoblast cells, methylation levels at 74 CpG sites across 3 disparate regions within D4Z4 were measured in FSHD2 patients and controls. RESULTS: We found that rates of demethylation caused by FSHD2 are not consistent across D4Z4. We identified a focal region of extreme demethylation within a 5' domain, which we named DR1. Other D4Z4 regions, including the DUX4 ORF, were hypomethylated but to a much lesser extent. CONCLUSIONS: These data challenge the simple view that FSHD is caused by a broad "opening" of D4Z4 and lead us to postulate that the region of focal demethylation is the site of action of the key D4Z4 chromatin regulatory factors that go awry in FSHD.

Motor and cognitive assessment of infants and young boys with Duchenne Muscular Dystrophy: results from the Muscular Dystrophy Association DMD Clinical Research Network.

Therapeutic trials in Duchenne Muscular Dystrophy (DMD) exclude young boys because traditional outcome measures rely on cooperation. The Bayley III Scales of Infant and Toddler Development (Bayley III) have been validated in developing children and those with developmental disorders but have not been studied in DMD. Expanded Hammersmith Functional Motor Scale (HFMSE) and North Star Ambulatory Assessment (NSAA) may also be useful in this young DMD population. Clinical evaluators from the MDA-DMD Clinical Research Network were trained in these assessment tools. Infants and boys with DMD (n = 24; 1.9 ± 0.7 years) were assessed. The mean Bayley III motor composite score was low (82.8 ± 8; p ≤ .0001) (normal = 100 ± 15). Mean gross motor and fine motor function scaled scores were low (both p ≤ .0001). The mean cognitive comprehensive (p=.0002), receptive language (p ≤ .0001), and expressive language (p = .0001) were also low compared to normal children. Age was negatively associated with Bayley III gross motor (r = -0.44; p = .02) but not with fine motor, cognitive, or language scores. HFMSE (n=23) showed a mean score of 31 ± 13. NSAA (n = 18 boys; 2.2 ± 0.4 years) showed a mean score of 12 ± 5. Outcome assessments of young boys with DMD are feasible and in this multicenter study were best demonstrated using the Bayley III.

Cerebral and muscle MRI abnormalities in myotonic dystrophy.

Pathophysiological mechanisms underlying the clinically devastating CNS features of myotonic dystrophy (DM) remain more enigmatic and controversial than do the muscle abnormalities of this common form of muscular dystrophy. To better define CNS and cranial muscle changes in DM, we used quantitative volumetric and diffusion tensor MRI methods to measure cerebral and masticatory muscle differences between controls (n=5) and adults with either congenital (n=5) or adult onset (n=5) myotonic dystrophy type 1 and myotonic dystrophy type 2 (n=5). Muscle volumes were diminished in DM1 and strongly correlated with reduced white matter integrity and gray matter volume. Moreover, correlation of reduced fractional anisotropy (white matter integrity) and gray matter volume in both DM1 and DM2 suggests that these abnormalities may share a common underlying pathophysiological mechanism. Further quantitative temporal and spatial characterization of these features will help delineate developmental and progressive neurological components of DM, and help determine the causative molecular and cellular mechanisms.

SNP haplotype mapping in a small ALS family.

The identification of genes for monogenic disorders has proven to be highly effective for understanding disease mechanisms, pathways and gene function in humans. Nevertheless, while thousands of Mendelian disorders have not yet been mapped there has been a trend away from studying single-gene disorders. In part, this is due to the fact that many of the remaining single-gene families are not large enough to map the disease locus to a single site in the genome. New tools and approaches are needed to allow researchers to effectively tap into this genetic gold-mine. Towards this goal, we have used haploid cell lines to experimentally validate the use of high-density single nucleotide polymorphism (SNP) arrays to define genome-wide haplotypes and candidate regions, using a small amyotrophic lateral sclerosis (ALS) family as a prototype. Specifically, we used haploid-cell lines to determine if high-density SNP arrays accurately predict haplotypes across entire chromosomes and show that haplotype information significantly enhances the genetic information in small families. Panels of haploid-cell lines were generated and a 5 centimorgan (cM) short tandem repeat polymorphism (STRP) genome scan was performed. Experimentally derived haplotypes for entire chromosomes were used to directly identify regions of the genome identical-by-descent in 5 affected individuals. Comparisons between experimentally determined and in silico haplotypes predicted from SNP arrays demonstrate that SNP analysis of diploid DNA accurately predicted chromosomal haplotypes. These methods precisely identified 12 candidate intervals, which are shared by all 5 affected individuals. Our study illustrates how genetic information can be maximized using readily available tools as a first step in mapping single-gene disorders in small families.

Myotonic dystrophy type 2: human founder haplotype and evolutionary conservation of the repeat tract.

Myotonic dystrophy (DM), the most common form of muscular dystrophy in adults, can be caused by a mutation on either chromosome 19 (DM1) or 3 (DM2). In 2001, we demonstrated that DM2 is caused by a CCTG expansion in intron 1 of the zinc finger protein 9 (ZNF9) gene. To investigate the ancestral origins of the DM2 expansion, we compared haplotypes for 71 families with genetically confirmed DM2, using 19 short tandem repeat markers that we developed that flank the repeat tract. All of the families are white, with the majority of Northern European/German descent and a single family from Afghanistan. Several conserved haplotypes spanning >700 kb appear to converge into a single haplotype near the repeat tract. The common interval that is shared by all families with DM2 immediately flanks the repeat, extending up to 216 kb telomeric and 119 kb centromeric of the CCTG expansion. The DM2 repeat tract contains the complex repeat motif (TG)(n)(TCTG)(n)(CCTG)(n). The CCTG portion of the repeat tract is interrupted on normal alleles, but, as in other expansion disorders, these interruptions are lost on affected alleles. We examined haplotypes of 228 control chromosomes and identified a potential premutation allele with an uninterrupted (CCTG)(20) on a haplotype that was identical to the most common affected haplotype. Our data suggest that the predominant Northern European ancestry of families with DM2 resulted from a common founder and that the loss of interruptions within the CCTG portion of the repeat tract may predispose alleles to further expansion. To gain insight into possible function of the repeat tract, we looked for evolutionary conservation. The complex repeat motif and flanking sequences within intron 1 are conserved among human, chimpanzee, gorilla, mouse, and rat, suggesting a conserved biological function.

Spinocerebellar ataxia type 8: molecular genetic comparisons and haplotype analysis of 37 families with ataxia.

We reported elsewhere that an untranslated CTG expansion causes the dominantly inherited neurodegenerative disorder spinocerebellar ataxia type 8 (SCA8). SCA8 shows a complex inheritance pattern with extremes of incomplete penetrance, in which often only one or two affected individuals are found in a given family. SCA8 expansions have also been found in control chromosomes, indicating that separate genetic or environmental factors increase disease penetrance among SCA8-expansion-carrying patients with ataxia. We describe the molecular genetic features and disease penetrance of 37 different families with SCA8 ataxia from the United States, Canada, Japan, and Mexico. Haplotype analysis using 17 STR markers spanning an approximately 1-Mb region was performed on the families with ataxia, on a group of expansion carriers in the general population, and on psychiatric patients, to clarify the genetic basis of the reduced penetrance and to investigate whether CTG expansions among different populations share a common ancestral background. Two major ancestrally related haplotypes (A and A') were found among white families with ataxia, normal controls, and patients with major psychosis, indicating a common ancestral origin of both pathogenic and nonpathogenic SCA8 expansions among whites. Two additional and distinct haplotypes were found among a group of Japanese families with ataxia (haplotype B) and a Mexican family with ataxia (haplotype C). Our finding that SCA8 expansions on three independently arising haplotypes are found among patients with ataxia and cosegregate with ataxia when multiple family members are affected further supports the direct role of the CTG expansion in disease pathogenesis.