Nonsense-mediated mRNA decay: terminating erroneous gene expression

Unraveling novel mRNA transcripts of the human DNA N-glycosylase 1 (NTHL1) gene with the implementation of an innovative targeted DNA-seq assay

The human NTHL1 gene encodes a DNA glycosylase that plays a key role in the base excision repair (BER) pathway, repairing oxidative DNA damage and maintaining genome integrity. The physiological activity of NTHL1 is crucial in preventing genetic alterations that can lead to cancer. In this study, we employed an innovative targeted DNA sequencing (DNA-seq) methodology to explore the transcriptional landscape of the NTHL1 gene, revealing previously uncharacterized alternative splicing events and novel exons. Our designed approach provided significantly improved sequencing depth and coverage, enabling the identification of novel NTHL1 mRNA transcripts. Bioinformatics analysis confirmed all annotated splice junctions of the main NTHL1 transcripts (v.1 - v.3) and revealed novel mRNA transcripts (NTHL1 v.4 - v.9) derived from splicing events between annotated exons as well as mRNAs containing previously uncharacterized exons (NTHL1 v.10 - v.14). Quantitative PCR analysis highlighted a diverse expression pattern of these novel transcripts across different human cell lines, suggesting cell-specific roles and regulatory mechanisms. Notably, NTHL1 v.5 was overexpressed in luminal A breast cancer cells (MCF-7), while v.13 was prominent in triple negative (BT-20), HER2 + breast cancer (SK-BR-3), prostate, colorectal cancer cells and HEK-293 cells. Our findings suggest that specific novel NTHL1 transcripts may encode protein isoforms with distinct structural features, as indicated by ribosome profiling datasets, while others containing premature termination codons could function as long non-coding RNAs. These insights enhance our understanding of NTHL1 regulatory role and its potential as a biomarker and therapeutic target in human malignancies. This study underscores the importance of exploring the transcriptional diversity of NTHL1 to fully elucidate its role in cancer pathobiology.

Preclinical alternative drug discovery programs for monogenic rare diseases. Should small molecules or gene therapy be used? The case of hereditary spastic paraplegias

Patients diagnosed with rare diseases and their and families search desperately to organize drug discovery campaigns. Alternative models that differ from default paradigms offer real opportunities. There are, however, no clear guidelines for the development of such models, which reduces success rates and raises costs. We address the main challenges in making the discovery of new preclinical treatments more accessible, using rare hereditary paraplegia as a paradigmatic case. First, we discuss the necessary expertise, and the patients' clinical and genetic data. Then, we revisit gene therapy, de novo drug development, and drug repurposing, discussing their applicability. Moreover, we explore a pool of recommended in silico tools for pathogenic variant and protein structure prediction, virtual screening, and experimental validation methods, discussing their strengths and weaknesses. Finally, we focus on successful case applications.

A PNPLA3-Deficient iPSC-Derived Hepatocyte Screen Identifies Pathways to Potentially Reduce Steatosis in Metabolic Dysfunction-Associated Fatty Liver Disease

The incidence of nonalcoholic fatty liver disease (NAFLD), or metabolic dysfunction-associated fatty liver disease (MAFLD), is increasing in adults and children. Unfortunately, effective pharmacological treatments remain unavailable. Single nucleotide polymorphisms (SNPs) in the patatin-like phospholipase domain-containing protein (PNPLA3 I148M) have the most significant genetic association with the disease at all stages of its progression. A roadblock to identifying potential treatments for PNPLA3-induced NAFLD is the lack of a human cell platform that recapitulates the PNPLA3 I148M-mediated onset of lipid accumulation. Hepatocyte-like cells were generated from PNPLA3-/- and PNPLA3I148M/M-induced pluripotent stem cells (iPSCs). Lipid levels were measured by staining with BODIPY 493/503 and were found to increase in PNPLA3 variant iPSC-derived hepatocytes. A small-molecule screen identified multiple compounds that target Src/PI3K/Akt signaling and could eradicate lipid accumulation in these cells. We found that drugs currently in clinical trials for cancer treatment that target the same pathways also reduced lipid accumulation in PNPLA3 variant cells.

Foxp2 deficiency impairs reproduction by modulating the hypothalamic-pituitary-gonadal axis in zebrafish†

FOXP2 was initially characterized as a transcription factor linked to speech and language disorders. Single-cell RNA sequencing reveals that Foxp2 is enriched in the gonadotrope cluster of the pituitary gland and colocalized with the hormones LHB and FSHB in chickens and mice, implying that FOXP2 might be associated with reproduction in vertebrates. Herein, we investigated the roles of foxp2 in reproduction in a Foxp2-deficient zebrafish model. The results indicated that the loss of Foxp2 inhibits courtship behavior in adult male zebrafish. Notably, Foxp2 deficiency disrupts gonad development, leading to retardation of follicle development and a decrease in oocytes in females at the full-growth stage, among other phenotypes. The transcriptome analysis (RNA-seq) also revealed that differentially expressed genes clustered into the estrogen signaling and ovarian steroidogenesis-related signaling pathways. In addition, we found that Foxp2 deficiency could modulate the hypothalamic-pituitary-gonadal axis, especially the regulation of lhb and fshb expression, in zebrafish. In contrast, the injection of human chorionic gonadotropin, a specific LH agonist, partially rescues Foxp2-impaired reproduction in zebrafish, suggesting that Foxp2 plays an important role in the regulation of reproduction via the hypothalamic-pituitary-gonadal axis in zebrafish. Thus, our findings reveal a new role for Foxp2 in the regulation of reproduction in vertebrates.

Identification of diverse RNA viruses in Obscuromonas flagellates (Euglenozoa: Trypanosomatidae: Blastocrithidiinae)

Trypanosomatids (Euglenozoa) are a diverse group of unicellular flagellates predominately infecting insects (monoxenous species) or circulating between insects and vertebrates or plants (dixenous species). Monoxenous trypanosomatids harbor a wide range of RNA viruses belonging to the families Narnaviridae, Totiviridae, Qinviridae, Leishbuviridae, and a putative group of tombus-like viruses. Here, we focus on the subfamily Blastocrithidiinae, a previously unexplored divergent group of monoxenous trypanosomatids comprising two related genera: Obscuromonas and Blastocrithidia. Members of the genus Blastocrithidia employ a unique genetic code, in which all three stop codons are repurposed to encode amino acids, with TAA also used to terminate translation. Obscuromonas isolates studied here bear viruses of three families: Narnaviridae, Qinviridae, and Mitoviridae. The latter viral group is documented in trypanosomatid flagellates for the first time. While other known mitoviruses replicate in the mitochondria, those of trypanosomatids appear to reside in the cytoplasm. Although no RNA viruses were detected in Blastocrithidia spp., we identified an endogenous viral element in the genome of B. triatomae indicating its past encounter(s) with tombus-like viruses.

Use of steric blocking antisense oligonucleotides for the targeted inhibition of junction containing precursor microRNAs

Antisense oligonucleotides (ASOs) are widely used as therapeutics for neurodegenerative diseases, cancers, and virus infections. One class of ASOs functions to enhance protein expression by sequestering the mature microRNA (miRNA) in a double-stranded structure within the RNA-induced silencing complex (RISC). An alternative approach for the targeted control of gene expression is to use ASOs that bind to the pre-elements of miRNAs (pre-miRNAs) and modulate their enzymatic processing. Here, we demonstrate that ASOs can be used to disrupt a specific structural feature, "junction," within pre-miR-31 that is important in directing efficient processing by the Dicer/TRBP complex. Furthermore, we extend and validate this strategy to pre-miR-144, which has a similar junction-dependent structure-function relationship. We found that a significant number of human pre-miRNAs are predicted to contain junctions, and validated our ASO approach on several members of this group. Importantly, we also verified the application of junction-targeting ASOs for the specific inhibition of pre-miRNA processing in cell . Our study reemphasizes the important roles of RNA structure in regulating Dicer/TRBP processing of pre-miRNAs and provides the framework to develop structure-informed ASOs that serve to inhibit miRNA production.

A premature stop codon in the CYP2C19 gene may explain the unexpected sensitivity of vultures to diclofenac toxicity

The unintended environmental exposure of vultures to diclofenac has resulted in the deaths of millions of old-world vultures on the Asian subcontinent. While toxicity has been since associated with a long half-life of elimination and zero order metabolism, the actual constraint in biotransformation is yet to be clarified. For this study we evaluated if the evident zero order metabolism could be due to defects in the CYP2C9/2C19 enzyme system. For this, using whole genome sequencing and de-novo transcriptome alignment, the vulture CYP2C19 open reading frame was identified through Splign analysis. The result sequence analysis revealed the presence of a premature stop codon on intron 7 of the identified open reading frame. Even if the stop codon was not present, amino acid residue analysis tended to suggest that the enzyme would be lower in activity than the equivalent human enzyme, with differences present at sites 105, 286 and 289. The defect was also conserved across the eight non-related vultures tested. From these results, we conclude that the sensitivity of the old-world vultures to diclofenac is due to the non-expression of a viable CYP2C19 enzyme system. This is not too dissimilar to the effects seen in certain people with a similar defective enzyme.

A frameshift deletion in F8 associated with hemophilia A in Labrador Retriever dogs

Hemophilia A is the most common inherited coagulation factor disorder in dogs. It manifests as excessive bleeding resulting from pathogenic variants in the X-chromosomal F8 gene encoding coagulation factor VIII (FVIII) protein. In this study, we performed careful clinical phenotyping to confirm hemophilia A in two distinct Labrador Retriever (LR) pedigrees. Whole-genome sequencing on an affected dog from litter 1 identified a case-specific frameshift deletion variant in F8 predicted to cause a premature stop codon (c.2923_2924del, p.(E975Kfs*8)). This variant was hemizygous in all the affected males from litter 1 (n = 3), while all the unaffected LRs in the pedigree were heterozygous or wild-type (n = 22). Additionally, screened samples from 199 LRs were all found to be wild-type. As a result of this study, a gene test can now be developed to screen dogs before breeding to prevent further cases. However, it is important to note that the affected LR with decreased FVIII activity from litter 2 was wild-type for the identified deletion variant, and no segregating F8 variants were detected when this dog's DNA sample was whole-genome sequenced. Thus, the cause of decreased FVIII activity in this dog remains to be unraveled in future studies.

A single amino acid change led to structural and functional differentiation of PvHd1 to control flowering in switchgrass

Switchgrass, a forage and bioenergy crop, occurs as two main ecotypes with different but overlapping ranges of adaptation. The two ecotypes differ in a range of characteristics, including flowering time. Flowering time determines the duration of vegetative development and therefore biomass accumulation, a key trait in bioenergy crops. No causal variants for flowering time differences between switchgrass ecotypes have, as yet, been identified. In this study, we mapped a robust flowering time quantitative trait locus (QTL) on chromosome 4K in a biparental F2 population and characterized the flowering-associated transcription factor gene PvHd1, an ortholog of CONSTANS in Arabidopsis and Heading date 1 in rice, as the underlying causal gene. Protein modeling predicted that a serine to glycine substitution at position 35 (p.S35G) in B-Box domain 1 greatly altered the global structure of the PvHd1 protein. The predicted variation in protein compactness was supported in vitro by a 4 °C shift in denaturation temperature. Overexpressing the PvHd1-p.35S allele in a late-flowering CONSTANS-null Arabidopsis mutant rescued earlier flowering, whereas PvHd1-p.35G had a reduced ability to promote flowering, demonstrating that the structural variation led to functional divergence. Our findings provide us with a tool to manipulate the timing of floral transition in switchgrass cultivars and, potentially, expand their cultivation range.

A temperature-tolerant CRISPR base editor mediates highly efficient and precise gene editing in Drosophila

CRISPR nucleases generate a broad spectrum of mutations that includes undesired editing outcomes. Here, we develop optimized C-to-T base editing systems for the generation of precise loss- or gain-of-function alleles in Drosophila and identify temperature as a crucial parameter for efficiency. We find that a variant of the widely used APOBEC1 deaminase has attenuated activity at 18° to 29°C and shows considerable dose-dependent toxicity. In contrast, the temperature-tolerant evoCDA1 domain mediates editing of typically more than 90% of alleles and is substantially better tolerated. Furthermore, formation of undesired mutations is exceptionally rare in Drosophila compared to other species. The predictable editing outcome, high efficiency, and product purity enables near homogeneous induction of STOP codons or alleles encoding protein variants in vivo. Last, we demonstrate how optimized expression enables conditional base editing in marked cell populations. This work substantially facilitates creation of precise alleles in Drosophila and provides key design parameters for developing efficient base editing systems in other ectothermic species.

GABRA1-Related Disorders: From Genetic to Functional Pathways

OBJECTIVE

Variants in GABRA1 have been associated with a broad epilepsy spectrum, ranging from genetic generalized epilepsies to developmental and epileptic encephalopathies. However, our understanding of what determines the phenotype severity and best treatment options remains inadequate. We therefore aimed to analyze the electroclinical features and the functional effects of GABRA1 variants to establish genotype-phenotype correlations.

METHODS

Genetic and electroclinical data of 27 individuals (22 unrelated and 2 families) harboring 20 different GABRA1 variants were collected and accompanied by functional analysis of 19 variants.

RESULTS

Individuals in this cohort could be assigned into different clinical subgroups based on the functional effect of their variant and its structural position within the GABRA1 subunit. A homogenous phenotype with mild cognitive impairment and infantile onset epilepsy (focal seizures, fever sensitivity, and electroencephalographic posterior epileptiform discharges) was described for variants in the extracellular domain and the small transmembrane loops. These variants displayed loss-of-function (LoF) effects, and the patients generally had a favorable outcome. A more severe phenotype was associated with variants in the pore-forming transmembrane helices. These variants displayed either gain-of-function (GoF) or LoF effects. GoF variants were associated with severe early onset neurodevelopmental disorders, including early infantile developmental and epileptic encephalopathy.

INTERPRETATION

Our data expand the genetic and phenotypic spectrum of GABRA1 epilepsies and permit delineation of specific subphenotypes for LoF and GoF variants, through the heterogeneity of phenotypes and variants. Generally, variants in the transmembrane helices cause more severe phenotypes, in particular GoF variants. These findings establish the basis for a better understanding of the pathomechanism and a precision medicine approach in GABRA1-related disorders. Further studies in larger populations are needed to provide a conclusive genotype-phenotype correlation. ANN NEUROL 2023.

Salt-inducible kinase inhibition promotes the adipocyte thermogenic program and adipose tissue browning

OBJECTIVE

Norepinephrine stimulates the adipose tissue thermogenic program through a β-adrenergic receptor (βAR)-cyclic adenosine monophosphate (cAMP)-protein kinase A (PKA) signaling cascade. We discovered that a noncanonical activation of the mechanistic target of rapamycin complex 1 (mTORC1) by PKA is required for the βAR-stimulation of adipose tissue browning. However, the downstream events triggered by PKA-phosphorylated mTORC1 activation that drive this thermogenic response are not well understood.

METHODS

We used a proteomic approach of Stable Isotope Labeling by/with Amino acids in Cell culture (SILAC) to characterize the global protein phosphorylation profile in brown adipocytes treated with the βAR agonist. We identified salt-inducible kinase 3 (SIK3) as a candidate mTORC1 substrate and further tested the effect of SIK3 deficiency or SIK inhibition on the thermogenic gene expression program in brown adipocytes and in mouse adipose tissue.

RESULTS

SIK3 interacts with RAPTOR, the defining component of the mTORC1 complex, and is phosphorylated at Ser884 in a rapamycin-sensitive manner. Pharmacological SIK inhibition by a pan-SIK inhibitor (HG-9-91-01) in brown adipocytes increases basal Ucp1 gene expression and restores its expression upon blockade of either mTORC1 or PKA. Short-hairpin RNA (shRNA) knockdown of Sik3 augments, while overexpression of SIK3 suppresses, Ucp1 gene expression in brown adipocytes. The regulatory PKA phosphorylation domain of SIK3 is essential for its inhibition. CRISPR-mediated Sik3 deletion in brown adipocytes increases type IIa histone deacetylase (HDAC) activity and enhances the expression of genes involved in thermogenesis such as Ucp1, Pgc1α, and mitochondrial OXPHOS complex protein. We further show that HDAC4 interacts with PGC1α after βAR stimulation and reduces lysine acetylation in PGC1α. Finally, a SIK inhibitor well-tolerated in vivo (YKL-05-099) can stimulate the expression of thermogenesis-related genes and browning of mouse subcutaneous adipose tissue.

CONCLUSIONS

Taken together, our data reveal that SIK3, with the possible contribution of other SIKs, functions as a phosphorylation switch for β-adrenergic activation to drive the adipose tissue thermogenic program and indicates that more work to understand the role of the SIKs is warranted. Our findings also suggest that maneuvers targeting SIKs could be beneficial for obesity and related cardiometabolic disease.

Rapid Identification by Resequencing-Based QTL Mapping of a Novel Allele RGA1-FH Decreasing Grain Length in a Rice Restorer Line 'Fuhui212'

Grain size is one of the most frequently selected traits during domestication and modern breeding. The continued discovery and characterization of new genes and alleles in controlling grain size are important in safeguarding the food supply for the world's growing population. Previously, a small grain size was observed in a rice restorer line 'Fuhui212', while the underlying genetic factors controlling this trait were unknown. In this study, by combining QTL mapping, variant effect prediction, and complementation experiments, we recovered a novel allele RGA1-FH that explains most of the phenotypic changes. The RGA1-FH allele contains an A-to-T splicing site variant that disrupts the normal function of RGA1. While population analysis suggests extremely strong artificial selection in maintaining a functional allele of RGA1, our study is the first, to the best of our knowledge, to prove that a dysfunctional RGA1 allele can also be beneficial in real agricultural production. Future breeding programs would benefit from paying more attention to the rational utilization of those overlooked 'unfavored' alleles.

Genome analyses reveal population structure and a purple stigma color gene candidate in finger millet

Finger millet is a key food security crop widely grown in eastern Africa, India and Nepal. Long considered a 'poor man's crop', finger millet has regained attention over the past decade for its climate resilience and the nutritional qualities of its grain. To bring finger millet breeding into the 21^st century, here we present the assembly and annotation of a chromosome-scale reference genome. We show that this ~1.3 million years old allotetraploid has a high level of homoeologous gene retention and lacks subgenome dominance. Population structure is mainly driven by the differential presence of large wild segments in the pericentromeric regions of several chromosomes. Trait mapping, followed by variant analysis of gene candidates, reveals that loss of purple coloration of anthers and stigma is associated with loss-of-function mutations in the finger millet orthologs of the maize R1/B1 and Arabidopsis GL3/EGL3 anthocyanin regulatory genes. Proanthocyanidin production in seed is not affected by these gene knockouts.

Identification of neoantigens and immunological subtypes in clear cell renal cell carcinoma for mRNA vaccine development and patient selection

Clear cell renal cell carcinoma (ccRCC) is a common urological malignancy with diverse histological types. This study aimed to detect neoantigens in ccRCC to develop mRNA vaccines and distinguish between ccRCC immunological subtypes for construction of an immune landscape to select patients suitable for vaccination. Using The Cancer Genome Atlas SpliceSeq database, The Cancer Genome Atlas, and the International Cancer Genome Consortium cohorts, we comprehensively analysed potential tumour antigens of ccRCC associated with aberrant alternative splicing, somatic mutation, nonsense-mediated mRNA decay factors, antigen-presenting cells, and overall survival. Immune subtypes (C1/C2) and nine immune gene modules of ccRCC were identified by consistency clustering and weighted correlation network analysis. The immune landscape as well as molecular and cellular characteristics of immunotypes were assessed. Rho-guanine nucleotide exchange factor 3 (ARHGEF3) was identified as a new ccRCC antigen for development of an mRNA vaccine. A higher tumour mutation burden, differential expression of immune checkpoints, and immunogenic cell death were observed in cases with the C2 immunotype. Cellular characteristics increased the complexity of the immune environment, and worse outcomes were observed in ccRCC cases with the C2 immunotype. We constructed the immune landscape for selecting patients with the C2 immunotype suitable for vaccination.

The catalytic-dead Pcif1 regulates gene expression and fertility in Drosophila

Eukaryotic mRNAs are modified at the 5' end with a methylated guanosine (m7G) that is attached to the transcription start site (TSS) nucleotide. The TSS nucleotide is 2'-O-methylated (Nm) by CMTR1 in organisms ranging from insects to human. In mammals, the TSS adenosine can be further N 6 -methylated by RNA polymerase II phosphorylated CTD-interacting factor 1 (PCIF1) to create m6Am. Curiously, the fly ortholog of mammalian PCIF1 is demonstrated to be catalytic-dead, and its functions are not known. Here, we show that Pcif1 mutant flies display a reduced fertility which is particularly marked in females. Deep sequencing analysis of Pcif1 mutant ovaries revealed transcriptome changes with a notable increase in expression of genes belonging to the mitochondrial ATP synthetase complex. Furthermore, the Pcif1 protein is distributed along euchromatic regions of polytene chromosomes, and the Pcif1 mutation behaved as a modifier of position-effect-variegation (PEV) suppressing the heterochromatin-dependent silencing of the white gene. Similar or stronger changes in the transcriptome and PEV phenotype were observed in flies that expressed a cytosolic version of Pcif1. These results point to a nuclear cotranscriptional gene regulatory role for the catalytic-dead fly Pcif1 that is probably based on its conserved ability to interact with the RNA polymerase II carboxy-terminal domain.

A frame-shift mutation in COMTD1 is associated with impaired pheomelanin pigmentation in chicken

The biochemical pathway regulating the synthesis of yellow/red pheomelanin is less well characterized than the synthesis of black/brown eumelanin. Inhibitor of gold (IG phenotype) is a plumage colour variant in chicken that provides an opportunity to further explore this pathway since the recessive allele (IG) at this locus is associated with a defect in the production of pheomelanin. IG/IG homozygotes display a marked dilution of red pheomelanin pigmentation, whilst black pigmentation (eumelanin) is only slightly affected. Here we show that a 2-base pair insertion (frame-shift mutation) in the 5th exon of the Catechol-O-methyltransferase containing domain 1 gene (COMTD1), expected to cause a complete or partial loss-of-function of the COMTD1 enzyme, shows complete concordance with the IG phenotype within and across breeds. We show that the COMTD1 protein is localized to mitochondria in pigment cells. Knockout of Comtd1 in a mouse melanocytic cell line results in a reduction in pheomelanin metabolites and significant alterations in metabolites of glutamate/glutathione, riboflavin, and the tricarboxylic acid cycle. Furthermore, COMTD1 overexpression enhanced cellular proliferation following chemical-induced transfection, a potential inducer of oxidative stress. These observations suggest that COMTD1 plays a protective role for melanocytes against oxidative stress and that this supports their ability to produce pheomelanin.

EPS8 variant causes deafness, autosomal recessive 102 (DFNB102) and literature review

Pathogenic variants in the EPS8 gene result in nonsyndromic hearing loss. This gene encodes the EPS8 protein in cochlear inner hair cells and performs critical roles in stimulating actin polymerization and bundling. Thus far, only four pathogenic variations in EPS8 have been described. In this study, we report the fifth pathogenic variant in the EPS8 gene in an Iranian patient with DFNB102. Furthermore, we review literature cases with EPS8 mutations.

Genomic analysis of a novel pathogenic variant in the gene LMNA associated with cardiac laminopathies found in Ecuadorian siblings: A case report

Introduction

Cardiac laminopathies are caused by mutations in the LMNA gene and include a wide range of clinical manifestations involving electrical and mechanical changes in cardiomyocytes. In Ecuador, cardiovascular diseases were the primary cause of death in 2019, accounting for 26.5% of total deaths. Cardiac laminopathy-associated mutations involve genes coding for structural proteins with functions related to heart development and physiology.

Family description

Two Ecuadorian siblings, self-identified as mestizos, were diagnosed with cardiac laminopathies and suffered embolic strokes. Moreover, by performing Next-Generation Sequencing, a pathogenic variant (NM_170707.3:c.1526del) was found in the gene LMNA.

Discussion and conclusion

Currently, genetic tests are an essential step for disease genetic counseling, including cardiovascular disease diagnosis. Identification of a genetic cause that may explain the risk of cardiac laminopathies in a family can help the post-test counseling and recommendations from the cardiologist. In the present report, a pathogenic variant ((NM_170707.3:c.1526del) has been identified in two Ecuadorian siblings with cardiac laminopathies. The LMNA gene codes for A-type laminar proteins that are associated with gene transcription regulation. Mutations in the LMNA gene cause laminopathies, disorders with diverse phenotypic manifestations. Moreover, understanding the molecular biology of the disease-causing mutations is essential in deciding the correct type of treatment.

Characterization of cellular phenotypes in neurons derived from induced pluripotent stem cells of male patients with Fabry disease

Fabry disease (FD) is an X-linked inherited lysosomal metabolism disorder in which globotriaosylceramide (Gb3) accumulates in various organs resulting from a deficiency in alpha-galactosidase A. The clinical features of FD include progressive impairments of the renal, cardiac, and peripheral nervous systems. In addition, patients with FD often develop neuropsychiatric symptoms, such as depression and dementia, which are believed to be induced by the cellular injury of cerebrovascular and partially neuronal cells due to Gb3 accumulation. Although the analysis of autopsy brain tissue from patients with FD showed no accumulation of Gb3, abnormal deposits of Gb3 were found in the neurons of several brain areas, including the hippocampus. Therefore, in this study, we generated induced pluripotent stem cells (iPSCs) from patients with FD and differentiated them into neuronal cells to investigate pathological and biological changes in the neurons of FD. Neural stem cells (NSCs) and neurons were successfully differentiated from the iPSCs we generated; however, cellular damage and morphological changes were not found in these cells. Immunostaining revealed no Gb3 accumulation in NSCs and neurons. Transmission electron microscopy did not reveal any zebra body-like structures or inclusion bodies, which are characteristic of FD. These results indicated that neuronal cells derived from FD-iPSCs exhibited normal morphology and no Gb3 accumulation. It is likely that more in vivo environment-like cultures are needed for iPSC-derived neurons to reproduce disease-specific features.

Novel clinical presentation and PAX6 mutation in families with congenital aniridia

Purpose

To explore the clinical phenotype and genetic defects of families with congenital aniridia.

Methods

Four Chinese families with aniridia were enrolled in this study. The detailed ocular presentations of the patients were recorded. Whole exome sequencing (BGI MGIEasy V4 chip) was used to detect the gene mutation. Sanger sequencing was performed to validate the potential pathogenic variants, and segregation analysis was performed on all available family members.

Results

By whole exome sequencing and Sanger sequencing, three recurrent mutations (c.112del, p.Arg38Glyfs*16; c.299G > A, p.Trp100* and c.718C > T, p.Arg240*) and one novel mutation (c.278_281del, p.Glu93Alafs*30) of PAX6 were identified. All the mutations were co-segregated with the phenotype in the families. We also observed spontaneous anterior lens capsule rupture in aniridia for the first time.

Conclusion

We report spontaneous anterior lens capsule rupture as a novel phenotype of aniridia and three recurrent mutations and one novel mutation of PAX6 in families with aniridia. Our results expanded the phenotype and genotype spectra of aniridia and can help us better understand the disease.

Genetic landscape in coagulation factor XIII associated defects – Advances in coagulation and beyond

Coagulation factor XIII (FXIII) acts as a fine fulcrum in blood plasma that maintains the balance between bleeding and thrombosis by covalently crosslinking the pre-formed fibrin clot into an insoluble one that is resistant to premature fibrinolysis. In plasma, FXIII circulates as a pro-transglutaminase complex composed of the dimeric catalytic FXIII-A encoded by the F13A1 gene and dimeric carrier/regulatory FXIII-B subunits encoded by the F13B gene. Growing evidence accumulated over decades of exhaustive research shows that not only does FXIII play major roles in both pathological extremes of hemostasis i.e. bleeding and thrombosis, but that it is, in fact, a pleiotropic protein with physiological roles beyond coagulation. However, the current FXIII genetic-epidemiological literature is overwhelmingly derived from the bleeding pathology associated with its deficiency. In this article we review the current clinical, functional, and molecular understanding of this fascinating multifaceted protein, especially putting into the same perspective its genetic landscape.

Computational biology insights into genotype-clinical phenotype-protein phenotype relationships between novel SLC26A2 variants identified in inherited skeletal dysplasias

BACKGROUND

Pathogenic variants in the transmembrane sulfate transporter protein SLC26A2 are associated with different phenotypes of inherited chondrodysplasias. As limited data is published from India, in this study we sought to elucidate the molecular basis of inherited chondrodysplasias in an Indian cohort.

METHODS

Molecular screening of 32 fetuses with antenatally diagnosed lethal skeletal dysplasia was performed by next generation sequencing and Sanger sequencing. The genotype-protein phenotype characterization was done using computational biology techniques like homology modelling, stability and pathogenicity predictions.

RESULTS

We identified five rare autosomal recessive SLC26A2 [NM_000112.4] variants, including three homozygous c.796dupA(p.Thr266Asnfs*12), c.1724delA(p.Lys575Serfs*10), and c.1375_1377dup(p.Val459dup) and two heterozygous variants (c.532C > T(p.Arg178*)) and (c.1382C > T(p.Ala461Val)) in compound heterozygous form in a total of four foetuses. Genotype-protein phenotype annotations highlighted that the clinically severe achondrogenesis 1B causative c.796dupA(p.Thr266Asnfs*12) and c.1724delA(p.Lys575Serfs*10)variants impact SLC26A2 protein structure by deletion of the protein core and transmembrane STAS domains, respectively. In clinically moderate atelosteogenesis type 2 phenotype, the c.1382C > T(p.Ala461Val) variant is predicted to distort alpha helix conformation and alter the bonding properties and free energy dynamics of transmembrane domains and the c.532C > T(p.Arg178*) variant results in loss of both core transmembrane and STAS domains of the SLC26A2 protein. The c.1375_1377dup(p.Val459dup) variant identified in clinically milder atelosteogenesis type II-diastrophic dysplasia spectrum lethal phenotype is predicted to decrease the Qualitative Model Energy Analysis (QMean), which affects major geometrical aspects of the SLC26A2 protein structure.

CONCLUSION

We expand the spectrum of SLC26A2 related lethal chondrodysplasia and report three novel variants correlating clinical severity and protein phenotype within the lethal spectrum of this rare dysplasia. We demonstrate the relevance of structural characterization to aid novel variant reclassification to provide better prenatal management and reproductive options to families with lethal antenatal skeletal disorder.

USP37 Deubiquitinates CDC73 in HPT-JT Syndrome

The CDC73/HRPT2 gene, a defect which causes hyperparathyroidism–jaw tumor (HPT-JT) syndrome, encodes CDC73/parafibromin. We aimed to investigate whether CDC73 would be a target for ubiquitin–proteasome degradation. We cloned full-length cDNAs encoding a family of 58 ubiquitin-specific deubiquitinating enzymes (DUBs), also known as ubiquitin-specific proteases (USPs). Use of the yeast two-hybrid system then enabled us to identify USP37 as interacting with CDC73. The biochemical interaction between the USP37 and CDC73 and their reciprocal binding domains were studied. Co-localization of CDC73 and USP37 was observed in cells. CDC73 was found to be polyubiquitinated, and polyubiquitination of CDC73 was prominent in mutants. CDC73 was deubiquitinated via K48-specific ubiquitin chains by USP37, but not by the catalytically inactive USP37^C350S mutant. Observation of the binding between deletion mutants of CDC73 and USP37 revealed that the β-catenin binding site of CDC73 and the ubiquitin-interacting motifs 2 and 3 (UIM2 and 3) of USP37 were responsible for the interaction between the two proteins. Moreover, these two enzymes co-existed within the nucleus of COS7 cells. We conclude that USP37 is a DUB for CDC73 and that the two proteins interact through specific domains, suggesting that USP37 is responsible for the stability of CDC73 in HPT-JT syndrome.

A non-canonical Raf function is required for dorsal-ventral patterning during Drosophila embryogenesis

Proper embryonic development requires directional axes to pattern cells into embryonic structures. In Drosophila, spatially discrete expression of transcription factors determines the anterior to posterior organization of the early embryo, while the Toll and TGFβ signalling pathways determine the early dorsal to ventral pattern. Embryonic MAPK/ERK signaling contributes to both anterior to posterior patterning in the terminal regions and to dorsal to ventral patterning during oogenesis and embryonic stages. Here we describe a novel loss of function mutation in the Raf kinase gene, which leads to loss of ventral cell fates as seen through the loss of the ventral furrow, the absence of Dorsal/NFκB nuclear localization, the absence of mesoderm determinants Twist and Snail, and the expansion of TGFβ. Gene expression analysis showed cells adopting ectodermal fates much like loss of Toll signaling. Our results combine novel mutants, live imaging, optogenetics and transcriptomics to establish a novel role for Raf, that appears to be independent of the MAPK cascade, in embryonic patterning.

Functional validation of variants of unknown significance using CRISPR gene editing and transcriptomics: A Kleefstra syndrome case study

There are an estimated > 400 million people living with a rare disease globally, with genetic variants the cause of approximately 80% of cases. Next Generation Sequencing (NGS) rapidly identifies genetic variants however they are often of unknown significance. Low throughput functional validation in specialist laboratories is the current ad hoc approach for functional validation of genetic variants, which creating major bottlenecks in patient diagnosis. This study investigates the application of CRISPR gene editing followed by genome wide transcriptomic profiling to facilitate patient diagnosis. As proof-of-concept, we introduced a variant in the Euchromatin histone methyl transferase (EHMT1) gene into HEK293T cells. We identified changes in the regulation of the cell cycle, neural gene expression and suppression of gene expression changes on chromosome 19 and chromosome X, that are in keeping with Kleefstra syndrome clinical phenotype and/or provide insight into disease mechanism. This study demonstrates the utility of genome editing followed by functional readouts to rapidly and systematically validating the function of variants of unknown significance in patients suffering from rare diseases.

Does size matter? Two new deletions in the HBB gene cause β0-thalassemia

Most β-thalassemias are caused by mutations involving one or a limited number of nucleotides within the gene or its adjacent regions. They can be substitutions or deletions; in these cases, the loss ranges from a single nucleotide to even the entire HBB gene, so we wonder if the phenotype is due to the size of the deletion or the location of the mutation. To clarify this, we present two new deletions in the β-globin gene that cause β⁰-thalassemia. The hematological parameters were determined with an automated cell counter; the Hb A2 and Hb F levels were measured by performance liquid chromatography. Hemoglobins were analyzed by capillary zone electrophoresis (Sebia Capillarys Flex system) and ion-exchange HPLC (BioRad Variant II β-thalassemia Short Program). Molecular characterization was performed by automatic Sanger sequencing. The screening of common α-thalassemia point mutations and deletions in the world (21 in total) were carried out using multiplex PCR followed by reverse-hybridization with a commercial Alpha-Globin StripAssay kit. We have characterized two new mutations-(1) 1-bp deletion [CD61/62(-G)] [HBB:c.186_187delG], (2) 105-bp deletion [IVS-2-nt767-CD111] [HBB:c.316-84_333del]-and we have described, for first time in Spain, the 25-bp deletion [β nts 252 - 276 deleted] [HBB:c.93-22_95del] mutation. These mutations were classified as pathogenic by UniProt Variants confirmed according to the American College of Medical Genetics and Genomics guidelines. These mutations present a phenotype compatible with β⁰-thalassemia, supported by hematological parameters that correlate the degree of reduction in the synthesis of the β-globin chain. Identification of this type of mutation is important for genetic counselling of partners where both are carriers, so that they are aware of the genetic risk of having affected children, allowing them to take an informed decision about their reproductive choices.

CRISPR single base editing, neuronal disease modelling and functional genomics for genetic variant analysis: pipeline validation using Kleefstra syndrome EHMT1 haploinsufficiency

Background

Over 400 million people worldwide are living with a rare disease. Next Generation Sequencing (NGS) identifies potential disease causative genetic variants. However, many are identified as variants of uncertain significance (VUS) and require functional laboratory validation to determine pathogenicity, and this creates major diagnostic delays.

Methods

In this study we test a rapid genetic variant assessment pipeline using CRISPR homology directed repair to introduce single nucleotide variants into inducible pluripotent stem cells (iPSCs), followed by neuronal disease modelling, and functional genomics on amplicon and RNA sequencing, to determine cellular changes to support patient diagnosis and identify disease mechanism.

Results

As proof-of-principle, we investigated an EHMT1 (Euchromatin histone methyltransferase 1; EHMT1 c.3430C > T; p.Gln1144*) genetic variant pathogenic for Kleefstra syndrome and determined changes in gene expression during neuronal progenitor cell differentiation. This pipeline rapidly identified Kleefstra syndrome in genetic variant cells compared to healthy cells, and revealed novel findings potentially implicating the key transcription factors REST and SP1 in disease pathogenesis.

Conclusion

The study pipeline is a rapid, robust method for genetic variant assessment that will support rare diseases patient diagnosis. The results also provide valuable information on genome wide perturbations key to disease mechanism that can be targeted for drug treatments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02740-3.

KIAA0319 influences cilia length, cell migration and mechanical cell-substrate interaction

Following its association with dyslexia in multiple genetic studies, the KIAA0319 gene has been extensively investigated in different animal models but its function in neurodevelopment remains poorly understood. We developed the first human cellular knockout model for KIAA0319 in RPE1 retinal pigment epithelia cells via CRISPR-Cas9n to investigate its role in processes suggested but not confirmed in previous studies, including cilia formation and cell migration. We observed in the KIAA0319 knockout increased cilia length and accelerated cell migration. Using Elastic Resonator Interference Stress Microscopy (ERISM), we detected an increase in cellular force for the knockout cells that was restored by a rescue experiment. Combining ERISM and immunostaining we show that RPE1 cells exert highly dynamic, piconewton vertical pushing forces through actin-rich protrusions that are surrounded by vinculin-rich pulling sites. This protein arrangement and force pattern has previously been associated to podosomes in other cells. KIAA0319 depletion reduces the fraction of cells forming these actin-rich protrusions. Our results suggest an involvement of KIAA0319 in cilia biology and cell-substrate force regulation.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Efficient knockout of the phytoene desaturase gene in a hybrid poplar (Populus alba × Populus glandulosa) using the CRISPR/Cas9 system with a single gRNA

The CRISPR/Cas9 system has been used for genome editing in several plant species; however, there are few reports on its use in trees. Here, CRISPR/Cas9 was used to mutate a target gene in Populus alba × Populus glandulosa hybrid poplars. The hybrid poplar is routinely used in molecular biological studies due to the well-established Agrobacterium-mediated transformation method. A single guide RNA (sgRNA) with reported high mutation efficiency in other popular species was designed with a protospacer adjacent motif sequence for the phytoene desaturase 1 (PagPDS1) gene. The pHSE/Cas9-PagPDS1 sgRNA vector was delivered into hybrid poplar cells using Agrobacterium-mediated transformation. The transgenic plants were propagated and classified them into three groups according to their phenotypes. Among a total of 110 lines of transgenic hybrid poplars, 82 lines showed either an albino or a pale green phenotype, indicating around 74.5% phenotypic mutation efficiency of the PagPDS1 gene. The albino phenotypes were observed when the CRISPR/Cas9-mediated mutations in both PagPDS1 alleles in the transgenic plants. There was no off-target modification of the PagPDS2 gene, which has a potential sgRNA target sequence with two mismatches. The results confirmed that the sgRNA can specifically edit PagPDS1 rather than PagPDS2, indicating that CRISPR/Cas9-mediated genome editing can effectively induce target mutations in the hybrid poplar. This technique will be useful to improve tree quality in hybrid poplars (P. alba × P. glandulosa); for example, by enhancing biomass or stress tolerance.

A novel CEP290 disease-causing variant identified in a patient with leber congenital amaurosis using a medical diagnostic panel sequencing

BACKGROUND

This study aims to identify the underlying genetic cause of a Chinese patient with Leber congenital amaurosis (LCA).

METHODS

Detailed clinical data and family history were collected. A medical diagnostic panel sequencing covering 4450 genes was conducted. Two candidate disease-causing mutations detected in CEP290 were then validated with Sanger sequencing and bioinformatic analysis. Reverse transcription polymerase chain reaction (RT-PCR) and cDNA sequencing were performed to understand the effect of the novel CEP290 mutation on CEP290 mRNA splicing.

RESULTS

A five-month-old LCA patient with both parents was enrolled. Medical diagnostic panel sequencing revealed that the patient is a compound heterozygote for two potentially pathogenic CEP290 mutations. Among them, c.1666dupA (p.I556NfsX20) was previously reported and has a significant association with LCA phenotype. A novel CEP290 mutation (c.3310-1_3313delGCTTA) was confirmed in both the patient and her father. RT-PCR and cDNA sequencing confirmed that the novel mutation affects the CEP290 mRNA splicing and results in a complete skipping of exon 29. The frameshift resulted in an early stop codon and truncation of the mutant protein by 1371 amino acid residues (p.L1104fsX6).

CONCLUSIONS

Our report provided a new mutation to the spectrum of CEP290-related diseases. The data suggested that the c.3310-1_3313delGCTTA mutation affects the CEP290 mRNA splicing and the CEP290 protein function. This valuable information is important for future studies on the mRNA splicing of CEP290 and the pathogenesis of CEP290-related diseases.

Parental Origin of the RB1 Gene Mutations in Families with Low Penetrance Hereditary Retinoblastoma

Simple Summary

Some families with hereditary retinoblastoma exhibit mild phenotype with low penetrance and variable expressivity, including complete absence of clinical signs of the disease in some carriers of the germline RB1 mutation. The identification of low-penetrance mutations in the RB1 gene and the study of their inheritance in pedigrees is contributing to understanding the mechanisms underlying the development of retinoblastoma with low penetrance. It is important both for further expansion of knowledge in the field of molecular genetics of retinoblastoma, and for competent genetic counseling and subsequent clinical management of families with this form of the disease. Our results support an assumption that parental origin of an RB1 mutation influences the likelihood of developing retinoblastoma. We also revealed a relatively high frequency of asymptomatic carriage of the RB1 mutations among the parents of retinoblastoma patients, highlighting the utmost necessity for molecular analysis among the probands’ relatives irrespective of their clinical status and family history of retinoblastoma.

Abstract

Our aim was to identify RB1 alterations causing hereditary low penetrance retinoblastoma and to evaluate how the parental origin of an RB1 mutation affects its phenotypic expression. By NGS and MLPA, RB1 mutations were found in 191 from 332 unrelated retinoblastoma patients. Among patients with identified RB1 mutations but without clinical family history of retinoblastoma, 7% (12/175) were found to have hereditary disease with one of the parents being an asymptomatic carrier of an RB1 mutation. Additionally, in two families with retinoblastoma history, mutations were inherited by probands from unaffected parents. Overall, nine probands inherited RB1 mutations from clinically unaffected fathers and five, from mothers. Yet, we gained explanations of maternal “unaffectedness” in most cases, either as somatic mosaicism or as clinical presentation of retinomas in involution, rendering the proportion of paternal to maternal truly asymptomatic mutation carriers as 9:1 (p = 0.005). This observation supports an assumption that parental origin of an RB1 mutation influences the likelihood of developing retinoblastoma. Additionally, our study revealed a relatively high frequency of asymptomatic carriage of the RB1 mutations among the parents of retinoblastoma patients, highlighting the utmost necessity of molecular analysis among the probands’ relatives irrespective of their clinical status and family history of retinoblastoma.

CRISPR/Cas9-based knockout reveals that the clock gene timeless is indispensable for regulating circadian behavioral rhythms in Bombyx mori

Circadian rhythms, which are ubiquitous and adaptive, occur across all species, from microbes to humans, in which they organize and modify behavior and physiology. timeless (tim) is a canonical clock gene. The core composition of the Drosophila melanogaster endogenous circadian clock has been extensively investigated; however, in lepidopteran insects, including Bombyx mori, the mechanism is complicated and little is known regarding the participation of tim in the negative feedback loop responsible for behavioral activities. To arrive at a comprehensive understanding of the role of tim in the B. mori endogenous circadian clock, we exploited the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 gene editing system. We attempted to elucidate the functions of tim in the circadian clock of B. mori using Bmtim mutants. The knockouts affected two circadian behavioral activities: adult emergence and embryo hatching rhythms. Quantitative real-time polymerase chain reaction results confirmed that tim-knockouts induced relative reductions in the expression levels, and thereby the oscillation amplitudes, of Bmper and Bmclk messenger RNAs during both the photophase and scotophase. Additionally, the daily rhythmic expression of Bmdbt was upregulated in the photophase and downregulated in the scotophase in a tim-knockout. Our study reveals that tim is integral to the B. mori circadian clock and may be involved in regulating eclosion and hatching rhythms.

CRISPR-Knockout of CSE Gene Improves Saccharification Efficiency by Reducing Lignin Content in Hybrid Poplar

Caffeoyl shikimate esterase (CSE) has been shown to play an important role in lignin biosynthesis in plants and is, therefore, a promising target for generating improved lignocellulosic biomass crops for sustainable biofuel production. Populus spp. has two CSE genes (CSE1 and CSE2) and, thus, the hybrid poplar (Populus alba × P. glandulosa) investigated in this study has four CSE genes. Here, we present transgenic hybrid poplars with knockouts of each CSE gene achieved by CRISPR/Cas9. To knockout the CSE genes of the hybrid poplar, we designed three single guide RNAs (sg1-sg3), and produced three different transgenic poplars with either CSE1 (CSE1-sg2), CSE2 (CSE2-sg3), or both genes (CSE1/2-sg1) mutated. CSE1-sg2 and CSE2-sg3 poplars showed up to 29.1% reduction in lignin deposition with irregularly shaped xylem vessels. However, CSE1-sg2 and CSE2-sg3 poplars were morphologically indistinguishable from WT and showed no significant differences in growth in a long-term living modified organism (LMO) field-test covering four seasons. Gene expression analysis revealed that many lignin biosynthetic genes were downregulated in CSE1-sg2 and CSE2-sg3 poplars. Indeed, the CSE1-sg2 and CSE2-sg3 poplars had up to 25% higher saccharification efficiency than the WT control. Our results demonstrate that precise editing of CSE by CRISPR/Cas9 technology can improve lignocellulosic biomass without a growth penalty.

Rare deleterious BUB1B variants induce premature ovarian insufficiency and early menopause

Losing of ovarian functions prior to natural menopause age causes female infertility and early menopause. Premature ovarian insufficiency (POI) is defined as the loss of ovarian activity before 40 years of age. Known genetic causes account for 25-30% of POI cases, demonstrating the high genetic heterogeneity of POI and the necessity for further genetic explorations. Here we conducted genetic analyses using whole-exome sequencing in a Chinese non-syndromic POI family with the affected mother and at least four affected daughters. Intriguingly, a rare missense variant of BUB1B c.273A>T (p.Gln91His) was shared by all the cases in this family. Furthermore, our replication study using targeted sequencing revealed a novel stop-gain variant of BUB1B c.1509T>A (p.Cys503*) in one of 200 sporadic POI cases. Both heterozygous BUB1B variants were evaluated to be deleterious by multiple in silico tools. BUB1B encodes BUBR1, a crucial spindle assembly checkpoint component involved in cell division. BUBR1 insufficiency may induce vulnerability to oxidative stress. Therefore, we generated a mouse model with a loss-of-function mutant of Bub1b, and also employed D-galactose-induced aging assays for functional investigations. Notably, Bub1b+/- female mice presented late-onset subfertility, and they were more sensitive to oxidative stress than wild-type female controls, mimicking the clinical phenotypes of POI cases affected by deleterious BUB1B variants. Our findings in human cases and mouse models consistently suggest, for the first time, that heterozygous deleterious variants of BUB1B are involved in late-onset POI and related disorders.

Mirtron-mediated RNA knockdown/replacement therapy for the treatment of dominant retinitis pigmentosa

Rhodopsin (RHO) gene mutations are a common cause of autosomal dominant retinitis pigmentosa (ADRP). The need to suppress toxic protein expression together with mutational heterogeneity pose challenges for treatment development. Mirtrons are atypical RNA interference effectors that are spliced from transcripts as short introns. Here, we develop a novel mirtron-based knockdown/replacement gene therapy for the mutation-independent treatment of RHO-related ADRP, and demonstrate efficacy in a relevant mammalian model. Splicing and potency of rhodopsin-targeting candidate mirtrons are initially determined, and a mirtron-resistant codon-modified version of the rhodopsin coding sequence is validated in vitro. These elements are then combined within a single adeno-associated virus (AAV) and delivered subretinally in a Rho^P23H knock-in mouse model of ADRP. This results in significant mouse-to-human rhodopsin RNA replacement and is associated with a slowing of retinal degeneration. This provides proof of principle that synthetic mirtrons delivered by AAV are capable of reducing disease severity in vivo.

TF-RBP-AS Triplet Analysis Reveals the Mechanisms of Aberrant Alternative Splicing Events in Kidney Cancer: Implications for Their Possible Clinical Use as Prognostic and Therapeutic Biomarkers

Aberrant alternative splicing (AS) is increasingly linked to cancer; however, how AS contributes to cancer development still remains largely unknown. AS events (ASEs) are largely regulated by RNA-binding proteins (RBPs) whose ability can be modulated by a variety of genetic and epigenetic mechanisms. In this study, we used a computational framework to investigate the roles of transcription factors (TFs) on regulating RBP-AS interactions. A total of 6519 TF–RBP–AS triplets were identified, including 290 TFs, 175 RBPs, and 16 ASEs from TCGA–KIRC RNA sequencing data. TF function categories were defined according to correlation changes between RBP expression and their targeted ASEs. The results suggested that most TFs affected multiple targets, and six different classes of TF-mediated transcriptional dysregulations were identified. Then, regulatory networks were constructed for TF–RBP–AS triplets. Further pathway-enrichment analysis showed that these TFs and RBPs involved in triplets were enriched in a variety of pathways that were associated with cancer development and progression. Survival analysis showed that some triplets were highly associated with survival rates. These findings demonstrated that the integration of TFs into alternative splicing regulatory networks can help us in understanding the roles of alternative splicing in cancer.

Bi-allelic variants in the ER quality-control mannosidase gene EDEM3 cause a congenital disorder of glycosylation

EDEM3 encodes a protein that converts Man8GlcNAc2 isomer B to Man7-5GlcNAc2. It is involved in the endoplasmic reticulum-associated degradation pathway, responsible for the recognition of misfolded proteins that will be targeted and translocated to the cytosol and degraded by the proteasome. In this study, through a combination of exome sequencing and gene matching, we have identified seven independent families with 11 individuals with bi-allelic protein-truncating variants and one individual with a compound heterozygous missense variant in EDEM3. The affected individuals present with an inherited congenital disorder of glycosylation (CDG) consisting of neurodevelopmental delay and variable facial dysmorphisms. Experiments in human fibroblast cell lines, human plasma, and mouse plasma and brain tissue demonstrated decreased trimming of Man8GlcNAc2 isomer B to Man7GlcNAc2, consistent with loss of EDEM3 enzymatic activity. In human cells, Man5GlcNAc2 to Man4GlcNAc2 conversion is also diminished with an increase of Glc1Man5GlcNAc2. Furthermore, analysis of the unfolded protein response showed a reduced increase in EIF2AK3 (PERK) expression upon stimulation with tunicamycin as compared to controls, suggesting an impaired unfolded protein response. The aberrant plasma N-glycan profile provides a quick, clinically available test for validating variants of uncertain significance that may be identified by molecular genetic testing. We propose to call this deficiency EDEM3-CDG.

Novel rare variants in FGFR1 and clinical characteristics analysis in a series of congenital hypogonadotropic hypogonadism patients

OBJECTIVE

We aimed to analyse FGFR1 rare variants in a series of Chinese congenital hypogonadotropic hypogonadism (CHH) patients. In addition, we intended to understand the clinical characteristics and the response to treatment of CHH patients with FGFR1 rare variants.

PATIENTS AND METHODS

A total of 357 CHH patients were recruited at Peking Union Medical College Hospital. We used Sanger sequencing to analyse FGFR1 gene. In silico analysis was carried out to study the pathogenicity of novel missense variants. The clinical, endocrinological and therapeutic effects from patients carrying FGFR1 rare variants were analysed retrospectively.

RESULTS

Thimissense mutations.rty patients in this series were found to harbour 29 FGFR1 rare variants, with 8 recurrent and 21 novel variants. After comprehensive analysis, 18 out of 21 novel variants were classified as likely pathogenic (LP) ones. These variants are widely spread throughout the FGFR1 gene and almost all FGFR1 functional domains, which exhibited no hot spot. Cryptorchidism, cleft palate and dental abnormality incidence in this CHH series that possessed FGFR1 LP variants were approximately 38.5%, 7.6% and 3.8%, respectively. Among patients who accepted the fertility-promoting treatment, 8 out of 10 patients succeeded in spermatogenesis.

CONCLUSIONS

Eighteen novel LP variants were found to expand the spectrum of FGFR1 rare variants. In CHH patients possessing FGFR1 variants, we found that the rate of spermatogenesis was high following fertility-promoting therapy and the existence of cryptorchidism may represent the underlying factors which affect spermatogenesis.

Differential expression of Dp71 and Dp40 isoforms in proliferating and differentiated neural stem cells: Identification of Dp40 splicing variants

Dp71 and Dp40 are the main products of the DMD gene in the central nervous system, and they are developmentally regulated from the early stages of embryonic development to adulthood. To further study the roles of Dp71 and Dp40 during cell proliferation and neural differentiation, we analyzed Dp71/Dp40 isoform expression at the mRNA level by RT-PCR assays to identify alternative splicing (AS) in the isoforms expressed in rat neural stem/progenitor cells (NSPCs) and in differentiated cells (neurons and glia). We found that proliferating NSPCs expressed Dp71d, Dp71dΔ71, Dp71f, Dp71fΔ71, Dp71dΔ74 and Dp40, as well as two Dp40 isoforms: Dp40Δ63,64 and Dp40Δ64-67. In differentiated cells we also found the expression of Dp71d, Dp71dΔ71, Dp71f, Dp71fΔ71 and Dp40. However, the expression frequencies were different in both stages. In addition, in differentiated cells, we found Dp71fΔ71-74, and interestingly, we did not find the expression of Dp71dΔ74 or the newly identified Dp40 isoforms. In this work we show that NSPC differentiation is accompanied by changes in Dp71/Dp40 isoform expression, suggesting different roles for these isoforms in NSPCs proliferation and neuronal differentiation, and we describe, for the first time, alternative splicing of Dp40.

Rcl1 depletion impairs 18S pre-rRNA processing at the A1-site and up-regulates a cohort of ribosome biogenesis genes in zebrafish

Yeast Rcl1 is a potential endonuclease that mediates pre-RNA cleavage at the A₂-site to separate 18S rRNA from 5.8S and 25S rRNAs. However, the biological function of Rcl1 in opisthokonta is poorly defined. Moreover, there is no information regarding the exact positions of 18S pre-rRNA processing in zebrafish. Here, we report that zebrafish pre-rRNA harbours three major cleavage sites in the 5′ETS, namely –477nt (A′-site), –97nt (A₀-site) and the 5′ETS and 18S rRNA link (A₁-site), as well as two major cleavage regions within the ITS1, namely 208–218nt (site 2) and 20–33nt (site E). We also demonstrate that depletion of zebrafish Rcl1 mainly impairs cleavage at the A₁-site. Phenotypically, rcl1^–/– mutants exhibit a small liver and exocrine pancreas and die before 15 days post-fertilization. RNA-seq analysis revealed that the most significant event in rcl1^–/– mutants is the up-regulated expression of a cohort of genes related to ribosome biogenesis and tRNA production. Our data demonstrate that Rcl1 is essential for 18S rRNA maturation at the A₁-site and for digestive organogenesis in zebrafish. Rcl1 deficiency, similar to deficiencies in other ribosome biogenesis factors, might trigger a common mechanism to upregulate the expression of genes responsible for ribosome biogenesis.

Structure and genetic variants of thyroglobulin: Pathophysiological implications

Thyroglobulin (TG) plays a main role in the biosynthesis of thyroid hormones (TH), and, thus, it is involved in a wide range of vital functions throughout the life cycle of all vertebrates. Deficiency of TH production due to TG genetic variants causes congenital hypothyroidism (CH), with devastating consequences such as intellectual disability and impaired growth if untreated. To this day, 229 variations in the human TG gene have been identified while the 3D structure of TG has recently appeared. Although TG deficiency is thought to be of autosomal recessive inheritance, the introduction of massive sequencing platforms led to the identification of a variety of monoallelic TG variants (combined with mutations in other thyroid gene products) opening new questions regarding the possibility of oligogenic inheritance of the disease. In this review we discuss remarkable advances in the understanding of the TG architecture and the pathophysiology of CH associated with TG defects, providing new insights for the management of congenital disorders as well as counseling benefits for families with a history of TG abnormalities. Moreover, we summarize relevant aspects of TH synthesis within TG and offer an updated analysis of animal and cellular models of TG deficiency for pathophysiological studies of thyroid dyshormonogenesis while highlighting perspectives for new investigations. All in all, even though there has been sustained progress in understanding the role of TG in thyroid pathophysiology during the past 50 years, functional characterization of TG variants remains an important area of study for future advancement in the field.

Next-Generation Sequencing Identifies Pathogenic Variants in HGF, POU3F4, TECTA, and MYO7A in Consanguineous Pakistani Deaf Families

Background

Approximately 70% of congenital deafness is attributable to genetic causes. Incidence of congenital deafness is known to be higher in families with consanguineous marriage. In this study, we investigated the genetic causes in three consanguineous Pakistani families segregating with prelingual, severe-to-profound deafness.

Results

Through targeted next-generation sequencing of 414 genes known to be associated with deafness, homozygous variants c.536del (p. Leu180Serfs∗20) in TECTA, c.3719 G>A (p. Arg1240Gln) in MYO7A, and c.482+1986_1988del in HGF were identified as the pathogenic causes of enrolled families. Interestingly, in one large consanguineous family, an additional c.706G>A (p. Glu236Lys) variant in the X-linked POU3F4 gene was also identified in multiple affected family members causing deafness. Genotype-phenotype cosegregation was confirmed in all participating family members by Sanger sequencing.

Conclusions

Our results showed that the genetic causes of deafness are highly heterogeneous. Even within a single family, the affected members with apparently indistinguishable clinical phenotypes may have different pathogenic variants.

Sycp1 Is Not Required for Subtelomeric DNA Double-Strand Breaks but Is Required for Homologous Alignment in Zebrafish Spermatocytes

In meiotic prophase I, homologous chromosomes are bound together by the synaptonemal complex, in which two axial elements are connected by transverse filaments and central element proteins. In human and zebrafish spermatocytes, homologous recombination and assembly of the synaptonemal complex initiate predominantly near telomeres. In mice, synapsis is not required for meiotic double-strand breaks (DSBs) and homolog alignment but is required for DSB repair; however, the interplay of these meiotic events in the context of peritelomeric bias remains unclear. In this study, we identified a premature stop mutation in the zebrafish gene encoding the transverse filament protein Sycp1. In sycp1 mutant zebrafish spermatocytes, axial elements were formed and paired at chromosome ends between homologs during early to mid-zygonema. However, they did not synapse, and their associations were mostly lost in late zygotene- or pachytene-like stages. In sycp1 mutant spermatocytes, γH2AX signals were observed, and Dmc1/Rad51 and RPA signals appeared predominantly near telomeres, resembling wild-type phenotypes. We observed persistent localization of Hormad1 along the axis in sycp1 mutant spermatocytes, while the majority of Iho1 signals appeared and disappeared with kinetics similar to those in wild-type spermatocytes. Notably, persistent Iho1 foci were observed in spo11 mutant spermatocytes, suggesting that Iho1 dissociation from axes occurs in a DSB-dependent manner. Our results demonstrated that Sycp1 is not required for peritelomeric DSB formation but is necessary for complete pairing of homologs in zebrafish meiosis.

Loss-of-Function Variants in the SYNPO2L Gene Are Associated With Atrial Fibrillation

Multiple genome-wide association studies (GWAS) have identified numerous loci associated with atrial fibrillation (AF). However, the genes driving these associations and how they contribute to the AF pathogenesis remains poorly understood. To identify genes likely to be driving the observed association, we searched the FinnGen study consisting of 12,859 AF cases and 73,341 controls for rare genetic variants predicted to cause loss-of-function. A specific splice site variant was found in the SYNPO2L gene, located in an AF associated locus on chromosome 10. This variant was associated with an increased risk of AF with a relatively high odds ratio of 3.5 (p = 9.9 × 10-8). SYNPO2L is an important gene involved in the structural development and function of the cardiac myocyte and our findings thus support the recent suggestions that AF can present as atrial cardiomyopathy.

A Novel Titin Truncation Variant Linked to Familial Dilated Cardiomyopathy Found in a Japanese Family and Its Functional Analysis in Genome-Edited Model Cells

Dilated cardiomyopathy (DCM) is a common cause of heart failure. TTN, which encodes titin protein, is a representative causative gene of DCM, and is presented mainly as a truncation variant. However, TTN truncation variants are also found in healthy individuals, and it is therefore important to evaluate the pathogenicity of each variant. In this study, we analyzed 67 cardiomyopathy-associated genes in a male Japanese patient who was hospitalized for recurrent severe heart failure and identified a novel truncation variant, TTN Ser17456Arg fs*14. This TTN truncation variant was located in the A-band region. Moreover, the patient's mother with heart failure harbored the same variant, whereas the father and brother without heart failure did not harbor the variant. To examine the functional changes associated with the truncation variant, H9c2 cells were subjected to genome editing to generate cells with a homologous truncation variant. The cells were differentiated using all-trans-retinoic acid, and the mRNA expression of skeletal actin and cardiac actin were found to be increased and decreased, respectively, consistent with known changes in patients with DCM or heart failure. In contrast, another cell with the titin truncation variant used as a control showed no changes in heart failure-related genes. In summary, we found a novel TTN truncation variant in familial DCM patients and confirmed its functional changes using a relatively simple cell model. The novel truncation variant was identified as a pathogenic and disease-causing mutation.