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.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease.

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

Phenotypic and mutational spectrum of ROR2-related Robinow syndrome.

Robinow syndrome is characterized by a triad of craniofacial dysmorphisms, disproportionate-limb short stature, and genital hypoplasia. A significant degree of phenotypic variability seems to correlate with different genes/loci. Disturbances of the noncanonical WNT-pathway have been identified as the main cause of the syndrome. Biallelic variants in ROR2 cause an autosomal recessive form of the syndrome with distinctive skeletal findings. Twenty-two patients with a clinical diagnosis of autosomal recessive Robinow syndrome were screened for variants in ROR2 using multiple molecular approaches. We identified 25 putatively pathogenic ROR2 variants, 16 novel, including single nucleotide variants and exonic deletions. Detailed phenotypic analyses revealed that all subjects presented with a prominent forehead, hypertelorism, short nose, abnormality of the nasal tip, brachydactyly, mesomelic limb shortening, short stature, and genital hypoplasia in male patients. A total of 19 clinical features were present in more than 75% of the subjects, thus pointing to an overall uniformity of the phenotype. Disease-causing variants in ROR2, contribute to a clinically recognizable autosomal recessive trait phenotype with multiple skeletal defects. A comprehensive quantitative clinical evaluation of this cohort delineated the phenotypic spectrum of ROR2-related Robinow syndrome. The identification of exonic deletion variant alleles further supports the contention of a loss-of-function mechanism in the etiology of the syndrome.

Novel pathogenic variants and quantitative phenotypic analyses of Robinow syndrome: WNT signaling perturbation and phenotypic variability.

Robinow syndrome (RS) is a genetically heterogeneous disorder with six genes that converge on the WNT/planar cell polarity (PCP) signaling pathway implicated (DVL1, DVL3, FZD2, NXN, ROR2, and WNT5A). RS is characterized by skeletal dysplasia and distinctive facial and physical characteristics. To further explore the genetic heterogeneity, paralog contribution, and phenotypic variability of RS, we investigated a cohort of 22 individuals clinically diagnosed with RS from 18 unrelated families. Pathogenic or likely pathogenic variants in genes associated with RS or RS phenocopies were identified in all 22 individuals, including the first variant to be reported in DVL2. We retrospectively collected medical records of 16 individuals from this cohort and extracted clinical descriptions from 52 previously published cases. We performed Human Phenotype Ontology (HPO) based quantitative phenotypic analyses to dissect allele-specific phenotypic differences. Individuals with FZD2 variants clustered into two groups with demonstrable phenotypic differences between those with missense and truncating alleles. Probands with biallelic NXN variants clustered together with the majority of probands carrying DVL1, DVL2, and DVL3 variants, demonstrating no phenotypic distinction between the NXN-autosomal recessive and dominant forms of RS. While phenotypically similar diseases on the RS differential matched through HPO analysis, clustering using phenotype similarity score placed RS-associated phenotypes in a unique cluster containing WNT5A, FZD2, and ROR2 apart from non-RS-associated paralogs. Through human phenotype analyses of this RS cohort and OMIM clinical synopses of Mendelian disease, this study begins to tease apart specific biologic roles for non-canonical WNT-pathway proteins.

RNA Sequencing in B-Cell Lymphomas.

High-throughput mRNA sequencing (RNA-Seq) provides both qualitative and quantitative evaluation of the transcriptome. This method uses complementary DNA (cDNA) to generate several millions of short sequence reads that are aligned to a reference genome allowing the comprehensive characterization of the transcripts in a cell. RNA-Seq has a wide variety of applications which lead to a pervasive adoption of this method well beyond the genomics community and a deployment of this technique as a standard part of the toolkit applied in life sciences. This chapter describes a protocol to perform mRNA sequencing using the Illumina NextSeq or MiSeq platforms, presents sequencing data quality metrics, and outlines a bioinformatic pipeline for sequence alignment, digital gene expression, identification of gene fusions, detection of transcript isoforms, description and annotation of genetic variants, and de novo immunoglobulin gene assembly.

Survival of human lymphoma cells requires B-cell receptor engagement by self-antigens.

The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) relies on chronic active B-cell receptor (BCR) signaling. BCR pathway inhibitors induce remissions in a subset of ABC DLBCL patients. BCR microclusters on the surface of ABC cells resemble those generated following antigen engagement of normal B cells. We speculated that binding of lymphoma BCRs to self-antigens initiates and maintains chronic active BCR signaling in ABC DLBCL. To assess whether antigenic engagement of the BCR is required for the ongoing survival of ABC cells, we developed isogenic ABC cells that differed solely with respect to the IgH V region of their BCRs. In competitive assays with wild-type cells, substitution of a heterologous V region impaired the survival of three ABC lines. The viability of one VH4-34(+) ABC line and the ability of its BCR to bind to its own cell surface depended on V region residues that mediate the intrinsic autoreactivity of VH4-34 to self-glycoproteins. The BCR of another ABC line reacted with self-antigens in apoptotic debris, and the survival of a third ABC line was sustained by reactivity of its BCR to an idiotypic epitope in its own V region. Hence, a diverse set of self-antigens is responsible for maintaining the malignant survival of ABC DLBCL cells. IgH V regions used by the BCRs of ABC DLBCL biopsy samples varied in their ability to sustain survival of these ABC lines, suggesting a screening procedure to identify patients who might benefit from BCR pathway inhibition.

Determination of molecular structures of HIV envelope glycoproteins using cryo-electron tomography and automated sub-tomogram averaging.

Since its discovery nearly 30 years ago, more than 60 million people have been infected with the human immunodeficiency virus (HIV) (www.usaid.gov). The virus infects and destroys CD4+ T-cells thereby crippling the immune system, and causing an acquired immunodeficiency syndrome (AIDS) (2). Infection begins when the HIV Envelope glycoprotein "spike" makes contact with the CD4 receptor on the surface of the CD4+ T-cell. This interaction induces a conformational change in the spike, which promotes interaction with a second cell surface co-receptor (5,9). The significance of these protein interactions in the HIV infection pathway makes them of profound importance in fundamental HIV research, and in the pursuit of an HIV vaccine. The need to better understand the molecular-scale interactions of HIV cell contact and neutralization motivated the development of a technique to determine the structures of the HIV spike interacting with cell surface receptor proteins and molecules that block infection. Using cryo-electron tomography and 3D image processing, we recently demonstrated the ability to determine such structures on the surface of native virus, at ˜20 Šresolution (9,14). This approach is not limited to resolving HIV Envelope structures, and can be extended to other viral membrane proteins and proteins reconstituted on a liposome. In this protocol, we describe how to obtain structures of HIV envelope glycoproteins starting from purified HIV virions and proceeding stepwise through preparing vitrified samples, collecting, cryo-electron microscopy data, reconstituting and processing 3D data volumes, averaging and classifying 3D protein subvolumes, and interpreting results to produce a protein model. The computational aspects of our approach were adapted into modules that can be accessed and executed remotely using the Biowulf GNU/Linux parallel processing cluster at the NIH (http://biowulf.nih.gov). This remote access, combined with low-cost computer hardware and high-speed network access, has made possible the involvement of researchers and students working from school or home.