Genome Sequencing for Diagnosing Rare Diseases.

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) - 8.2% of families in the initial cohort - had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.).

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly.

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

Critical assessment of variant prioritization methods for rare disease diagnosis within the rare genomes project.

BACKGROUND: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting. METHODS: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values. RESULTS: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency. CONCLUSIONS: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed.

Unique variants in CLCN3, encoding an endosomal anion/proton exchanger, underlie a spectrum of neurodevelopmental disorders.

The genetic causes of global developmental delay (GDD) and intellectual disability (ID) are diverse and include variants in numerous ion channels and transporters. Loss-of-function variants in all five endosomal/lysosomal members of the CLC family of Cl- channels and Cl-/H+ exchangers lead to pathology in mice, humans, or both. We have identified nine variants in CLCN3, the gene encoding CIC-3, in 11 individuals with GDD/ID and neurodevelopmental disorders of varying severity. In addition to a homozygous frameshift variant in two siblings, we identified eight different heterozygous de novo missense variants. All have GDD/ID, mood or behavioral disorders, and dysmorphic features; 9/11 have structural brain abnormalities; and 6/11 have seizures. The homozygous variants are predicted to cause loss of ClC-3 function, resulting in severe neurological disease similar to the phenotype observed in Clcn3-/- mice. Their MRIs show possible neurodegeneration with thin corpora callosa and decreased white matter volumes. Individuals with heterozygous variants had a range of neurodevelopmental anomalies including agenesis of the corpus callosum, pons hypoplasia, and increased gyral folding. To characterize the altered function of the exchanger, electrophysiological analyses were performed in Xenopus oocytes and mammalian cells. Two variants, p.Ile607Thr and p.Thr570Ile, had increased currents at negative cytoplasmic voltages and loss of inhibition by luminal acidic pH. In contrast, two other variants showed no significant difference in the current properties. Overall, our work establishes a role for CLCN3 in human neurodevelopment and shows that both homozygous loss of ClC-3 and heterozygous variants can lead to GDD/ID and neuroanatomical abnormalities.

KCND2 variants associated with global developmental delay differentially impair Kv4.2 channel gating.

Here, we report on six unrelated individuals, all presenting with early-onset global developmental delay, associated with impaired motor, speech and cognitive development, partly with developmental epileptic encephalopathy and physical dysmorphisms. All individuals carry heterozygous missense variants of KCND2, which encodes the voltage-gated potassium (Kv) channel α-subunit Kv4.2. The amino acid substitutions associated with the variants, p.(Glu323Lys) (E323K), p.(Pro403Ala) (P403A), p.(Val404Leu) (V404L) and p.(Val404Met) (V404M), affect sites known to be critical for channel gating. To unravel their likely pathogenicity, recombinant mutant channels were studied in the absence and presence of auxiliary ß-subunits under two-electrode voltage clamp in Xenopus oocytes. All channel mutants exhibited slowed and incomplete macroscopic inactivation, and the P403A variant in addition slowed activation. Co-expression of KChIP2 or DPP6 augmented the functional expression of both wild-type and mutant channels; however, the auxiliary ß-subunit-mediated gating modifications differed from wild type and among mutants. To simulate the putative setting in the affected individuals, heteromeric Kv4.2 channels (wild type + mutant) were studied as ternary complexes (containing both KChIP2 and DPP6). In the heteromeric ternary configuration, the E323K variant exhibited only marginal functional alterations compared to homomeric wild-type ternary, compatible with mild loss-of-function. By contrast, the P403A, V404L and V404M variants displayed strong gating impairment in the heteromeric ternary configuration, compatible with loss-of-function or gain-of-function. Our results support the etiological involvement of Kv4.2 channel gating impairment in early-onset monogenic global developmental delay. In addition, they suggest that gain-of-function mechanisms associated with a substitution of V404 increase epileptic seizure susceptibility.

DNM1L variant presenting as adolescent-onset sensory neuronopathy, spasticity, dystonia, and ataxia.

DMN1L encodes for dynamin-like protein 1 (DLP1) which plays a key role in perixosomal and mitochondrial fission. Individuals with heterozygous variants in DNM1L present with a wide range of neurologic symptoms, including encephalopathy, epilepsy, and motor deficits. Here we report on a woman presenting with adolescence onset of sensory neuronopathy, spasticity, dystonia, and ataxia. Trio genome sequencing identified a heterozygous variant in DNM1L (NM_012062.3 c.121G>A/p.Val41Met) which was thought to be pathogenic. This case describes the latest known symptomatic onset of DMN1L-related disease described in literature. We highlight our approach to a challenging diagnostic workup and interpretation of a specific variant that has not been previously reported. Furthermore, the case highlights the diagnostic importance of utilizing genomic sequencing and research studies for patients with rare disease.

Delineation of a novel neurodevelopmental syndrome associated with PAX5 haploinsufficiency.

PAX5 is a transcription factor associated with abnormal posterior midbrain and cerebellum development in mice. PAX5 is highly loss-of-function intolerant and missense constrained, and has been identified as a candidate gene for autism spectrum disorder (ASD). We describe 16 individuals from 12 families who carry deletions involving PAX5 and surrounding genes, de novo frameshift variants that are likely to trigger nonsense-mediated mRNA decay, a rare stop-gain variant, or missense variants that affect conserved amino acid residues. Four of these individuals were published previously but without detailed clinical descriptions. All these individuals have been diagnosed with one or more neurodevelopmental phenotypes including delayed developmental milestones (DD), intellectual disability (ID), and/or ASD. Seizures were documented in four individuals. No recurrent patterns of brain magnetic resonance imaging (MRI) findings, structural birth defects, or dysmorphic features were observed. Our findings suggest that PAX5 haploinsufficiency causes a neurodevelopmental disorder whose cardinal features include DD, variable ID, and/or ASD.

seqr: A web-based analysis and collaboration tool for rare disease genomics.

Exome and genome sequencing have become the tools of choice for rare disease diagnosis, leading to large amounts of data available for analyses. To identify causal variants in these datasets, powerful filtering and decision support tools that can be efficiently used by clinicians and researchers are required. To address this need, we developed seqr - an open-source, web-based tool for family-based monogenic disease analysis that allows researchers to work collaboratively to search and annotate genomic callsets. To date, seqr is being used in several research pipelines and one clinical diagnostic lab. In our own experience through the Broad Institute Center for Mendelian Genomics, seqr has enabled analyses of over 10,000 families, supporting the diagnosis of more than 3,800 individuals with rare disease and discovery of over 300 novel disease genes. Here, we describe a framework for genomic analysis in rare disease that leverages seqr's capabilities for variant filtration, annotation, and causal variant identification, as well as support for research collaboration and data sharing. The seqr platform is available as open source software, allowing low-cost participation in rare disease research, and a community effort to support diagnosis and gene discovery in rare disease.

De novo variants in MPP5 cause global developmental delay and behavioral changes.

Membrane Protein Palmitoylated 5 (MPP5) is a highly conserved apical complex protein essential for cell polarity, fate and survival. Defects in cell polarity are associated with neurologic disorders including autism and microcephaly. MPP5 is essential for neurogenesis in animal models, but human variants leading to neurologic impairment have not been described. We identified three patients with heterozygous MPP5 de novo variants (DNV) and global developmental delay (GDD) and compared their phenotypes and magnetic resonance imaging (MRI) to ascertain how MPP5 DNV leads to GDD. All three patients with MPP5 DNV experienced GDD with language delay/regression and behavioral changes. MRI ranged from normal to decreased gyral folding and microcephaly. The effects of MPP5 depletion on the developing brain were assessed by creating a heterozygous conditional knock out (het CKO) murine model with central nervous system (CNS)-specific Nestin-Cre drivers. In the het CKO model, Mpp5 depletion led to microcephaly, decreased cerebellar volume and cortical thickness. Het CKO mice had decreased ependymal cells and Mpp5 at the apical surface of cortical ventricular zone compared with wild type. Het CKO mice also failed to maintain progenitor pools essential for neurogenesis. The proportion of cortical cells undergoing apoptotic cell death increased, suggesting that cell death reduces progenitor population and neuron number. Het CKO mice also showed behavioral changes, similar to our patients. To our knowledge, this is the first report to show that variants in MPP5 are associated with GDD, behavioral abnormalities and language regression/delay. Murine modeling shows that neurogenesis is likely altered in these individuals, with cell death and skewed cellular composition playing significant roles.

Wide range of phenotypic severity in individuals with late truncations unique to the predominant CDKL5 transcript in the brain.

Cyclin-dependent kinase-like 5 (CDKL5) deficiency disorder (CDD) is caused by heterozygous or hemizygous variants in CDKL5 and is characterized by refractory epilepsy, cognitive and motor impairments, and cerebral visual impairment. CDKL5 has multiple transcripts, of which the longest transcripts, NM_003159 and NM_001037343, have been used historically in clinical laboratory testing. However, the transcript NM_001323289 is the most highly expressed in brain and contains 170 nucleotides at the 3' end of its last exon that are noncoding in other transcripts. Two truncating variants in this region have been reported in association with a CDD phenotype. To clarify the significance and range of phenotypes associated with late truncating variants in this region of the predominant transcript in the brain, we report detailed information on two individuals, updated clinical information on a third individual, and a summary of published and unpublished individuals reported in ClinVar. The two new individuals (one male and one female) each had a relatively mild clinical presentation including periods of pharmaco-responsive epilepsy, independent walking and limited purposeful communication skills. A previously reported male continued to have a severe phenotype. Overall, variants in this region demonstrate a range of clinical severity consistent with reports in CDD but with the potential for milder presentation.

Prospective, phenotype-driven selection of critically ill neonates for rapid exome sequencing is associated with high diagnostic yield.

PURPOSE: To investigate the impact of rapid-turnaround exome sequencing in critically ill neonates using phenotype-based subject selection criteria. METHODS: Intensive care unit babies aged <6 months with hypotonia, seizures, a complex metabolic phenotype, and/or multiple congenital malformations were prospectively enrolled for rapid (<7 day) trio-based exome sequencing. Genomic variants relevant to the presenting phenotype were returned to the medical team. RESULTS: A genetic diagnosis was attained in 29 of 50 (58%) sequenced cases. Twenty-seven (54%) patients received a molecular diagnosis involving known disease genes; two additional cases (4%) were solved with pathogenic variants found in novel disease genes. In 24 of the solved cases, diagnosis had impact on patient management and/or family members. Management changes included shift to palliative care, medication changes, involvement of additional specialties, and the consideration of new experimental therapies. CONCLUSION: Phenotype-based patient selection is effective at identifying critically ill neonates with a high likelihood of receiving a molecular diagnosis via rapid-turnaround exome sequencing, leading to faster and more accurate diagnoses, reducing unnecessary testing and procedures, and informing medical care.

Parental interest in genomic sequencing of newborns: enrollment experience from the BabySeq Project.

PURPOSE: Newborn genomic sequencing (nGS) has great potential to improve pediatric care. Parental interest and concerns about genomics are relatively unexplored. Understanding why parents decline research consent for nGS may reveal implementation barriers. METHODS: We evaluated parental interest in a randomized trial of nGS in well-baby and intensive care unit nursery settings. Interested families attended an informational enrollment session (ES) with a genetic counselor prior to consenting. Reason(s) for declining participation and sociodemographic associations were analyzed. RESULTS: Of 3860 eligible approached families, 10% attended ES, 67% of whom enrolled. Of 1760 families queried for decline reasons, 58% were uninterested in research. Among 499 families considering research, principal reasons for decline prior to ES included burdensome study logistics (48%), feeling overwhelmed postpartum (17%), and lack of interest/discomfort with genetic testing (17%). Decliners after ES more often cited concerns about privacy/insurability (41%) and uncertain/unfavorable results (23%). CONCLUSION: Low interest in research and study logistics were major initial barriers to postpartum enrollment and are likely generic to many postpartum research efforts. Concerns over privacy and result implications were most commonly cited in decliners after ES. Understanding parental concerns around research nGS may inform future integration of nGS into newborn screening, predictive testing, and pediatric diagnostics.

De novo variant of TRRAP in a patient with very early onset psychosis in the context of non-verbal learning disability and obsessive-compulsive disorder: a case report.

BACKGROUND: TRRAP encodes a multidomain protein kinase that works as a genetic cofactor to influence DNA methylation patterns, DNA damage repair, and chromatin remodeling. TRRAP protein is vital to early neural developmental processes, and variants in this gene have been associated with schizophrenia and childhood disintegrative disorder. CASE PRESENTATION: Here, we report on a patient with a de novo nonsynonymous TRRAP single-nucleotide variant (EST00000355540.3:c.5957G > A, p.Arg1986Gln) and early onset major depression accompanied by a psychotic episode (before age 10) that occurred in the context of longer standing nonverbal learning disability and a past history of obsessions and compulsions. CONCLUSIONS: The de novo variant and presentation of very early onset psychosis indicate a rare Mendelian disorder inheritance model. The genotype and behavioral abnormalities of this patient are reviewed.

Monogenic Hashimoto thyroiditis associated with a variant in the thyroglobulin (TG) gene.

BACKGROUND: Risk of autoimmune thyroid disease (AITD) is strongly heritable. Multiple genes confer increased risk for AITD, but a monogenic origin has not yet been described. We studied a family with apparent autosomal dominant, early onset Hashimoto thyroiditis. METHODS: The family was enrolled in an IRB-approved protocol. Whole exome sequencing was used to study the proband and an affected sibling. The identified variant was studied in other family members by Sanger sequencing. RESULTS: We identified a previously unreported splice site variant in the thyroglobulin gene (TG c.1076-1G > C). This variant was confirmed in all affected family members who underwent testing, and also noted in one unaffected child. The variant is associated with exon 9 skipping, resulting in a novel in-frame variant transcript of TG. CONCLUSION: We discovered a monogenic form of AITD associated with a splice site variant in the thyroglobulin gene. This finding raises questions about the origins of thyroid autoimmunity; possible explanations include increased immunogenicity of the mutated protein or thyroid toxicity with secondary development of anti-thyroid antibodies. Further study into the effects of this variant on thyroid function and thyroid autoimmunity are warranted.

Three-generation family with novel contiguous gene deletion on chromosome 2p22 associated with thoracic aortic aneurysm syndrome.

Latent transforming growth factor binding proteins (LTBP) are a family of extracellular matrix glycoproteins that play an important role in the regulation of transforming growth factor beta (TGF-ß) activation. Dysregulation of the TGF-ß pathway has been implicated in the pathogenesis of inherited disorders predisposing to thoracic aortic aneurysms syndromes (TAAS) including Marfan syndrome (MFS; FBN1) and Loeys-Dietz syndrome (LDS; TGFBR1, TGFBR2, TGFB2, TGFB3, SMAD2, SMAD3). While these syndromes have distinct clinical criteria, they share clinical features including aortic root dilation and musculoskeletal findings. LTBP1 is a component of the TGF-ß pathway that binds to fibrillin-1 in the extracellular matrix rendering TGF-ß inactive. We describe a three-generation family case series with a heterozygous ∼5.1 Mb novel contiguous gene deletion of chromosome 2p22.3-p22.2 involving 11 genes, including LTBP1. The deletion has been identified in the proband, father and grandfather, who all have a phenotype consistent with a TAAS. Findings include thoracic aortic dilation, ptosis, malar hypoplasia, high arched palate, retrognathia, pes planus, hindfoot deformity, obstructive sleep apnea, and low truncal tone during childhood with joint laxity that progressed to reduced joint mobility over time. While the three affected individuals did not meet criteria for either MFS or LDS, they shared features of both. Although the deletion includes 11 genes, given the relationship between LTBP1, TGF-ß, and fibrillin-1, LTBP1 stands out as one of the possible candidate genes for the clinical syndrome observed in this family. More studies are necessary to evaluate the potential role of LTBP1 in the pathophysiology of TAAS.

STRchive: a dynamic resource detailing population-level and locus-specific insights at tandem repeat disease loci.

Approximately 3% of the human genome consists of repetitive elements called tandem repeats (TRs), which include short tandem repeats (STRs) of 1-6bp motifs and variable number tandem repeats (VNTRs) of 7+bp motifs. TR variants contribute to several dozen mono- and polygenic diseases but remain understudied and "enigmatic," particularly relative to single nucleotide variants. It remains comparatively challenging to interpret the clinical significance of TR variants. Although existing resources provide portions of necessary data for interpretation at disease-associated loci, it is currently difficult or impossible to efficiently invoke the additional details critical to proper interpretation, such as motif pathogenicity, disease penetrance, and age of onset distributions. It is also often unclear how to apply population information to analyses. We present STRchive (S-T-archive, http://strchive.org/ ), a dynamic resource consolidating information on TR disease loci in humans from research literature, up-to-date clinical resources, and large-scale genomic databases, with the goal of streamlining TR variant interpretation at disease-associated loci. We apply STRchive -including pathogenic thresholds, motif classification, and clinical phenotypes-to a gnomAD cohort of ∼18.5k individuals genotyped at 60 disease-associated loci. Through detailed literature curation, we demonstrate that the majority of TR diseases affect children despite being thought of as adult diseases. Additionally, we show that pathogenic genotypes can be found within gnomAD which do not necessarily overlap with known disease prevalence, and leverage STRchive to interpret locus-specific findings therein. We apply a diagnostic blueprint empowered by STRchive to relevant clinical vignettes, highlighting possible pitfalls in TR variant interpretation. As a living resource, STRchive is maintained by experts, takes community contributions, and will evolve as understanding of TR diseases progresses.

Mono and biallelic variants in HCN2 cause severe neurodevelopmental disorders.

Hyperpolarization activated Cyclic Nucleotide (HCN) gated channels are crucial for various neurophysiological functions, including learning and sensory functions, and their dysfunction are responsible for brain disorders, such as epilepsy. To date, HCN2 variants have only been associated with mild epilepsy and recently, one monoallelic missense variant has been linked to developmental and epileptic encephalopathy. Here, we expand the phenotypic spectrum of HCN2- related disorders by describing twenty-one additional individuals from fifteen unrelated families carrying HCN2 variants. Seventeen individuals had developmental delay/intellectual disability (DD/ID), two had borderline DD/ID, and one had borderline DD. Ten individuals had epilepsy with DD/ID, with median age of onset of 10 months, and one had epilepsy with normal development. Molecular diagnosis identified thirteen different pathogenic HCN2 variants, including eleven missense variants affecting highly conserved amino acids, one frameshift variant, and one in-frame deletion. Seven variants were monoallelic of which five occurred de novo, one was not maternally inherited, one was inherited from a father with mild learning disabilities, and one was of unknown inheritance. The remaining six variants were biallelic, with four homozygous and two compound heterozygous variants. Functional studies using two-electrode voltage-clamp recordings in Xenopus laevis oocytes were performed on three monoallelic variants, p.(Arg324His), p.(Ala363Val), and p.(Met374Leu), and three biallelic variants, p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp). The p.(Arg324His) variant induced a strong increase of HCN2 conductance, while p.(Ala363Val) and p.(Met374Leu) displayed dominant negative effects, leading to a partial loss of HCN2 channel function. By confocal imaging, we found that the p.(Leu377His), p.(Pro493Leu) and p.(Gly587Asp) pathogenic variants impaired membrane trafficking, resulting in a complete loss of HCN2 elicited currents in Xenopus oocytes. Structural 3D-analysis in depolarized and hyperpolarized states of HCN2 channels, revealed that the pathogenic variants p.(His205Gln), p.(Ser409Leu), p.(Arg324Cys), p.(Asn369Ser) and p.(Gly460Asp) modify molecular interactions altering HCN2 function. Taken together, our data broadens the clinical spectrum associated with HCN2 variants, and disclose that HCN2 is involved in developmental encephalopathy with or without epilepsy.

Emerging roles and opportunities for rare disease patient advocacy groups.

Background: Patient advocacy groups (PAGs) serve a vital role for rare disease patients and families by providing educational resources, support, and a sense of community. Motivated by patient need, PAGs are increasingly at the forefront of policy, research, and drug development for their disease of interest. Objectives: The study explored the current landscape of PAGs in order to guide new and existing PAGs on available resources and challenges to research engagement. We aim to inform industry, advocates, and healthcare personnel about PAG achievements and ways they are increasingly involved in research. Design: We chose PAGs from the Rare Diseases Clinical Research Network (RDCRN) Coalition for Patient Advocacy Groups (CPAG) listserv and the National Organization for Rare Disorders (NORD) 'Find a patient organization'. Methods: We surveyed eligible PAG leaders about the demographics, goals, and research activities of their organization. For analysis, PAGs were bucketed by size, age, prevalence of disease, and budget. Data were de-identified for cross-tabulation and multinomial logistic regression analysis with R. Results: Research engagement was an extremely important goal for most PAGs (81%), though ultra-rare disease and high-budget PAGs were most likely to cite it as the top priority. In total, 79% reported research engagement in some capacity, including registries, translational research, and clinical trials. 'Ultra-rare' PAGs were less likely than 'rare' PAGs to have an ongoing clinical trial. Conclusion: While PAGs of varying sizes, budgets, and maturity levels reported an interest in research, limited funding and lack of disease awareness continue to create barriers to achieving their goals. While support tools exist to make research more accessible, often their utility depends on the funding, sustainability, maturity of the PAG itself, and the level of investment of collaborators. Despite the availability of current support systems, there are challenges related to both the start-up and sustainability of patient-centric research efforts.

Advancing Understanding of Inequities in Rare Disease Genomics.

PURPOSE: Advances in genomic research have facilitated rare disease diagnosis for thousands of individuals. Unfortunately, the benefits of advanced genetic diagnostic technology are not distributed equitably among the population, as has been seen in many other health care contexts. Quantifying and describing inequities in genetic diagnostic yield is inherently challenging due to barriers to both clinical and research genetic testing. We therefore present an implementation protocol developed to expand access to our rare disease genomic research study and to further understand existing inequities. METHODS AND FINDINGS: The Rare Genomes Project (RGP) at the Broad Institute of MIT and Harvard offers research genome sequencing to individuals with rare disease who remain genetically undiagnosed through direct interaction with the individual or family. This presents an opportunity for diagnosis beyond the clinical context, thus eliminating many barriers to access. An initial goal of RGP was to equalize access to genomic sequencing by decoupling testing access from proximity to a major medical center and physician referral. However, study participants over the initial 3 years of this project were predominantly white and well resourced. To further understand and address the lack of diversity within RGP, we developed a novel protocol embedded within the larger RGP study, in an approach informed by an implementation science framework. The aims of this protocol were: (1) to diversify recruitment and enrollment within RGP; (2) understand the process and context of implementing genomic medicine for rare disease diagnosis; and (3) investigate the value of a diagnosis for underserved populations. IMPLICATIONS: Improved understanding of existing inequities and potential strategies to address them are needed to advance equity in rare disease genetic diagnosis and research. In addition to the moral imperative of equity in genomic medicine, this approach is critical in order to fully understand the genomic underpinnings of rare disease.

Critical assessment of variant prioritization methods for rare disease diagnosis within the Rare Genomes Project.

Background: A major obstacle faced by rare disease families is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years, and causal variants are identified in under 50%. The Rare Genomes Project (RGP) is a direct-to-participant research study on the utility of genome sequencing (GS) for diagnosis and gene discovery. Families are consented for sharing of sequence and phenotype data with researchers, allowing development of a Critical Assessment of Genome Interpretation (CAGI) community challenge, placing variant prioritization models head-to-head in a real-life clinical diagnostic setting. Methods: Predictors were provided a dataset of phenotype terms and variant calls from GS of 175 RGP individuals (65 families), including 35 solved training set families, with causal variants specified, and 30 test set families (14 solved, 16 unsolved). The challenge tasked teams with identifying the causal variants in as many test set families as possible. Ranked variant predictions were submitted with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on rank position of true positive causal variants and maximum F-measure, based on precision and recall of causal variants across EPCR thresholds. Results: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performing teams recalled the causal variants in up to 13 of 14 solved families by prioritizing high quality variant calls that were rare, predicted deleterious, segregating correctly, and consistent with reported phenotype. In unsolved families, newly discovered diagnostic variants were returned to two families following confirmatory RNA sequencing, and two prioritized novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant, in an unsolved proband with phenotype overlap with asparagine synthetase deficiency. Conclusions: By objective assessment of variant predictions, we provide insights into current state-of-the-art algorithms and platforms for genome sequencing analysis for rare disease diagnosis and explore areas for future optimization. Identification of diagnostic variants in unsolved families promotes synergy between researchers with clinical and computational expertise as a means of advancing the field of clinical genome interpretation.