Mitochondrial and Nuclear DNA Variants in Amyotrophic Lateral Sclerosis: Enrichment in the Mitochondrial Control Region and Sirtuin Pathway Genes in Spinal Cord Tissue.

Amyotrophic Lateral Sclerosis (ALS) is a progressive disease with prevalent mitochondrial dysfunctions affecting both upper and lower motor neurons in the motor cortex, brainstem, and spinal cord. Despite mitochondria having their own genome (mtDNA), in humans, most mitochondrial genes are encoded by the nuclear genome (nDNA). Our study aimed to simultaneously screen for nDNA and mtDNA genomes to assess for specific variant enrichment in ALS compared to control tissues. Here, we analysed whole exome (WES) and whole genome (WGS) sequencing data from spinal cord tissues, respectively, of 6 and 12 human donors. A total of 31,257 and 301,241 variants in nuclear-encoded mitochondrial genes were identified from WES and WGS, respectively, while mtDNA reads accounted for 73 and 332 variants. Despite technical differences, both datasets consistently revealed a specific enrichment of variants in the mitochondrial Control Region (CR) and in several of these genes directly associated with mitochondrial dynamics or with Sirtuin pathway genes within ALS tissues. Overall, our data support the hypothesis of a variant burden in specific genes, highlighting potential actionable targets for therapeutic interventions in ALS.

Detection of ribonucleotides embedded in DNA by Nanopore sequencing.

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes.

Macroalgal microbiomes unveil a valuable genetic resource for halogen metabolism.

BACKGROUND: Macroalgae, especially reds (Rhodophyta Division) and browns (Phaeophyta Division), are known for producing various halogenated compounds. Yet, the reasons underlying their production and the fate of these metabolites remain largely unknown. Some theories suggest their potential antimicrobial activity and involvement in interactions between macroalgae and prokaryotes. However, detailed investigations are currently missing on how the genetic information of prokaryotic communities associated with macroalgae may influence the fate of organohalogenated molecules. RESULTS: To address this challenge, we created a specialized dataset containing 161 enzymes, each with a complete enzyme commission number, known to be involved in halogen metabolism. This dataset served as a reference to annotate the corresponding genes encoded in both the metagenomic contigs and 98 metagenome-assembled genomes (MAGs) obtained from the microbiome of 2 red (Sphaerococcus coronopifolius and Asparagopsis taxiformis) and 1 brown (Halopteris scoparia) macroalgae. We detected many dehalogenation-related genes, particularly those with hydrolytic functions, suggesting their potential involvement in the degradation of a wide spectrum of halocarbons and haloaromatic molecules, including anthropogenic compounds. We uncovered an array of degradative gene functions within MAGs, spanning various bacterial orders such as Rhodobacterales, Rhizobiales, Caulobacterales, Geminicoccales, Sphingomonadales, Granulosicoccales, Microtrichales, and Pseudomonadales. Less abundant than degradative functions, we also uncovered genes associated with the biosynthesis of halogenated antimicrobial compounds and metabolites. CONCLUSION: The functional data provided here contribute to understanding the still largely unexplored role of unknown prokaryotes. These findings support the hypothesis that macroalgae function as holobionts, where the metabolism of halogenated compounds might play a role in symbiogenesis and act as a possible defense mechanism against environmental chemical stressors. Furthermore, bacterial groups, previously never connected with organohalogen metabolism, e.g., Caulobacterales, Geminicoccales, Granulosicoccales, and Microtrichales, functionally characterized through MAGs reconstruction, revealed a biotechnologically relevant gene content, useful in synthetic biology, and bioprospecting applications. Video Abstract.

OPA1 mutation affects autophagy and triggers senescence in autosomal dominant optic atrophy plus fibroblasts.

In several cases of mitochondrial diseases, the underlying genetic and bioenergetic causes of reduced oxidative phosphorylation (OxPhos) in mitochondrial dysfunction are well understood. However, there is still limited knowledge about the specific cellular outcomes and factors involved for each gene and mutation, which contributes to the lack of effective treatments for these disorders. This study focused on fibroblasts from a patient with Autosomal Dominant Optic Atrophy (ADOA) plus syndrome harboring a mutation in the Optic Atrophy 1 (OPA1) gene. By combining functional and transcriptomic approaches, we investigated the mitochondrial function and identified cellular phenotypes associated with the disease. Our findings revealed that fibroblasts with the OPA1 mutation exhibited a disrupted mitochondrial network and function, leading to altered mitochondrial dynamics and reduced autophagic response. Additionally, we observed a premature senescence phenotype in these cells, suggesting a previously unexplored role of the OPA1 gene in inducing senescence in ADOA plus patients. This study provides novel insights into the mechanisms underlying mitochondrial dysfunction in ADOA plus and highlights the potential importance of senescence in disease progression.

Unraveling C-to-U RNA editing events from direct RNA sequencing.

In mammals, RNA editing events involve the conversion of adenosine (A) in inosine (I) by ADAR enzymes or the hydrolytic deamination of cytosine (C) in uracil (U) by the APOBEC family of enzymes, mostly APOBEC1. RNA editing has a plethora of biological functions, and its deregulation has been associated with various human disorders. While the large-scale detection of A-to-I is quite straightforward using the Illumina RNAseq technology, the identification of C-to-U events is a non-trivial task. This difficulty arises from the rarity of such events in eukaryotic genomes and the challenge of distinguishing them from background noise. Direct RNA sequencing by Oxford Nanopore Technology (ONT) permits the direct detection of Us on sequenced RNA reads. Surprisingly, using ONT reads from wild-type (WT) and APOBEC1-knock-out (KO) murine cell lines as well as in vitro synthesized RNA without any modification, we identified a systematic error affecting the accuracy of the Cs call, thereby leading to incorrect identifications of C-to-U events. To overcome this issue in direct RNA reads, here we introduce a novel machine learning strategy based on the isolation Forest (iForest) algorithm in which C-to-U editing events are considered as sequencing anomalies. Using in vitro synthesized and human ONT reads, our model optimizes the signal-to-noise ratio improving the detection of C-to-U editing sites with high accuracy, over 90% in all samples tested. Our results suggest that iForest, known for its rapid implementation and minimal memory requirements, is a promising tool to denoise ONT reads and reliably identify RNA modifications.

Cyclic AMP induces reversible EPAC1 condensates that regulate histone transcription.

The second messenger cyclic AMP regulates many nuclear processes including transcription, pre-mRNA splicing and mitosis. While most functions are attributed to protein kinase A, accumulating evidence suggests that not all nuclear cyclic AMP-dependent effects are mediated by this kinase, implying that other effectors may be involved. Here we explore the nuclear roles of Exchange Protein Activated by cyclic AMP 1. We find that it enters the nucleus where forms reversible biomolecular condensates in response to cyclic AMP. This phenomenon depends on intrinsically disordered regions present at its amino-terminus and is independent of protein kinase A. Finally, we demonstrate that nuclear Exchange Protein Activated by cyclic AMP 1 condensates assemble at genomic loci on chromosome 6 in the proximity of Histone Locus Bodies and promote the transcription of a histone gene cluster. Collectively, our data reveal an unexpected mechanism through which cyclic AMP contributes to nuclear spatial compartmentalization and promotes the transcription of specific genes.

Structural and Comparative Analyses of Insects Suggest the Presence of an Ultra-Conserved Regulatory Element of the Genes Encoding Vacuolar-Type ATPase Subunits and Assembly Factors.

Gene and genome comparison represent an invaluable tool to identify evolutionarily conserved sequences with possible functional significance. In this work, we have analyzed orthologous genes encoding subunits and assembly factors of the V-ATPase complex, an important enzymatic complex of the vacuolar and lysosomal compartments of the eukaryotic cell with storage and recycling functions, respectively, as well as the main pump in the plasma membrane that energizes the epithelial transport in insects. This study involves 70 insect species belonging to eight insect orders. We highlighted the conservation of a short sequence in the genes encoding subunits of the V-ATPase complex and their assembly factors analyzed with respect to their exon-intron organization of those genes. This study offers the possibility to study ultra-conserved regulatory elements under an evolutionary perspective, with the aim of expanding our knowledge on the regulation of complex gene networks at the basis of organellar biogenesis and cellular organization.

Dysregulation of testis mRNA expression levels in hatchery-produced vs wild greater amberjack Seriola dumerili.

Reproductive dysfunctions have been recently documented in male greater amberjack Seriola dumerili caught from the wild and reared in captivity. In the present study, we compared testis transcriptome in wild fish (WILD), hatchery-produced fish with apparently normal spermatogenesis (Normal Farmed; NormalF) and hatchery-produced fish with evident reproductive dysfunction (Dysfunctional Farmed; DysF). Gene expression analysis identified 2157, 1985 and 74 differentially expressed genes (DEGs) in DysF vs WILD, NormalF vs DysF and NormalF vs WILD comparisons, respectively. In DysF, a dysregulation of several interconnected biological processes, including cell assembly, steroidogenesis and apoptosis was found. Gene enrichment of progesterone-mediated oocyte maturation, oocyte meiosis and cell cycle pathways were identified in the DysF vs NormalF comparison. Most of the DEGs involved in the enriched pathways were downregulated in DysF. The comparison of NormalF vs WILD showed that most of the DEGs were downregulated in NormalF, including a gene that encodes for a regulatory protein with a protective role in apoptosis regulation (ptpn6), indicating that spermatogenesis was dysfunctional also in the apparently "normal" hatchery-produced fish. Hence, rearing of male greater amberjack in captivity, from eggs produced by captive breeders, did not prevent the appearance of reproductive dysfunctions, and these dysfunctions involved several biological processes and metabolic pathways.

The ADAR1 editome reveals drivers of editing-specificity for ADAR1-isoforms.

Adenosine deaminase acting on RNA ADAR1 promotes A-to-I conversion in double-stranded and structured RNAs. ADAR1 has two isoforms transcribed from different promoters: cytoplasmic ADAR1p150 is interferon-inducible while ADAR1p110 is constitutively expressed and primarily localized in the nucleus. Mutations in ADAR1 cause Aicardi - Goutières syndrome (AGS), a severe autoinflammatory disease associated with aberrant IFN production. In mice, deletion of ADAR1 or the p150 isoform leads to embryonic lethality driven by overexpression of interferon-stimulated genes. This phenotype is rescued by deletion of the cytoplasmic dsRNA-sensor MDA5 indicating that the p150 isoform is indispensable and cannot be rescued by ADAR1p110. Nevertheless, editing sites uniquely targeted by ADAR1p150 remain elusive. Here, by transfection of ADAR1 isoforms into ADAR-less mouse cells we detect isoform-specific editing patterns. Using mutated ADAR variants, we test how intracellular localization and the presence of a Z-DNA binding domain-α affect editing preferences. These data show that ZBDα only minimally contributes to p150 editing-specificity while isoform-specific editing is primarily directed by the intracellular localization of ADAR1 isoforms. Our study is complemented by RIP-seq on human cells ectopically expressing tagged-ADAR1 isoforms. Both datasets reveal enrichment of intronic editing and binding by ADAR1p110 while ADAR1p150 preferentially binds and edits 3'UTRs.

Cumulus Cell Transcriptome after Cumulus-Oocyte Complex Exposure to Nanomolar Cadmium in an In Vitro Animal Model of Prepubertal and Adult Age.

Cadmium (Cd), a highly toxic pollutant, impairs oocyte fertilization, through oxidative damage on cumulus cells (CCs). This study analysed the transcriptomic profile of CCs of cumulus-oocyte complexes (COCs) from adult and prepubertal sheep, exposed to Cd nanomolar concentration during in vitro maturation. In both age-groups, CCs of matured oocytes underwent RNA-seq, data analysis and validation. Differentially expressed genes (DEGs) were identified in adult (n = 99 DEGs) and prepubertal (n = 18 DEGs) CCs upon Cd exposure. Transcriptomes of adult CCs clustered separately between Cd-exposed and control samples, whereas prepubertal ones did not as observed by Principal Component Analysis. The transcriptomic signature of Cd-induced CC toxicity was identified by gene annotation and literature search. Genes associated with previous studies on ovarian functions and/or Cd effects were confirmed and new genes were identified, thus implementing the knowledge on their involvement in such processes. Enrichment and validation analysis showed that, in adult CCs, Cd acted as endocrine disruptor on DEGs involved in hormone biosynthesis, cumulus expansion, regulation of cell signalling, growth and differentiation and oocyte maturation, whereas in prepubertal CCs, Cd affected DEGs involved in CC development and viability and CC-oocyte communications. In conclusion, these DEGs could be used as valuable non-invasive biomarkers for oocyte competence.

Profiling RNA Editing in Single Cells.

RNA editing is a widespread molecular phenomenon occurring in a variety of organisms. In humans, it mainly involves the deamination of adenosine to inosine (A-to-I) in double-stranded RNAs by ADAR enzymes. A-to-I RNA editing has been investigated in different tissues as well as in diverse experimental and pathological conditions. By contrast, its biological role in single cells has not been explored in depth. Recent methodologies for cell sorting in combination with deep sequencing technologies have enabled the study of eukaryotic transcriptomes at single cell resolution, paving the way to the profiling of their epitranscriptomic dynamics.Here we describe a step-by-step protocol to detect and characterize A-to-I events occurring in publicly available single-cell RNAseq experiments from human alpha and beta pancreatic cells.

UTRdb 2.0: a comprehensive, expert curated catalog of eukaryotic mRNAs untranslated regions.

The 5' and 3' untranslated regions of eukaryotic mRNAs (UTRs) play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization, and message stability. Since 1996, we have developed and maintained UTRdb, a specialized database of UTR sequences. Here we present UTRdb 2.0, a major update of UTRdb featuring an extensive collection of eukaryotic 5' and 3' UTR sequences, including over 26 million entries from over 6 million genes and 573 species, enriched with a curated set of functional annotations. Annotations include CAGE tags and polyA signals to label the completeness of 5' and 3'UTRs, respectively. In addition, uORFs and IRES are annotated in 5'UTRs as well as experimentally validated miRNA targets in 3'UTRs. Further annotations include evolutionarily conserved blocks, Rfam motifs, ADAR-mediated RNA editing events, and m6A modifications. A web interface allowing a flexible selection and retrieval of specific subsets of UTRs, selected according to a combination of criteria, has been implemented which also provides comprehensive download facilities. UTRdb 2.0 is accessible at http://utrdb.cloud.ba.infn.it/utrdb/.

Mitofusin 2 mutation drives cell proliferation in Charcot-Marie-Tooth 2A fibroblasts.

Dominant mutations in ubiquitously expressed mitofusin 2 gene (MFN2) cause Charcot-Marie-Tooth type 2A (CMT2A; OMIM 609260), an inherited sensory-motor neuropathy that affects peripheral nerve axons. Mitofusin 2 protein has been found to take part in mitochondrial fusion, mitochondria-endoplasmic reticulum tethering, mitochondrial trafficking along axons, mitochondrial quality control and various types of cancer, in which MFN2 has been indicated as a tumor suppressor gene. Discordant data on the mitochondrial altered phenotypes in patient-derived fibroblasts harboring MFN2 mutations and in animal models have been reported. We addressed some of these issues by focusing on mitochondria behavior during autophagy and mitophagy in fibroblasts derived from a CMT2AMFN2 patient with an MFN2650G > T/C217F mutation in the GTPase domain. This study investigated mitochondrial dynamics, respiratory capacity and autophagy/mitophagy, to tackle the multifaceted MFN2 contribution to CMT2A pathogenesis. We found that MFN2 mutated fibroblasts showed impairment of mitochondrial morphology, bioenergetics capacity, and impairment of the early stages of autophagy, but not mitophagy. Unexpectedly, transcriptomic analysis of mutated fibroblasts highlighted marked differentially expressed pathways related to cell population proliferation and extracellular matrix organization. We consistently found the activation of mTORC2/AKT signaling and accelerated proliferation in the CMT2AMFN2 fibroblasts. In conclusion, our evidence indicates that MFN2 mutation can positively drive cell proliferation in CMT2AMFN2 fibroblasts.

Comparative transcriptomic profiling of peach and nectarine cultivars reveals cultivar-specific responses to chilled postharvest storage.

Introduction: Peach (Prunus persica (L.) Batsch,) and nectarine fruits (Prunus persica (L.) Batsch, var nectarine), are characterized by a rapid deterioration at room temperature. Therefore, cold storage is widely used to delay fruit post-harvest ripening and extend fruit commercial life. Physiological disorders, collectively known as chilling injury, can develop typically after 3 weeks of low-temperature storage and affect fruit quality. Methods: A comparative transcriptomic analysis was performed to identify regulatory pathways that develop before chilling injury symptoms are detectable using next generation sequencing on the fruits of two contrasting cultivars, one peach (Sagittaria) and one nectarine, (Big Top), over 14 days of postharvest cold storage. Results: There was a progressive increase in the number of differentially expressed genes between time points (DEGs) in both cultivars. More (1264) time point DEGs were identified in 'Big Top' compared to 'Sagittaria' (746 DEGs). Both cultivars showed a downregulation of pathways related to photosynthesis, and an upregulation of pathways related to amino sugars, nucleotide sugar metabolism and plant hormone signal transduction with ethylene pathways being most affected. Expression patterns of ethylene related genes (including biosynthesis, signaling and ERF transcription factors) correlated with genes involved in cell wall modification, membrane composition, pathogen and stress response, which are all involved later during storage in development of chilling injury. Discussion: Overall, the results show that common pathways are activated in the fruit of 'Big Top' nectarine and 'Sagittaria' peach in response to cold storage but include also differences that are cultivar-specific responses.

Correction: Tusup et al. Evaluation of the Interplay between the ADAR Editome and Immunotherapy in Melanoma. Non-Coding RNA 2021, 7, 5.

In the original article [...].

The A to I editing landscape in melanoma and its relation to clinical outcome.

RNA editing refers to non-transient RNA modifications that occur after transcription and prior to translation by the ribosomes. RNA editing is more widespread in cancer cells than in non-transformed cells and is associated with tumorigenesis of various cancer tissues. However, RNA editing can also generate neo-antigens that expose tumour cells to host immunosurveillance. Global RNA editing in melanoma and its relevance to clinical outcome currently remain poorly characterized. The present study compared RNA editing as well as gene expression in tumour cell lines from melanoma patients of short or long metastasis-free survival, patients relapsing or not after immuno- and targeted therapy and tumours harbouring BRAF or NRAS mutations. Overall, our results showed that NTRK gene expression can be a marker of resistance to BRAF and MEK inhibition and gives some insights of candidate genes as potential biomarkers. In addition, this study revealed an increase in Adenosine-to-Inosine editing in Alu regions and in non-repetitive regions, including the hyperediting of the MOK and DZIP3 genes in relapsed tumour samples during targeted therapy and of the ZBTB11 gene in NRAS mutated melanoma cells. Therefore, RNA editing could be a promising tool for identifying predictive markers, tumour neoantigens and targetable pathways that could help in preventing relapses during immuno- or targeted therapies.

Yeast as a Model to Unravel New BRCA2 Functions in Cell Metabolism.

Mutations in BRCA2 gene increase the risk for breast cancer and for other cancer types, including pancreatic and prostate cancer. Since its first identification as an oncosupressor in 1995, the best-characterized function of BRCA2 is in the repair of DNA double-strand breaks (DSBs) by homologous recombination. BRCA2 directly interacts with both RAD51 and single-stranded DNA, mediating loading of RAD51 recombinase to sites of single-stranded DNA. In the absence of an efficient homologous recombination pathway, DSBs accumulate resulting in genome instability, thus supporting tumorigenesis. Yet the precise mechanism by which BRCA2 exerts its tumor suppressor function remains unclear. BRCA2 has also been involved in other biological functions including protection of telomere integrity and stalled replication forks, cell cycle progression, transcriptional control and mitophagy. Recently, we and others have reported a role of BRCA2 in modulating cell death programs through a molecular mechanism conserved in yeast and mammals. Here we hypothesize that BRCA2 is a multifunctional protein which exerts specific functions depending on cell stress response pathway. Based on a differential RNA sequencing analysis carried out on yeast cells either growing or undergoing a regulated cell death process, either in the absence or in the presence of BRCA2, we suggest that BRCA2 causes central carbon metabolism reprogramming in response to death stimuli and encourage further investigation on the role of metabolic reprogramming in BRCA2 oncosuppressive function.

A Machine Learning Approach to Parkinson's Disease Blood Transcriptomics.

The increased incidence and the significant health burden associated with Parkinson's disease (PD) have stimulated substantial research efforts towards the identification of effective treatments and diagnostic procedures. Despite technological advancements, a cure is still not available and PD is often diagnosed a long time after onset when irreversible damage has already occurred. Blood transcriptomics represents a potentially disruptive technology for the early diagnosis of PD. We used transcriptome data from the PPMI study, a large cohort study with early PD subjects and age matched controls (HC), to perform the classification of PD vs. HC in around 550 samples. Using a nested feature selection procedure based on Random Forests and XGBoost we reached an AUC of 72% and found 493 candidate genes. We further discussed the importance of the selected genes through a functional analysis based on GOs and KEGG pathways.