Detection of elusive DNA copy-number variations in hereditary disease and cancer through the use of noncoding and off-target sequencing reads.

Copy-number variants (CNVs) play a substantial role in the molecular pathogenesis of hereditary disease and cancer, as well as in normal human interindividual variation. However, they are still rather difficult to identify in mainstream sequencing projects, especially involving exome sequencing, because they often occur in DNA regions that are not targeted for analysis. To overcome this problem, we developed OFF-PEAK, a user-friendly CNV detection tool that builds on a denoising approach and the use of "off-target" DNA reads, which are usually discarded by sequencing pipelines. We benchmarked OFF-PEAK on data from targeted sequencing of 96 cancer samples, as well as 130 exomes of individuals with inherited retinal disease from three different populations. For both sets of data, OFF-PEAK demonstrated excellent performance (>95% sensitivity and >80% specificity vs. experimental validation) in detecting CNVs from in silico data alone, indicating its immediate applicability to molecular diagnosis and genetic research.

Long-read genome sequencing identifies cryptic structural variants in congenital aniridia cases.

BACKGROUND: Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved "PAX6-negative" cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches. RESULTS: Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1. CONCLUSIONS: In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

PRPH2-Related Retinal Dystrophies: Mutational Spectrum in 103 Families from a Spanish Cohort.

PRPH2, one of the most frequently inherited retinal dystrophy (IRD)-causing genes, implies a high phenotypic variability. This study aims to analyze the PRPH2 mutational spectrum in one of the largest cohorts worldwide, and to describe novel pathogenic variants and genotype-phenotype correlations. A study of 220 patients from 103 families recruited from a database of 5000 families. A molecular diagnosis was performed using classical molecular approaches and next-generation sequencing. Common haplotypes were ascertained by analyzing single-nucleotide polymorphisms. We identified 56 variants, including 11 novel variants. Most of them were missense variants (64%) and were located in the D2-loop protein domain (77%). The most frequently occurring variants were p.Gly167Ser, p.Gly208Asp and p.Pro221_Cys222del. Haplotype analysis revealed a shared region in families carrying p.Leu41Pro or p.Pro221_Cys222del. Patients with retinitis pigmentosa presented an earlier disease onset. We describe the largest cohort of IRD families associated with PRPH2 from a single center. Most variants were located in the D2-loop domain, highlighting its importance in interacting with other proteins. Our work suggests a likely founder effect for the variants p.Leu41Pro and p.Pro221_Cys222del in our Spanish cohort. Phenotypes with a primary rod alteration presented more severe affectation. Finally, the high phenotypic variability in PRPH2 hinders the possibility of drawing genotype-phenotype correlations.

Atypical Mutational Spectrum of SARS-CoV-2 Replicating in the Presence of Ribavirin.

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

VIsoQLR: an interactive tool for the detection, quantification and fine-tuning of isoforms in selected genes using long-read sequencing.

DNA variants altering the pre-mRNA splicing process represent an underestimated cause of human genetic diseases. Their association with disease traits should be confirmed using functional assays from patient cell lines or alternative models to detect aberrant mRNAs. Long-read sequencing is a suitable technique to identify and quantify mRNA isoforms. Available isoform detection and/or quantification tools are generally designed for the whole transcriptome analysis. However experiments focusing on genes of interest need more precise data fine-tuning and visualization tools.Here we describe VIsoQLR, an interactive analyzer, viewer and editor for the semi-automated identification and quantification of known and novel isoforms using long-read sequencing data. VIsoQLR is tailored to thoroughly analyze mRNA expression in splicing assays of selected genes. Our tool takes sequences aligned to a reference, and for each gene, it defines consensus splice sites and quantifies isoforms. VIsoQLR introduces features to edit the splice sites through dynamic and interactive graphics and tables, allowing accurate manual curation. Known isoforms detected by other methods can also be imported as references for comparison. A benchmark against two other popular transcriptome-based tools shows VIsoQLR accurate performance on both detection and quantification of isoforms. Here, we present VIsoQLR principles and features and its applicability in a case study example using nanopore-based long-read sequencing. VIsoQLR is available at https://github.com/TBLabFJD/VIsoQLR .

Biallelic intragenic tandem duplication of CPLANE1 in Joubert syndrome: A case report.

Joubert syndrome (JS) is a clinically and genetically heterogeneous genetic disorder. To date, 40 JS-causing genes have been reported and CPLANE1 is one of the most frequently mutated, with biallelic pathogenic missense and truncating variants explaining up to 14% of JS cases. We present a case of JS diagnosed after the identification of a novel biallelic intragenic duplication of exons 20-46 of CPLANE1. The quadruplication was identified by short-read sequencing and copy number variant analysis and confirmed in tandem by long PCR with the breakpoints defined by a nanopore-based long-read sequencing approach. Based on the genetic findings and the clinical presentation of the patient, a brain MRI was ordered, evidencing the molar tooth sign, which confirmed the diagnosis of JS in the patient. This is, to the best of our knowledge, the first report of an intragenic duplication in this gene as the potential molecular mechanism of JS.

CSVS, a crowdsourcing database of the Spanish population genetic variability.

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

Minigene Splicing Assays and Long-Read Sequencing to Unravel Pathogenic Deep-Intronic Variants in PAX6 in Congenital Aniridia.

PAX6 haploinsufficiency causes aniridia, a congenital eye disorder that involves the iris, and foveal hypoplasia. Comprehensive screening of the PAX6 locus, including the non-coding regions, by next-generation sequencing revealed four deep-intronic variants with potential effects on pre-RNA splicing. Nevertheless, without a functional analysis, their pathogenicity could not be established. We aimed to decipher their impact on the canonical PAX6 splicing using in vitro minigene splicing assays and nanopore-based long-read sequencing. Two multi-exonic PAX6 constructs were generated, and minigene assays were carried out. An aberrant splicing pattern was observed for two variants in intron 6, c.357+136G>A and c.357+334G>A. In both cases, several exonization events, such as pseudoexon inclusions and partial intronic retention, were observed due to the creation or activation of new/cryptic non-canonical splicing sites, including a shared intronic donor site. In contrast, two variants identified in intron 11, c.1032+170A>T and c.1033-275A>C, seemed not to affect splicing processes. We confirmed the high complexity of alternative splicing of PAX6 exon 6, which also involves unreported cryptic intronic sites. Our study highlights the importance of integrating functional studies into diagnostic algorithms to decipher the potential implication of non-coding variants, usually classified as variants of unknown significance, thus allowing variant reclassification to achieve a conclusive genetic diagnosis.

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.

Comparison of Extracellular Vesicle Isolation Methods for miRNA Sequencing.

MicroRNAs (miRNAs) encapsulated in extracellular vesicles (EVs) are potential diagnostic and prognostic biomarkers. However, discrepancies in miRNA patterns and their validation are still frequent due to differences in sample origin, EV isolation, and miRNA sequencing methods. The aim of the present study is to find a reliable EV isolation method for miRNA sequencing, adequate for clinical application. To this aim, two comparative studies were performed in parallel with the same human plasma sample: (i) isolation and characterization of EVs obtained using three procedures: size exclusion chromatography (SEC), iodixanol gradient (GRAD), and its combination (SEC+GRAD) and (ii) evaluation of the yield of miRNA sequences obtained using NextSeq 500 (Illumina) and three miRNA library preparation protocols: NEBNext, NEXTFlex, and SMARTer smRNA-seq. The conclusion of comparison (i) is that recovery of the largest amount of EVs and reproducibility were attained with SEC, but GRAD and SEC+GRAD yielded purer EV preparations. The conclusion of (ii) is that the NEBNext library showed the highest reproducibility in the number of miRNAs recovered and the highest diversity of miRNAs. These results render the combination of GRAD EV isolation and NEBNext library preparation for miRNA retrieval as adequate for clinical applications using plasma samples.

IL-6-based mortality prediction model for COVID-19: Validation and update in multicenter and second wave cohorts.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is a highly variable condition. Validated tools to assist in the early detection of patients at high risk of mortality can help guide medical decisions. OBJECTIVE: We sought to validate externally, as well as in patients from the second pandemic wave in Europe, our previously developed mortality prediction model for hospitalized COVID-19 patients. METHODS: Three validation cohorts were generated: 2 external with 185 and 730 patients from the first wave and 1 internal with 119 patients from the second wave. The probability of death was calculated for all subjects using our prediction model, which includes peripheral blood oxygen saturation/fraction of inspired oxygen ratio, neutrophil-to-lymphocyte ratio, lactate dehydrogenase, IL-6, and age. Discrimination and calibration were evaluated in the validation cohorts. The prediction model was updated by reestimating individual risk factor effects in the overall cohort (N = 1477). RESULTS: The mortality prediction model showed good performance in the external validation cohorts 1 and 2, and in the second wave validation cohort 3 (area under the receiver-operating characteristic curve, 0.94, 0.86, and 0.86, respectively), with excellent calibration (calibration slope, 0.86, 0.94, and 0.79; intercept, 0.05, 0.03, and 0.10, respectively). The updated model accurately predicted mortality in the overall cohort (area under the receiver-operating characteristic curve, 0.91), which included patients from both the first and second COVID-19 waves. The updated model was also useful to predict fatal outcome in patients without respiratory distress at the time of evaluation. CONCLUSIONS: This is the first COVID-19 mortality prediction model validated in patients from the first and second pandemic waves. The COR+12 online calculator is freely available to facilitate its implementation (https://utrero-rico.shinyapps.io/COR12_Score/).

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.

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.

Asthma diagnosis using integrated analysis of eosinophil microRNAs.

BACKGROUND: Asthma is a syndrome characterized by airway inflammation and obstruction. Due to its heterogeneity, the difficulties in asthma diagnosis and treatment make the discovery of new biomarkers a focus of research. So, we determined the differential miRNA expression of eosinophils between healthy and asthmatic patients and to establish a differentially expressed miRNA profile detectable in sera for use as biomarker. METHODS: MicroRNAs from peripheral eosinophils from healthy and asthmatic subjects were isolated and analyzed by next-generation sequencing and confirmed by quantitative PCR in 29 asthmatics and 10 healthy individuals. The levels of serum miRNAs were performed by quantitative PCR in 138 asthmatics and 39 healthy subjects. Regression analysis and Random Forest models were performed. RESULTS: We found a set of miRNAs whose expression differs between eosinophils from asthmatics and healthy subjects. These miRNAs can classify asthmatics into two clusters that differed in the number of eosinophils and periostin concentration in serum. Some of these miRNAs were also confirmed in sera, as miR-185-5p which discriminates asthmatics from healthy subjects. Together with other two miRNAs, miR-185-5p allowed us to create a logistic regression model to discriminate better both conditions and a Random Forest model that can even sort the asthmatics into intermittent, mild persistent, moderate persistent, and severe persistent asthma. CONCLUSION: Our data show that miRNAs profile in eosinophils can be used as asthma diagnosis biomarker in serum and that this profile is able to rank asthma severity.

Dosimetry-based high-activity therapy with 131I-metaiodobenzylguanidine (131I-mIBG) and topotecan for the treatment of high-risk refractory neuroblastoma.

PURPOSE: Patients with high-risk neuroblastoma have an increased risk of recurrence and relapse of disease and a very poor prognosis. 131I-metaiodobenzylguanidine (131I-mIBG) in combination with topotecan as a radiosensitizer can be an effective and relatively well-tolerated agent for the treatment of refractory neuroblastoma. The aim of this retrospective study was to evaluate response and outcome of combined therapy with 131I-mIBG and topotecan. METHODS: Ten patients, between 3 and 20 years of age, were included. Nine patients had been refractory to several lines of chemotherapy and radiotherapy. One patient with a very high-risk neuroblastoma had received only induction therapy. Response was graded according to the International Neuroblastoma Staging System. RESULTS: Regarding treatment response, two patients achieved complete remission, one with relapse at 16 months, five achieved a partial remission, four showed progression at between 1 and 18 months; two showed stable disease with progression at between 1 and 5 months, and one showed progressive disease. Eight of the ten patients died with overall survival between 4 and 63 months, and two patients were still alive without disease at the time of this report: 52 and 32 months (patient had received only induction therapy). Acute and subacute adverse effects were mainly haematological, and one patient developed a differentiated thyroid cancer. CONCLUSION: In patients with high-risk refractory neuroblastoma, administration of high activities of 131I-mIBG in combination with topotecan was found to be an effective therapy, increasing overall survival and progression-free survival. Further studies including a larger number of patients and using 131I-mIBG for first-line up-front therapy are warranted.

Systematic identification of phosphorylation-mediated protein interaction switches.

Proteomics techniques can identify thousands of phosphorylation sites in a single experiment, the majority of which are new and lack precise information about function or molecular mechanism. Here we present a fast method to predict potential phosphorylation switches by mapping phosphorylation sites to protein-protein interactions of known structure and analysing the properties of the protein interface. We predict 1024 sites that could potentially enable or disable particular interactions. We tested a selection of these switches and showed that phosphomimetic mutations indeed affect interactions. We estimate that there are likely thousands of phosphorylation mediated switches yet to be uncovered, even among existing phosphorylation datasets. The results suggest that phosphorylation sites on globular, as distinct from disordered, parts of the proteome frequently function as switches, which might be one of the ancient roles for kinase phosphorylation.