Closing the gap: Systematic integration of multiplexed functional data resolves variants of uncertain significance in BRCA1, TP53, and PTEN.

Clinical interpretation of missense variants is challenging because the majority identified by genetic testing are rare and their functional effects are unknown. Consequently, most variants are of uncertain significance and cannot be used for clinical diagnosis or management. Although not much can be done to ameliorate variant rarity, multiplexed assays of variant effect (MAVEs), where thousands of single-nucleotide variant effects are simultaneously measured experimentally, provide functional evidence that can help resolve variants of unknown significance (VUSs). However, a rigorous assessment of the clinical value of multiplexed functional data for variant interpretation is lacking. Thus, we systematically combined previously published BRCA1, TP53, and PTEN multiplexed functional data with phenotype and family history data for 324 VUSs identified by a single diagnostic testing laboratory. We curated 49,281 variant functional scores from MAVEs for these three genes and integrated four different TP53 multiplexed functional datasets into a single functional prediction for each variant by using machine learning. We then determined the strength of evidence provided by each multiplexed functional dataset and reevaluated 324 VUSs. Multiplexed functional data were effective in driving variant reclassification when combined with clinical data, eliminating 49% of VUSs for BRCA1, 69% for TP53, and 15% for PTEN. Thus, multiplexed functional data, which are being generated for numerous genes, are poised to have a major impact on clinical variant interpretation.

CRaTER enrichment for on-target gene editing enables generation of variant libraries in hiPSCs.

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CRISPRa On-Target Editing Retrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of variants in MYH7, a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different variants. We differentiated these hiPSCs to cardiomyocytes and show MHC-ß fusion proteins can localize as expected. Additionally, single-cell contractility analyses revealed cardiomyocytes with a pathogenic, hypertrophic cardiomyopathy-associated MYH7 variant exhibit salient HCM physiology relative to isogenic controls. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of functional transgenic cell lines at unprecedented scale.

Interpretation of Genomic Sequencing Results in Healthy and Ill Newborns: Results from the BabySeq Project.

Genomic sequencing provides many opportunities in newborn clinical care, but the challenges of interpreting and reporting newborn genomic sequencing (nGS) results need to be addressed for its broader and effective application. The BabySeq Project is a pilot randomized clinical trial that explores the medical, behavioral, and economic impacts of nGS in well newborns and those admitted to a neonatal intensive care unit (NICU). Here we present childhood-onset and actionable adult-onset disease risk, carrier status, and pharmacogenomics findings from nGS of 159 newborns in the BabySeq Project. nGS revealed a risk of childhood-onset disease in 15/159 (9.4%) newborns; none of the disease risks were anticipated based on the infants' known clinical or family histories. nGS also revealed actionable adult-onset disease risk in 3/85 (3.5%) newborns whose parents consented to receive this information. Carrier status for recessive diseases and pharmacogenomics variants were reported in 88% and 5% of newborns, respectively. Additional indication-based analyses were performed in 29/32 (91%) NICU newborns and 6/127 (5%) healthy newborns who later had presentations that prompted a diagnostic analysis. No variants that sufficiently explained the reason for the indications were identified; however, suspicious but uncertain results were reported in five newborns. Testing parental samples contributed to the interpretation and reporting of results in 13/159 (8%) newborns. Our results suggest that nGS can effectively detect risk and carrier status for a wide range of disorders that are not detectable by current newborn screening assays or predicted based on the infant's known clinical or family history, and the interpretation of results can substantially benefit from parental testing.

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.

The BabySeq project: implementing genomic sequencing in newborns.

BACKGROUND: The greatest opportunity for lifelong impact of genomic sequencing is during the newborn period. The "BabySeq Project" is a randomized trial that explores the medical, behavioral, and economic impacts of integrating genomic sequencing into the care of healthy and sick newborns. METHODS: Families of newborns are enrolled from Boston Children's Hospital and Brigham and Women's Hospital nurseries, and half are randomized to receive genomic sequencing and a report that includes monogenic disease variants, recessive carrier variants for childhood onset or actionable disorders, and pharmacogenomic variants. All families participate in a disclosure session, which includes the return of results for those in the sequencing arm. Outcomes are collected through review of medical records and surveys of parents and health care providers and include the rationale for choice of genes and variants to report; what genomic data adds to the medical management of sick and healthy babies; and the medical, behavioral, and economic impacts of integrating genomic sequencing into the care of healthy and sick newborns. DISCUSSION: The BabySeq Project will provide empirical data about the risks, benefits and costs of newborn genomic sequencing and will inform policy decisions related to universal genomic screening of newborns. TRIAL REGISTRATION: The study is registered in ClinicalTrials.gov Identifier: NCT02422511 . Registration date: 10 April 2015.

Robertsonian translocations modify genomic distribution of γH2AFX and H3.3 in mouse germ cells.

Heterozygosity for Robertsonian translocations hampers pairing and synapsis between the translocated chromosome and its normal homologs during meiotic prophase I. This causes meiotic silencing of unsynapsed chromatin in pericentromeric regions. Several lines of evidence suggest that autosomal asynapsis leads to meiotic arrest in males and two underlying mechanisms have been proposed: (1) reactivation of the X and Y chromosomes due to competition for silencing factors and (2) meiotic silencing of genes that are located in the unsynapsed regions and are essential for meiotic progression. The latter mechanism requires that asynapsis and meiotic silencing spread beyond the p-arms of the normal homologs into gene-rich regions. We used chromatin immunoprecipitation assays to determine whether histones γH2AFX and H3.3, both marks of asynapsis and meiotic silencing, are enriched in gene-rich regions of the translocated chromosomes and their homologs in the spermatocytes of heterozygous carriers of Robertsonian translocations. We also asked if γH2AFX and H3.3 enrichment was reduced at the X chromosome and if γH2AFX and H3.3 enrichment was higher on the normal homolog. Our data show that γH2AFX enrichment extends as far as 9-15 Mb of the annotated genomic sequence of the q-arms of the translocated chromosomal trivalents and that both γH2AFX and H3.3 levels are reduced over the X chromosome. Our data are also suggestive of an asymmetry in γH2AFX and H3.3 enrichment with a bias toward the non-translocated homolog.

Multiplexed Functional Assessments of MYH7 Variants in Human Cardiomyocytes.

BACKGROUND: Pathogenic autosomal-dominant missense variants in MYH7 (myosin heavy chain 7), which encodes the sarcomeric protein (ß-MHC [beta myosin heavy chain]) expressed in cardiac and skeletal myocytes, are a leading cause of hypertrophic cardiomyopathy and are clinically actionable. However, ≈75% of MYH7 missense variants are of unknown significance. While human-induced pluripotent stem cells (hiPSCs) can be differentiated into cardiomyocytes to enable the interrogation of MYH7 variant effect in a disease-relevant context, deep mutational scanning has not been executed using diploid hiPSC derivates due to low hiPSC gene-editing efficiency. Moreover, multiplexable phenotypes enabling deep mutational scanning of MYH7 variant hiPSC-derived cardiomyocytes are unknown. METHODS: To overcome these obstacles, we used CRISPRa On-Target Editing Retrieval enrichment to generate an hiPSC library containing 113 MYH7 codon variants suitable for deep mutational scanning. We first established that ß-MHC protein loss occurs in a hypertrophic cardiomyopathy human heart with a pathogenic MYH7 variant. We then differentiated the MYH7 missense variant hiPSC library to cardiomyocytes for multiplexed assessment of ß-MHC variant abundance by massively parallel sequencing and hiPSC-derived cardiomyocyte survival. RESULTS: Both the multiplexed assessment of ß-MHC abundance and hiPSC-derived cardiomyocyte survival accurately segregated all known pathogenic variants from synonymous variants. Functional data were generated for 4 variants of unknown significance and 58 additional MYH7 missense variants not yet detected in patients. CONCLUSIONS: This study leveraged hiPSC differentiation into disease-relevant cardiomyocytes to enable multiplexed assessments of MYH7 missense variants for the first time. Phenotyping strategies used here enable the application of deep mutational scanning to clinically actionable genes, which should reduce the burden of variants of unknown significance on patients and clinicians.

Defining and Reducing Variant Classification Disparities.

Background: Multiplexed Assays of Variant Effects (MAVEs) can test all possible single variants in a gene of interest. The resulting saturation-style data may help resolve variant classification disparities between populations, especially for variants of uncertain significance (VUS). Methods: We analyzed clinical significance classifications in 213,663 individuals of European-like genetic ancestry versus 206,975 individuals of non-European-like genetic ancestry from All of Us and the Genome Aggregation Database. Then, we incorporated clinically calibrated MAVE data into the Clinical Genome Resource's Variant Curation Expert Panel rules to automate VUS reclassification for BRCA1, TP53, and PTEN . Results: Using two orthogonal statistical approaches, we show a higher prevalence ( p ≤5.95e-06) of VUS in individuals of non-European-like genetic ancestry across all medical specialties assessed in all three databases. Further, in the non-European-like genetic ancestry group, higher rates of Benign or Likely Benign and variants with no clinical designation ( p ≤2.5e-05) were found across many medical specialties, whereas Pathogenic or Likely Pathogenic assignments were higher in individuals of European-like genetic ancestry ( p ≤2.5e-05). Using MAVE data, we reclassified VUS in individuals of non-European-like genetic ancestry at a significantly higher rate in comparison to reclassified VUS from European-like genetic ancestry ( p =9.1e-03) effectively compensating for the VUS disparity. Further, essential code analysis showed equitable impact of MAVE evidence codes but inequitable impact of allele frequency ( p =7.47e-06) and computational predictor ( p =6.92e-05) evidence codes for individuals of non-European-like genetic ancestry. Conclusions: Generation of saturation-style MAVE data should be a priority to reduce VUS disparities and produce equitable training data for future computational predictors.

Assigning credit where it's due: An information content score to capture the clinical value of Multiplexed Assays of Variant Effect.

Background: A variant can be pathogenic or benign with relation to a human disease. Current classification categories from benign to pathogenic reflect a probabilistic summary of current understanding. A primary metric of clinical utility for multiplexed assays of variant effect (MAVE) is the number of variants that can be reclassified from uncertain significance (VUS). However, we hypothesized that this measure of utility underrepresents the information gained from MAVEs and that an information theory approach which includes data that does not reclassify variants will better reflect true information gain. We used this information theory approach to evaluate the information gain, in bits, for MAVEs of BRCA1, PTEN, and TP53. Here, one bit represents the amount of information required to completely classify a single variant starting from no information. Results: BRCA1 MAVEs produced a total of 831.2 bits of information, 6.58% of the total missense information in BRCA1 and a 22-fold increase over the information that only contributed to VUS reclassification. PTEN MAVEs produced 2059.6 bits of information which represents 32.8% of the total missense information in PTEN and an 85-fold increase over the information that contributed to VUS reclassification. TP53 MAVEs produced 277.8 bits of information which represents 6.22% of the total missense information in TP53 and a 3.5-fold increase over the information that contributed to VUS reclassification. Conclusions: An information content approach will more accurately portray information gained through MAVE mapping efforts than counting the number of variants reclassified. This information content approach may also help define the impact of modifying information definitions used to classify many variants, such as guideline rule changes.

CRaTER enrichment for on-target gene-editing enables generation of variant libraries in hiPSCs.

Standard transgenic cell line generation requires screening 100-1000s of colonies to isolate correctly edited cells. We describe CR ISPR a On- T arget E diting R etrieval (CRaTER) which enriches for cells with on-target knock-in of a cDNA-fluorescent reporter transgene by transient activation of the targeted locus followed by flow sorting to recover edited cells. We show CRaTER recovers rare cells with heterozygous, biallelic-editing of the transcriptionally-inactive MYH7 locus in human induced pluripotent stem cells (hiPSCs), enriching on average 25-fold compared to standard antibiotic selection. We leveraged CRaTER to enrich for heterozygous knock-in of a library of single nucleotide variants (SNVs) in MYH7 , a gene in which missense mutations cause cardiomyopathies, and recovered hiPSCs with 113 different MYH7 SNVs. We differentiated these hiPSCs to cardiomyocytes and show MYH7 fusion proteins can localize as expected. Thus, CRaTER substantially reduces screening required for isolation of gene-edited cells, enabling generation of transgenic cell lines at unprecedented scale.

Reconciling newborn screening and a novel splice variant in BTD associated with partial biotinidase deficiency: a BabySeq Project case report.

Here, we report a newborn female infant from the well-baby cohort of the BabySeq Project who was identified with compound heterozygous BTD gene variants. The two identified variants included a well-established pathogenic variant (c.1612C>T, p.Arg538Cys) that causes profound biotinidase deficiency (BTD) in homozygosity. In addition, a novel splice variant (c.44+1G>A, p.?) was identified in the invariant splice donor region of intron 1, potentially predictive of loss of function. The novel variant was predicted to impact splicing of exon 1; however, given the absence of any reported pathogenic variants in exon 1 and the presence of alternative splicing with exon 1 absent in most tissues in the GTEx database, we assigned an initial classification of uncertain significance. Follow-up medical record review of state-mandated newborn screen (NBS) results revealed an initial out-of-range biotinidase activity level. Levels from a repeat NBS sample barely passed cutoff into the normal range. To determine whether the infant was biotinidase-deficient, subsequent diagnostic enzyme activity testing was performed, confirming partial BTD, and resulted in a change of management for this patient. This led to reclassification of the novel splice variant based on these results. In conclusion, combining the genetic and NBS results together prompted clinical follow-up that confirmed partial BTD and informed this novel splice site's reclassification, emphasizing the importance of combining iterative genetic and phenotypic evaluations.

Dynamics of response to asynapsis and meiotic silencing in spermatocytes from Robertsonian translocation carriers.

Failure of homologous synapsis during meiotic prophase triggers transcriptional repression. Asynapsis of the X and Y chromosomes and their consequent silencing is essential for spermatogenesis. However, asynapsis of portions of autosomes in heterozygous translocation carriers may be detrimental for meiotic progression. In fact, a wide range of phenotypic outcomes from meiotic arrest to normal spermatogenesis have been described and the causes of such a variation remain elusive. To better understand the consequences of asynapsis in male carriers of Robertsonian translocations, we focused on the dynamics of recruitment of markers of asynapsis and meiotic silencing at unsynapsed autosomal trivalents in the spermatocytes of Robertsonian translocation carrier mice. Here we report that the enrichment of breast cancer 1 (BRCA1) and histone γH2AX at unsynapsed trivalents declines during the pachytene stage of meiosis and differs from that observed in the sex body. Furthermore, histone variant H3.3S31, which associates with the sex chromosomes in metaphase I/anaphase I spermatocytes, localizes to autosomes in 12% and 31% of nuclei from carriers of one and three translocations, respectively. These data suggest that the proportion of spermatocytes with markers of meiotic silencing of unsynapsed chromatin (MSUC) at trivalents depends on both, the stage of meiosis and the number of translocations. This may explain some of the variability in phenotypic outcomes associated with Robertsonian translocations. In addition our data suggest that the dynamics of response to asynapsis in Robertsonian translocations differs from the response to sex chromosomal asynapsis in the male germ line.