Structural analysis of pathogenic missense mutations in GABRA2 and identification of a novel de novo variant in the desensitization gate

Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

Mapping the Constrained Coding Regions in the Human Genome to Their Corresponding Proteins

Constrained Coding Regions (CCRs) in the human genome have been derived from DNA sequencing data of large cohorts of healthy control populations, available in the Genome Aggregation Database (gnomAD) [1]. They identify regions depleted of protein-changing variants and thus identify segments of the genome that have been constrained during human evolution. By mapping these DNA-defined regions from genomic coordinates onto the corresponding protein positions and combining this information with protein annotations, we have explored the distribution of CCRs and compared their co-occurrence with different protein functional features, previously annotated at the amino acid level in public databases. As expected, our results reveal that functional amino acids involved in interactions with DNA/RNA, protein-protein contacts and catalytic sites are the protein features most likely to be highly constrained for variation in the control population. More surprisingly, we also found that linear motifs, linear interacting peptides (LIPs), disorder-order transitions upon binding with other protein partners and liquid-liquid phase separating (LLPS) regions are also strongly associated with high constraint for variability. We also compared intra-species constraints in the human CCRs with inter-species conservation and functional residues to explore how such CCRs may contribute to the analysis of protein variants. As has been previously observed, CCRs are only weakly correlated with conservation, suggesting that intraspecies constraints complement interspecies conservation and can provide more information to interpret variant effects.

Assessing the clinical utility of protein structural analysis in genomic variant classification: experiences from a diagnostic laboratory

BACKGROUND

The widespread clinical application of genome-wide sequencing has resulted in many new diagnoses for rare genetic conditions, but testing regularly identifies variants of uncertain significance (VUS). The remarkable rise in the amount of genomic data has been paralleled by a rise in the number of protein structures that are now publicly available, which may have clinical utility for the interpretation of missense and in-frame insertions or deletions.

METHODS

Within a UK National Health Service genomic medicine diagnostic laboratory, we investigated the number of VUS over a 5-year period that were evaluated using protein structural analysis and how often this analysis aided variant classification.

RESULTS

We found 99 novel missense and in-frame variants across 67 genes that were initially classified as VUS by our diagnostic laboratory using standard variant classification guidelines and for which further analysis of protein structure was requested. Evidence from protein structural analysis was used in the re-assessment of 64 variants, of which 47 were subsequently reclassified as pathogenic or likely pathogenic and 17 remained as VUS. We identified several case studies where protein structural analysis aided variant interpretation by predicting disease mechanisms that were consistent with the observed phenotypes, including loss-of-function through thermodynamic destabilisation or disruption of ligand binding, and gain-of-function through de-repression or escape from proteasomal degradation.

CONCLUSIONS

We have shown that using in silico protein structural analysis can aid classification of VUS and give insights into the mechanisms of pathogenicity. Based on our experience, we propose a generic evidence-based workflow for incorporating protein structural information into diagnostic practice to facilitate variant classification.

Refinements and considerations for trio whole genome sequence analysis when investigating Mendelian Diseases presenting in early childhood

To facilitate early deployment of whole-genome sequencing (WGS) for severely ill children, a standardized pipeline for WGS analysis with timely turnaround and primary care pediatric uptake is needed. We developed a bioinformatics pipeline for comprehensive gene-agnostic trio WGS analysis of children suspected of having an undiagnosed monogenic disease that included detection and interpretation of primary genetic mechanisms of disease, including SNVs/indels, CNVs/SVs, uniparental disomy (UPD), imprinted genes, short tandem repeat expansions, mobile element insertions, SMN1/2 copy number calling, and mitochondrial genome variants. We assessed primary care practitioner experience and competence in a large cohort of 521 families (comprising 90% WGS trios). Children were identified by primary practitioners for recruitment, and we used the UK index of multiple deprivation to confirm lack of patient socio-economic status ascertainment bias. Of the 521 children sequenced, 176 (34%) received molecular diagnoses, with rates as high as 45% for neurology clinics. Twenty-three of the diagnosed cases (13%) required bespoke methods beyond routine SNV/CNV analysis. In our multidisciplinary clinician user experience assessment, both pediatricians and clinical geneticists expressed strong support for rapid WGS early in the care pathway, but requested further training in determining patient selection, consenting, and variant interpretation. Rapid trio WGS provides an efficacious single-pass screening test for children when deployed by primary practitioners in clinical settings that carry high a priori risk for rare pediatric disease presentations.

De Novo VPS4A Mutations Cause Multisystem Disease with Abnormal Neurodevelopment

The endosomal sorting complexes required for transport (ESCRTs) are essential for multiple membrane modeling and membrane-independent cellular processes. Here we describe six unrelated individuals with de novo missense variants affecting the ATPase domain of VPS4A, a critical enzyme regulating ESCRT function. Probands had structural brain abnormalities, severe neurodevelopmental delay, cataracts, growth impairment, and anemia. In cultured cells, overexpression of VPS4A mutants caused enlarged endosomal vacuoles resembling those induced by expression of known dominant-negative ATPase-defective forms of VPS4A. Proband-derived fibroblasts had enlarged endosomal structures with abnormal accumulation of the ESCRT protein IST1 on the limiting membrane. VPS4A function was also required for normal endosomal morphology and IST1 localization in iPSC-derived human neurons. Mutations affected other ESCRT-dependent cellular processes, including regulation of centrosome number, primary cilium morphology, nuclear membrane morphology, chromosome segregation, mitotic spindle formation, and cell cycle progression. We thus characterize a distinct multisystem disorder caused by mutations affecting VPS4A and demonstrate that its normal function is required for multiple human developmental and cellular processes.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.

Structural analysis of pathogenic missense mutations in GABRA2 and identification of a novel de novo variant in the desensitization gate

Background

Cys‐loop receptors control neuronal excitability in the brain and their dysfunction results in numerous neurological disorders. Recently, six missense variants in GABRA2, a member of this family, have been associated with early infantile epileptic encephalopathy (EIEE). We identified a novel de novo missense variant in GABRA2 in a patient with EIEE and performed protein structural analysis of the seven variants.

Methods

The novel variant was identified by trio whole‐genome sequencing. We performed protein structural analysis of the seven variants, and compared them to previously reported pathogenic mutations at equivalent positions in other Cys‐loop receptors. Additionally, we studied the distribution of disease‐associated variants in the transmembrane helices of these proteins.

Results

The seven variants are in the transmembrane domain, either close to the desensitization gate, the activation gate, or in inter‐subunit interfaces. Six of them have pathogenic mutations at equivalent positions in other Cys‐loop receptors, emphasizing the importance of these residues. Also, pathogenic mutations are more common in the pore‐lining helix, consistent with this region being highly constrained for variation in control populations.

Conclusion

Our study reports a novel pathogenic variant in GABRA2, characterizes the regions where pathogenic mutations are in the transmembrane helices, and underscores the value of considering sequence, evolutionary, and structural information as a strategy for variant interpretation of novel missense mutations.