Optimizing in silico drug discovery: simulation of connected differential expression signatures and applications to benchmarking.

BACKGROUND: We present a novel simulation method for generating connected differential expression signatures. Traditional methods have struggled with the lack of reliable benchmarking data and biases in drug-disease pair labeling, limiting the rigorous benchmarking of connectivity-based approaches. OBJECTIVE: Our aim is to develop a simulation method based on a statistical framework that allows for adjustable levels of parametrization, especially the connectivity, to generate a pair of interconnected differential signatures. This could help to address the issue of benchmarking data availability for connectivity-based drug repurposing approaches. METHODS: We first detailed the simulation process and how it reflected real biological variability and the interconnectedness of gene expression signatures. Then, we generated several datasets to enable the evaluation of different existing algorithms that compare differential expression signatures, providing insights into their performance and limitations. RESULTS: Our findings demonstrate the ability of our simulation to produce realistic data, as evidenced by correlation analyses and the log2 fold-change distribution of deregulated genes. Benchmarking reveals that methods like extreme cosine similarity and Pearson correlation outperform others in identifying connected signatures. CONCLUSION: Overall, our method provides a reliable tool for simulating differential expression signatures. The data simulated by our tool encompass a wide spectrum of possibilities to challenge and evaluate existing methods to estimate connectivity scores. This may represent a critical gap in connectivity-based drug repurposing research because reliable benchmarking data are essential for assessing and advancing in the development of new algorithms. The simulation tool is available as a R package (General Public License (GPL) license) at https://github.com/cgonzalez-gomez/cosimu.

Mucosal bivalent live attenuated vaccine protects against human metapneumovirus and respiratory syncytial virus in mice.

Live-Attenuated Vaccines (LAVs) stimulate robust mucosal and cellular responses and have the potential to protect against Respiratory Syncytial Virus (RSV) and Human Metapneumovirus (HMPV), the main etiologic agents of viral bronchiolitis and pneumonia in children. We inserted the RSV-F gene into an HMPV-based LAV (Metavac®) we previously validated for the protection of mice against HMPV challenge, and rescued a replicative recombinant virus (Metavac®-RSV), exposing both RSV- and HMPV-F proteins at the virion surface and expressing them in reconstructed human airway epithelium models. When administered to BALB/c mice by the intranasal route, bivalent Metavac®-RSV demonstrated its capacity to replicate with reduced lung inflammatory score and to protect against both RSV and lethal HMPV challenges in vaccinated mice while inducing strong IgG and broad RSV and HMPV neutralizing antibody responses. Altogether, our results showed the versatility of the Metavac® platform and suggested that Metavac®-RSV is a promising mucosal bivalent LAV candidate to prevent pneumovirus-induced diseases.

Transcriptional Profiling of Immune and Inflammatory Responses in the Context of SARS-CoV-2 Fungal Superinfection in a Human Airway Epithelial Model.

An increasing amount of evidence indicates a relatively high prevalence of superinfections associated with coronavirus disease 2019 (COVID-19), including invasive aspergillosis, but the underlying mechanisms remain to be characterized. In the present study, to better understand the biological impact of superinfection, we determine and compare the host transcriptional response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) versus Aspergillus superinfection, using a model of reconstituted human airway epithelium. Our analyses reveal that both simple infection and superinfection induce strong deregulation of core components of innate immune and inflammatory responses, with a stronger response to superinfection in the bronchial epithelial model compared to its nasal counterpart. Our results also highlight unique transcriptional footprints of SARS-CoV-2 Aspergillus superinfection, such as an imbalanced type I/type III IFN, and an induction of several monocyte and neutrophil associated chemokines, that could be useful for the understanding of Aspergillus-associated COVID-19 and also the management of severe forms of aspergillosis in this specific context.

Sequencing the Genome of Indian Flying Fox, Natural Reservoir of Nipah Virus, Using Hybrid Assembly and Conservative Secondary Scaffolding.

Indian fruit bats, flying fox Pteropus medius was identified as an asymptomatic natural host of recently emerged Nipah virus, which is known to induce a severe infectious disease in humans. The absence of P. medius genome sequence presents an important obstacle for further studies of virus-host interactions and better understanding of mechanisms of zoonotic viral emergence. Generation of the high-quality genome sequence is often linked to a considerable effort associated to elevated costs. Although secondary scaffolding methods have reduced sequencing expenses, they imply the development of new tools for the integration of different data sources to achieve more reliable sequencing results. We initially sequenced the P. medius genome using the combination of Illumina paired-end and Nanopore sequencing, with a depth of 57.4x and 6.1x, respectively. Then, we introduced the novel scaff2link software to integrate multiple sources of information for secondary scaffolding, allowing to remove the association with discordant information among two sources. Different quality metrics were next produced to validate the benefits from secondary scaffolding. The P. medius genome, assembled by this method, has a length of 1,985 Mb and consists of 33,613 contigs and 16,113 scaffolds with an NG50 of 19 Mb. At least 22.5% of the assembled sequences is covered by interspersed repeats already described in other species and 19,823 coding genes are annotated. Phylogenetic analysis demonstrated the clustering of P. medius genome with two other Pteropus bat species, P. alecto and P. vampyrus, for which genome sequences are currently available. SARS-CoV entry receptor ACE2 sequence of P. medius was 82.7% identical with ACE2 of Rhinolophus sinicus bats, thought to be the natural host of SARS-CoV. Altogether, our results confirm that a lower depth of sequencing is enough to obtain a valuable genome sequence, using secondary scaffolding approaches and demonstrate the benefits of the scaff2link application. The genome sequence is now available to the scientific community to (i) proceed with further genomic analysis of P. medius, (ii) to characterize the underlying mechanism allowing Nipah virus maintenance and perpetuation in its bat host, and (iii) to monitor their evolutionary pathways toward a better understanding of bats' ability to control viral infections.

Characterization and Treatment of SARS-CoV-2 in Nasal and Bronchial Human Airway Epithelia.

In the current COVID-19 pandemic context, proposing and validating effective treatments represents a major challenge. However, the scarcity of biologically relevant pre-clinical models of SARS-CoV-2 infection imposes a significant barrier for scientific and medical progress, including the rapid transition of potentially effective treatments to the clinical setting. We use reconstituted human airway epithelia to isolate and then characterize the viral infection kinetics, tissue-level remodeling of the cellular ultrastructure, and transcriptional early immune signatures induced by SARS-CoV-2 in a physiologically relevant model. Our results emphasize distinctive transcriptional immune signatures between nasal and bronchial HAE, both in terms of kinetics and intensity, hence suggesting putative intrinsic differences in the early response to SARS-CoV-2 infection. Most important, we provide evidence in human-derived tissues on the antiviral efficacy of remdesivir monotherapy and explore the potential of the remdesivir-diltiazem combination as an option worthy of further investigation to respond to the still-unmet COVID-19 medical need.

High Pathogenicity of Nipah Virus from Pteropus lylei Fruit Bats, Cambodia.

We conducted an in-depth characterization of the Nipah virus (NiV) isolate previously obtained from a Pteropus lylei bat in Cambodia in 2003 (CSUR381). We performed full-genome sequencing and phylogenetic analyses and confirmed CSUR381 is part of the NiV-Malaysia genotype. In vitro studies revealed similar cell permissiveness and replication of CSUR381 (compared with 2 other NiV isolates) in both bat and human cell lines. Sequence alignments indicated conservation of the ephrin-B2 and ephrin-B3 receptor binding sites, the glycosylation site on the G attachment protein, as well as the editing site in phosphoprotein, suggesting production of nonstructural proteins V and W, known to counteract the host innate immunity. In the hamster animal model, CSUR381 induced lethal infections. Altogether, these data suggest that the Cambodia bat-derived NiV isolate has high pathogenic potential and, thus, provide insight for further studies and better risk assessment for future NiV outbreaks in Southeast Asia.

Massively parallel and multiplex blood group genotyping using next-generation-sequencing.

OBJECTIVES: Thirty-six blood group systems are listed by the International Society of Blood Transfusion, containing almost 350 antigens. Most of these result from a single nucleotide polymorphism (SNP). Serology is the standard method for blood group typing. However, this technique has some limitations and cannot respond to the growing demand of blood product typing for a large number of antigens. Here we describe a blood group genotyping assay directly from whole blood samples using Next-Generation Sequencing (NGS), allowing the simultaneous identification of 15 SNPs associated with the blood group systems of 95 patients in a single run. DESIGN AND METHOD: After an automated DNA extraction, targets are amplified by multiplex polymerase chain reaction (PCRm). Two panels addressing 9 groups have been developed (MNS, Lutheran, Kell, Duffy, Kidd, Diego, Yt, Dombrock, and Colton), one for 8 SNPs, the other for 7 SNPs. For each sample, both panels corresponding to 14 amplicons (1 amplicon containing 2 SNPs) are pooled. Then a dual-indexed library is generated from each pool by linking Illumina adaptors directly onto amplicons, followed by sequencing using the MiSeq platform (Illumina). RESULTS: In a single experiment, 95 blood donor samples have been sequenced for the genes of interest. Among the 1425 targeted single nucleotide polymorphisms, 1420 were identified by sequencing, reflecting a coverage of 99.65%. The obtained data shows a good correlation (99% for all SNPs) with other blood group typing methods. Depending on the allele pairs analyzed, correlations vary between 97.12 and 100%. CONCLUSION: Next-Generation sequencing would supplement serological and molecular techniques and, in the near future, could replace it with complete and fast results acquisition for pre-screening and identification of rare blood bags.