Prediction of human mRNA donor and acceptor sites from the DNA sequence

Occurrence of cancer in Marfan syndrome: Report of two patients with neuroblastoma and review of the literature

Marfan syndrome (MFS) is an autosomal dominant connective tissue disorder caused by pathogenic variants in FBN1, with a hitherto unknown association with cancer. Here, we present two females with MFS who developed pediatric neuroblastoma. Patient 1 presented with neonatal MFS and developed an adrenal neuroblastoma with unfavorable tumor genetics at 10 months of age. Whole genome sequencing revealed a germline de novo missense FBN1 variant (NP_000129.3:p.(Asp1322Asn)), resulting in intron 32 inclusion and exon 32 retention. Patient 2 was diagnosed with classic MFS, caused by a germline de novo frameshift variant in FBN1 (NP_000129.3:p.(Cys805Ter)). At 18 years, she developed high-risk neuroblastoma with a somatic ALK pathogenic variant (NP_004295.2:p.(Arg1275Gln)). We identified 32 reported cases of MFS with cancer in PubMed, yet none with neuroblastoma. Among patients, we observed an early cancer onset and high frequency of MFS complications. We also queried cancer databases for somatic FBN1 variants, finding 49 alterations reported in PeCan, and variants in 2% of patients in cBioPortal. In conclusion, we report the first two patients with MFS and neuroblastoma and highlight an early age at cancer diagnosis in reported patients with MFS. Further epidemiological and functional studies are needed to clarify the growing evidence linking MFS and cancer.

Colocalised Genetic Associations Reveal Alternative Splicing Variants as Candidate Causal Links for Breast Cancer Risk in 10 Loci

Genome-wide association studies (GWASs) have revealed numerous loci associated with breast cancer risk, yet the precise causal variants, their impact on molecular mechanisms, and the affected genes often remain elusive. We hypothesised that specific variants exert their influence by affecting cis-regulatory alternative splice elements. An analysis of splicing quantitative trait loci (sQTL) in healthy breast tissue from female individuals identified multiple variants linked to alterations in splicing ratios. Through colocalisation analysis, we pinpointed 43 variants within twelve genes that serve as candidate causal links between sQTL and GWAS findings. In silico splice analysis highlighted a potential mechanism for three genes-FDPS, SGCE, and MRPL11-where variants in proximity to or on the splice site modulate usage, resulting in alternative splice transcripts. Further in vitro/vivo studies are imperative to fully understand how these identified changes contribute to breast oncogenesis. Moreover, investigating their potential as biomarkers for breast cancer risk could enhance screening strategies and early detection methods for breast cancer.

Genetic factors in the pathogenesis of cardio-oncology

In recent years, with advancements in medicine, the survival period of patients with tumours has significantly increased. The adverse effects of tumour treatment on patients, especially cardiac toxicity, have become increasingly prominent. In elderly patients with breast cancer, treatment-related cardiovascular toxicity has surpassed cancer itself as the leading cause of death. Moreover, in recent years, an increasing number of novel antitumour drugs, such as multitargeted agents, antibody‒drug conjugates (ADCs), and immunotherapies, have been applied in clinical practice. The cardiotoxicity induced by these drugs has become more pronounced, leading to a complex and diverse mechanism of cardiac damage. The risks of unintended cardiovascular toxicity are increased by high-dose anthracyclines, immunotherapies, and concurrent radiation, in addition to traditional cardiovascular risk factors such as smoking, hypertension, diabetes, hyperlipidaemia, and obesity. However, these factors do not fully explain why only a subset of individuals experience treatment-related cardiac toxicity, whereas others with similar clinical features do not. Recent studies indicate that genetics play a significant role in susceptibility to the development of cardiovascular toxicity from cancer therapies. These genes are involved in drug metabolism, oxidative damage, cardiac dysfunction, and other processes. Moreover, emerging evidence suggests that epigenetics also plays a role in drug-induced cardiovascular toxicity. We conducted a review focusing on breast cancer as an example to help oncologists and cardiologists better understand the mechanisms and effects of genetic factors on cardiac toxicity. In this review, we specifically address the relationship between genetic alterations and cardiac toxicity, including chemotherapy-related genetic changes, targeted therapy-related genetic changes, and immune therapy-related genetic changes. We also discuss the role of epigenetic factors in cardiac toxicity. We hope that this review will improve the risk stratification of patients and enable therapeutic interventions that mitigate these unintended adverse consequences of life-saving cancer treatments.

Investigating Splice Defects in USH2A Using Targeted Long-Read Sequencing

Biallelic variants in USH2A are associated with retinitis pigmentosa (RP) and Type 2 Usher Syndrome (USH2), leading to impaired vision and, additionally, hearing loss in the latter. Although the introduction of next-generation sequencing into clinical diagnostics has led to a significant uplift in molecular diagnostic rates, many patients remain molecularly unsolved. It is thought that non-coding variants or variants of uncertain significance contribute significantly to this diagnostic gap. This study aims to demonstrate the clinical utility of the reverse transcription-polymerase chain reaction (RT-PCR)-Oxford Nanopore Technology (ONT) sequencing of USH2A mRNA transcripts from nasal epithelial cells to determine the splice-altering effect of candidate variants. Five affected individuals with USH2 or non-syndromic RP who had undergone whole genome sequencing were recruited for further investigation. All individuals had uncertain genotypes in USH2A, including deep intronic rare variants, c.8682-654C>G, c.9055+389G>A, and c.9959-2971C>T; a synonymous variant of uncertain significance, c.2139C>T; p.(Gly713=); and a predicted loss of function duplication spanning an intron/exon boundary, c.3812-3_3837dup p.(Met1280Ter). In silico assessment using SpliceAI provided splice-altering predictions for all candidate variants which were investigated using ONT sequencing. All predictions were found to be accurate; however, in the case of c.3812-3_3837dup, the outcome was a complex cryptic splicing pattern with predominant in-frame exon 18 skipping and a low level of exon 18 inclusion leading to the predicted stop gain. This study detected and functionally characterised simple and complex mis-splicing patterns in USH2A arising from previously unknown deep intronic variants and previously reported variants of uncertain significance, confirming the pathogenicity of the variants.

Comparison of human leukocyte antigen in patients with paroxysmal nocturnal hemoglobinuria of different clone sizes

Glycosylphosphatidylinositol-anchored protein-deficient hematopoietic stem and progenitor cell development caused by PIGA mutations cannot fully explain the pathogenesis of paroxysmal nocturnal hemoglobinuria (PNH). Herein, patients newly diagnosed with PNH at our hospital between April 2019 and April 2021 were recruited. The human leukocyte antigen (HLA) class I and II loci were analyzed, and patients were stratified by PNH clone sizes: small (< 50%) and large (≥ 50%). In 40 patients (29 males; 72.5%), the median PNH clone size was 72%. Thirteen (32.5%) and twenty-seven (67.5%) patients harbored small and large PNH clones, respectively. DRB1*15:01 and DQB1*06:02 had higher frequencies in patients with PNH than in healthy controls (adjusted P-value = 4.10 × 10-4 and 4.10 × 10-4, respectively). Whole HLA class I and II allele contributions differed (P = 0.046 and 0.065, not significant difference) when comparing patients with small and large PNH clones. B*13:01 and C*04:01 allelic frequencies were significantly higher in patients with small clones (P = 0.032 and P = 0.032, respectively). Patients with small clones had higher class II HLA evolutionary divergence (HED) (P = 0.041) and global class I and II HED (P = 0.019). In the entire cohort, 17 HLA aberrations were found in 11 (27.5%) patients. No significant differences in HLA aberrations were found between patients with small or large clones. In conclusion, patients with small clones tended to have a higher frequency of immune attack-associated alleles. A higher HED in patients with small clones may reflect a propensity for T cell-mediated autoimmunity. HLA aberrations were similar between patients with small and large clones.

Identification of an ultra-rare Alu insertion in the CFTR gene: Pitfalls and challenges in genetic test interpretation

Cystic fibrosis (CF) is a life-limiting genetic disorder characterized by defective chloride ion transport due to mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. Early detection through newborn screening programs significantly improves outcomes for individuals with CF by enabling timely intervention. Here, we report the identification of an Alu element insertion within the exon 15 of CFTR gene, initially overlooked in standard next-generation sequencing analyses. However, using traditional molecular techniques, based on polymerase chain reaction and Sanger sequencing, allowed the identification of the Alu element and the reporting of a correct diagnosis. Our analysis, based on bioinformatics tools and molecular techniques, revealed that the Alu element insertion severely affects the gene expression, splicing patterns, and structure of CFTR protein. In conclusion, this study emphasizes the importance of how the integration of human expertise and modern technologies represents a pivotal step forward in genomic medicine, ensuring the delivery of precision healthcare to individuals affected by genetic diseases.

The Influence of a Genetic Variant in CCDC78 on LMNA-Associated Skeletal Muscle Disease

Mutations in the LMNA gene-encoding A-type lamins can cause Limb-Girdle muscular dystrophy Type 1B (LGMD1B). This disease presents with weakness and wasting of the proximal skeletal muscles and has a variable age of onset and disease severity. This variability has been attributed to genetic background differences among individuals; however, such variants have not been well characterized. To identify such variants, we investigated a multigeneration family in which affected individuals are diagnosed with LGMD1B. The primary genetic cause of LGMD1B in this family is a dominant mutation that activates a cryptic splice site, leading to a five-nucleotide deletion in the mature mRNA. This results in a frame shift and a premature stop in translation. Skeletal muscle biopsies from the family members showed dystrophic features of variable severity, with the muscle fibers of some family members possessing cores, regions of sarcomeric disruption, and a paucity of mitochondria, not commonly associated with LGMD1B. Using whole genome sequencing (WGS), we identified 21 DNA sequence variants that segregate with the family members possessing more profound dystrophic features and muscle cores. These include a relatively common variant in coiled-coil domain containing protein 78 (CCDC78). This variant was given priority because another mutation in CCDC78 causes autosomal dominant centronuclear myopathy-4, which causes cores in addition to centrally positioned nuclei. Therefore, we analyzed muscle biopsies from family members and discovered that those with both the LMNA mutation and the CCDC78 variant contain muscle cores that accumulated both CCDC78 and RyR1. Muscle cores containing mislocalized CCDC78 and RyR1 were absent in the less profoundly affected family members possessing only the LMNA mutation. Taken together, our findings suggest that a relatively common variant in CCDC78 can impart profound muscle pathology in combination with a LMNA mutation and accounts for variability in skeletal muscle disease phenotypes.

DRANetSplicer: A Splice Site Prediction Model Based on Deep Residual Attention Networks

The precise identification of splice sites is essential for unraveling the structure and function of genes, constituting a pivotal step in the gene annotation process. In this study, we developed a novel deep learning model, DRANetSplicer, that integrates residual learning and attention mechanisms for enhanced accuracy in capturing the intricate features of splice sites. We constructed multiple datasets using the most recent versions of genomic data from three different organisms, Oryza sativa japonica, Arabidopsis thaliana and Homo sapiens. This approach allows us to train models with a richer set of high-quality data. DRANetSplicer outperformed benchmark methods on donor and acceptor splice site datasets, achieving an average accuracy of (96.57%, 95.82%) across the three organisms. Comparative analyses with benchmark methods, including SpliceFinder, Splice2Deep, Deep Splicer, EnsembleSplice, and DNABERT, revealed DRANetSplicer's superior predictive performance, resulting in at least a (4.2%, 11.6%) relative reduction in average error rate. We utilized the DRANetSplicer model trained on O. sativa japonica data to predict splice sites in A. thaliana, achieving accuracies for donor and acceptor sites of (94.89%, 94.25%). These results indicate that DRANetSplicer possesses excellent cross-organism predictive capabilities, with its performance in cross-organism predictions even surpassing that of benchmark methods in non-cross-organism predictions. Cross-organism validation showcased DRANetSplicer's excellence in predicting splice sites across similar organisms, supporting its applicability in gene annotation for understudied organisms. We employed multiple methods to visualize the decision-making process of the model. The visualization results indicate that DRANetSplicer can learn and interpret well-known biological features, further validating its overall performance. Our study systematically examined and confirmed the predictive ability of DRANetSplicer from various levels and perspectives, indicating that its practical application in gene annotation is justified.

Evaluating the efficacy of a long-read sequencing-based approach in the clinical diagnosis of neonatal congenital adrenocortical hyperplasia

Congenital adrenal hyperplasia (CAH) is a group of autosomal recessive disorders predominantly characterized by impaired corticosteroid synthesis. Clinical phenotypes include hypoadrenocorticism, electrolyte disturbances, abnormal gonadal development, and short stature, of which severe hyponadrenocorticism and salt wasting can be life-threatening. Genetic analysis can help in the clinical diagnosis of CAH. However, the 21-OHD-causing gene CYP21A2 is arranged in tandem with the highly homologous CYP21A1P pseudogene, making it difficult to determine the exact genotypes using the traditional method of multiplex ligation-dependent probe amplification (MLPA) plus Sanger sequencing or next-generation sequencing (NGS). We applied a long-read sequencing-based approach termed comprehensive analysis of CAH (CACAH) to 48 newborns with CAH that were diagnosed by clinical features and the traditional MLPA plus Sanger sequencing method for retrospective analysis, to evaluate its efficacy in the clinical diagnosis of neonatal CAH. Compared with the MLPA plus Sanger sequencing method, CACAH showed 100 % consistency in detecting SNV/indel variants located in exons and exon-intron boundary regions of CAH-related genes. It can directly determine the cis-trans relationship without the need to analyze parental genotypes, which reduces the time to diagnosis. Moreover, CACAH was able to distinguish different CYP21A1P/CYP21A2 and TNXA/TNXB chimeras, and detect additional variants (CYP21A2 variants c.-121C > T, c.*13G > A, c.*52C > T, c.*440C > T, c.*443 T > C, and TNXB variants c.12463 + 2 T > C, c.12204 + 5G > A). We also identified the TNXB variant c.11435_11524 + 30del alone instead of as a part of the TNXA/TNXB-CH-1 chimera in two newborns, which might be introduced by gene conversion. All of these characteristics enabled clinicians to better explain the phenotype of subjects and manage them more effectively. CACAH has a great advantage over the traditional MLPA and Sanger sequencing methods, showing substantial potential in the genetic diagnosis and screening of neonatal CAH.

Structural similarity of human papillomavirus E4 and polyomaviral VP4 exhibited by genomic analysis of the common kestrel (Falco tinnunculus) polyomavirus

Polyomaviruses are widely distributed viruses of birds that may induce developmental deformities and internal organ disorders primarily in nestlings. In this study, polyomavirus sequence was detected in kidney and liver samples of a common kestrel (Falco tinnunculus) that succumbed at a rescue station in Hungary. The amplified 5025 nucleotide (nt) long genome contained the early (large and small T antigen, LTA and STA) and late (viral proteins, VP1, VP2, VP3) open reading frames (ORFs) typical for polyomaviruses. One of the additional putative ORFs (named VP4) showed identical localization with the VP4 and ORF-X of gammapolyomaviruses, but putative splicing sites could not be found in its sequence. Interestingly, the predicted 123 amino acid (aa) long protein sequence showed the highest similarity with human papillomavirus E4 early proteins in respect of the aa distribution and motif arrangement implying similar functions. The LTA of the kestrel polyomavirus shared <59.2% nt and aa pairwise identity with the LTA sequence of other polyomaviruses and formed a separated branch in the phylogenetic tree among gammapolyomaviruses. Accordingly, the kestrel polyomavirus may be the first member of a novel species within the Gammapolyomavirus genus, tentatively named Gammapolyomavirus faltin.

Novel compound heterozygous MYO15A splicing variants in autosomal recessive non-syndromic hearing loss

BACKGROUND

Hereditary hearing loss is a highly heterogeneous disorder. This study aimed to identify the genetic cause of a Chinese family with autosomal recessive non-syndromic sensorineural hearing loss (ARNSHL).

METHODS

Clinical information and peripheral blood samples were collected from the proband and its parents. Two-step high-throughput next-generation sequencing on the Ion Torrent platform was applied to detect variants as follows. First, long-range PCR was performed to amplify all the regions of the GJB2, GJB3, SLC26A4, and MT-RNR1 genes, followed by next-generation sequencing. If no candidate pathogenetic variants were found, the targeted exon sequencing with AmpliSeq technology was employed to examine another 64 deafness-associated genes. Sanger sequencing was used to identify variants and the lineage co-segregation. The splicing of the MYO15A gene was assessed by in silico bioinformatics prediction and minigene assays.

RESULTS

Two candidate MYO15A gene (OMIM, #602,666) heterozygous splicing variants, NG_011634.2 (NM_016239.3): c.6177 + 1G > T and c.9690 + 1G > A, were identified in the proband, and these two variants were both annotated as pathogenic according to the American College of Medical Genetics and Genomics (ACMG) guidelines. Further bioinformatic analysis predicted that the c.6177 + 1G > T variant might cause exon skipping and that the c.9690 + 1G > A variant might activate a cryptic splicing donor site in the downstream intronic region. An in vitro minigene assay confirmed the above predictions.

CONCLUSIONS

We identified a compound heterozygous splicing variant in the MYO15A gene in a Han Chinese family with ARNSHL. Our results broaden the spectrum of MYO15A variants, potentially benefiting the early diagnosis, prevention, and treatment of the disease.

SignalP: The Evolution of a Web Server

SignalP ( https://services.healthtech.dtu.dk/services/SignalP-6.0/ ) is a very popular prediction method for signal peptides, the intrinsic signals that make proteins secretory. The SignalP web server has existed since 1995 and is now in its sixth major version. In this historical account, we (three authors who have taken part in the entire journey plus the first author of the latest version) describe the differences between the versions and discuss the various decisions taken along the way.

Comparison of Endogenous Alpharetroviruses (ALV-like) across Galliform Species: New Distant Proviruses

The Genus Alpharetrovirus contains viruses pathogenic mainly for chickens, forming the Avian Sarcoma and Leukosis Virus group (ASLV). Cells of most Galliform species, besides chickens, contain genetic elements (endogenous retroviruses, ERVs) that could recombine with other alpharetroviruses or express proteins, complementing defective ASLV, which may successfully replicate and cause disease. However, they are quite unknown, and only ALV-F, from ring-necked pheasants, has been partially published. Upon scrutiny of 53 genomes of different avian species, we found Alpharetrovirus-like sequences only in 12 different Galliformes, including six full-length (7.4-7.6 Kbp) and 27 partial sequences. Phylogenetic studies of the regions studied (LTR, gag, pol, and env) consistently resulted in five almost identical clades containing the same ERVs: Clade I (presently known ASLVs); Clade II (Callipepla spp. ERVs); Clade IIIa (Phasianus colchicus ERVs); Clade IIIb (Alectoris spp. ERVs); and Clade IV (Centrocercus spp. ERVs). The low pol identity scores suggested that each of these Clades may be considered a different species. ORF analysis revealed that putatively encoded proteins would be very similar in length and domains to those of other alpharetroviruses and thus potentially functional. This will undoubtedly contribute to better understanding the biology of defective viruses, especially in wild Galliformes, their evolution, and the danger they may represent for other wild species and the poultry industry.

Molecular Challenges in the Diagnosis of X-Linked Chronic Granulomatous Disease: CNVs, Intronic Variants, Skewed X-Chromosome Inactivation, and Gonosomal Mosaicism

Chronic granulomatous disease (CGD) is a prototypical inborn error of immunity affecting phagocytes, in which these cells are unable to produce reactive oxygen species. CGD is caused by defects in genes encoding subunits of the NADPH oxidase enzyme complex (CYBA, CYBB, CYBC1, NCF1, NCF2, NCF4); inflammatory responses are dysregulated, and patients are highly susceptible to recurrent severe bacterial and fungal infections. X-linked CGD (XL-CGD), caused by mutations in the CYBB gene, is the most common and severe form of CGD. In this study, we describe the analytical processes undertaken in 3 families affected with XL-CGD to illustrate several molecular challenges in the genetic diagnosis of this condition: in family 1, a girl with a heterozygous deletion of CYBB exon 13 and skewed X-chromosome inactivation (XCI); in family 2, a boy with a hemizygous deletion of CYBB exon 7, defining its consequences at the mRNA level; and in family 3, 2 boys with the same novel intronic variant in CYBB (c.1151 + 6 T > A). The variant affected the splicing process, although a small fraction of wild-type mRNA was produced. Their mother was a heterozygous carrier, while their maternal grandmother was a carrier in form of gonosomal mosaicism. In summary, using a variety of techniques, including an NGS-based targeted gene panel and deep amplicon sequencing, copy number variation calling strategies, microarray-based comparative genomic hybridization, and cDNA analysis to define splicing defects and skewed XCI, we show how to face and solve some uncommon genetic mechanisms in the diagnosis of XL-CGD.

Identification of Protein-Coding Gene Structure and Protein-Related Genes and Their Splicing Sites in Kidney Stone Disease: A Protein Big Data Analysis

The study of protein-coding gene structure and protein-related genes in kidney stone disease is used for accurate identification of splicing sites and accurate location of gene exon boundaries, which is one of the difficulties and key problems in understanding the genome and discovering new genes. Prediction techniques based on signal characteristics of conserved sequences around splicing sites, such as the weighted array model (WAM), are widely used. On this basis, several other features that can be used for splicing site recognition (such as the base composition of splicing site upstream and downstream sequences, the change of signal and base composition of upstream and downstream sequences with the C + G content of adjacent sequences) were mined further, and a model was developed to describe these features. In this study, a log-linear model that can effectively integrate these features for splicing site recognition was designed, and a SpliceKey programme was developed. The findings reveal that SpliceKey's splicing site identification accuracy is not only much better than the WAM approach, but also better than DGSplice.

Noncanonical splice-site variant in peripheral myelin protein 22 gene (PMP22) in a patient with hereditary neuropathy with liability to pressure palsies

AIM

Hereditary neuropathy with liability to pressure palsies (HNPP) is a peripheral neuropathy with autosomal dominant inheritance. Diagnosis can be made from the characteristic abnormalities determined by nerve conduction studies (NCS), including subclinical deficits at physiological compression sites. Heterozygous deletion of the chromosome 17p11.2-p12 region including the peripheral myelin protein 22 gene (PMP22) is the cause in the majority of cases. However, the loss of function of PMP22 due to frameshift-causing insertion/deletion, missense, nonsense, or splice-site disrupting variants cause HNPP in some patients. We report a case of a patient diagnosed with HNPP on the basis of clinical features and the results of NCS. No deletions of PMP22 were detected by fluorescence in situ hybridization.

METHODS

We performed direct nucleotide sequence analysis and identified a heterozygous variant, c.78 + 3G > T, in PMP22. Since this variant is located outside the canonical splice site at the exon 2-intron 2 junction, we investigated whether the variant causes aberrant splicing and leads to the skipping of exon 2 of PMP22 by in vitro minigene splicing assay.

RESULTS

We demonstrated that the c.78 + 3G > T variant causes the skipping of exon 2 and leads to loss of function of the mutant allele.

CONCLUSION

Searching for sequence variants located outside the canonical splice sites should also be considered even when deletion of PMP22 is not found in a patient with a clinical diagnosis suggesting HNPP.

Genetic Investigation of Consanguineous Pakistani Families Segregating Rare Spinocerebellar Disorders

Spinocerebellar disorders are a vast group of rare neurogenetic conditions, generally characterized by overlapping clinical symptoms including progressive cerebellar ataxia, spastic paraparesis, cognitive deficiencies, skeletal/muscular and ocular abnormalities. The objective of the present study is to identify the underlying genetic causes of the rare spinocerebellar disorders in the Pakistani population. Herein, nine consanguineous families presenting different spinocerebellar phenotypes have been investigated using whole exome sequencing. Sanger sequencing was performed for segregation analysis in all the available individuals of each family. The molecular analysis of these families identified six novel pathogenic/likely pathogenic variants; ZFYVE26: c.1093del, SACS: c.1201C>T, BICD2: c.2156A>T, ALS2: c.2171-3T>G, ALS2: c.3145T>A, and B4GALNT1: c.334_335dup, and three already reported pathogenic variants; FA2H: c.159_176del, APTX: c.689T>G, and SETX: c.5308_5311del. The clinical features of all patients in each family are concurrent with the already reported cases. Hence, the current study expands the mutation spectrum of rare spinocerebellar disorders and implies the usefulness of next-generation sequencing in combination with clinical investigation for better diagnosis of these overlapping phenotypes.

Impact of genetic structural variants in factor XI deficiency: identification, accurate characterization, and inferred mechanism by long-read sequencing

BACKGROUND

Congenital factor XI (FXI) deficiency is a probably underestimated coagulopathy that confers antithrombotic protection. Characterization of genetic defects in F11 is mainly focused on the identification of single-nucleotide variants and small insertion/deletions because they represent up to 99% of the alterations accounting for factor deficiency, with only 3 gross gene defects of structural variants (SVs) having been described.

OBJECTIVES

To identify and characterize the SVs affecting F11.

METHODS

The study was performed in 93 unrelated subjects with FXI deficiency recruited in Spanish hospitals over a period of 25 years (1997-2022). F11 was analyzed by next-generation sequencing, multiplex ligand probe amplification, and long-read sequencing.

RESULTS

Our study identified 30 different genetic variants. Interestingly, we found 3 SVs, all heterozygous: a complex duplication affecting exons 8 and 9, a tandem duplication of exon 14, and a large deletion affecting the whole gene. Nucleotide resolution obtained by long-read sequencing revealed Alu repetitive elements involved in all breakpoints. The large deletion was probably generated de novo in the paternal allele during gametogenesis, and despite affecting 30 additional genes, no syndromic features were described.

CONCLUSION

SVs may account for a high proportion of F11 genetic defects implicated in the molecular pathology of congenital FXI deficiency. These SVs, likely caused by a nonallelic homologous recombination involving repetitive elements, are heterogeneous in both type and length and may be de novo. These data support the inclusion of methods to detect SVs in this disorder, with long-read-based methods being the most appropriate because they detect all SVs and achieve adequate nucleotide resolution.

Functional characterization of human genomic variation linked to polygenic diseases

The burden of human disease lies predominantly in polygenic diseases. Since the early 2000s, genome-wide association studies (GWAS) have identified genetic variants and loci associated with complex traits. These have ranged from variants in coding sequences to mutations in regulatory regions, such as promoters and enhancers, as well as mutations affecting mediators of mRNA stability and other downstream regulators, such as 5' and 3'-untranslated regions (UTRs), long noncoding RNA (lncRNA), and miRNA. Recent research advances in genetics have utilized a combination of computational techniques, high-throughput in vitro and in vivo screening modalities, and precise genome editing to impute the function of diverse classes of genetic variants identified through GWAS. In this review, we highlight the vastness of genomic variants associated with polygenic disease risk and address recent advances in how genetic tools can be used to functionally characterize them.

Molecular Basis of Unequal Alternative Splicing of Human SCD5 and Its Alteration by Natural Genetic Variations

Alternative splicing (AS) is a major means of post-transcriptional control of gene expression, and provides a dynamic versatility of protein isoforms. Cancer-related AS disorders have diagnostic, prognostic and therapeutic values. Changes in the expression and AS of human stearoyl-CoA desaturase-5 (SCD5) are promising specific tumor markers, although the transcript variants (TVs) of the gene have not yet been confirmed. Our in silico, in vitro and in vivo study focuses on the distribution of SCD5 TVs (A and B) in human tissues, the functionality of the relevant splice sites, and their modulation by certain single-nucleotide variations (SNVs). An order of magnitude higher SCD5A expression was found compared with SCD5B. This unequal splicing is attributed to a weaker recognition of the SCD5B-specific splicing acceptor site, based on predictions confirmed by an optimized minigene assay. The pronounced dominance of SCD5A was largely modified (rs1430176385_A, rs1011850309_A) or even inverted (rs1011850309_C) by natural SNVs at the TV-specific splice sites. Our results provide long missing data on the proportion of SCD5 TVs in human tissues and reveal mutation-driven changes in SCD5 AS, potentially affecting tumor-associated reprogramming of lipid metabolism, thus having prognostic significance, which may be utilized for novel and personalized therapeutic approaches.

Phylogenetic and functional characterization of water bears (Tardigrada) tubulins

Tardigrades are microscopic ecdysozoans that can withstand extreme environmental conditions. Several tardigrade species undergo reversible morphological transformations and enter into cryptobiosis, which helps them to survive periods of unfavorable environmental conditions. However, the underlying molecular mechanisms of cryptobiosis are mostly unknown. Tubulins are evolutionarily conserved components of the microtubule cytoskeleton that are crucial in many cellular processes. We hypothesize that microtubules are necessary for the morphological changes associated with successful cryptobiosis. The molecular composition of the microtubule cytoskeleton in tardigrades is unknown. Therefore, we analyzed and characterized tardigrade tubulins and identified 79 tardigrade tubulin sequences in eight taxa. We found three α-, seven β-, one γ-, and one ε-tubulin isoform. To verify in silico identified tardigrade tubulins, we also isolated and sequenced nine out of ten predicted Hypsibius exemplaris tubulins. All tardigrade tubulins were localized as expected when overexpressed in mammalian cultured cells: to the microtubules or to the centrosomes. The presence of a functional ε-tubulin, clearly localized to centrioles, is attractive from a phylogenetic point of view. Although the phylogenetically close Nematoda lost their δ- and ε-tubulins, some groups of Arthropoda still possess them. Thus, our data support the current placement of tardigrades into the Panarthropoda clade.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

Species-specific dynamics of specialized metabolism in germinating sorghum grain revealed by temporal and tissue-resolved transcriptomics and metabolomics

Seed germination is crucial for plant productivity, and the biochemical changes during germination affect seedling survival, plant health and yield. While the general metabolism of germination is extensively studied, the role of specialized metabolism is less investigated. We therefore analyzed the metabolism of the defense compound dhurrin during sorghum (Sorghum bicolor) grain germination and early seedling development. Dhurrin is a cyanogenic glucoside, which is catabolized into different bioactive compounds at other stages of plant development, but its fate and role during germination is unknown. We dissected sorghum grain into three different tissues and investigated dhurrin biosynthesis and catabolism at the transcriptomic, metabolomic and biochemical level. We further analyzed transcriptional signature differences of cyanogenic glucoside metabolism between sorghum and barley (Hordeum vulgare), which produces similar specialized metabolites. We found that dhurrin is de novo biosynthesized and catabolized in the growing embryonic axis as well as the scutellum and aleurone layer, two tissues otherwise mainly acknowledged for their involvement in release and transport of general metabolites from the endosperm to the embryonic axis. In contrast, genes encoding cyanogenic glucoside biosynthesis in barley are exclusively expressed in the embryonic axis. Glutathione transferase enzymes (GSTs) are involved in dhurrin catabolism and the tissue-resolved analysis of GST expression identified new pathway candidate genes and conserved GSTs as potentially important in cereal germination. Our study demonstrates a highly dynamic tissue- and species-specific specialized metabolism during cereal grain germination, highlighting the importance of tissue-resolved analyses and identification of specific roles of specialized metabolites in fundamental plant processes.

Use of expanded carrier screening for retrospective diagnosis of two deceased siblings with Van Maldergem syndrome 2: case report

Van Maldergem syndrome (VMLDS) is a recessive disease which affects multiple organs including the face, ear, and limb extremities. It can be caused by pathogenic variants in either the gene DCHS1 or FAT4. Diagnosis of VMLDS is complicated, especially regarding its similarity of symptoms to Hennekam syndrome, another disorder caused by FAT4 variants. Reported patients are two infantile siblings with multiple congenital anomalies, who deceased without clinical diagnosis. Whole exome sequencing was exploited for expanded carrier screening (ECS) of their parents, which revealed a novel splicing variant in the gene FAT4, NM_024582.6: c.7018+1G>A. In silico analysis of the variant indicates loss of canonical donor splice site of intron 6. This variant is classified as pathogenic based on ACMG criteria. Reverse phenotyping of patients resulted in likely diagnosis of VMLDS2. This study reaffirms the possibility of using ECS, leading to the genetic diagnosis of a rare disease with complicated clinical features.

Phenotypic screening models for rapid diagnosis of genetic variants and discovery of personalized therapeutics

Precision medicine strives for highly individualized treatments for disease under the notion that each individual's unique genetic makeup and environmental exposures imprints upon them not only a disposition to illness, but also an optimal therapeutic approach. In the realm of rare disorders, genetic predisposition is often the predominant mechanism driving disease presentation. For such, mostly, monogenic disorders, a causal gene to phenotype association is likely. As a result, it becomes important to query the patient's genome for the presence of pathogenic variations that are likely to cause the disease. Determining whether a variant is pathogenic or not is critical to these analyses and can be challenging, as many disease-causing variants are novel and, ergo, have no available functional data to help categorize them. This problem is exacerbated by the need for rapid evaluation of pathogenicity, since many genetic diseases present in young children who will experience increased morbidity and mortality without rapid diagnosis and therapeutics. Here, we discuss the utility of animal models, with a focus mainly on C. elegans, as a contrast to tissue culture and in silico approaches, with emphasis on how these systems are used in determining pathogenicity of variants with uncertain significance and then used to screen for novel therapeutics.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

BACKGROUND

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

RESULTS

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in "cryptic" splicing regions beyond canonical donor or acceptor splice sites.

CONCLUSION

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

Generation of the First Human In Vitro Model for McArdle Disease Based on iPSC Technology

McArdle disease is a rare autosomal recessive disorder caused by mutations in the PYGM gene. This gene encodes for the skeletal muscle isoform of glycogen phosphorylase (myophosphorylase), the first enzyme in glycogenolysis. Patients with this disorder are unable to obtain energy from their glycogen stored in skeletal muscle, prompting an exercise intolerance. Currently, there is no treatment for this disease, and the lack of suitable in vitro human models has prevented the search for therapies against it. In this article, we have established the first human iPSC-based model for McArdle disease. For the generation of this model, induced pluripotent stem cells (iPSCs) from a patient with McArdle disease (harbouring the homozygous mutation c.148C>T; p.R50* in the PYGM gene) were differentiated into myogenic cells able to contract spontaneously in the presence of motor neurons and generate calcium transients, a proof of their maturity and functionality. Additionally, an isogenic skeletal muscle model of McArdle disease was created. As a proof-of-concept, we have tested in this model the rescue of PYGM expression by two different read-through compounds (PTC124 and RTC13). The developed model will be very useful as a platform for testing drugs or compounds with potential pharmacological activity.

Palmoplantar Keratoderma: A Molecular Genetic Analysis of Family Cases

Palmoplantar keratoderma is a clinically polymorphic disorder with a heterogeneous etiology characterized by marked hyperkeratotic lesions on the surface of palms and soles. Hereditary forms of palmoplantar keratoderma usually have autosomal dominant inheritance and are caused by mutations in dozens of genes, most of which belong to the keratin family. We carried out clinical and molecular genetic analysis of the affected and healthy members of four families with autosomal dominant palmoplantar keratoderma. In three out of four family cases of autosomal dominant palmoplantar keratoderma, the following molecular genetic causes were established: in two families—previously non-described missense mutations in the AQP5 gene (NM_001651.4): c.369C>G (p.(Asn123Lys)) and c.103T>G (p.(Trp35Gly)); in one family—a described splice site mutation in the KRT9 gene (NM_000226.4): c.31T>G. In one family, the possible cause of palmoplantar keratoderma was detected—a variant in the KRT1 gene (NM_006121.4): c.931G>A (p.(Glu311Lys)).

Novel phenotype and genotype spectrum of NARS2 and literature review of previous mutations

BACKGROUND

Mutations in NARS2 (MIM: 612803) are associated with combined oxidative phosphorylation deficiency 24 (COXPD24; MIM: 616239) that is a rare mitochondrial and a multisystem autosomal recessive disorder.

AIMS

We aimed to detect the underlying genetic factors in two siblings with progressive ataxia, epilepsy, and severe-to-profound hearing impairment.

METHODS

After doing medical assessments and pertinent tests (i.e., auditory brainstem responses, pure tone otoacoustic emission test, cardiac examinations, computed tomography, and electroencephalogram), because of the clinical and probable genetic heterogeneity, whole-exome sequencing was performed, and co-segregation analysis was confirmed by Sanger sequencing. Biological impacts of the novel variant were evaluated using sequence-to-function bioinformatics tools.

RESULTS

A novel homozygous missense variant, NM_024678.6:c.545 T > A; p.(Ile182Lys), in exon 5 of NARS2 was identified in both patients and verified by Sanger sequencing. In silico analyses introduced this variant as pathogenic. Mitral valve prolapses with mild regurgitation, brachymetatarsia, severe hallux valgus, and clubbed fingers were reported as novel manifestations in association with NARS2 gene. By doing a literature review, we also underscored the high heterogeneity of disease phenotype.

CONCLUSIONS

Herein, we report some novel phenotype and genotype features of two female patients in an Iranian consanguineous family with COXPD24, caused by a variant in NARS2-NM_024678.6: c.545 T > A; p.(Ile182Lys). Moreover, our data expanded the phenotype and genotype spectrum of NARS2-related disorder and confirmed an unpredictable nature of genotype-phenotype correlation in COXPD24.

Characterisation of a novel OPA1 splice variant resulting in cryptic splice site activation and mitochondrial dysfunction

Autosomal dominant optic atrophy (DOA) is an inherited optic neuropathy that results in progressive, bilateral visual acuity loss and field defects. OPA1 is the causative gene in around 60% of cases of DOA. The majority of patients have a pure ocular phenotype, but 20% have extra-ocular features (DOA +). We report on a patient with DOA + manifesting as bilateral optic atrophy, spastic paraparesis, urinary incontinence and white matter changes in the central nervous system associated with a novel heterozygous splice variant NM_015560.2(OPA1):c.2356-1 G > T. Further characterisation, which was performed using fibroblasts obtained from a skin biopsy, demonstrated that this variant altered mRNA splicing of the OPA1 transcript, specifically a 21 base pair deletion at the start of exon 24, NM_015560.2(OPA1):p.Cys786_Lys792del. The majority of variant transcripts were shown to escape nonsense-mediated decay and modelling of the predicted protein structure suggests that the in-frame 7 amino acid deletion may affect OPA1 oligomerisation. Fibroblasts carrying the c.2356-1 G > T variant demonstrated impaired mitochondrial bioenergetics, membrane potential, increased cell death, and disrupted and fragmented mitochondrial networks in comparison to WT cells. This study suggests that the c.2356-1 G > T OPA1 splice site variant leads to a cryptic splice site activation and may manifest in a dominant-negative manner, which could account for the patient's severe syndromic phenotype.

Analysis of TMIE gene mutations including the first large deletion of exon 1 with autosomal recessive non-syndromic deafness

Background

Transmembrane inner ear (TMIE) protein is an essential component of the mechanotransduction complex. In collaboration with other components, TMIE aids the maintenance and function of the sensory hair cells. Autosomal recessive deafness-6 (DFNB6) is caused by mutated TMIE, a gene in the high genetic heterogeneity spectrum of deafness. Hearing loss has a significant impact on the global economy and the quality of life of affected persons, their families, and society. Here, three unrelated families with TMIE variants are presented. All three cases were found while studying the genetic causes of an Iranian cohort of subjects with cochlear implants.

Methods

Whole exome sequencing was performed to find possible genetic etiology in probands of families after a comprehensive medical evaluation for hearing loss. Co-segregation analysis in probands and other family members was performed by Sanger sequencing. The variants were interpreted per the American College of Medical Genetics and Genomics guidelines.

Results

Three different variants associated with TMIE were confirmed as reasons for autosomal recessive non-syndromic deafness. The first novel ~ 10-kb deletion surrounding exon 1 of TMIE along with two previously reported variants co-segregated with families including a frameshift variant c.122_125dup (p.Pro43fs) and a missense variant c.250 C > T; p.(Arg84Trp) in exons 2, and 3, respectively.

Conclusion

This study increases the mutational spectrum of the TMIE gene and highlights the importance of the large deletion of this gene as a reason for hearing loss. Moreover, an efficient and simple multiplex PCR assay was developed to determine the exact breakpoints of the TMIE deletion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01287-9.

Novel Loss-of-Function Variants in CHD2 Cause Childhood-Onset Epileptic Encephalopathy in Chinese Patients

Developmental and epileptic encephalopathy-94 (DEE94) is a severe form of epilepsy characterized by a broad spectrum of neurodevelopmental disorders. It is caused by pathogenic CHD2 variants. While only a few pathogenic CHD2 variants have been reported with detailed clinical phenotypes, most of which lack molecular analysis. In this study, next-generation sequencing (NGS) was performed to identify likely pathogenic CHD2 variants in patients with epilepsy. Three likely pathogenic variants were finally identified in different patients. The seizure onset ages were from two years to six years. Patients 1 and 2 had developmental delays before epilepsy, while patient 3 had intellectual regression after the first seizure onset. The observed seizures were myoclonic, febrile, and generalized tonic-clonic, which had been controlled by different combinations of antiepileptic drugs. Two de novo (c.1809_1809+1delGGinsTT, p.? and c.3455+2_3455+3insTG, p.?) and one maternal (c.3783G>A, p.W1261*) variant were identified, which were all predicted to be pathogenic/likely pathogenic. Molecular analysis was performed in patient 1, and we detected aberrantly spliced products, proving the pathogenicity of this CHD2 variant. New cases with novel variants, along with a detailed clinical and molecular analysis, are important for a better understanding of CHD2-related epileptic encephalopathy.

A Severe Dementia Syndrome Caused by Intron Retention and Cryptic Splice Site Activation in STUB1 and Exacerbated by TBP Repeat Expansions

Heterozygous pathogenic variants in the STIP1 homologous and U-box containing protein 1 (STUB1) gene have been identified as causes of autosomal dominant inherited spinocerebellar ataxia type 48 (SCA48). SCA48 is characterized by an ataxic movement disorder that is often, but not always, accompanied by a cognitive affective syndrome. We report a severe early onset dementia syndrome that mimics frontotemporal dementia and is caused by the intronic splice donor variant c.524+1G>A in STUB1. Impaired splicing was demonstrated by RNA analysis and in minigene assays of mutated and wild-type constructs of STUB1. The most striking consequence of this splicing impairment was retention of intron 3 in STUB1, which led to an in-frame insertion of 63 amino acids (aa) (p.Arg175_Glu176ins63) into the highly conserved coiled-coil domain of its encoded protein, C-terminus of HSP70-interacting protein (CHIP). To a lesser extent, activation of two cryptic splice sites in intron 3 was observed. The almost exclusively used one, c.524+86, was not predicted by in silico programs. Variant c.524+86 caused a frameshift (p.Arg175fs*93) that resulted in a truncated protein and presumably impairs the C-terminal U-box of CHIP, which normally functions as an E3 ubiquitin ligase. The cryptic splice site c.524+99 was rarely used and led to an in-frame insertion of 33 aa (p.Arg175_Glu176ins33) that resulted in disruption of the coiled-coil domain, as has been previously postulated for complete intron 3 retention. We additionally detected repeat expansions in the range of reduced penetrance in the TATA box-binding protein (TBP) gene by excluding other genes associated with dementia syndromes. The repeat expansion was heterozygous in one patient but compound heterozygous in the more severely affected patient. Therefore, we concluded that the observed severe dementia syndrome has a digenic background, making STUB1 and TBP important candidate genes responsible for early onset dementia syndromes.

Germline variants in tumor suppressor FBXW7 lead to impaired ubiquitination and a neurodevelopmental syndrome

Neurodevelopmental disorders are highly heterogenous conditions resulting from abnormalities of brain architecture and/or function. FBXW7 (F-box and WD-repeat-domain-containing 7), a recognized developmental regulator and tumor suppressor, has been shown to regulate cell-cycle progression and cell growth and survival by targeting substrates including CYCLIN E1/2 and NOTCH for degradation via the ubiquitin proteasome system. We used a genotype-first approach and global data-sharing platforms to identify 35 individuals harboring de novo and inherited FBXW7 germline monoallelic chromosomal deletions and nonsense, frameshift, splice-site, and missense variants associated with a neurodevelopmental syndrome. The FBXW7 neurodevelopmental syndrome is distinguished by global developmental delay, borderline to severe intellectual disability, hypotonia, and gastrointestinal issues. Brain imaging detailed variable underlying structural abnormalities affecting the cerebellum, corpus collosum, and white matter. A crystal-structure model of FBXW7 predicted that missense variants were clustered at the substrate-binding surface of the WD40 domain and that these might reduce FBXW7 substrate binding affinity. Expression of recombinant FBXW7 missense variants in cultured cells demonstrated impaired CYCLIN E1 and CYCLIN E2 turnover. Pan-neuronal knockdown of the Drosophila ortholog, archipelago, impaired learning and neuronal function. Collectively, the data presented herein provide compelling evidence of an F-Box protein-related, phenotypically variable neurodevelopmental disorder associated with monoallelic variants in FBXW7.

Whole-exome sequencing in a Japanese multiplex family identifies new susceptibility genes for intracranial aneurysms

Background

Intracranial aneurysms (IAs) cause subarachnoid hemorrhage, which has high rates of mortality and morbidity when ruptured. Recently, the role of rare variants in the genetic background of complex diseases has been increasingly recognized. The aim of this study was to identify rare variants for susceptibility to IA.

Methods

Whole-exome sequencing was performed on seven members of a Japanese pedigree with highly aggregated IA. Candidate genes harboring co-segregating rare variants with IA were re-sequenced and tested for association with IA using additional 500 probands and 323 non-IA controls. Functional analysis of rare variants detected in the pedigree was also conducted.

Results

We identified two gene variants shared among all four affected participants in the pedigree. One was the splicing donor c.1515+1G>A variant in NPNT (Nephronectin), which was confirmed to cause aberrant splicing by a minigene assay. The other was the missense p.P83T variant in CBY2 (Chibby family member 2). Overexpression of p.P83T CBY2 fused with red fluorescent protein tended to aggregate in the cytoplasm. Although Nephronectin has been previously reported to be involved in endothelial angiogenic functions, CBY2 is a novel molecule in terms of vascular pathophysiology. We confirmed that CBY2 was expressed in cerebrovascular smooth muscle cells in an isoform2-specific manner. Targeted CBY2 re-sequencing in additional case-control samples identified three deleterious rare variants (p.R46H, p.P83T, and p.L183R) in seven probands, showing a significant enrichment in the overall probands (8/501) compared to the controls (0/323) (p = 0.026, Fisher’s extract test).

Conclusions

NPNT and CBY2 were identified as novel susceptibility genes for IA. The highly heterogeneous and polygenic architecture of IA susceptibility can be uncovered by accumulating extensive analyses that focus on each pedigree with a high incidence of IA.

In silico profiling of non-synonymous SNPs in IDS gene for early diagnosis of Hunter syndrome

Background

Single amino acid substitutions in the Iduronate-2-sulfatase enzyme result in destabilization of the protein and cause a genetic disorder called Hunter syndrome. To gain functional insight into the mutations causing Hunter syndrome, various bioinformatics tools were employed, and special significance is given to molecular docking.

Results

In-silico tools available online for preliminary analysis including SIFT, PolyPhen 2.0, etc., were primarily employed and have identified 51 Non-synonymous Single Nucleotide Polymorphisms (ns-SNPs) as possibly deleterious. Further, modelling and energy minimization followed by Root Mean Square Deviation (RMSD) calculation has labelled 42 mutations as probably deleterious ns-SNPs. Later, trajectory analysis was performed using online tools like PSIPRED, SRide, etc., and has predicted six ns-SNPs as potentially deleterious. Additionally, docking was performed, and three candidate ns-SNPs were identified. Finally, these three ns-SNPs were confirmed to play a significant role in causing syndrome through root mean square fluctuation (RMSF) calculations.

Conclusion

From the observed results, G134E, V503D, and E521D were predicted to be candidate ns-SNPs in comparison with other in-silico tools and confirmed by RMSF calculations. Thus, the identified candidate ns-SNPs can be employed as a potential genetic marker in the early diagnosis of Hunter syndrome after clinical validation.

CYP2D6 gene resequencing in the Malagasy, a population at the crossroads between Asia and Africa: a pilot study

Background:

Plasmodium vivax malaria is endemic in Madagascar, where populations have genetic inheritance from Southeast Asia and East Africa. Primaquine, a drug of choice for vivax malaria, is metabolized principally via CYP2D6. CYP2D6 variation was characterized by locus-specific gene sequencing and was compared with TaqMan™ genotype data.

Materials & methods:

Long-range PCR amplicons were generated from 96 Malagasy samples and subjected to next-generation sequencing.

Results:

The authors observed high concordance between TaqMan™-based CYP2D6 genotype calls and the base calls from sequencing. In addition, there are new variants and haplotypes present in the Malagasy.

Conclusion:

Sequencing unique admixed populations provides more detailed and accurate insights regarding CYP2D6 variability, which may help optimize primaquine treatment across human genetic diversity.

A functional cellular framework for sex and estrous cycle-dependent gene expression and behavior

Sex hormones exert a profound influence on gendered behaviors. How individual sex hormone-responsive neuronal populations regulate diverse sex-typical behaviors is unclear. We performed orthogonal, genetically targeted sequencing of four estrogen receptor 1-expressing (Esr1+) populations and identified 1,415 genes expressed differentially between sexes or estrous states. Unique subsets of these genes were distributed across all 137 transcriptomically defined Esr1+ cell types, including estrous stage-specific ones, that comprise the four populations. We used differentially expressed genes labeling single Esr1+ cell types as entry points to functionally characterize two such cell types, BNSTprTac1/Esr1 and VMHvlCckar/Esr1. We observed that these two cell types, but not the other Esr1+ cell types in these populations, are essential for sex recognition in males and mating in females, respectively. Furthermore, VMHvlCckar/Esr1 cell type projections are distinct from those of other VMHvlEsr1 cell types. Together, projection and functional specialization of dimorphic cell types enables sex hormone-responsive populations to regulate diverse social behaviors.

Dyserythropoietic anaemia with an intronic GATA1 splicing mutation in patients suspected to have Diamond-Blackfan anaemia

Diamond-Blackfan anaemia (DBA) shares clinical features with two recently reported sporadic cases of dyserythropoietic anaemia with a cryptic GATA1 splicing mutation (c.871-24 C>T). We hypothesized that some patients clinically diagnosed with DBA but whose causative genes were unknown may carry the intronic GATA1 mutation. Here, we examined 79 patients in our DBA cohort, who had no detectable causative genes. The intronic GATA1 mutation was identified in two male patients sharing the same pedigree that included multiple cases with anaemia. Cosegregation of this mutation and disease in multiple family members provide evidence to support the pathogenicity of the intronic GATA1 mutation.

Tra2beta-Dependent Regulation of RIO Kinase 3 Splicing During Rift Valley Fever Virus Infection Underscores the Links Between Alternative Splicing and Innate Antiviral Immunity

Rift Valley fever virus (RVFV) is an emerging pathogen that has potential to cause severe disease in humans and domestic livestock. Propagation of RVFV strain MP-12 is negatively impacted by the actions of RIOK3, a protein involved in the cellular immune response to viral infection. During RVFV infection, RIOK3 mRNA is alternatively spliced to produce an isoform that correlates with the inhibition of interferon β signaling. Here, we identify splicing factor TRA2-β (also known as TRA2beta and hTRA2-β) as a key regulator governing the relative abundance of RIOK3 splicing isoforms. Using RT-PCR and minigenes, we determined that TRA2-β interaction with RIOK3 pre-mRNA was necessary for constitutive splicing of RIOK3 mRNA, and conversely, lack of TRA2-β engagement led to increased alternative splicing. Expression of TRA2-β was found to be necessary for RIOK3's antiviral effect against RVFV. Intriguingly, TRA2-β mRNA is also alternatively spliced during RVFV infection, leading to a decrease in cellular TRA2-β protein levels. These results suggest that splicing modulation serves as an immune evasion strategy by RVFV and/or is a cellular mechanism to prevent excessive immune response. Furthermore, the results suggest that TRA2-β can act as a key regulator of additional steps of the innate immune response to viral infection.

Crimean-Congo Hemorrhagic Fever Virus Past Infections Are Associated with Two Innate Immune Response Candidate Genes in Dromedaries

Dromedaries are an important livestock, used as beasts of burden and for meat and milk production. However, they can act as an intermediate source or vector for transmitting zoonotic viruses to humans, such as the Middle East respiratory syndrome coronavirus (MERS-CoV) or Crimean-Congo hemorrhagic fever virus (CCHFV). After several outbreaks of CCHFV in the Arabian Peninsula, recent studies have demonstrated that CCHFV is endemic in dromedaries and camel ticks in the United Arab Emirates (UAE). There is no apparent disease in dromedaries after the bite of infected ticks; in contrast, fever, myalgia, lymphadenopathy, and petechial hemorrhaging are common symptoms in humans, with a case fatality ratio of up to 40%. We used the in-solution hybridization capture of 100 annotated immune genes to genotype 121 dromedaries from the UAE tested for seropositivity to CCHFV. Through univariate linear regression analysis, we identified two candidate genes belonging to the innate immune system: FCAR and CLEC2B. These genes have important functions in the host defense against viral infections and in stimulating natural killer cells, respectively. This study opens doors for future research into immune defense mechanisms in an enzootic host against an important zoonotic disease.

Six splice site variations, three of them novel, in the ABO gene occurring in nine individuals with ABO subtypes

Background

Nucleotide mutations in the ABO gene may reduce the activity of glycosyltransferase, resulting in lower levels of A or B antigen expression in red blood cells. Six known splice sites have been identified according to the database of red cell immunogenetics and the blood group terminology of the International Society of Blood Transfusion. Here, we describe six distinct splice site variants in individuals with ABO subtypes.

Methods

The ABO phenotype was examined using a conventional serological method. A polymerase chain reaction sequence-based typing method was used to examine the whole coding sequence of the ABO gene. The ABO gene haplotypes were studied using allele-specific primer amplification or cloning technology. In silico analytic tools were used to assess the functional effect of splice site variations.

Results

Six distinct variants in the ABO gene splice sites were identified in nine individuals with ABO subtypes, including c.28 + 1_2delGT, c.28 + 5G > A, c.28 + 5G > C, c.155 + 5G > A, c.204-1G > A and c.374 + 5G > A. c.28 + 1_2delGT was detected in an A_w individual, while c.28 + 5G > A, c.28 + 5G > C, and c.204-1G > A were detected in B_el individuals. c.155 + 5G > A was detected in one B₃ and two AB₃ individuals, whereas c.374 + 5G > A was identified in two A_el individuals. Three novel splice site variants (c.28 + 1_2delGT, c.28 + 5G > A and c.28 + 5G > C) in the ABO gene were discovered, all of which resulted in low antigen expression. In silico analysis revealed that all variants had the potential to alter splice transcripts.

Conclusions

Three novel splice site variations in the ABO gene were identified in Chinese individuals, resulting in decreased A or B antigen expression and the formation of ABO subtypes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-03141-5.

Spinal cord-predominant neuropathology in an adult-onset case of POLR3A-related spastic ataxia

Biallelic mutations in POLR3A have been associated with childhood‐onset hypomyelinating leukodystrophies and adolescent‐to‐adult‐onset spastic ataxia, the latter of which has been linked to the intronic variant c.1909 + 22G>A. We report a case of adult‐onset spastic ataxia in a 75‐year‐old man, being a compound heterozygous carrier of this variant, whose brain and spinal cord were for the first time investigated by neuropathological examination. We describe prominent degeneration of the posterior columns, spinocerebellar tracts, and anterior corticospinal tracts of the spinal cord in a pattern resembling Friedreich's ataxia, with a notable lack of significant white matter pathology throughout the brain, in marked contrast with childhood‐onset cases. Immunohistochemical examination for the POLR3A protein demonstrated no apparent differences in localization or staining intensity between the proband and an age‐matched control subject. We demonstrate the clinicopathologic description of POLR3A‐related neurodegenerative disease and also mention the differential diagnosis of the childhood‐onset hypomyelinating leukodystrophy and late‐onset spastic ataxia phenotypes.

Diagnostic Yield of Targeted Hearing Loss Gene Panel Sequencing in a Large German Cohort With a Balanced Age Distribution from a Single Diagnostic Center: An Eight-year Study

Objectives:

Hereditary hearing loss exhibits high degrees of genetic and clinical heterogeneity. To elucidate the population-specific and age-related genetic and clinical spectra of hereditary hearing loss, we investigated the sequencing data of causally associated hearing loss genes in a large cohort of hearing-impaired probands with a balanced age distribution from a single center in Southwest Germany.

Design:

Genetic testing was applied to 305 hearing-impaired probands/families with a suspected genetic hearing loss etiology and a balanced age distribution over a period of 8 years (2011–2018). These individuals were representative of the regional population according to age and sex distributions. The genetic testing workflow consisted of single-gene screening (n = 21) and custom-designed hearing loss gene panel sequencing (n = 284) targeting known nonsyndromic and syndromic hearing loss genes in a diagnostic setup. Retrospective reanalysis of sequencing data was conducted by applying the current American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines.

Results:

A genetic diagnosis was established for 75 (25%) of the probands that involved 75 causal variants in 35 genes, including 16 novel causal variants and 9 medically significant variant reclassifications. Nearly half of the solved cases (47%; n = 35) were related to variants in the five most frequently affected genes: GJB2 (25%), MYO15A, WFS1, SLC26A4, and COL11A1 (all 5%). Nearly one-quarter of the cases (23%; n = 17) were associated with variants in seven additional genes (TMPRSS3, COL4A3, LOXHD1, EDNRB, MYO6, TECTA, and USH2A). The remaining one-third of single cases (33%; n = 25) were linked to variants in 25 distinct genes. Diagnostic rates and gene distribution were highly dependent on phenotypic characteristics. A positive family history of autosomal-recessive inheritance in combination with early onset and higher grades of hearing loss significantly increased the solve rate up to 60%, while late onset and lower grades of hearing loss yielded significantly fewer diagnoses. Regarding genetic diagnoses, autosomal-dominant genes accounted for 37%, autosomal-recessive genes for 60%, and X-linked genes for 3% of the solved cases. Syndromic/nonsyndromic hearing loss mimic genes were affected in 27% of the genetic diagnoses.

Conclusions:

The genetic epidemiology of the largest German cohort subjected to comprehensive targeted sequencing for hereditary hearing loss to date revealed broad causal gene and variant spectra in this population. Targeted hearing loss gene panel analysis proved to be an effective tool for ensuring an appropriate diagnostic yield in a routine clinical setting including the identification of novel variants and medically significant reclassifications. Solve rates were highly sensitive to phenotypic characteristics. The unique population-adapted and balanced age distribution of the cohort favoring late hearing loss onset uncovered a markedly large contribution of autosomal-dominant genes to the diagnoses which may be a representative for other age balanced cohorts in other populations.

Genetic Variants in the Protein S ( PROS1 ) Gene and Protein S Deficiency in a Danish Population

Protein S (PS) deficiency is a risk factor for venous thromboembolism (VTE) and can be caused by variants of the gene encoding PS ( PROS1 ). This study aimed to evaluate the clinical value of molecular analysis of the PROS1 gene in PS-deficient participants. We performed Sanger sequencing of the coding region of the PROS1 gene and multiplex ligation-dependent probe amplification to exclude large structural rearrangements. Free PS was measured by a particle-enhanced immunoassay, while PS activity was assessed by a clotting method.

A total of 87 PS-deficient participants and family members were included. In 22 index participants, we identified 13 PROS1 coding variants. Five variants were novel. In 21 index participants, no coding sequence variants or structural rearrangements were identified. The free PS level was lower in index participants carrying a PROS1 variant compared with index participants with no variant (0.51 [0.32–0.61] vs. 0.62 [0.57–0.73] × 10 ³ IU/L; p  < 0.05). The p.(Thr78Met) variant was associated with only slightly decreased free PS levels (0.59 [0.53–0.66] × 10 ³ IU/L) compared with the p.(Glu390Lys) variant (0.27 [0.24–0.37] × 10 ³ IU/L, p  < 0.01). The frequency of VTE in participants with a coding PROS1 variant was 43 and 17% in the group with normal PROS1 gene ( p  = 0.05).

In conclusion, we report 13 PROS1 coding variants including five novel variants. PS levels differ by PROS1 variant and the frequency of VTE was higher when a coding PROS1 variant was present. Hence, molecular analysis of the PROS1 gene may add clinical value in the diagnostic work-up of PS deficiency.

Effects of SLCO1B1 and SLCO1B3 Genetic Polymorphisms on Valsartan Pharmacokinetics in Healthy Korean Volunteers

Purpose: This study aimed to examine OATP1B1 (SLCO1B1) and OATP1B3 (SLCO1B3) on the pharmacokinetics of valsartan. Twenty-five subjects were genotyped for 16 single-nucleotide polymorphisms of the SLCO1B1 and SLCO1B3 genes. Methods: After a single dose of 160 mg of valsartan was orally administered to healthy male volunteers, drug concentrations were assayed up to 48 h. The 25 subjects were genotyped for 16 single-nucleotide polymorphisms (SNPs) of the SLCO1B1 and SLCO1B3 genes. Subjects were classified into groups according to their SLCO1B1*1B haplotype; 23 subjects were carriers of SLCO1B1*1B and two subjects were included in the reference group with SLCO1B1*1A/*1A. Alternations of the splicing factor-binding site pattern caused by the given mutation were evaluated with the Human Splicing Finder (HSF) 3.1. Results: The subjects who carried SLCO1B1*1B showed a 2.3-fold higher clearance than those without the *1B haplotype. Mean C_max and AUC_inf were reduced by 45% and 54%, respectively, in the SLCO1B1*1B genotype group compared to the reference group with the *1A/*1A genotype (p < 0.01). The carriers of the rs4149153 T allele of SLCO1B3 had a 27% lower mean C_max and a 1.5-fold higher Vd compared to homozygotic CC carriers (p < 0.05). In a combined analysis of SLCO1B1 and SLCO1B3, subjects not carrying SLCO1B1 *1B and carrying SLCO1B3 rs4149153 T allele showed a 1.6-fold higher clearance than those with the other genotypes, whereas mean C_max and AUC_last were reduced by 35% and 42%, respectively (p < 0.05), in the subjects. HSF 3.1 analysis showed that rs4149153 could cause alterations of the acceptor splice site (TAAATACTAAAGAC to TAAATATTAAAGAC) with scoring change (from 72.57 to 71.92, difference = −0.9). Conclusion: It was found that plasma exposure to valsartan is significantly decreased in SLCO1B1*1B carriers and carriers of the rs4149153 T allele of SLCO1B3, possibly as a result of increased hepatic uptake.

Cancer predisposition and germline CTNNA1 variants

Hereditary Diffuse Gastric Cancer (HDGC) is a cancer predisposing syndrome mainly caused by germline inactivating variants in CDH1, encoding E-cadherin. Early-onset diffuse gastric cancer (DGC) and/or invasive lobular breast cancer (LBC) are the main phenotypes in CDH1-associated HDGC. CTNNA1, encoding for α-E-catenin, and E-cadherin-partner in the adherens junction complex, has been recently classified as a HDGC predisposing gene. Nevertheless, little is known about CTNNA1 tumor spectrum in variant carriers and variant-type associated causality. Herein, we systematically reviewed the literature searching for CTNNA1 germline variants carriers, further categorized them according to HDGC clinical criteria (HDGC vs non-HDGC), collected phenotypes, classified variants molecularly and according to CDH1 ACMG/AMP guidelines and performed genotype-phenotype analysis. We found 41 families carrying CTNNA1 germline variants encompassing in total 105 probands and relatives. All probands from 13 HDGC families presented DGC and their average age of onset was 40 ± 17 years; 10/13 (77%) HDGC families carried a pathogenic (P) variant. Most probands from 28 non-HDGC families developed unspecified-BC, as well as most of their relatives; 4/28 (14%) carried a P variant, 16/28 (57%) carried a likely pathogenic (LP) variant, 7/28 (25%) carried variants of unknown significance (VUS) and 1/28 (4%) carried a likely benign variant. Regardless of clinical criteria, 97% (32/33) of probands and relatives from P variant-carrier families had DGC/unspecified-GC. In LP variant-carrier families, 82% (28/34) of probands and relatives had unspecified-BC. Only 2/105 individuals had LBC. A cluster of frameshift and nonsense variants was found in CTNNA1 last exon of non-HDGC families and classified as VUS. In conclusion, current available data confirms an association of CTNNA1 P variants with early-onset DGC, but not with LBC. We demonstrate that in ascertained cohorts, CTNNA1 P variants explain <2% of HDGC families and support the use of ACMG/AMP CDH1 specific variant curation guidelines, while no specific guidelines are developed for CTNNA1 variant classification. Moreover, we demonstrated that truncating variants at the CTNNA1 NMD-incompetent last exon have limited deleteriousness, and that CTNNA1 LP variants have lower actionability than CDH1 LP variants. Current knowledge supports considering only CTNNA1 P variants as clinically actionable in HDGC carrying families.

New Parvoviruses and Picornavirus in Tissues and Feces of Foals with Interstitial Pneumonia

Six foals with interstitial pneumonia of undetermined etiology from Southern California were analyzed by viral metagenomics. Spleen, lung, and colon content samples obtained during necropsy from each animal were pooled, and nucleic acids from virus-like particles enriched for deep sequencing. The recently described equine copiparvovirus named eqcopivirus, as well as three previously uncharacterized viruses, were identified. The complete ORFs genomes of two closely related protoparvoviruses, and of a bocaparvovirus, plus the partial genome of a picornavirus were assembled. The parvoviruses were classified as members of new ungulate protoparvovirus and bocaparvovirus species in the Parvoviridae family. The picornavirus was classified as a new species in the Salivirus genus of the Picornaviridae family. Spleen, lung, and colon content samples from each foal were then tested for these viral genomes by nested PCR and RT-PCR. When present, parvoviruses were detected in both feces and spleen. The picornavirus, protoparvovirus, and eqcopivirus genomes were detected in the lungs of one animal each. Three foals were co-infected with the picornavirus and either a protoparvovirus, bocaparvovirus, or eqcopivirus. Two other foals were infected with a protoparvovirus only. No viral infection was detected in one animal. The complete ORFs of the first equine protoparvoviruses and bocaparvovirus, the partial ORF of the third equine picornavirus, and their detection in tissues of foals with interstitial pneumonia are described here. Testing the involvement of these viruses in fatal interstitial pneumonia or other equine diseases will require larger epidemiological and/or inoculation studies.