A novel method for the multiplexed target enrichment of MinION next generation sequencing libraries using PCR-generated baits

The enrichment of targeted regions within complex next generation sequencing libraries commonly uses biotinylated baits to capture the desired sequences. This method results in high read coverage over the targets and their flanking regions. Oxford Nanopore Technologies recently released an USB3.0-interfaced sequencer, the MinION. To date no particular method for enriching MinION libraries has been standardized. Here, using biotinylated PCR-generated baits in a novel approach, we describe a simple and efficient way for multiplexed enrichment of MinION libraries, overcoming technical limitations related with the chemistry of the sequencing-adapters and the length of the DNA fragments. Using Phage Lambda and Escherichia coli as models we selectively enrich for specific targets, significantly increasing the corresponding read-coverage, eliminating unwanted regions. We show that by capturing genomic fragments, which contain the target sequences, we recover reads extending targeted regions and thus can be used for the determination of potentially unknown flanking sequences. By pooling enriched libraries derived from two distinct E. coli strains and analyzing them in parallel, we demonstrate the efficiency of this method in multiplexed format. Crucially we evaluated the optimal bait size for large fragment libraries and we describe for the first time a standardized method for target enrichment in MinION platform

Nanopore sequencing and full genome de novo assembly of human cytomegalovirus TB40/E reveals clonal diversity and structural variations.

BACKGROUND: Human cytomegalovirus (HCMV) has a double-stranded DNA genome of approximately 235 Kbp that is structurally complex including extended GC-rich repeated regions. Genomic recombination events are frequent in HCMV cultures but have also been observed in vivo. Thus, the assembly of HCMV whole genomes from technologies producing shorter than 500 bp sequences is technically challenging. Here we improved the reconstruction of HCMV full genomes by means of a hybrid, de novo genome-assembly bioinformatics pipeline upon data generated from the recently released MinION MkI B sequencer from Oxford Nanopore Technologies. RESULTS: The MinION run of the HCMV (strain TB40/E) library resulted in ~ 47,000 reads from a single R9 flowcell and in ~ 100× average read depth across the virus genome. We developed a novel, self-correcting bioinformatics algorithm to assemble the pooled HCMV genomes in three stages. In the first stage of the bioinformatics algorithm, long contigs (N50 = 21,892) of lower accuracy were reconstructed. In the second stage, short contigs (N50 = 5686) of higher accuracy were assembled, while in the final stage the high quality contigs served as template for the correction of the longer contigs resulting in a high-accuracy, full genome assembly (N50 = 41,056). We were able to reconstruct a single representative haplotype without employing any scaffolding steps. The majority (98.8%) of the genomic features from the reference strain were accurately annotated on this full genome construct. Our method also allowed the detection of multiple alternative sub-genomic fragments and non-canonical structures suggesting rearrangement events between the unique (UL /US) and the repeated (T/IRL/S) genomic regions. CONCLUSIONS: Third generation high-throughput sequencing technologies can accurately reconstruct full-length HCMV genomes including their low-complexity and highly repetitive regions. Full-length HCMV genomes could prove crucial in understanding the genetic determinants and viral evolution underpinning drug resistance, virulence and pathogenesis

The interferon receptor-1 promoter polymorphisms affect the outcome of Caucasians with HBeAg-negative chronic HBV infection

The outcome of HBeAg-negative chronic hepatitis B virus (HBV) patients who may remain in the inactive carrier state (IC) or progress to HBeAg-negative chronic hepatitis B may be affected by the host genetic profile. Genetic polymorphisms within not only the promoter but also the coding sequence of the interferon receptor 1 (INFAR1) gene have been associated with susceptibility to chronic HBV infection, but their role on the outcomes of HBeAg-negative patients has not been evaluated. We examined the association of INFAR1 promoter polymorphisms with the phase of chronic HBV infection in a demographically characterized Caucasian cohort of 183 consecutive HBeAg-negative chronic HBV patients.Using a combination of conventional and allele-specific polymerase chain reactions, bidirectional sequencing and DNA-fragment analysis, we performed typing of three Single Nucleotide Polymorphisms (SNPs -568G/C, -408C/T, -3C/T) and one Variable Number Tandem Repeat [VNTR -77(GT)n] within the INFR1 promoter sequence.The genetic polymorphisms examined were found to be associated with the phase of HBeAg-negative chronic HBV patients. Using a multiple logistic regression model adjusting for age, gender and origin of the individuals, we found that patients with linked genotypes -408CT_-3CT were more likely to be ICs (OR = 2.42 vs. CC, P = 0.036). Also, given the partial linkage between SNP -568G/C and VNTR -77(GT)n, we found that linked genotypes -77(GT)n ≤ 8/≤8_-568GC and -77(GT)n ≤ 8/≤8_-568CC were detected more frequently among ICs (OR = 11.69, P = 0.005 and OR = 7.56, P = 0.001 vs. -77(GT)n >8/>8_-568GG respectively).These findings suggest that these genetic variations represent important factors associated with the clinical phase of HBeAg-negative chronic HBV infection

Molecular and Clinical Prognostic Biomarkers of COVID-19 Severity and Persistence

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), poses several challenges to clinicians, due to its unpredictable clinical course. The identification of laboratory biomarkers, specific cellular, and molecular mediators of immune response could contribute to the prognosis and management of COVID-19 patients. Of utmost importance is also the detection of differentially expressed genes, which can serve as transcriptomic signatures, providing information valuable to stratify patients into groups, based on the severity of the disease. The role of biomarkers such as IL-6, procalcitonin, neutrophil–lymphocyte ratio, white blood cell counts, etc. has already been highlighted in recently published studies; however, there is a notable amount of new evidence that has not been summarized yet, especially regarding transcriptomic signatures. Hence, in this review, we assess the latest cellular and molecular data and determine the significance of abnormalities in potential biomarkers for COVID-19 severity and persistence. Furthermore, we applied Gene Ontology (GO) enrichment analysis using the genes reported as differentially expressed in the literature in order to investigate which biological pathways are significantly enriched. The analysis revealed a number of processes, such as inflammatory response, and monocyte and neutrophil chemotaxis, which occur as part of the complex immune response to SARS-CoV-2

Human endogenous retrovirus-K HML-2 integration within RASGRF2 is associated with intravenous drug abuse and modulates transcription in a cell-line model

HERV-K HML-2 (HK2) has been proliferating in the germ line of humans at least as recently as 250,000 years ago, with some integrations that remain polymorphic in the modern human population. One of the solitary HK2 LTR polymorphic integrations lies between exons 17 and 18 of RASGRF2, a gene that affects dopaminergic activity and is thus related to addiction. Here we show that this antisense HK2 integration (namely RASGRF2-int) is found more frequently in persons who inject drugs compared with the general population. In a Greek HIV-1–positive population (n = 202), we found RASGRF2-int 2.5 times (14 versus 6%) more frequently in patients infected through i.v. drug use compared with other transmission route controls (P = 0.03). Independently, in a United Kingdom-based hepatitis C virus-positive population (n = 184), we found RASGRF2-int 3.6 times (34 versus 9.5%) more frequently in patients infected during chronic drug abuse compared with controls (P ≺ 0.001). We then tested whether RASGRF2-int could be mechanistically responsible for this association by modulating transcription of RASGRF2. We show that the CRISPR/Cas9-mediated insertion of HK2 in HEK293 cells in the exact RASGRF2 intronic position found in the population resulted in significant transcriptional and phenotypic changes. We also explored mechanistic features of other intronic HK2 integrations and show that HK2 LTRs can be responsible for generation of cis-natural antisense transcripts, which could interfere with the transcription of nearby genes. Our findings suggest that RASGRF2-int is a strong candidate for dopaminergic manipulation, and emphasize the importance of accurate mapping of neglected HERV polymorphisms in human genomic studies

Spatiotemporal Characteristics of the Largest HIV-1 CRF02_AG Outbreak in Spain: Evidence for Onward Transmissions.

Background and Aim: The circulating recombinant form 02_AG (CRF02_AG) is the predominant clade among the human immunodeficiency virus type-1 (HIV-1) non-Bs with a prevalence of 5.97% (95% Confidence Interval-CI: 5.41-6.57%) across Spain. Our aim was to estimate the levels of regional clustering for CRF02_AG and the spatiotemporal characteristics of the largest CRF02_AG subepidemic in Spain. Methods: We studied 396 CRF02_AG sequences obtained from HIV-1 diagnosed patients during 2000-2014 from 10 autonomous communities of Spain. Phylogenetic analysis was performed on the 391 CRF02_AG sequences along with all globally sampled CRF02_AG sequences (N = 3,302) as references. Phylodynamic and phylogeographic analysis was performed to the largest CRF02_AG monophyletic cluster by a Bayesian method in BEAST v1.8.0 and by reconstructing ancestral states using the criterion of parsimony in Mesquite v3.4, respectively. Results: The HIV-1 CRF02_AG prevalence differed across Spanish autonomous communities we sampled from (

Upregulation of Human Endogenous Retroviruses in Bronchoalveolar Lavage Fluid of COVID-19 Patients

Severe COVID-19 pneumonia has been associated with the development of intense inflammatory responses during the course of infections with SARS-CoV-2. Given that human endogenous retroviruses (HERVs) are known to be activated during and participate in inflammatory processes, we examined whether HERV dysregulation signatures are present in COVID-19 patients. By comparing transcriptomes of bronchoalveolar lavage fluid (BALF) of COVID-19 patients and healthy controls, and peripheral blood monocytes (PBMCs) from patients and controls, we have shown that HERVs are intensely dysregulated in BALF of COVID-19 patients compared to those in BALF of healthy control patients but not in PBMCs. In particular, upregulation in the expression of specific HERV families was detected in BALF samples of COVID-19 patients, with HERV-FRD being the most highly upregulated family among the families analyzed. In addition, we compared the expression of HERVs in human bronchial epithelial cells (HBECs) without and after senescence induction in an oncogene-induced senescence model in order to quantitatively measure changes in the expression of HERVs in bronchial cells during the process of cellular senescence. This apparent difference of HERV dysregulation between PBMCs and BALF warrants further studies in the involvement of HERVs in inflammatory pathogenetic mechanisms as well as exploration of HERVs as potential biomarkers for disease progression. Furthermore, the increase in the expression of HERVs in senescent HBECs in comparison to that in noninduced HBECs provides a potential link for increased COVID-19 severity and mortality in aged populations

HCV Defective Genomes Promote Persistent Infection by Modulating the Viral Life Cycle

Defective interfering (DI) RNAs have been detected in several human viruses. HCV in-frame deletions mutants (IFDMs), missing mainly the envelope proteins, have been found in patient sera and liver tissues. IFDMs replicate independently and can be trans-packaged into infectious virions in the presence of full length viral genome. So far, their biological role is unclear. In this study, we have isolated and cloned IFDMs from sera samples and liver tissues of patients infected with HCV genotypes 1b, 2a, and 3a. IFDMs were present in up to 26% of samples tested. Using the in vitro HCV cell culture system, co-expression of the wild type (wt) HCV replicon with HCV IFDMs RNA resulted in increased HCV replication. Additionally, co-transfection of the HCV full length genome RNA and a defective mutant missing the envelope region led to increased viral release, collectively suggesting an important biological role for IFDMs in the virus life cycle. Recently, exosomes, masters of intercellular communication, have been implicated in the transport of HCV viral genomes. We report for the first time that exosomal RNA isolated from HCV sera samples contains HCV defective genomes. We also demonstrate that inhibition of exosomal biogenesis and release influences HCV viral replication. Overall, we provide evidence that the presence of HCV IFDMs affects both viral replication and release. IFDMs exploit exosomes as means of transport, a way to evade the immune system, to spread more efficiently and possibly maintain persistent infection

Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations

Erratum in: Nat Commun. 2023 Mar 20;14(1):1539. doi: 10.1038/s41467-023-37302-5.Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori fromEurope and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks bymigration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.This work was supported by Sequencing Grants-in-aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of Japan (221S0002, 18KK0266, 19H03473, 21H00346 and 22H02871) to Y.Y. F.F.V. is financed by FCT through Assistant Researcher grant CEECIND/03023/2017 and a project grant PTDC/BTM-TEC/3238/ 2020. I.K. studentship was funded by the National Strategic Reference Framework Operational Program “Competitiveness, Entrepreneurship and Innovation” (NSRF 2014-2020, project No. MIS5002486) and sequencing of strains was supported by the InfeNeutra Project (NSRF 2007-2013, project no. MIS450598) of the Ministry of Culture and Edu- cation, Greece. K.T. and the sequencing of KI isolates was supported by Erik Philip-Sörensen Foundation grant G2016-08, and Swedish Society for Medical research (SSMF). All primary bioinformatics and parts of the comparative genomics were performed on resources provided by Swedish National Infrastructure for Computing (SNIC) through Uppsala Multidisciplinary Center for Advanced Computational Science (UPPMAX) under projects snic2018-8-24 and uppstore2017270. Work by S.S. was supported by the German Research Foundation (DFG, project number 158 989 968–SFB 900/A1) and by the Bavarian Ministry of Sci- ence and the Arts in the framework of the Bavarian Research Network “New Strategies Against Multi-Resistant Pathogens by Means of Digital Networking—bayresq.net”. D.F. was supported by Shanghai Municipal Science and Technology Major Project No. 2019SHZDZX02.info:eu-repo/semantics/publishedVersio

SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave

International audienceBackground: The spatiotemporal profiling of molecular transmission clusters (MTCs) using viral genomic data can effectively identify transmission networks in order to inform public health actions targeting SARS-CoV-2 spread. Methods: We used whole genome SARS-CoV-2 sequences derived from ten European regions belonging to eight countries to perform phylogenetic and phylodynamic analysis. We developed dedicated bioinformatics pipelines to identify regional MTCs and to assess demographic factors potentially associated with their formation. Results: The total number and the scale of MTCs varied from small household clusters identified in all regions, to a super-spreading event found in Uusimaa-FI. Specific age groups were more likely to belong to MTCs in different regions. The clustered sequences referring to the age groups 50–100 years old (y.o.) were increased in all regions two weeks after the establishment of the lockdown, while those referring to the age group 0–19 y.o. decreased only in those regions where schools’ closure was combined with a lockdown. Conclusions: The spatiotemporal profiling of the SARS-CoV-2 MTCs can be a useful tool to monitor the effectiveness of the interventions and to reveal cryptic transmissions that have not been identified through contact tracing