(Re)-definition of the holo- and apo-Fur direct regulons of Helicobacter pylori.

Iron homeostasis is a critical process for living organisms because this metal is an essential co-factor for fundamental biochemical activities, like energy production and detoxification, albeit its excess quickly leads to cell intoxication. The protein Fur (ferric uptake regulator) controls iron homeostasis in bacteria by switching from its apo- to holo-form as a function of the cytoplasmic level of ferrous ions, thereby modulating gene expression. The Helicobacter pylori HpFur protein has the rare ability to operate as a transcriptional commutator; apo- and holo-HpFur function as two different repressors with distinct DNA binding recognition properties for specific sets of target genes. Although the regulation of apo- and holo-HpFur in this bacterium has been extensively investigated, we propose a genome-wide redefinition of holo-HpFur direct regulon in H. pylori by integration of RNA-seq and ChIP-seq data, and a large extension of the apo-HpFur direct regulon. We show that in response to iron availability, new coding sequences, non-coding RNAs, toxin-antitoxin systems, and transcripts within open reading frames are directly regulated by apo- or holo-HpFur. These new targets and the more thorough validation and deeper characterization of those already known provide a complete and updated picture of the direct regulons of this two-faced transcriptional regulator.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility.

BACKGROUND: Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the "missing heritability" phenomenon. OBJECTIVE: To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. METHODS: We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. RESULTS: Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. CONCLUSIONS: No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

COVID-19 mortality in Italy varies by patient age, sex and pandemic wave.

SARS-CoV-2 has caused a worldwide epidemic of enormous proportions, which resulted in different mortality rates in different countries for unknown reasons. We analyzed factors associated with mortality using data from the Italian national database of more than 4 million SARS-CoV-2-positive cases diagnosed between January 2020 and July 2021, including > 415 thousand hospitalized for coronavirus disease-19 (COVID-19) and > 127 thousand deceased. For patients for whom age, sex and date of infection detection were available, we determined the impact of these variables on mortality 30 days after the date of diagnosis or hospitalization. Multivariable weighted Cox analysis showed that each of the analyzed variables independently affected COVID-19 mortality. Specifically, in the overall series, age was the main risk factor for mortality, with HR > 100 in the age groups older than 65 years compared with a reference group of 15-44 years. Male sex presented a two-fold higher risk of death than female sex. Patients infected after the first pandemic wave (i.e. after 30 June 2020) had an approximately threefold lower risk of death than those infected during the first wave. Thus, in a series of all confirmed SARS-CoV-2-infected cases in an entire European nation, elderly age was by far the most significant risk factor for COVID-19 mortality, confirming that protecting the elderly should be a priority in pandemic management. Male sex and being infected during the first wave were additional risk factors associated with COVID-19 mortality.

Interplay between non-coding RNA transcription, stringent phenotype and antibiotic production in Streptomyces.

While in recent years the key role of non-coding RNAs (ncRNAs) in regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces, and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: i.) the wild type strain; ii.) an isogenic pirA-defective mutant with central carbon metabolism imbalance, "relaxed" phenotype, and repression of antibiotic production; iii.) a pirA-derivative strain harboring a "stringent" RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Genomic and functional evaluation of TNFSF14 in multiple sclerosis susceptibility.

Among multiple sclerosis (MS) susceptibility genes, the strongest non-human leukocyte antigen (HLA) signal in the Italian population maps to the TNFSF14 gene encoding LIGHT, a glycoprotein involved in dendritic cell (DC) maturation. Through fine-mapping in a large Italian dataset (4,198 patients with MS and 3,903 controls), we show that the TNFSF14 intronic SNP rs1077667 is the primarily MS-associated variant in the region. Expression quantitative trait locus (eQTL) analysis indicates that the MS risk allele is significantly associated with reduced TNFSF14 messenger RNA levels in blood cells, which is consistent with the allelic imbalance in RNA-Seq reads (P < 0.0001). The MS risk allele is associated with reduced levels of TNFSF14 gene expression (P < 0.01) in blood cells from 84 Italian patients with MS and 80 healthy controls (HCs). Interestingly, patients with MS are lower expressors of TNFSF14 compared to HC (P < 0.007). Individuals homozygous for the MS risk allele display an increased percentage of LIGHT-positive peripheral blood myeloid DCs (CD11c+, P = 0.035) in 37 HCs, as well as in in vitro monocyte-derived DCs from 22 HCs (P = 0.04). Our findings suggest that the intronic variant rs1077667 alters the expression of TNFSF14 in immune cells, which may play a role in MS pathogenesis.

Interplay between Non-Coding RNA Transcription, Stringent/Relaxed Phenotype and Antibiotic Production in Streptomyces ambofaciens.

While in recent years the key role of non-coding RNAs (ncRNAs) in the regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: (i) the wild-type strain; (ii) an isogenic pirA-defective mutant with central carbon metabolism imbalance, "relaxed" phenotype, and repression of antibiotic production; and (iii) a pirA-derivative strain harboring a "stringent" RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that the expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Contribution of Rare and Low-Frequency Variants to Multiple Sclerosis Susceptibility in the Italian Continental Population.

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

Author Correction: Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A novel homozygous mutation in the TRDN gene causes a severe form of pediatric malignant ventricular arrhythmia.

BACKGROUND: Triadin is a protein expressed in cardiac and skeletal muscle that has an essential role in the structure and functional regulation of calcium release units and excitation-contraction coupling. Mutations in the triadin gene (TRDN) have been described in different forms of human arrhythmia syndromes with early onset and severe arrhythmogenic phenotype, including triadin knockout syndrome. OBJECTIVE: The purpose of this study was to characterize the pathogenetic mechanism underlying a case of severe pediatric malignant arrhythmia associated with a defect in the TRDN gene. METHODS: We used a trio whole exome sequencing approach to identify the genetic defect in a 2-year-old boy who had been resuscitated from sudden cardiac arrest and had frequent episodes of ventricular fibrillation and a family history positive for sudden death. We then performed in vitro functional analysis to investigate possible pathogenic mechanisms underlying this severe phenotype. RESULTS: We identified a novel homozygous missense variant (p.L56P) in the TRDN gene in the proband that was inherited from the heterozygous unaffected parents. Expression of a green fluorescent protein (GFP)-tagged mutant human cardiac triadin isoform (TRISK32-L56P-GFP) in heterologous systems revealed that the mutation alters protein dynamics. Furthermore, when co-expressed with the type 2 ryanodine receptor, caffeine-induced calcium release from TRISK32-L56P-GFP was relatively lower compared to that observed with the wild-type construct. CONCLUSION: The results of this study allowed us to hypothesize a pathogenic mechanism underlying this rare arrhythmogenic recessive form, suggesting that the mutant protein potentially can trigger arrhythmias by altering calcium homeostasis.

Candidate Genes and MiRNAs Linked to the Inverse Relationship Between Cancer and Alzheimer's Disease: Insights From Data Mining and Enrichment Analysis.

The incidence of cancer and Alzheimer's disease (AD) increases exponentially with age. A growing body of epidemiological evidence and molecular investigations inspired the hypothesis of an inverse relationship between these two pathologies. It has been proposed that the two diseases might utilize the same proteins and pathways that are, however, modulated differently and sometimes in opposite directions. Investigation of the common processes underlying these diseases may enhance the understanding of their pathogenesis and may also guide novel therapeutic strategies. Starting from a text-mining approach, our in silico study integrated the dispersed biological evidence by combining data mining, gene set enrichment, and protein-protein interaction (PPI) analyses while searching for common biological hallmarks linked to AD and cancer. We retrieved 138 genes (ALZCAN gene set), computed a significant number of enriched gene ontology clusters, and identified four PPI modules. The investigation confirmed the relevance of autophagy, ubiquitin proteasome system, and cell death as common biological hallmarks shared by cancer and AD. Then, from a closer investigation of the PPI modules and of the miRNAs enrichment data, several genes (SQSTM1, UCHL1, STUB1, BECN1, CDKN2A, TP53, EGFR, GSK3B, and HSPA9) and miRNAs (miR-146a-5p, MiR-34a-5p, miR-21-5p, miR-9-5p, and miR-16-5p) emerged as promising candidates. The integrative approach uncovered novel miRNA-gene networks (e.g., miR-146 and miR-34 regulating p62 and Beclin1 in autophagy) that might give new insights into the complex regulatory mechanisms of gene expression in AD and cancer.

Unraveling gut microbiota in Parkinson's disease and atypical parkinsonism.

BACKGROUND: Although several studies have suggested that abnormalities in gut microbiota may play a critical role in the pathogenesis of PD, data are still extremely heterogeneous. METHODS: 16S gene ribosomal RNA sequencing was performed on fecal samples of 350 individuals, subdivided into idiopathic PD (n = 193, of whom 39 were drug naïve) stratified by disease duration, PSP (n = 22), MSA (n = 22), and healthy controls (HC; n = 113). Several confounders were taken into account, including dietary habits. RESULTS: Despite the fact that unadjusted comparison of PD and HC showed several differences in relative taxa abundances, the significant results were greatly reduced after adjusting for confounders. Although most of these differences were associated with disease duration, lower abundance in Lachnospiraceae was the only difference between de novo PD and HC (remaining lower across almost all PD duration strata). Decreased Lachnospiraceae and increased Lactobacillaceae and Christensenellaceae were associated with a worse clinical profile, including higher frequencies of cognitive impairment, gait disturbances, and postural instability. When compared with HC, MSA and PSP patients shared the changes in PD, with a few exceptions: in MSA, Lachnospiraceae were not lower, and Prevotellaceae were reduced; in PSP, Lactobacillaceae were similar, and Streptococcaceae were reduced. CONCLUSIONS: Gut microbiota may be an environmental modulator of the pathogenesis of PD and contribute to the interindividual variability of clinical features. Data are influenced by PD duration and several confounders that need to be taken into account in future studies. Prospective studies in de novo PD patients are needed to elucidate the net effect of dysbiosis on the progression of the disease. © 2018 International Parkinson and Movement Disorder Society.

Correction: MYC-containing amplicons in acute myeloid leukemia: genomic structures, evolution, and transcriptional consequences.

In the original version of this Article, the affiliation details for Giovanni Martinelli were incorrectly given as 'Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, 40138, Bologna, Italy' and it should have been given as 'Istituto Scientifico Romagnolo per lo Studio e la Cura dei Tumori (IRST) IRCCS, Meldola, Italy and not Department of Experimental, Diagnostic and Specialty Medicine, University of Bologna, 40138, Bologna, Italy.'Furthermore, the original version of this Article contained an error in the spelling of the authors Alberto L'Abbate and Pietro D'Addabbo, an acute accent was used instead of an apostrophe for these authors names.These errors have now been corrected in both the PDF and HTML versions of the Article.

Pirin: A novel redox-sensitive modulator of primary and secondary metabolism in Streptomyces.

Pirins are evolutionarily conserved iron-containing proteins that are found in all kingdoms of life, and have been implicated in diverse molecular processes, mostly associated with cellular stress. In the present study, we started from the evidence that the insertional inactivation of pirin-like gene SAM23877_RS18305 (pirA) by ΦC31 Att/Int system-based vectors in spiramycin-producing strain Streptomyces ambofaciens ATCC 23877 resulted in marked effects on central carbon and energy metabolism gene expression, high sensitivity to oxidative injury and repression of polyketide antibiotic production. By using integrated transcriptomic, proteomic and metabolite profiling, together with genetic complementation, we here show that most of these effects could be traced to the inability of the pirA-defective strain to modulate beta-oxidation pathway, leading to an unbalanced supply of precursor monomers for polyketide biosynthesis. Indeed, in silico protein-protein interaction modeling and in vitro experimental validation allowed us to demonstrate that PirA is a novel redox-sensitive negative modulator of very long-chain acyl-CoA dehydrogenase, which catalyzes the first committed step of the beta-oxidation pathway.

A Census of Tandemly Repeated Polymorphic Loci in Genic Regions Through the Comparative Integration of Human Genome Assemblies.

Polymorphic Tandem Repeat (PTR) is a common form of polymorphism in the human genome. A PTR consists in a variation found in an individual (or in a population) of the number of repeating units of a Tandem Repeat (TR) locus of the genome with respect to the reference genome. Several phenotypic traits and diseases have been discovered to be strongly associated with or caused by specific PTR loci. PTR are further distinguished in two main classes: Short Tandem Repeats (STR) when the repeating unit has size up to 6 base pairs, and Variable Number Tandem Repeats (VNTR) for repeating units of size above 6 base pairs. As larger and larger populations are screened via high throughput sequencing projects, it becomes technically feasible and desirable to explore the association between PTR and a panoply of such traits and conditions. In order to facilitate these studies, we have devised a method for compiling catalogs of PTR from assembled genomes, and we have produced a catalog of PTR for genic regions (exons, introns, UTR and adjacent regions) of the human genome (GRCh38). We applied four different TR discovery software tools to uncover in the first phase 55,223,485 TR (after duplicate removal) in GRCh38, of which 373,173 were determined to be PTR in the second phase by comparison with five assembled human genomes. Of these, 263,266 are not included by state-of-the-art PTR catalogs. The new methodology is mainly based on a hierarchical and systematic application of alignment-based sequence comparisons to identify and measure the polymorphism of TR. While previous catalogs focus on the class of STR of small total size, we remove any size restrictions, aiming at the more general class of PTR, and we also target fuzzy TR by using specific detection tools. Similarly to other previous catalogs of human polymorphic loci, we focus our catalog toward applications in the discovery of disease-associated loci. Validation by cross-referencing with existing catalogs on common clinically-relevant loci shows good concordance. Overall, this proposed census of human PTR in genic regions is a shared resource (web accessible), complementary to existing catalogs, facilitating future genome-wide studies involving PTR.

MYC-containing amplicons in acute myeloid leukemia: genomic structures, evolution, and transcriptional consequences.

Double minutes (dmin), homogeneously staining regions, and ring chromosomes are vehicles of gene amplification in cancer. The underlying mechanism leading to their formation as well as their structure and function in acute myeloid leukemia (AML) remain mysterious. We combined a range of high-resolution genomic methods to investigate the architecture and expression pattern of amplicons involving chromosome band 8q24 in 23 cases of AML (AML-amp). This revealed that different MYC-dmin architectures can coexist within the same leukemic cell population, indicating a step-wise evolution rather than a single event origin, such as through chromothripsis. This was supported also by the analysis of the chromothripsis criteria, that poorly matched the model in our samples. Furthermore, we found that dmin could evolve toward ring chromosomes stabilized by neocentromeres. Surprisingly, amplified genes (mainly PVT1) frequently participated in fusion transcripts lacking a corresponding DNA template. We also detected a significant overexpression of the circular RNA of PVT1 (circPVT1) in AML-amp cases versus AML with a normal karyotype. Our results show that 8q24 amplicons in AML are surprisingly plastic DNA structures with an unexpected association to novel fusion transcripts and circular RNAs.

The Hidden Genomic and Transcriptomic Plasticity of Giant Marker Chromosomes in Cancer.

Genome amplification in the form of rings or giant rod-shaped marker chromosomes (RGMs) is a common genetic alteration in soft tissue tumors. The mitotic stability of these structures is often rescued by perfectly functioning analphoid neocentromeres, which therefore significantly contribute to cancer progression. Here, we disentangled the genomic architecture of many neocentromeres stabilizing marker chromosomes in well-differentiated liposarcoma and lung sarcomatoid carcinoma samples. In cells carrying heavily rearranged RGMs, these structures were assembled as patchworks of multiple short amplified sequences, disclosing an extremely high level of complexity and definitely ruling out the existence of regions prone to neocentromere seeding. Moreover, by studying two well-differentiated liposarcoma samples derived from the onset and the recurrence of the same tumor, we documented an expansion of the neocentromeric domain that occurred during tumor progression, which reflects a strong selective pressure acting toward the improvement of the neocentromeric functionality in cancer. In lung sarcomatoid carcinoma cells we documented, extensive "centromere sliding" phenomena giving rise to multiple, closely mapping neocentromeric epialleles on separate coexisting markers occur, likely due to the instability of neocentromeres arising in cancer cells. Finally, by investigating the transcriptional activity of neocentromeres, we came across a burst of chimeric transcripts, both by extremely complex genomic rearrangements, and cis/trans-splicing events. Post-transcriptional editing events have been reported to expand and variegate the genetic repertoire of higher eukaryotes, so they might have a determining role in cancer. The increased incidence of fusion transcripts, might act as a driving force for the genomic amplification process, together with the increased transcription of oncogenes.

Retroviral Scanning: Mapping MLV Integration Sites to Define Cell-specific Regulatory Regions.

Moloney murine leukemia (MLV) virus-based retroviral vectors integrate predominantly in acetylated enhancers and promoters. For this reason, mLV integration sites can be used as functional markers of active regulatory elements. Here, we present a retroviral scanning tool, which allows the genome-wide identification of cell-specific enhancers and promoters. Briefly, the target cell population is transduced with an mLV-derived vector and genomic DNA is digested with a frequently cutting restriction enzyme. After ligation of genomic fragments with a compatible DNA linker, linker-mediated polymerase chain reaction (LM-PCR) allows the amplification of the virus-host genome junctions. Massive sequencing of the amplicons is used to define the mLV integration profile genome-wide. Finally, clusters of recurrent integrations are defined to identify cell-specific regulatory regions, responsible for the activation of cell-type specific transcriptional programs. The retroviral scanning tool allows the genome-wide identification of cell-specific promoters and enhancers in prospectively isolated target cell populations. Notably, retroviral scanning represents an instrumental technique for the retrospective identification of rare populations (e.g. somatic stem cells) that lack robust markers for prospective isolation.

WoPPER: Web server for Position Related data analysis of gene Expression in Prokaryotes.

The structural and conformational organization of chromosomes is crucial for gene expression regulation in eukaryotes and prokaryotes as well. Up to date, gene expression data generated using either microarray or RNA-sequencing are available for many bacterial genomes. However, differential gene expression is usually investigated with methods considering each gene independently, thus not taking into account the physical localization of genes along a bacterial chromosome. Here, we present WoPPER, a web tool integrating gene expression and genomic annotations to identify differentially expressed chromosomal regions in bacteria. RNA-sequencing or microarray-based gene expression data are provided as input, along with gene annotations. The user can select genomic annotations from an internal database including 2780 bacterial strains, or provide custom genomic annotations. The analysis produces as output the lists of positionally related genes showing a coordinated trend of differential expression. Graphical representations, including a circular plot of the analyzed chromosome, allow intuitive browsing of the results. The analysis procedure is based on our previously published R-package PREDA. The release of this tool is timely and relevant for the scientific community, as WoPPER will fill an existing gap in prokaryotic gene expression data analysis and visualization tools. WoPPER is open to all users and can be reached at the following URL: https://WoPPER.ba.itb.cnr.it.

Time-Resolved Transcriptomics and Constraint-Based Modeling Identify System-Level Metabolic Features and Overexpression Targets to Increase Spiramycin Production in Streptomyces ambofaciens.

In this study we have applied an integrated system biology approach to characterize the metabolic landscape of Streptomyces ambofaciens and to identify a list of potential metabolic engineering targets for the overproduction of the secondary metabolites in this microorganism. We focused on an often overlooked growth period (i.e., post-first rapid growth phase) and, by integrating constraint-based metabolic modeling with time resolved RNA-seq data, we depicted the main effects of changes in gene expression on the overall metabolic reprogramming occurring in S. ambofaciens. Moreover, through metabolic modeling, we unraveled a set of candidate overexpression gene targets hypothetically leading to spiramycin overproduction. Model predictions were experimentally validated by genetic manipulation of the recently described ethylmalonyl-CoA metabolic node, providing evidence that spiramycin productivity may be increased by enhancing the carbon flow through this pathway. The goal was achieved by over-expressing the ccr paralog srm4 in an ad hoc engineered plasmid. This work embeds the first metabolic reconstruction of S. ambofaciens and the successful experimental validation of model predictions and demonstrates the validity and the importance of in silico modeling tools for the overproduction of molecules with a biotechnological interest. Finally, the proposed metabolic reconstruction, which includes manually refined pathways for several secondary metabolites with antimicrobial activity, represents a solid platform for the future exploitation of S. ambofaciens biotechnological potential.