New Candidates for Autism/Intellectual Disability Identified by Whole-Exome Sequencing

Intellectual disability (ID) is characterized by impairments in the cognitive processes and in the tasks of daily life. It encompasses a clinically and genetically heterogeneous group of neurodevelopmental disorders often associated with autism spectrum disorder (ASD). Social and communication abilities are strongly compromised in ASD. The prevalence of ID/ASD is 1–3%, and approximately 30% of the patients remain without a molecular diagnosis. Considering the extreme genetic locus heterogeneity, next-generation sequencing approaches have provided powerful tools for candidate gene identification. Molecular diagnosis is crucial to improve outcome, prevent complications, and hopefully start a therapeutic approach. Here, we performed parent–offspring trio whole-exome sequencing (WES) in a cohort of 60 mostly syndromic ID/ASD patients and we detected 8 pathogenic variants in genes already known to be associated with ID/ASD (SYNGAP1, SMAD6, PACS1, SHANK3, KMT2A, KCNQ2, ACTB, and POGZ). We found four de novo disruptive variants of four novel candidate ASD/ID genes: MBP, PCDHA1, PCDH15, PDPR. We additionally selected via bioinformatic tools many variants in unknown genes that alone or in combination can contribute to the phenotype. In conclusion, our data confirm the efficacy of WES in detecting pathogenic variants of known and novel ID/ASD genes.     

Variant Enrichment Analysis to Explore Pathways Functionality in Complex Autoinflammatory Skin Disorders through Whole Exome Sequencing Analysis

The challenge of unravelling the molecular basis of multifactorial disorders nowadays cannot rely just on association studies searching for potential causative variants shared by groups of patients and not present in healthy individuals; indeed, association studies have as a main limitation the lack of information on the interactions between the disease-causing variants. Thus, new genomic analysis tools focusing on disrupted pathways rather than associated gene variants are required to better understand the complexity of a disease. Therefore, we developed the Variant Enrichment Analysis (VEA) workflow, a tool applicable for whole exome sequencing data, able to find differences between the numbers of genetic variants in a given pathway in comparison with a reference dataset. In this study, we applied VEA to discover novel pathways altered in patients with complex autoinflammatory skin disorders, namely PASH (n = 9), 3 of whom are overlapping with SAPHO) and PAPASH (n = 3). With this approach we have been able to identify pathways related to neutrophil and endothelial cells homeostasis/activations, as disrupted in our patients. We hypothesized that unregulated neutrophil transendothelial migration could elicit increased neutrophil infiltration and tissue damage. Based on our findings, VEA, in our experimental dataset, allowed us to predict novel pathways impaired in subjects with autoinflammatory skin disorders.     

Deep Resequencing of 9 Candidate Genes Identifies a Role for ARAP1 and IGF2BP2 in Modulating Insulin Secretion Adjusted for Insulin Resistance in Obese Southern Europeans

Type 2 diabetes is characterized by impairment in insulin secretion, with an established genetic contribution. We aimed to evaluate common and low-frequency (1–5%) variants in nine genes strongly associated with insulin secretion by targeted sequencing in subjects selected from the extremes of insulin release measured by the disposition index. Collapsing data by gene and/or function, the association between disposition index and nonsense variants were significant, also after adjustment for confounding factors (OR = 0.25, 95% CI = 0.11–0.59, p = 0.001). Evaluating variants individually, three novel variants in ARAP1, IGF2BP2 and GCK, out of eight reaching significance singularly, remained associated after adjustment. Constructing a genetic risk model combining the effects of the three variants, only carriers of the ARAP1 and IGF2BP2 variants were significantly associated with a reduced probability to be in the lower, worst, extreme of insulin secretion (OR = 0.223, 95% CI = 0.105–0.473, p < 0.001). Observing a high number of normal glucose tolerance between carriers, a regression posthoc analysis was performed. Carriers of genetic risk model variants had higher probability to be normoglycemic, also after adjustment (OR = 2.411, 95% CI = 1.136–5.116, p = 0.022). Thus, in our southern European cohort, nonsense variants in all nine candidate genes showed association with better insulin secretion adjusted for insulin resistance, and we established the role of ARAP1 and IGF2BP2 in modulating insulin secretion.     

Identification of Rare LRP5 Variants in a Cohort of Males with Impaired Bone Mass

Osteoporosis is the most common bone disease characterized by reduced bone mass and increased bone fragility. Genetic contribution is one of the main causes of primary osteoporosis; therefore, both genders are affected by this skeletal disorder. Nonetheless, osteoporosis in men has received little attention, thus being underestimated and undertreated. The aim of this study was to identify novel genetic variants in a cohort of 128 males with idiopathic low bone mass using a next-generation sequencing (NGS) panel including genes whose mutations could result in reduced bone mineral density (BMD). Genetic analysis detected in eleven patients ten rare heterozygous variants within the LRP5 gene, which were categorized as VUS (variant of uncertain significance), likely pathogenic and benign variants according to American College of Medical Genetics and Genomics (ACMG) guidelines. Protein structural and Bayesian analysis performed on identified LRP5 variants pointed out p.R1036Q and p.R1135C as pathogenic, therefore suggesting the likely association of these two variants with the low bone mass phenotype. In conclusion, this study expands our understanding on the importance of a functional LRP5 protein in bone formation and highlights the necessity to sequence this gene in subjects with idiopathic low BMD.     

Complex Transposon Insertion as a Novel Cause of Pompe Disease

Pompe disease (OMIM#232300) is an autosomal recessive lysosomal storage disorder caused by mutations in the GAA gene. According to public mutation databases, more than 679 pathogenic variants have been described in GAA, none of which are associated with mobile genetic elements. In this article, we report a novel molecular genetic cause of Pompe disease, which could be hardly detected using routine molecular genetic analysis. Whole genome sequencing followed by comprehensive functional analysis allowed us to discover and characterize a complex mobile genetic element insertion deep in the intron 15 of the GAA gene in a patient with infantile onset Pompe disease.     

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study’s novel candidate genes and allele combinations will provide a valuable reference for future breeding work.     

Clinical Phenotype of PDE6B-Associated Retinitis Pigmentosa

In this retrospective, longitudinal, observational cohort study, we investigated the phenotypic and genotypic features of retinitis pigmentosa associated with variants in the PDE6B gene. Patients underwent clinical examination and genetic testing at a single tertiary referral center, including best-corrected visual acuity (BCVA), kinetic visual field (VF), full-field electroretinography, full-field stimulus threshold, spectral domain optical coherence tomography, and fundus autofluorescence imaging. The genetic testing comprised candidate gene sequencing, inherited retinal disease gene panel sequencing, whole-genome sequencing, and testing for familial variants by Sanger sequencing. Twenty-four patients with mutations in PDE6B from 21 families were included in the study (mean age at the first visit: 32.1 ± 13.5 years). The majority of variants were putative splicing defects (8/23) and missense (7/23) mutations. Seventy-nine percent (38/48) of eyes had no visual acuity impairment at the first visit. Visual acuity impairment was mild in 4% (2/48), moderate in 13% (6/48), and severe in 4% (2/48). BCVA was symmetrical in the right and left eyes. The kinetic VF measurements were highly symmetrical in the right and left eyes, as was the horizontal ellipsoid zone (EZ) width. Regarding the genetic findings, 43% of the PDE6B variants found in our patients were novel. Thus, this study contributed substantially to the PDE6B mutation spectrum. The visual acuity impairment was mild in 83% of eyes, providing a window of opportunity for investigational new drugs. The EZ width was reduced in all patients and was highly symmetric between the eyes, making it a promising outcome measure. We expect these findings to have implications on the design of future PDE6B-related retinitis pigmentosa (RP) clinical trials.     

Molecular Inversion Probe-Based Sequencing of USH2A Exons and Splice Sites as a Cost-Effective Screening Tool in USH2 and arRP Cases

A substantial proportion of subjects with autosomal recessive retinitis pigmentosa (arRP) or Usher syndrome type II (USH2) lacks a genetic diagnosis due to incomplete USH2A screening in the early days of genetic testing. These cases lack eligibility for optimal genetic counseling and future therapy. USH2A defects are the most frequent cause of USH2 and are also causative in individuals with arRP. Therefore, USH2A is an important target for genetic screening. The aim of this study was to assess unscreened or incompletely screened and unexplained USH2 and arRP cases for (likely) pathogenic USH2A variants. Molecular inversion probe (MIP)-based sequencing was performed for the USH2A exons and their flanking regions, as well as published deep-intronic variants. This was done to identify single nucleotide variants (SNVs) and copy number variants (CNVs) in 29 unscreened or partially pre-screened USH2 and 11 partially pre-screened arRP subjects. In 29 out of these 40 cases, two (likely) pathogenic variants were successfully identified. Four of the identified SNVs and one CNV were novel. One previously identified synonymous variant was demonstrated to affect pre-mRNA splicing. In conclusion, genetic diagnoses were obtained for a majority of cases, which confirms that MIP-based sequencing is an effective screening tool for USH2A. Seven unexplained cases were selected for future analysis with whole genome sequencing.     

Identification of Novel Genomic-Variant Patterns of OR56A5, OR52L1, and CTSD in Retinitis Pigmentosa Patients by Whole-Exome Sequencing

Inherited retinal dystrophies (IRDs) are rare but highly heterogeneous genetic disorders that affect individuals and families worldwide. However, given its wide variability, its analysis of the driver genes for over 50% of the cases remains unexplored. The present study aims to identify novel driver genes, disease-causing variants, and retinitis pigmentosa (RP)-associated pathways. Using family-based whole-exome sequencing (WES) to identify putative RP-causing rare variants, we identified a total of five potentially pathogenic variants located in genes OR56A5, OR52L1, CTSD, PRF1, KBTBD13, and ATP2B4. Of the variants present in all affected individuals, genes OR56A5, OR52L1, CTSD, KBTBD13, and ATP2B4 present as missense mutations, while PRF1 and CTSD present as frameshift variants. Sanger sequencing confirmed the presence of the novel pathogenic variant PRF1 (c.124_128del) that has not been reported previously. More causal-effect or evidence-based studies will be required to elucidate the precise roles of these SNPs in the RP pathogenesis. Taken together, our findings may allow us to explore the risk variants based on the sequencing data and upgrade the existing variant annotation database in Taiwan. It may help detect specific eye diseases such as retinitis pigmentosa in East Asia.     

Accurate Molecular Diagnosis of Gaucher Disease Using Clinical Exome Sequencing as a First-Tier Test

Gaucher disease (GD) is an autosomal recessive lysosomal disorder due to beta-glucosidase gene (GBA) mutations. The molecular diagnosis of GD is complicated by the presence of recombinant alleles originating from a highly homologous pseudogene. Clinical exome sequencing (CES) is a rapid genetic approach for identifying disease-causing mutations. However, copy number variation and recombination events are poorly detected, and further investigations are required to avoid mis-genotyping. The aim of this work was to set-up an integrated strategy for GD patients genotyping using CES as a first-line test. Eight patients diagnosed with GD were analyzed by CES. Five patients were fully genotyped, while three were revealed to be homozygous for mutations that were not confirmed in the parents. Therefore, MLPA (multiplex ligation-dependent probe amplification) and specific long-range PCR were performed, and two recombinant alleles, one of them novel, and one large deletion were identified. Furthermore, an MLPA assay performed in one family resulted in the identification of an additional novel mutation (p.M124V) in a relative, in trans with the known p.N409S mutation. In conclusion, even though CES has become extensively used in clinical practice, our study emphasizes the importance of a comprehensive molecular strategy to provide proper GBA genotyping and genetic counseling.     

Rare Titin (TTN) Variants in Diseases Associated with Sudden Cardiac Death

A leading cause of death in western countries is sudden cardiac death, and can be associated with genetic disease. Next-generation sequencing has allowed thorough analysis of genes associated with this entity, including, most recently, titin. We aimed to identify potentially pathogenic genetic variants in titin. A total of 1126 samples were analyzed using a custom sequencing panel including major genes related to sudden cardiac death. Our cohort was divided into three groups: 432 cases from patients with cardiomyopathies, 130 cases from patients with channelopathies, and 564 post-mortem samples from individuals showing anatomical healthy hearts and non-conclusive causes of death after comprehensive autopsy. None of the patients included had definite pathogenic variants in the genes analyzed by our custom cardio-panel. Retrospective analysis comparing the in-house database and available public databases also was performed. We identified 554 rare variants in titin, 282 of which were novel. Seven were previously reported as pathogenic. Of these 554 variants, 493 were missense variants, 233 of which were novel. Of all variants identified, 399 were unique and 155 were identified at least twice. No definite pathogenic variants were identified in any of genes analyzed. We identified rare, mostly novel, titin variants that seem to play a potentially pathogenic role in sudden cardiac death. Additional studies should be performed to clarify the role of these variants in sudden cardiac death.     

Succinic Semialdehyde Dehydrogenase Deficiency: In Vitro and In Silico Characterization of a Novel Pathogenic Missense Variant and Analysis of the Mutational Spectrum of ALDH5A1

Succinic semialdehyde dehydrogenase deficiency (SSADHD) is a rare, monogenic disorder affecting the degradation of the main inhibitory neurotransmitter γ-amino butyric acid (GABA). Pathogenic variants in the ALDH5A1 gene that cause an enzymatic dysfunction of succinic semialdehyde dehydrogenase (SSADH) lead to an accumulation of potentially toxic metabolites, including γ-hydroxybutyrate (GHB). Here, we present a patient with a severe phenotype of SSADHD caused by a novel genetic variant c.728T > C that leads to an exchange of leucine to proline at residue 243, located within the highly conserved nicotinamide adenine dinucleotide (NAD)⁺ binding domain of SSADH. Proline harbors a pyrrolidine within its side chain known for its conformational rigidity and disruption of protein secondary structures. We investigate the effect of this novel variant in vivo, in vitro, and in silico. We furthermore examine the mutational spectrum of all previously described disease-causing variants and computationally assess all biologically possible missense variants of ALDH5A1 to identify mutational hotspots.     

Identification of Novel Genes Associated with Partial Resistance to Aphanomyces Root Rot in Field Pea by BSR-Seq Analysis

Aphanomyces root rot, caused by Aphanomyces euteiches, causes severe yield loss in field pea (Pisum sativum). The identification of a pea germplasm resistant to this disease is an important breeding objective. Polygenetic resistance has been reported in the field pea cultivar ‘00-2067’. To facilitate marker-assisted selection (MAS), bulked segregant RNA-seq (BSR-seq) analysis was conducted using an F₈ RIL population derived from the cross of ‘Carman’ × ‘00-2067’. Root rot development was assessed under controlled conditions in replicated experiments. Resistant (R) and susceptible (S) bulks were constructed based on the root rot severity in a greenhouse study. The BSR-seq analysis of the R bulks generated 44,595,510~51,658,688 reads, of which the aligned sequences were linked to 44,757 genes in a reference genome. In total, 2356 differentially expressed genes were identified, of which 44 were used for gene annotation, including defense-related pathways (jasmonate, ethylene and salicylate) and the GO biological process. A total of 344.1 K SNPs were identified between the R and S bulks, of which 395 variants were located in 31 candidate genes. The identification of novel genes associated with partial resistance to Aphanomyces root rot in field pea by BSR-seq may facilitate efforts to improve management of this important disease.     

Lynch Syndrome Identification in Saudi Cohort of Endometrial Cancer Patients Screened by Universal Approach

Lynch syndrome (LS) is the most common cause of inherited endometrial cancer (EC). The prevalence and molecular characteristic of LS in Middle Eastern women with EC have been underexplored. To evaluate the frequency of LS in a cohort of EC patients from Saudi Arabia, a total of 436 EC cases were screened utilizing immunohistochemistry (IHC), MLH1 promoter methylation analysis and next-generation sequencing technology. A total of 53 of 436 (12.2%) ECs were classified as DNA mismatch repair-deficient (dMMR). MLH1 promoter hypermethylation was detected in 30 ECs (6.9%). Three ECs (0.7%) were found to be LS harboring germline pathogenic variants (PVs)/likely pathogenic variants (LPVs): two in the MSH2 gene and one in the MSH6 gene. Three ECs (0.7%) were Lynch-like syndrome (LLS) carrying double somatic MSH2 PVs/LPVs. Seven cases were found to have variants of uncertain significance in cancer-related genes other than MMR genes. Our results indicate that LS prevalence is low among Saudi EC patients and LLS is as common as LS in this ethnicity. Our findings could help in better understanding of the prevalence and mutational spectrum of this syndrome in Saudi Arabia, which may help in defining best strategies for LS identification, prevention and genetic counseling for EC patients.     

Cathepsin B p.Gly284Val Variant in Parkinson’s Disease Pathogenesis

Parkinson’s disease (PD) is generally considered a sporadic disorder, but a strong genetic background is often found. The aim of this study was to identify the underlying genetic cause of PD in two affected siblings and to subsequently assess the role of mutations in Cathepsin B (CTSB) in susceptibility to PD. A typical PD family was identified and whole-exome sequencing was performed in two affected siblings. Variants of interest were validated using Sanger sequencing. CTSB p.Gly284Val was genotyped in 2077 PD patients and 615 unrelated healthy controls from the Czech Republic, Ireland, Poland, Ukraine, and the USA. The gene burden analysis was conducted for the CTSB gene in an additional 769 PD probands from Mayo Clinic Florida familial PD cohort. CTSB expression and activity in patient-derived fibroblasts and controls were evaluated by qRT-PCR, western blot, immunocytochemistry, and enzymatic assay. The CTSB p.Gly284Val candidate variant was only identified in affected family members. Functional analysis of CTSB patient-derived fibroblasts under basal conditions did not reveal overt changes in endogenous expression, subcellular localization, or enzymatic activity in the heterozygous carrier of the CTSB variant. The identification of the CTSB p.Gly284Val may support the hypothesis that the CTSB locus harbors variants with differing penetrance that can determine the disease risk.     

Gene Variants Related to Cardiovascular and Pulmonary Diseases May Correlate with Severe Outcome of COVID-19

Background: Severe outcomes of COVID-19 account for up to 15% of all cases. The study aims to check if any gene variants related to cardiovascular (CVD) and pulmonary diseases (PD) are correlated with a severe outcome of COVID-19 in a Polish cohort of COVID-19 patients. Methods: In this study, a subset of 747 samples from unrelated individuals collected across Poland in 2020 and 2021 was used and whole-genome sequencing was performed. Results: The GWAS analysis of SNPs and short indels located in genes related to CVD identified one variant significant in COVID-19 severe outcome in the HADHA gene, while for the PD gene panel, we found two significant variants in the DRC1 gene. In this study, both potentially protective and risk variants were identified, of which variants in the HADHA gene deserve the most attention. Conclusions: This is the first study reporting the association between the HADHA and DRC1 genetic variants and COVID-19 severe outcome based on the cohort WGS analysis. Although all the identified variants are localised in introns, they may be correlated and therefore inherited along with other risk variants, potentially causative to severe outcome of COVID-19 but not discovered yet.