Natural proteome diversity links aneuploidy tolerance to protein turnover.

Accessing the natural genetic diversity of species unveils hidden genetic traits, clarifies gene functions and allows the generalizability of laboratory findings to be assessed. One notable discovery made in natural isolates of Saccharomyces cerevisiae is that aneuploidy-an imbalance in chromosome copy numbers-is frequent1,2 (around 20%), which seems to contradict the substantial fitness costs and transient nature of aneuploidy when it is engineered in the laboratory3-5. Here we generate a proteomic resource and merge it with genomic1 and transcriptomic6 data for 796 euploid and aneuploid natural isolates. We find that natural and lab-generated aneuploids differ specifically at the proteome. In lab-generated aneuploids, some proteins-especially subunits of protein complexes-show reduced expression, but the overall protein levels correspond to the aneuploid gene dosage. By contrast, in natural isolates, more than 70% of proteins encoded on aneuploid chromosomes are dosage compensated, and average protein levels are shifted towards the euploid state chromosome-wide. At the molecular level, we detect an induction of structural components of the proteasome, increased levels of ubiquitination, and reveal an interdependency of protein turnover rates and attenuation. Our study thus highlights the role of protein turnover in mediating aneuploidy tolerance, and shows the utility of exploiting the natural diversity of species to attain generalizable molecular insights into complex biological processes.

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.

Data-driven recombination detection in viral genomes.

Recombination is a key molecular mechanism for the evolution and adaptation of viruses. The first recombinant SARS-CoV-2 genomes were recognized in 2021; as of today, more than ninety SARS-CoV-2 lineages are designated as recombinant. In the wake of the COVID-19 pandemic, several methods for detecting recombination in SARS-CoV-2 have been proposed; however, none could faithfully confirm manual analyses by experts in the field. We hereby present RecombinHunt, an original data-driven method for the identification of recombinant genomes, capable of recognizing recombinant SARS-CoV-2 genomes (or lineages) with one or two breakpoints with high accuracy and within reduced turn-around times. ReconbinHunt shows high specificity and sensitivity, compares favorably with other state-of-the-art methods, and faithfully confirms manual analyses by experts. RecombinHunt identifies recombinant viral genomes from the recent monkeypox epidemic in high concordance with manually curated analyses by experts, suggesting that our approach is robust and can be applied to any epidemic/pandemic virus.

Ancient and recent origins of shared polymorphisms in yeast.

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

The pancancer overexpressed NFYC Antisense 1 controls cell cycle mitotic progression through in cis and in trans modes of action.

Antisense RNAs (asRNAs) represent an underappreciated yet crucial layer of gene expression regulation. Generally thought to modulate their sense genes in cis through sequence complementarity or their act of transcription, asRNAs can also regulate different molecular targets in trans, in the nucleus or in the cytoplasm. Here, we performed an in-depth molecular characterization of NFYC Antisense 1 (NFYC-AS1), the asRNA transcribed head-to-head to NFYC subunit of the proliferation-associated NF-Y transcription factor. Our results show that NFYC-AS1 is a prevalently nuclear asRNA peaking early in the cell cycle. Comparative genomics suggests a narrow phylogenetic distribution, with a probable origin in the common ancestor of mammalian lineages. NFYC-AS1 is overexpressed pancancer, preferentially in association with RB1 mutations. Knockdown of NFYC-AS1 by antisense oligonucleotides impairs cell growth in lung squamous cell carcinoma and small cell lung cancer cells, a phenotype recapitulated by CRISPR/Cas9-deletion of its transcription start site. Surprisingly, expression of the sense gene is affected only when endogenous transcription of NFYC-AS1 is manipulated. This suggests that regulation of cell proliferation is at least in part independent of the in cis transcription-mediated effect on NFYC and is possibly exerted by RNA-dependent in trans effects converging on the regulation of G2/M cell cycle phase genes. Accordingly, NFYC-AS1-depleted cells are stuck in mitosis, indicating defects in mitotic progression. Overall, NFYC-AS1 emerged as a cell cycle-regulating asRNA with dual action, holding therapeutic potential in different cancer types, including the very aggressive RB1-mutated tumors.

Effect of diet and genotype on the miRNome of mice with altered lipoprotein metabolism.

The molecular mechanism by which lipid/lipoprotein biosynthesis is regulated in mammals involves a very large number of genes that are subject to multiple levels of regulation. miRNAs are recognized contributors to lipid homeostasis at the post-transcriptional level, although the elucidation of their role is made difficult by the multiplicity of their targets and the ability of more miRNAs to affect the same mRNAs. In this study, an evaluation of how miRNA expression varies in organs playing a key role in lipid/lipoprotein metabolism was conducted in control mice and in two mouse models carrying genetic ablations which differently affect low-density lipoprotein metabolism. Mice were fed a lipid-poor standard diet and a diet enriched in cholesterol and saturated fat. The results obtained showed that there are no miRNAs whose expression constantly vary with dietary or genetic changes. Furthermore, it appears that diet, more than genotype, impacts on organ-specific miRNA expression profiles.

Whole-Exome and Transcriptome Sequencing Expands the Genotype of Majewski Osteodysplastic Primordial Dwarfism Type II.

Microcephalic Osteodysplastic Primordial Dwarfism type II (MOPDII) represents the most common form of primordial dwarfism. MOPD clinical features include severe prenatal and postnatal growth retardation, postnatal severe microcephaly, hypotonia, and an increased risk for cerebrovascular disease and insulin resistance. Autosomal recessive biallelic loss-of-function genomic variants in the centrosomal pericentrin (PCNT) gene on chromosome 21q22 cause MOPDII. Over the past decade, exome sequencing (ES) and massive RNA sequencing have been effectively employed for both the discovery of novel disease genes and to expand the genotypes of well-known diseases. In this paper we report the results both the RNA sequencing and ES of three patients affected by MOPDII with the aim of exploring whether differentially expressed genes and previously uncharacterized gene variants, in addition to PCNT pathogenic variants, could be associated with the complex phenotype of this disease. We discovered a downregulation of key factors involved in growth, such as IGF1R, IGF2R, and RAF1, in all three investigated patients. Moreover, ES identified a shortlist of genes associated with deleterious, rare variants in MOPDII patients. Our results suggest that Next Generation Sequencing (NGS) technologies can be successfully applied for the molecular characterization of the complex genotypic background of MOPDII.

VariantHunter: a method and tool for fast detection of emerging SARS-CoV-2 variants.

With the progression of the COVID-19 pandemic, large datasets of SARS-CoV-2 genome sequences were collected to closely monitor the evolution of the virus and identify the novel variants/strains. By analyzing genome sequencing data, health authorities can 'hunt' novel emerging variants of SARS-CoV-2 as early as possible, and then monitor their evolution and spread. We designed VariantHunter, a highly flexible and user-friendly tool for systematically monitoring the evolution of SARS-CoV-2 at global and regional levels. In VariantHunter, amino acid changes are analyzed over an interval of 4 weeks in an arbitrary geographical area (continent, country, or region); for every week in the interval, the prevalence is computed and changes are ranked based on their increase or decrease in prevalence. VariantHunter supports two main types of analysis: lineage-independent and lineage-specific. The former considers all the available data and aims to discover new viral variants. The latter evaluates specific lineages/viral variants to identify novel candidate designations (sub-lineages and sub-variants). Both analyses use simple statistics and visual representations (diffusion charts and heatmaps) to track viral evolution. A dataset explorer allows users to visualize available data and refine their selection. VariantHunter is a web application free to all users. The two types of supported analysis (lineage-independent and lineage-specific) allow user-friendly monitoring of the viral evolution, empowering genomic surveillance without requiring any computational background. Database URL http://gmql.eu/variant_hunter/.

An Intestinal Th17 Subset is Associated with Inflammation in Crohn's Disease and Activated by Adherent-invasive Escherichia coli.

IFNγ-producing ex-Th17 cells ['Th1/17'] were shown to play a key pathogenic role in experimental colitis and are abundant in the intestine. Here, we identified and characterised a novel, potentially colitogenic subset of Th17 cells in the intestine of patients with Crohn's disease [CD]. Human Th17 cells expressing CCR5 ['pTh17'] co-expressed T-bet and RORC/γt and produced very high levels of IL-17, together with IFN-γ. They had a gene signature of Th17 effector cells and were distinct from established Th1/17 cells. pTh17 cells, but not Th1/17 cells, were associated with intestinal inflammation in CD, and decreased upon successful anti-TNF therapy with infliximab. Conventional CCR5[-]Th17 cells differentiated to pTh17 cells with IL-23 in vitro. Moreover, anti-IL-23 therapy with risankizumab strongly reduced pTh17 cells in the intestine. Importantly, intestinal pTh17 cells were selectively activated by adherent-invasive Escherichia coli [AIEC], but not by a commensal/probiotic E. coli strain. AIEC induced high levels of IL-23 and RANTES from dendritic cells [DC]. Intestinal CCR5+Th1/17 cells responded instead to cytomegalovirus and were reduced in ulcerative colitis [UC], suggesting an unexpected protective role. In conclusion, we identified an IL-23-inducible subset of human intestinal Th17 cells. pTh17 cells produced high levels of pro-inflammatory cytokines, were selectively associated with intestinal inflammation in CD, and responded to CD-associated AIEC, suggesting a key colitogenic role.

Telomere-to-telomere assemblies of 142 strains characterize the genome structural landscape in Saccharomyces cerevisiae.

Pangenomes provide access to an accurate representation of the genetic diversity of species, both in terms of sequence polymorphisms and structural variants (SVs). Here we generated the Saccharomyces cerevisiae Reference Assembly Panel (ScRAP) comprising reference-quality genomes for 142 strains representing the species' phylogenetic and ecological diversity. The ScRAP includes phased haplotype assemblies for several heterozygous diploid and polyploid isolates. We identified circa (ca.) 4,800 nonredundant SVs that provide a broad view of the genomic diversity, including the dynamics of telomere length and transposable elements. We uncovered frequent cases of complex aneuploidies where large chromosomes underwent large deletions and translocations. We found that SVs can impact gene expression near the breakpoints and substantially contribute to gene repertoire evolution. We also discovered that horizontally acquired regions insert at chromosome ends and can generate new telomeres. Overall, the ScRAP demonstrates the benefit of a pangenome in understanding genome evolution at population scale.

Single nucleotide polymorphisms associated with wine fermentation and adaptation to nitrogen limitation in wild and domesticated yeast strains.

For more than 20 years, Saccharomyces cerevisiae has served as a model organism for genetic studies and molecular biology, as well as a platform for biotechnology (e.g., wine production). One of the important ecological niches of this yeast that has been extensively studied is wine fermentation, a complex microbiological process in which S. cerevisiae faces various stresses such as limited availability of nitrogen. Nitrogen deficiencies in grape juice impair fermentation rate and yeast biomass production, leading to sluggish or stuck fermentations, resulting in considerable economic losses for the wine industry. In the present work, we took advantage of the "1002 Yeast Genomes Project" population, the most complete catalogue of the genetic variation in the species and a powerful resource for genotype-phenotype correlations, to study the adaptation to nitrogen limitation in wild and domesticated yeast strains in the context of wine fermentation. We found that wild and domesticated yeast strains have different adaptations to nitrogen limitation, corroborating their different evolutionary trajectories. Using a combination of state-of-the-art bioinformatic (GWAS) and molecular biology (CRISPR-Cas9) methodologies, we validated that PNP1, RRT5 and PDR12 are implicated in wine fermentation, where RRT5 and PDR12 are also involved in yeast adaptation to nitrogen limitation. In addition, we validated SNPs in these genes leading to differences in fermentative capacities and adaptation to nitrogen limitation. Altogether, the mapped genetic variants have potential applications for the genetic improvement of industrial yeast strains.

Rice florigens control a common set of genes at the shoot apical meristem including the F-BOX BROADER TILLER ANGLE 1 that regulates tiller angle and spikelet development.

Rice flowering is triggered by transcriptional reprogramming at the shoot apical meristem (SAM) mediated by florigenic proteins produced in leaves in response to changes in photoperiod. Florigens are more rapidly expressed under short days (SDs) compared to long days (LDs) and include the HEADING DATE 3a (Hd3a) and RICE FLOWERING LOCUS T1 (RFT1) phosphatidylethanolamine binding proteins. Hd3a and RFT1 are largely redundant at converting the SAM into an inflorescence, but whether they activate the same target genes and convey all photoperiodic information that modifies gene expression at the SAM is currently unclear. We uncoupled the contribution of Hd3a and RFT1 to transcriptome reprogramming at the SAM by RNA sequencing of dexamethasone-inducible over-expressors of single florigens and wild-type plants exposed to photoperiodic induction. Fifteen highly differentially expressed genes common to Hd3a, RFT1, and SDs were retrieved, 10 of which still uncharacterized. Detailed functional studies on some candidates revealed a role for LOC_Os04g13150 in determining tiller angle and spikelet development and the gene was renamed BROADER TILLER ANGLE 1 (BRT1). We identified a core set of genes controlled by florigen-mediated photoperiodic induction and defined the function of a novel florigen target controlling tiller angle and spikelet development.

Chitinolytic Enzymes of the Hyperparasite Fungus Aphanocladium album: Genome-Wide Survey and Characterization of A Selected Enzyme.

The filamentous fungus Aphanocladium album is known as a hyperparasite of plant pathogenic fungi; hence, it has been studied as a possible agent for plant protection. Chitinases secreted by A. album have proven to be essential for its fungicidal activity. However, no complete analysis of the A. album chitinase assortment has been carried out, nor have any of its chitinases been characterized yet. In this study, we report the first draft assembly of the genome sequence of A. album (strain MX-95). The in silico functional annotation of the genome allowed the identification of 46 genes encoding chitinolytic enzymes of the GH18 (26 genes), GH20 (8 genes), GH75 (8 genes), and GH3 (4 genes) families. The encoded proteins were investigated by comparative and phylogenetic analysis, allowing clustering in different subgroups. A. album chitinases were also characterized according to the presence of different functional protein domains (carbohydrate-binding modules and catalytic domains) providing the first complete description of the chitinase repertoire of A. album. A single chitinase gene was then selected for complete functional characterization. The encoded protein was expressed in the yeast Pichia pastoris, and its activity was assayed under different conditions of temperature and pH and with different substrates. It was found that the enzyme acts mainly as a chitobiosidase, with higher activity in the 37-50 °C range.

HaploCoV: unsupervised classification and rapid detection of novel emerging variants of SARS-CoV-2.

Accurate and timely monitoring of the evolution of SARS-CoV-2 is crucial for identifying and tracking potentially more transmissible/virulent viral variants, and implement mitigation strategies to limit their spread. Here we introduce HaploCoV, a novel software framework that enables the exploration of SARS-CoV-2 genomic diversity through space and time, to identify novel emerging viral variants and prioritize variants of potential epidemiological interest in a rapid and unsupervised manner. HaploCoV can integrate with any classification/nomenclature and incorporates an effective scoring system for the prioritization of SARS-CoV-2 variants. By performing retrospective analyses of more than 11.5 M genome sequences we show that HaploCoV demonstrates high levels of accuracy and reproducibility and identifies the large majority of epidemiologically relevant viral variants - as flagged by international health authorities - automatically and with rapid turn-around times.Our results highlight the importance of the application of strategies based on the systematic analysis and integration of regional data for rapid identification of novel, emerging variants of SARS-CoV-2. We believe that the approach outlined in this study will contribute to relevant advances to current and future genomic surveillance methods.

Disadvantaged Economic Conditions and Stricter Border Rules Shape Afghan Refugees' Ethnobotany: Insights from Kohat District, NW Pakistan.

The study of migrants' ethnobotany can help to address the diverse socio-ecological factors affecting temporal and spatial changes in local ecological knowledge (LEK). Through semi-structured and in-depth conversations with ninety interviewees among local Pathans and Afghan refugees in Kohat District, NW Pakistan, one hundred and forty-five wild plant and mushroom folk taxa were recorded. The plants quoted by Afghan refugees living inside and outside the camps tend to converge, while the Afghan data showed significant differences with those collected by local Pakistani Pathans. Interviewees mentioned two main driving factors potentially eroding folk plant knowledge: (a) recent stricter border policies have made it more difficult for refugees to visit their home regions in Afghanistan and therefore to also procure plants in their native country; (b) their disadvantaged economic conditions have forced them to engage more and more in urban activities in the host country, leaving little time for farming and foraging practices. Stakeholders should foster the exposure that refugee communities have to their plant resources, try to increase their socio-economic status, and facilitate both their settling outside the camps and their transnational movement for enhancing their use of wild plants, ultimately leading to improvements in their food security and health status.

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/.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones.

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

GIGANTEA Is a Negative Regulator of Abscisic Acid Transcriptional Responses and Sensitivity in Arabidopsis.

Transcriptional reprogramming plays a key role in drought stress responses, preceding the onset of morphological and physiological acclimation. The best-characterized signal regulating gene expression in response to drought is the phytohormone abscisic acid (ABA). ABA-regulated gene expression, biosynthesis and signaling are highly organized in a diurnal cycle, so that ABA-regulated physiological traits occur at the appropriate time of day. The mechanisms that underpin such diel oscillations in ABA signals are poorly characterized. Here we uncover GIGANTEA (GI) as a key gatekeeper of ABA-regulated transcriptional and physiological responses. Time-resolved gene expression profiling by RNA sequencing under different irrigation scenarios indicates that gi mutants produce an exaggerated ABA response, despite accumulating wild-type levels of ABA. Comparisons with ABA-deficient mutants confirm the role of GI in controlling ABA-regulated genes, and the analysis of leaf temperature, a read-out for transpiration, supports a role for GI in the control of ABA-regulated physiological processes. Promoter regions of GI/ABA-regulated transcripts are directly targeted by different classes of transcription factors (TFs), especially PHYTOCHROME-INTERACTING FACTOR and -BINDING FACTOR, together with GI itself. We propose a model whereby diel changes in GI control oscillations in ABA responses. Peak GI accumulation at midday contributes to establishing a phase of reduced ABA sensitivity and related physiological responses, by gating DNA binding or function of different classes of TFs that cooperate or compete with GI at target regions.

miRNAs Copy Number Variations Repertoire as Hallmark Indicator of Cancer Species Predisposition.

Aging is one of the hallmarks of multiple human diseases, including cancer. We hypothesized that variations in the number of copies (CNVs) of specific genes may protect some long-living organisms theoretically more susceptible to tumorigenesis from the onset of cancer. Based on the statistical comparison of gene copy numbers within the genomes of both cancer-prone and -resistant species, we identified novel gene targets linked to tumor predisposition, such as CD52, SAT1 and SUMO. Moreover, considering their genome-wide copy number landscape, we discovered that microRNAs (miRNAs) are among the most significant gene families enriched for cancer progression and predisposition. Through bioinformatics analyses, we identified several alterations in miRNAs copy number patterns, involving miR-221, miR-222, miR-21, miR-372, miR-30b, miR-30d and miR-31, among others. Therefore, our analyses provide the first evidence that an altered miRNAs copy number signature can statistically discriminate species more susceptible to cancer from those that are tumor resistant, paving the way for further investigations.

SARS-CoV-2 variants: what have we learnt so far? Commentary.

Besides the timely detection of different SARS-CoV-2 variants through surveillance systems, functional and modelling studies are essential to better inform public health response and preparedness. Here, the knowledge available so far on SARS-CoV-2 variants is discussed from different perspectives, in order to highlight the relevance of a multidisciplinary approach in countering the threat posed by this insidious virus.