Prioritization of New Candidate Genes for Rare Genetic Diseases by a Disease-Aware Evaluation of Heterogeneous Molecular Networks.

Screening for pathogenic variants in the diagnosis of rare genetic diseases can now be performed on all genes thanks to the application of whole exome and genome sequencing (WES, WGS). Yet the repertoire of gene-disease associations is not complete. Several computer-based algorithms and databases integrate distinct gene-gene functional networks to accelerate the discovery of gene-disease associations. We hypothesize that the ability of every type of information to extract relevant insights is disease-dependent. We compiled 33 functional networks classified into 13 knowledge categories (KCs) and observed large variability in their ability to recover genes associated with 91 genetic diseases, as measured using efficiency and exclusivity. We developed GLOWgenes, a network-based algorithm that applies random walk with restart to evaluate KCs' ability to recover genes from a given list associated with a phenotype and modulates the prediction of new candidates accordingly. Comparison with other integration strategies and tools shows that our disease-aware approach can boost the discovery of new gene-disease associations, especially for the less obvious ones. KC contribution also varies if obtained using recently discovered genes. Applied to 15 unsolved WES, GLOWgenes proposed three new genes to be involved in the phenotypes of patients with syndromic inherited retinal dystrophies.

Analysis of Differentially Expressed MicroRNAs in Serum and Lung Tissues from Individuals with Severe Asthma Treated with Oral Glucocorticoids.

Nowadays, microRNAs (miRNAs) are increasingly used as biomarkers due to their potential contribution to the diagnosis and targeted treatment of a range of diseases. The aim of the study was to analyze the miRNA expression profiles in serum and lung tissue from patients with severe asthma treated with oral corticosteroids (OCS) and those without OCS treatment. For this purpose, serum and lung tissue miRNAs of OCS and non-OCS asthmatic individuals were evaluated by miRNAs-Seq, and subsequently miRNA validation was performed using RT-qPCR. Additionally, pathway enrichment analysis of deregulated miRNAs was conducted. We observed altered expression by the next-generation sequencing (NGS) of 11 miRNAs in serum, of which five (hsa-miR-148b-3p, hsa-miR-221-5p, hsa-miR-618, hsa-miR-941, and hsa-miR-769-5p) were validated by RT-qPCR, and three miRNAs in lung tissue (hsa-miR-144-3p, hsa-miR-144-5p, and hsa-miR-451a). The best multivariate logistic regression model to differentiate individuals with severe asthma, treated and untreated with OCS, was to combine the serum miRNAs hsa-miR-221-5p and hsa-miR-769-5p. Expression of hsa-miR-148b-3p and hsa-miR-221-5p correlated with FEV1/FVC (%) and these altered miRNAs act in key signaling pathways for asthma disease and the regulated expression of some genes (FOXO3, PTEN, and MAPK3) involved in these pathways. In conclusion, there are miRNA profiles differentially expressed in OCS-treated individuals with asthma and could be used as biomarkers of OCS treatment.

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies.

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort's IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

SQANTI: extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification.

High-throughput sequencing of full-length transcripts using long reads has paved the way for the discovery of thousands of novel transcripts, even in well-annotated mammalian species. The advances in sequencing technology have created a need for studies and tools that can characterize these novel variants. Here, we present SQANTI, an automated pipeline for the classification of long-read transcripts that can assess the quality of data and the preprocessing pipeline using 47 unique descriptors. We apply SQANTI to a neuronal mouse transcriptome using Pacific Biosciences (PacBio) long reads and illustrate how the tool is effective in characterizing and describing the composition of the full-length transcriptome. We perform extensive evaluation of ToFU PacBio transcripts by PCR to reveal that an important number of the novel transcripts are technical artifacts of the sequencing approach and that SQANTI quality descriptors can be used to engineer a filtering strategy to remove them. Most novel transcripts in this curated transcriptome are novel combinations of existing splice sites, resulting more frequently in novel ORFs than novel UTRs, and are enriched in both general metabolic and neural-specific functions. We show that these new transcripts have a major impact in the correct quantification of transcript levels by state-of-the-art short-read-based quantification algorithms. By comparing our iso-transcriptome with public proteomics databases, we find that alternative isoforms are elusive to proteogenomics detection. SQANTI allows the user to maximize the analytical outcome of long-read technologies by providing the tools to deliver quality-evaluated and curated full-length transcriptomes.

Mip6 binds directly to the Mex67 UBA domain to maintain low levels of Msn2/4 stress-dependent mRNAs.

RNA-binding proteins (RBPs) participate in all steps of gene expression, underscoring their potential as regulators of RNA homeostasis. We structurally and functionally characterize Mip6, a four-RNA recognition motif (RRM)-containing RBP, as a functional and physical interactor of the export factor Mex67. Mip6-RRM4 directly interacts with the ubiquitin-associated (UBA) domain of Mex67 through a loop containing tryptophan 442. Mip6 shuttles between the nucleus and the cytoplasm in a Mex67-dependent manner and concentrates in cytoplasmic foci under stress. Photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation experiments show preferential binding of Mip6 to mRNAs regulated by the stress-response Msn2/4 transcription factors. Consistent with this binding, MIP6 deletion affects their export and expression levels. Additionally, Mip6 interacts physically and/or functionally with proteins with a role in mRNA metabolism and transcription such as Rrp6, Xrn1, Sgf73, and Rpb1. These results reveal a novel role for Mip6 in the homeostasis of Msn2/4-dependent transcripts through its direct interaction with the Mex67 UBA domain.

Targeted Therapy Modulates the Secretome of Cancer-Associated Fibroblasts to Induce Resistance in HER2-Positive Breast Cancer.

The combination of trastuzumab plus pertuzumab plus docetaxel as a first-line therapy in patients with HER2-positive metastatic breast cancer has provided significant clinical benefits compared to trastuzumab plus docetaxel alone. However, despite the therapeutic success of existing therapies targeting HER2, tumours invariably relapse. Therefore, there is an urgent need to improve our understanding of the mechanisms governing resistance, so that specific therapeutic strategies can be developed to provide improved efficacy. It is well known that the tumour microenvironment (TME) has a significant impact on cancer behaviour. Cancer-associated fibroblasts (CAFs) are essential components of the tumour stroma that have been linked to acquired therapeutic resistance and poor prognosis in breast cancer. For this reason, it would be of interest to identify novel biomarkers in the tumour stroma that could emerge as therapeutic targets for the modulation of resistant phenotypes. Conditioned medium experiments carried out in our laboratory with CAFs derived from HER2-positive patients showed a significant capacity to promote resistance to trastuzumab plus pertuzumab therapies in two HER2-positive breast cancer cell lines (BCCLs), even in the presence of docetaxel. In order to elucidate the components of the CAF-conditioned medium that may be relevant in the promotion of BCCL resistance, we implemented a multiomics strategy to identify cytokines, transcription factors, kinases and miRNAs in the secretome that have specific targets in cancer cells. The combination of cytokine arrays, label-free LC-MS/MS quantification and miRNA analysis to explore the secretome of CAFs under treatment conditions revealed several up- and downregulated candidates. We discuss the potential role of some of the most interesting candidates in generating resistance in HER2-positive breast cancer.

Fine Breakpoint Mapping by Genome Sequencing Reveals the First Large X Inversion Disrupting the NHS Gene in a Patient with Syndromic Cataracts.

Inversions are structural variants that are generally balanced. However, they could lead to gene disruptions or have positional effects leading to diseases. Mutations in the NHS gene cause Nance-Horan syndrome, an X-linked disorder characterised by congenital cataracts and dental anomalies. Here, we aimed to characterise a balanced pericentric inversion X(p22q27), maternally inherited, in a child with syndromic bilateral cataracts by breakpoint mapping using whole-genome sequencing (WGS). 30× Illumina paired-end WGS was performed in the proband, and breakpoints were confirmed by Sanger sequencing. EdU assays and FISH analysis were used to assess skewed X-inactivation patterns. RNA expression of involved genes in the breakpoint boundaries was evaluated by droplet-digital PCR. We defined the breakpoint position of the inversion at Xp22.13, with a 15 bp deletion, disrupting the unusually large intron 1 of the canonical NHS isoform, and also perturbing topologically-associated domains (TADs). Moreover, a microhomology region of 5 bp was found on both sides. RNA analysis confirmed null and reduced NHS expression in the proband and his unaffected mother, respectively. In conclusion, we report the first chromosomal inversion disrupting NHS, fine-mapped by WGS. Our data expand the clinical spectrum and the pathogenic mechanisms underlying the NHS defects.

Association between neutrophil-lymphocyte ratio, platelet-lymphocyte ratio, and C-reactive protein levels and metabolic status in patients with a bipolar disorder.

OBJECTIVES: Neutrophil-lymphocyte ratio (NLR), platelet-lymphocyte ratio (PLR), and C-Reactive Protein (CRP) are markers of inflammation that are elevated in bipolar disorder (BD) and are also related to a higher risk of metabolic syndrome (MetS). This study aimed at investigating for the first time the association between NLR, PLR, and CRP and the metabolic status in BD. METHODS: We assessed the association between biomarkers and the metabolic status: number of metabolic risk factors, presence of MetS, insulin sensitivity (Quantitative Insulin Sensitivity Check Index, QUICKI) and insulin resistance (Homeostatic Model Assessment for Insulin Resistance, HOMA-IR index), in a sample of 219 outpatients with BD. RESULTS: 25.9% of the sample met the criteria for MetS. High levels of CRP were found in 12% of the sample. Older age, low PLR, high NLR, and high CRP levels significantly predicted a higher number of MetS risk factors (p < 0.001). Older age and low PLR were associated with a greater likelihood of developing MetS (p = 0.007). CONCLUSIONS: Although further studies are needed to replicate and validate these findings, inflammatory biomarkers as CRP, PLR and NLR could be useful tools to identify patients with a BD at risk for a metabolic adverse outcome.

New roles for AP-1/JUNB in cell cycle control and tumorigenic cell invasion via regulation of cyclin E1 and TGF-ß2.

BACKGROUND: JUNB transcription factor contributes to the formation of the ubiquitous transcriptional complex AP-1 involved in the control of many physiological and disease-associated functions. The roles of JUNB in the control of cell division and tumorigenic processes are acknowledged but still unclear. RESULTS: Here, we report the results of combined transcriptomic, genomic, and functional studies showing that JUNB promotes cell cycle progression via induction of cyclin E1 and repression of transforming growth factor (TGF)-ß2 genes. We also show that high levels of JUNB switch the response of TGF-ß2 stimulation from an antiproliferative to a pro-invasive one, induce endogenous TGF-ß2 production by promoting TGF-ß2 mRNA translation, and enhance tumor growth and metastasis in mice. Moreover, tumor genomic data indicate that JUNB amplification associates with poor prognosis in breast and ovarian cancer patients. CONCLUSIONS: Our results reveal novel functions for JUNB in cell proliferation and tumor aggressiveness through regulation of cyclin E1 and TGF-ß2 expression, which might be exploited for cancer prognosis and therapy.

An evaluation of pipelines for DNA variant detection can guide a reanalysis protocol to increase the diagnostic ratio of genetic diseases.

Clinical exome (CE) sequencing has become a first-tier diagnostic test for hereditary diseases; however, its diagnostic rate is around 30-50%. In this study, we aimed to increase the diagnostic yield of CE using a custom reanalysis algorithm. Sequencing data were available for three cohorts using two commercial protocols applied as part of the diagnostic process. Using these cohorts, we compared the performance of general and clinically relevant variant calling and the efficacy of an in-house bioinformatic protocol (FJD-pipeline) in detecting causal variants as compared to commercial protocols. On the whole, the FJD-pipeline detected 99.74% of the causal variants identified by the commercial protocol in previously solved cases. In the unsolved cases, FJD-pipeline detects more INDELs and non-exonic variants, and is able to increase the diagnostic yield in 2.5% and 3.2% in the re-analysis of 78 cancer and 62 cardiovascular cases. These results were considered to design a reanalysis, filtering and prioritization algorithm that was tested by reassessing 68 inconclusive cases of monoallelic autosomal recessive retinal dystrophies increasing the diagnosis by 4.4%. In conclusion, a guided NGS reanalysis of unsolved cases increases the diagnostic yield in genetic disorders, making it a useful diagnostic tool in medical genetics.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations.

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes.

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Peripheral T-cell lymphoma: molecular profiling recognizes subclasses and identifies prognostic markers.

Peripheral T-cell lymphoma (PTCL) is a clinically aggressive disease, with a poor response to therapy and a low overall survival rate of approximately 30% after 5 years. We have analyzed a series of 105 cases with a diagnosis of PTCL using a customized NanoString platform (NanoString Technologies, Seattle, WA) that includes 208 genes associated with T-cell differentiation, oncogenes and tumor suppressor genes, deregulated pathways, and stromal cell subpopulations. A comparative analysis of the various histological types of PTCL (angioimmunoblastic T-cell lymphoma [AITL]; PTCL with T follicular helper [TFH] phenotype; PTCL not otherwise specified [NOS]) showed that specific sets of genes were associated with each of the diagnoses. These included TFH markers, cytotoxic markers, and genes whose expression was a surrogate for specific cellular subpopulations, including follicular dendritic cells, mast cells, and genes belonging to precise survival (NF-κB) and other pathways. Furthermore, the mutational profile was analyzed using a custom panel that targeted 62 genes in 76 cases distributed in AITL, PTCL-TFH, and PTCL-NOS. The main differences among the 3 nodal PTCL classes involved the RHOAG17V mutations (P < .0001), which were approximately twice as frequent in AITL (34.09%) as in PTCL-TFH (16.66%) cases but were not detected in PTCL-NOS. A multivariate analysis identified gene sets that allowed the series of cases to be stratified into different risk groups. This study supports and validates the current division of PTCL into these 3 categories, identifies sets of markers that can be used for a more precise diagnosis, and recognizes the expression of B-cell genes as an IPI-independent prognostic factor for AITL.

tappAS: a comprehensive computational framework for the analysis of the functional impact of differential splicing.

Recent advances in long-read sequencing solve inaccuracies in alternative transcript identification of full-length transcripts in short-read RNA-Seq data, which encourages the development of methods for isoform-centered functional analysis. Here, we present tappAS, the first framework to enable a comprehensive Functional Iso-Transcriptomics (FIT) analysis, which is effective at revealing the functional impact of context-specific post-transcriptional regulation. tappAS uses isoform-resolved annotation of coding and non-coding functional domains, motifs, and sites, in combination with novel analysis methods to interrogate different aspects of the functional readout of transcript variants and isoform regulation. tappAS software and documentation are available at https://app.tappas.org.

Association between genetic variation in the myo-inositol monophosphatase 2 (IMPA2) gene and age at onset of bipolar disorder.

INTRODUCTION: The age at onset of bipolar disorder (BD) has significant implications for severity, duration of affective episodes, response to treatment, and psychiatric comorbidities. It has been suggested that early-onset BD (EO-BD) could represent a clinically distinct subtype with probable genetic risk factors different from those of late-onset BD (LO-BD). To date, several genes have been associated with BD risk but few studies have investigated the genetic differences between EO-BD and LO-BD. The aim of this study was to evaluate if variants of the gene coding for myo-inositol monophosphatase (IMPA2) are linked to age at onset of BD. METHOD: 235 bipolar patients were recruited and assessed. The final sample consisting of 192 euthymic individuals, was compared according to the age at onset. Polymorphisms were genotyped in the IMPA2 gene (rs669838, rs1020294, rs1250171, and rs630110). Early-onset was defined by the appearance of a first affective episode before the age of 18. RESULTS: The analyses showed that in the genotype distribution rs1020294 (p = .01) and rs1250171 (p = .01) were associated with the age at onset. The significant effect remained only in the rs1020294 SNP in which G carriers were more likely to debut later compared to patients presenting the AA genotype (p = .002; OR = 9.57, CI95%[2.37-38.64]). The results also showed that EO-BD tended to experience more alcohol misuse (p = .003; OR = .197, CI95%[.07-.58]) compared to LO-BD. CONCLUSIONS: Our results provide evidence for genetic differences between EO-BD and LO-BD at the IMPA2 gene as well as clinical differences between subgroups with therapeutic implications.

Notch-regulated miR-223 targets the aryl hydrocarbon receptor pathway and increases cytokine production in macrophages from rheumatoid arthritis patients.

Evidence links aryl hydrocarbon receptor (AHR) activation to rheumatoid arthritis (RA) pathogenesis, although results are inconsistent. AHR agonists inhibit pro-inflammatory cytokine expression in macrophages, pivotal cells in RA aetiopathogenesis, which hints at specific circuits that regulate the AHR pathway in RA macrophages. We compared microRNA (miR) expression in CD14(+) cells from patients with active RA or with osteoarthritis (OA). Seven miR were downregulated and one (miR-223) upregulated in RA compared to OA cells. miR-223 upregulation correlated with reduced Notch3 and Notch effector expression in RA patients. Overexpression of the Notch-induced repressor HEY-1 and co-culture of healthy donor monocytes with Notch ligand-expressing cells showed direct Notch-mediated downregulation of miR-223. Bioinformatics predicted the AHR regulator ARNT (AHR nuclear translocator) as a miR-223 target. Pre-miR-223 overexpression silenced ARNT 3'UTR-driven reporter expression, reduced ARNT (but not AHR) protein levels and prevented AHR/ARNT-mediated inhibition of pro-inflammatory cytokine expression. miR-223 counteracted AHR/ARNT-induced Notch3 upregulation in monocytes. Levels of ARNT and of CYP1B1, an AHR/ARNT signalling effector, were reduced in RA compared to OA synovial tissue, which correlated with miR-223 levels. Our results associate Notch signalling to miR-223 downregulation in RA macrophages, and identify miR-223 as a negative regulator of the AHR/ARNT pathway through ARNT targeting.