Several pathogenic viruses such as hepatitis B and human immunodeficiency viruses may integrate into the host genome. These virus/host integrations are detectable using paired-end next generation sequencing. However, the low number of expected true virus integrations may be difficult to distinguish from the noise of many false positive candidates. Here, we propose a novel filtering approach that increases specificity without compromising sensitivity for virus/host chimera detection. Our detection pipeline termed Vy-PER (Virus integration detection bY Paired End Reads) outperforms existing similar tools in speed and accuracy. We analysed whole genome data from childhood acute lymphoblastic leukemia (ALL), which is characterised by genomic rearrangements and usually associated with radiation exposure. This analysis was motivated by the recently reported virus integrations at genomic rearrangement sites and association with chromosomal instability in liver cancer. However, as expected, our analysis of 20 tumour and matched germline genomes from ALL patients finds no significant evidence for integrations by known viruses. Nevertheless, our method eliminates 12,800 false positives per genome (80× coverage) and only our method detects singleton human-phiX174-chimeras caused by optical errors of the Illumina HiSeq platform. This high accuracy is useful for detecting low virus integration levels as well as non-integrated viruses.

Post-mitotic cardiomyocytes have been considered to be non-permissive to precise targeted integration including homology-directed repair (HDR) after CRISPR/Cas9 genome editing. Here, we demonstrate that direct delivery of large amounts of transgene encoding guide RNA (gRNA) and repair template DNA via intra-ventricular injection of adeno-associated virus (AAV) promotes precise targeted genome replacement in adult murine cardiomyocytes expressing Cas9. Neither systemic injection of AAV nor direct injection of adenovirus promotes targeted integration, suggesting that high copy numbers of single-stranded transgenes are required in cardiomyocytes. Notably, AAV-mediated targeted integration in cardiomyocytes both in vitro and in vivo depends on the Fanconi anemia pathway, a key component of the single-strand template repair mechanism. In human cardiomyocytes differentiated from induced pluripotent stem cells, AAV-mediated targeted integration fluorescently labeled Mlc2v protein after differentiation, independently of DNA synthesis, and enabled real-time detection of sarcomere contraction in monolayered beating cardiomyocytes. Our findings provide a wide range of applications for targeted genome replacement in non-dividing cardiomyocytes.

Between October 2016 and December 2017, several European Countries had been involved in a massive Highly Pathogenic Avian Influenza (HPAI) epidemic sustained by H5N8 subtype virus. Starting on December 2016, also Italy was affected by H5N8 HPAI virus, with cases occurring in two epidemic waves: the first between December 2016 and May 2017, and the second in July-December 2017. Eighty-three outbreaks were recorded in poultry, 67 of which (80.72%) occurring in the second wave. A total of 14 cases were reported in wild birds. Epidemiological information and genetic analyses were conjointly used to get insight on the spread dynamics. Analyses indicated multiple introductions from wild birds to the poultry sector in the first epidemic wave, and noteworthy lateral spread from October 2017 in a limited geographical area with high poultry densities. Turkeys, layers and backyards were the mainly affected types of poultry production. Two genetic sub-groups were detected in the second wave in non-overlapping geographical areas, leading to speculate on the involvement of different wild bird populations. The integration of epidemiological data and genetic analyses allowed to unravel the transmission dynamics of H5N8 virus in Italy, and could be exploited to timely support in implementing tailored control measures.

Bluetongue virus (BTV), a complex double-stranded segmented RNA virus, has been found to initiate cellular autophagy for its own benefit. Here, with a view to understanding the underlying mechanisms, we first systematically dissected the exact signaling network in BTV-induced autophagy. We found that the activity of mTOR, a crucial pivot, was inhibited by BTV1 infection, subsequently leading to downstream p70S6K suppression and autophagy initiation. We then explored the upstream regulators of mTOR and analyzed their activities via a series of assays. We found BTV1-induced autophagy to be independent of the ERK1/2 signaling pathway. However, the BTV1-induced inhibition of PI3K/Akt was found to be partially responsible for mTOR inactivation and subsequent autophagy initiation. Furthermore, we found unexpectedly that AMPK seemed to play a more important role in BTV1-induced autophagy. Elevated [Ca(2+)]cyto-mediated activation of CaMKKβ exactly managed the activation of AMPK, which then positively regulated autophagy through suppressing mTOR. We must emphasize that TSC2 is a fatal mediator between upstream Akt or AMPK and downstream mTOR through its phosphorylation. Taken together, our data suggested that the BTV1-induced inhibition of the Akt-TSC2-mTOR pathway and the upregulation of the AMPK-TSC2-mTOR pathway both contributed to autophagy initiation and further favored virus replication.

The present study investigated the association between single nucleotide polymorphisms (SNPs) in immune- or DNA repair-related genes and the integration pattern of human papillomavirus (HPV), a promising prognostic marker in cervical cancer. The HPV integration patterns of cervical cancer patients were determined by polymerase chain reaction and in situ hybridization, and categorized as episomal (group A), single-copy or multi-copy tandem repetition integrated (group B), and undetectable HPV types (group C). After sample and SNP quality control, 166,505 SNPs in 161 samples (38, 111, and 12 patients in groups A, B, and C, respectively) were examined. None of the SNPs reached genome-wide significance, and several candidate SNPs for future study were selected, including rs10999435 on chromosome 10q22, rs1322054 on chromosome 9q32-33, and rs10902171 on chromosome 11p15. Luciferase assay identified rs1322054 as the primary functional variant to regulate gene expression in immune cell. Further studies are needed to determine the genetic background of different integration patterns of HPV in cervical cancer patients.

Circulating tumor DNA (ctDNA) in blood plasma is present at very low concentrations compared to cell-free DNA (cfDNA) of non-tumor origin. To enhance ctDNA detection, recent studies have been focused on understanding the non-random fragmentation pattern of cfDNA. These studies have investigated fragment sizes, genomic position of fragment end points, and fragment end motifs. Although these features have been described and shown to be aberrant in cancer patients, there is a lack of understanding of how the individual and integrated analysis of these features enrich ctDNA fraction and enhance ctDNA detection. Using whole genome sequencing and copy number analysis of plasma samples from 5 high grade serious ovarian cancer patients, we observed that (1) ctDNA is enriched not only in fragments shorter than mono-nucleosomes (~ 167 bp), but also in those shorter than di-nucleosomes (~ 240-330 bp) (28-159% enrichment). (2) fragments that start and end at the border or within the nucleosome core are enriched in ctDNA (5-46% enrichment). (3) certain DNA motifs conserved in regions 10 bp up- and down- stream of fragment ends (i.e. cleavage sites) could be used to detect tumor-derived fragments (10-44% enrichment). We further show that the integrated analysis of these three features resulted in a higher enrichment of ctDNA when compared to using fragment size alone (additional 7-25% enrichment after fragment size selection). We believe these genome wide features, which are independent of genetic mutational changes, could allow new ways to analyze and interpret cfDNA data, as significant aberrations of these features from a healthy state could improve its utility as a diagnostic biomarker.

The genetic modification of cattle has many agricultural and biomedical applications. However, random integration often results in the unstable expression of transgenes and unpredictable phenotypes. Targeting genes to the "safe locus" and stably expressing foreign genes at a high level are desirable methods for overcoming these hurdles. The Rosa26 locus has been widely used to produce genetically modified animals in some species expressing transgenes at high and consistent levels. For the first time, we identified a bovine orthologue of the mouse Rosa26 locus through a genomic sequence homology analysis. According to 5' rapid-amplification of cDNA ends (5'RACE), 3' rapid-amplification of cDNA ends (3'RACE), reverse transcription PCR (RT-PCR) and quantitative PCR (Q-PCR) experiments, this locus encodes a long noncoding RNA (lncRNA) comprising two exons that is expressed ubiquitously and stably in different tissues. The bovine Rosa26 (bRosa26) locus appears to be highly amenable to transcription activator-like effector nucleases (TALENs)-mediated knock-in, and ubiquitous expression of enhanced green fluorescent protein (EGFP) inserted in the bRosa26 locus was observed in various stages, including cells, embryos, fetus and cattle. Finally, we created a valuable master bRosa26-EGFP fetal fibroblast cell line in which any gene of interest can be efficiently introduced and stably expressed using recombinase-mediated cassette exchange (RMCE). The new tools described here will be useful for a variety of studies using cattle.

The ability to generate integration-free induced hepatocyte-like cells (iHeps) from somatic fibroblasts has the potential to advance their clinical application. Here, we have generated integration-free, functional, and expandable iHeps from mouse somatic fibroblasts. To elicit this direct conversion, we took advantage of an oriP/EBNA1-based episomal system to deliver a set of transcription factors, Gata4, Hnf1a, and Foxa3, to the fibroblasts. The established iHeps exhibit similar morphology, marker expression, and functional properties to primary hepatocytes. Furthermore, integration-free iHeps prolong the survival of fumarylacetoacetate-hydrolase-deficient (Fah(-/-)) mice after cell transplantation. Our study provides a novel concept for generating functional and expandable iHeps using a non-viral, non-integrating, plasmid-based system that could facilitate their pharmaceutical and biomedical application.

GCN2 is a serine/threonine kinase involved in cellular stress response related to amino acid starvation. Previously, we showed that GCN2 interacts with HIV-1 integrase and is activated during HIV-1 infection. Herein, we identified HIV-1 integrase as a previously unknown substrate of GCN2 in vitro with a major site of phosphorylation at residue S255 located in the C-terminal domain of HIV-1 integrase. The underlying mechanism was investigated and it appeared that the integrase active site was required in order for GCN2 to target the integrase residue S255. Moreover, various integrases from other retroviruses (e.g. MLV, ASV) were also recognized as a substrate by GCN2. In cells, HIV-1 lentiviral particles harboring mutation at integrase position 255 were affected in their replication. Preventing phosphorylation resulted in an increase in infectivity that correlated with an increase in viral DNA integration. Infectivity of MLV was also higher in cells knocked-out for GCN2 suggesting a conserved mechanism to control viral replication. Altogether, our data suggest that GCN2 may constitute a general guardian of genome stability by regulating foreign DNA integration and as such be part of the antiviral armamentarium of the cell.

Breast cancer (BC) has a higher incidence in young Lebanese woman as compared to the West. We assessed the microRNA (miRNA) microarray profile of tissues derived from Lebanese patients with early BC and performed mRNA-miRNA integration analysis. 173 miRNAs were significantly dysregulated in 45 BC versus 17 normal adjacent breast tissues, including 74 with a fold change more than two of which 17 were never reported before in cancer. Integration analysis of mRNA-miRNA microarray data revealed a potential role of 51 dysregulated miRNA regulating 719 tumor suppressive or oncogenic mRNA associated with increased proliferation and decreased migration and invasion. We then performed a comparative miRNA microarray profile analysis of BC tissue between these 45 Lebanese and 197 matched American BC patients. Notably, Lebanese BC patients had 21 exclusively dysregulated miRNA (e.g. miR-31, 362-3p, and 663) and 4 miRNA with different expression manner compared to American patients (e.g. miR-1288-star and 324-3p). Some of these differences could reflect variation in patient age at diagnosis or ethnic variation affecting miRNA epigenetic regulation or sequence of miRNA precursors. Our data provide a basis for genetic/epigenetic investigations to explore the role of miRNA in early stage BC in young women, including ethnic specific differences.

HPV genomic variability and chromosomal integration are important in the HPV-induced carcinogenic process. To uncover these genomic events in an HPV infection, we have developed an innovative and cost-effective sequencing approach named TaME-seq (tagmentation-assisted multiplex PCR enrichment sequencing). TaME-seq combines tagmentation and multiplex PCR enrichment for simultaneous analysis of HPV variation and chromosomal integration, and it can also be adapted to other viruses. For method validation, cell lines (n = 4), plasmids (n = 3), and HPV16, 18, 31, 33 and 45 positive clinical samples (n = 21) were analysed. Our results showed deep HPV genome-wide sequencing coverage. Chromosomal integration breakpoints and large deletions were identified in HPV positive cell lines and in one clinical sample. HPV genomic variability was observed in all samples allowing identification of low frequency variants. In contrast to other approaches, TaME-seq proved to be highly efficient in HPV target enrichment, leading to reduced sequencing costs. Comprehensive studies on HPV intra-host variability generated during a persistent infection will improve our understanding of viral carcinogenesis. Efficient identification of both HPV variability and integration sites will be important for the study of HPV evolution and adaptability and may be an important tool for use in cervical cancer diagnostics.

Antibiotics at sub-inhibitory concentrations are often found in the environment. Here they could impose selective pressure on bacteria, leading to the selection and dissemination of antibiotic resistance, despite being under the inhibitory threshold. The goal of this study was to evaluate the effects of sub-inhibitory concentrations of gentamicin on environmental class 1 integron cassettes in natural river microbial communities. Gentamicin at sub-inhibitory concentrations promoted the integration and selection of gentamicin resistance genes (GmRG) in class 1 integrons after only a one-day exposure. Therefore, sub-inhibitory concentrations of gentamicin induced integron rearrangements, increasing the mobilization potential of gentamicin resistance genes and potentially increasing their dissemination in the environment. This study demonstrates the effects of antibiotics at sub-inhibitory concentrations in the environment and supports concerns about antibiotics as emerging pollutants.

Bovine leukemia virus (BLV) is an oncogenic retrovirus which induces malignant lymphoma termed enzootic bovine leukosis (EBL) after a long incubation period. Insertion sites of the BLV proviral genome as well as the associations between disease progression and polymorphisms of the virus and host genome are not fully understood. To characterize the biological coherence between virus and host, we developed a DNA-capture-seq approach, in which DNA probes were used to efficiently enrich target sequence reads from the next-generation sequencing (NGS) library. In addition, enriched reads can also be analyzed for detection of proviral integration sites and clonal expansion of infected cells since the reads include chimeric reads of the host and proviral genomes. To validate this DNA-capture-seq approach, a persistently BLV-infected fetal lamb kidney cell line (FLK-BLV), four EBL tumor samples and four non-EBL blood samples were analyzed to identify BLV integration sites. The results showed efficient enrichment of target sequence reads and oligoclonal integrations of the BLV proviral genome in the FLK-BLV cell line. Moreover, three out of four EBL tumor samples displayed multiple integration sites of the BLV proviral genome, while one sample displayed a single integration site. In this study, we found the evidence for the first time that the integrated provirus defective at the 5' end was present in the persistent lymphocytosis cattle. The efficient and sensitive identification of BLV variability, integration sites and clonal expansion described in this study provide support for use of this innovative tool for understanding the detailed mechanisms of BLV infection during the course of disease progression.

Transgenic (Tg) mice containing bacterial artificial chromosome (BAC) DNA are widely used for gene expression analysis and gene therapy models because BAC transgenes provide gene expression at physiological levels with the same developmental timing as endogenous genes. To ensure correct interpretation of transgene functions, investigation of the genomic organisation and integration of the BAC transgene is required. Here, we describe a reliable method based on droplet digital PCR (ddPCR) and inverse PCR to estimate copy number, genomic organisation and insertion sites of BAC transgenes in the mouse genome. We generated BAC Tg mice containing fragments of BAC clone RP23-59P20. ddPCR and iPCR analysis showed that the transgene consisted of five fragments of the BAC clone containing the Mkrn3 gene region, and that the transgene was inserted into Bckdhb, homozygous deletion of which causes the maple syrup urine disease phenotype. The ddPCR method described here should prove useful for analysis of genomic organisation and integration of BAC transgenes.

Despite numerous studies done on understanding the role of DNA methylation, limited work has focused on systems integration of cell type-specific interplay between DNA methylation and gene transcription. Through a genome-wide analysis of DNA methylation across 19 cell types with T-47D as reference, we identified 106,252 cell type-specific differentially-methylated CpGs categorized into 7,537 differentially (46.6% hyper- and 53.4% hypo-) methylated regions. We found 44% promoter regions and 75% CpG islands were T-47D cell type-specific methylated. Pyrosequencing experiments validated the cell type-specific methylation across three benchmark cell lines. Interestingly, these DMRs overlapped with 1,145 known tumor suppressor genes. We then developed a Bayesian Gaussian Regression model to measure the relationship among DNA methylation, genomic segment distribution, differential gene expression and tumor suppressor gene status. The model uncovered that 3'UTR methylation has much less impact on transcriptional activity than other regions. Integration of DNA methylation and 82 transcription factor binding information across the 19 cell types suggested diverse interplay patterns between the two regulators. Our integrative analysis reveals cell type-specific and genomic region-dependent regulatory patterns and provides a perspective for integrating hundreds of various omics-seq data together.

Cassava is a major staple food for about 800 million people in the tropics and sub-tropical regions of the world. Production of cassava is significantly hampered by cassava brown streak disease (CBSD), caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). The disease is suppressing cassava yields in eastern Africa at an alarming rate. Previous studies have documented that CBSV is more devastating than UCBSV because it more readily infects both susceptible and tolerant cassava cultivars, resulting in greater yield losses. Using whole genome sequences from NGS data, we produced the first coalescent-based species tree estimate for CBSV and UCBSV. This species framework led to the finding that CBSV has a faster rate of evolution when compared with UCBSV. Furthermore, we have discovered that in CBSV, nonsynonymous substitutions are more predominant than synonymous substitution and occur across the entire genome. All comparative analyses between CBSV and UCBSV presented here suggest that CBSV may be outsmarting the cassava immune system, thus making it more devastating and harder to control.

Ebola virus (EBOV), isolate Makona, was the causative agent of the West African epidemic devastating predominantly Guinea, Liberia and Sierra Leone from 2013-2016. While several experimental vaccine and treatment approaches have been accelerated through human clinical trials, there is still no approved countermeasure available against this disease. Here, we report the construction and preclinical efficacy testing of a novel recombinant modified vaccinia Ankara (MVA)-based vaccine expressing the EBOV-Makona glycoprotein GP and matrix protein VP40 (MVA-EBOV). GP and VP40 form EBOV-like particles and elicit protective immune responses. In this study, we report 100% protection against lethal EBOV infection in guinea pigs after prime/boost vaccination with MVA-EBOV. Furthermore, this MVA-EBOV protected macaques from lethal disease after a single dose or prime/boost vaccination. The vaccine elicited a variety of antibody responses to both antigens, including neutralizing antibodies and antibodies with antibody-dependent cellular cytotoxic activity specific for GP. This is the first report that a replication-deficient MVA vector can confer full protection against lethal EBOV challenge after a single dose vaccination in macaques.

In situ detoxification of lignocellulose-derived microbial inhibitory compounds is an economical strategy for the fermentation of lignocellulose-derived sugars to fuels and chemicals. In this study, we investigated homologous integration and constitutive expression of Cbei_3974 and Cbei_3904, which encode aldo-keto reductase and previously annotated short chain dehydrogenase/reductase, respectively, in Clostridium beijerinckii NCIMB 8052 (Cb), resulting in two strains: Cb_3974 and Cb_3904. Expression of Cbei_3974 led to 2-fold increase in furfural detoxification relative to Cb_3904 and Cb_wild type. Correspondingly, butanol production was up to 1.2-fold greater in furfural-challenged cultures of Cb_3974 relative to Cb_3904 and Cb_wild type. With 4-hydroxybezaldehyde and syringaldehyde supplementation, Cb_3974 showed up to 2.4-fold increase in butanol concentration when compared to Cb_3904 and Cb_wild type. Syringic and vanillic acids were considerably less deleterious to all three strains of Cb tested. Overall, Cb_3974 showed greater tolerance to furfural, 4-hydroxybezaldehyde, and syringaldehyde with improved capacity for butanol production. Hence, development of Cb_3974 represents a significant progress towards engineering solventogenic Clostridium species that are tolerant to lignocellulosic biomass hydrolysates as substrates for ABE fermentation.

Retroviral gene therapy offers immense potential to treat many genetic diseases and has already shown efficacy in clinical trials. However, retroviral vector mediated genotoxicity remains a major challenge and clinically relevant approaches to reduce integration near genes and proto-oncogenes are needed. Foamy retroviral vectors have several advantages over gammaretroviral and lentiviral vectors including a potentially safer integration profile and a lower propensity to activate nearby genes. Here we successfully retargeted foamy retroviral vectors away from genes and into satellite regions enriched for trimethylated histone H3 at lysine 9 by modifying the foamy virus Gag and Pol proteins. Retargeted foamy retroviral vectors integrated near genes and proto-oncogenes less often (p < 0.001) than controls. Importantly, retargeted foamy retroviral vectors can be produced at high, clinically relevant titers (>107 transducing units/ml), and unlike other reported retargeting approaches engineered target cells are not needed to achieve retargeting. As proof of principle for use in the clinic we show efficient transduction and retargeting in human cord blood CD34+ cells. The modified Gag and Pol helper constructs we describe will allow any investigator to simply use these helper plasmids during vector production to retarget therapeutic foamy retroviral vectors.

The gut microbiota interacts with infectious diseases and affects host immunity. Liver disease is also reportedly associated with changes in the gut microbiota. To elucidate the changes in the gut microbiota before and after hepatitis C virus (HCV) eradication through direct-acting antiviral (DAA) treatment in patients with chronic hepatitis C (CHC), we investigated 42 samples from 14 patients who received DAA therapy for HCV. Fecal samples were obtained before treatment (Pre), when treatment ended (EOT), and 24 weeks after treatment ended (Post24). The target V3-4 region of the 16S rRNA gene from fecal samples was amplified using the Illumina Miseq sequencing platform. The diversity of the gut microbiota did not significantly differ between Pre, EOT, and Post24. Principal coordinates analysis showed that for each patient, the values at Pre, EOT, and Post24 were concentrated within a small area. The linear discriminant analysis of effect size showed that the relative abundances of Faecalibacterium and Bacillus increased at EOT, further increased at Post24, and were significantly increased at Post24 compared to Pre. These suggest that changes in the gut microbiota should be considered as among the various effects observed on living organisms after HCV eradication.