Genomics in the long-read sequencing era.

Long-read sequencing (LRS) technologies have provided extremely powerful tools to explore genomes. While in the early years these methods suffered technical limitations, they have recently made significant progress in terms of read length, throughput, and accuracy and bioinformatics tools have strongly improved. Here, we aim to review the current status of LRS technologies, the development of novel methods, and the impact on genomics research. We will explore the most impactful recent findings made possible by these technologies focusing on high-resolution sequencing of genomes and transcriptomes and the direct detection of DNA and RNA modifications. We will also discuss how LRS methods promise a more comprehensive understanding of human genetic variation, transcriptomics, and epigenetics for the coming years.

Deciphering and modelling the TGF-ß signalling interplays specifying the dorsal-ventral axis of the sea urchin embryo.

During sea urchin development, secretion of Nodal and BMP2/4 ligands and their antagonists Lefty and Chordin from a ventral organiser region specifies the ventral and dorsal territories. This process relies on a complex interplay between the Nodal and BMP pathways through numerous regulatory circuits. To decipher the interplay between these pathways, we used a combination of treatments with recombinant Nodal and BMP2/4 proteins and a computational modelling approach. We assembled a logical model focusing on cell responses to signalling inputs along the dorsal-ventral axis, which was extended to cover ligand diffusion and enable multicellular simulations. Our model simulations accurately recapitulate gene expression in wild-type embryos, accounting for the specification of ventral ectoderm, ciliary band and dorsal ectoderm. Our model simulations further recapitulate various morphant phenotypes, reveal a dominance of the BMP pathway over the Nodal pathway and stress the crucial impact of the rate of Smad activation in dorsal-ventral patterning. These results emphasise the key role of the mutual antagonism between the Nodal and BMP2/4 pathways in driving early dorsal-ventral patterning of the sea urchin embryo.

ClustMMRA v2: A Scalable Computational Pipeline for the Identification of MicroRNA Clusters Acting Cooperatively on Tumor Molecular Subgroups.

In recent cancer genomics programs, large-scale profiling of microRNAs has been routinely used in order to better understand the role of microRNAs in gene regulation and disease. To support the analysis of such amount of data, scalability of bioinformatics pipelines is increasingly important to handle larger datasets.Here, we describe a scalable implementation of the clustered miRNA Master Regulator Analysis (clustMMRA) pipeline, developed to search for genomic clusters of microRNAs potentially driving cancer molecular subtyping. Genomically clustered microRNAs can be simultaneously expressed to work in a combined manner and jointly regulate cell phenotypes. However, the majority of computational approaches for the identification of microRNA master regulators are typically designed to detect the regulatory effect of a single microRNA.We have applied the clustMMRA pipeline to multiple pediatric tumor datasets, up to a hundred samples in size, demonstrating very satisfying performances of the software on large datasets. Results have highlighted genomic clusters of microRNAs potentially involved in several subgroups of the different pediatric cancers or specifically involved in the phenotype of a subgroup. In particular, we confirmed the cluster of microRNAs at the 14q32 locus to be involved in multiple pediatric cancers, showing its specific downregulation in tumor subgroups with aggressive phenotype.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements.

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

ChIP-exo signal associated with DNA-binding motifs provides insight into the genomic binding of the glucocorticoid receptor and cooperating transcription factors.

The classical DNA recognition sequence of the glucocorticoid receptor (GR) appears to be present at only a fraction of bound genomic regions. To identify sequences responsible for recruitment of this transcription factor (TF) to individual loci, we turned to the high-resolution ChIP-exo approach. We exploited this signal by determining footprint profiles of TF binding at single-base-pair resolution using ExoProfiler, a computational pipeline based on DNA binding motifs. When applied to our GR and the few available public ChIP-exo data sets, we find that ChIP-exo footprints are protein- and recognition sequence-specific signatures of genomic TF association. Furthermore, we show that ChIP-exo captures information about TFs other than the one directly targeted by the antibody in the ChIP procedure. Consequently, the shape of the ChIP-exo footprint can be used to discriminate between direct and indirect (tethering to other DNA-bound proteins) DNA association of GR. Together, our findings indicate that the absence of classical recognition sequences can be explained by direct GR binding to a broader spectrum of sequences than previously known, either as a homodimer or as a heterodimer binding together with a member of the ETS or TEAD families of TFs, or alternatively by indirect recruitment via FOX or STAT proteins. ChIP-exo footprints also bring structural insights and locate DNA:protein cross-link points that are compatible with crystal structures of the studied TFs. Overall, our generically applicable footprint-based approach uncovers new structural and functional insights into the diverse ways of genomic cooperation and association of TFs.

Histone Chaperone SSRP1 is Essential for Wnt Signaling Pathway Activity During Osteoblast Differentiation.

Cellular differentiation is accompanied by dramatic changes in chromatin structure which direct the activation of lineage-specific transcriptional programs. Structure-specific recognition protein-1 (SSRP1) is a histone chaperone which is important for chromatin-associated processes such as transcription, DNA replication and repair. Since the function of SSRP1 during cell differentiation remains unclear, we investigated its potential role in controlling lineage determination. Depletion of SSRP1 in human mesenchymal stem cells elicited lineage-specific effects by increasing expression of adipocyte-specific genes and decreasing the expression of osteoblast-specific genes. Consistent with a role in controlling lineage specification, transcriptome-wide RNA-sequencing following SSRP1 depletion and the induction of osteoblast differentiation revealed a specific decrease in the expression of genes involved in biological processes related to osteoblast differentiation. Importantly, we observed a specific downregulation of target genes of the canonical Wnt signaling pathway, which was accompanied by decreased nuclear localization of active ß-catenin. Together our data uncover a previously unknown role for SSRP1 in promoting the activation of the Wnt signaling pathway activity during cellular differentiation. Stem Cells 2016;34:1369-1376.

RSAT 2015: Regulatory Sequence Analysis Tools.

RSAT (Regulatory Sequence Analysis Tools) is a modular software suite for the analysis of cis-regulatory elements in genome sequences. Its main applications are (i) motif discovery, appropriate to genome-wide data sets like ChIP-seq, (ii) transcription factor binding motif analysis (quality assessment, comparisons and clustering), (iii) comparative genomics and (iv) analysis of regulatory variations. Nine new programs have been added to the 43 described in the 2011 NAR Web Software Issue, including a tool to extract sequences from a list of coordinates (fetch-sequences from UCSC), novel programs dedicated to the analysis of regulatory variants from GWAS or population genomics (retrieve-variation-seq and variation-scan), a program to cluster motifs and visualize the similarities as trees (matrix-clustering). To deal with the drastic increase of sequenced genomes, RSAT public sites have been reorganized into taxon-specific servers. The suite is well-documented with tutorials and published protocols. The software suite is available through Web sites, SOAP/WSDL Web services, virtual machines and stand-alone programs at http://www.rsat.eu/.

Shotgun ecotoxicoproteomics of Daphnia pulex: biochemical effects of the anticancer drug tamoxifen.

Among pollutants released into the environment by human activities, residues of pharmaceuticals are an increasing matter of concern because of their potential impact on ecosystems. The aim of this study was to analyze differences of protein expression resulting from acute (2 days) and middle-term (7 days) exposure of aquatic microcrustacean Daphnia pulex to the anticancer drug tamoxifen. Using a liquid chromatography-mass spectrometry shotgun approach, about 4000 proteins could be identified, providing the largest proteomics data set of D. pulex published up to now. Considering both time points and tested concentrations, 189 proteins showed a significant fold change. The identity of regulated proteins suggested a decrease in translation, an increase in protein degradation and changes in carbohydrate and lipid metabolism as the major effects of the drug. Besides these impacted processes, which reflect a general stress response of the organism, some other regulated proteins play a role in Daphnia reproduction. These latter results are in accordance with our previous observations of the impact of tamoxifen on D. pulex reproduction and illustrate the potential of ecotoxicoproteomics to unravel links between xenobiotic effects at the biochemical and organismal levels. Data are available via ProteomeXchange with identifier PXD001257.

Quantitative proteomics identifies the membrane-associated peroxidase GPx8 as a cellular substrate of the hepatitis C virus NS3-4A protease.

UNLABELLED: The hepatitis C virus (HCV) NS3-4A protease is not only an essential component of the viral replication complex and a prime target for antiviral intervention but also a key player in the persistence and pathogenesis of HCV. It cleaves and thereby inactivates two crucial adaptor proteins in viral RNA sensing and innate immunity, mitochondrial antiviral signaling protein (MAVS) and TRIF, a phosphatase involved in growth factor signaling, T-cell protein tyrosine phosphatase (TC-PTP), and the E3 ubiquitin ligase component UV-damaged DNA-binding protein 1 (DDB1). Here we explored quantitative proteomics to identify novel cellular substrates of the NS3-4A protease. Cell lines inducibly expressing the NS3-4A protease were analyzed by stable isotopic labeling using amino acids in cell culture (SILAC) coupled with protein separation and mass spectrometry. This approach identified the membrane-associated peroxidase GPx8 as a bona fide cellular substrate of the HCV NS3-4A protease. Cleavage by NS3-4A occurs at Cys 11, removing the cytosolic tip of GPx8, and was observed in different experimental systems as well as in liver biopsies from patients with chronic HCV. Overexpression and RNA silencing studies revealed that GPx8 is involved in viral particle production but not in HCV entry or RNA replication. CONCLUSION: We provide proof-of-concept for the use of quantitative proteomics to identify cellular substrates of a viral protease and describe GPx8 as a novel proviral host factor targeted by the HCV NS3-4A protease.

Genome-wide analysis of host mRNA translation during hepatitis C virus infection.

In the model of Huh-7.5.1 hepatocyte cells infected by the JFH1 hepatitis C virus (HCV) strain, transcriptomic and proteomic studies have revealed modulations of pathways governing mainly apoptosis and cell cycling. Differences between transcriptomic and proteomic studies pointed to regulations occurring at the posttranscriptional level, including the control of mRNA translation. In this study, we investigated at the genome-wide level the translational regulation occurring during HCV infection. Sucrose gradient ultracentrifugation followed by microarray analysis was used to identify translationally regulated mRNAs (mRNAs associated with ribosomes) from JFH1-infected and uninfected Huh-7.5.1 cells. Translationally regulated mRNAs were found to correspond to genes enriched in specific pathways, including vesicular transport and posttranscriptional regulation. Interestingly, the strongest translational regulation was found for mRNAs encoding proteins involved in pre-mRNA splicing, mRNA translation, and protein folding. Strikingly, these pathways were not previously identified, through transcriptomic studies, as being modulated following HCV infection. Importantly, the observed changes in host mRNA translation were directly due to HCV replication rather than to HCV entry, since they were not observed in JFH1-infected Huh-7.5.1 cells treated with a potent HCV NS3 protease inhibitor. Overall, this study highlights the need to consider, beyond transcriptomic or proteomic studies, the modulation of host mRNA translation as an important aspect of HCV infection.

ExPASy: SIB bioinformatics resource portal.

ExPASy (http://www.expasy.org) has worldwide reputation as one of the main bioinformatics resources for proteomics. It has now evolved, becoming an extensible and integrative portal accessing many scientific resources, databases and software tools in different areas of life sciences. Scientists can henceforth access seamlessly a wide range of resources in many different domains, such as proteomics, genomics, phylogeny/evolution, systems biology, population genetics, transcriptomics, etc. The individual resources (databases, web-based and downloadable software tools) are hosted in a 'decentralized' way by different groups of the SIB Swiss Institute of Bioinformatics and partner institutions. Specifically, a single web portal provides a common entry point to a wide range of resources developed and operated by different SIB groups and external institutions. The portal features a search function across 'selected' resources. Additionally, the availability and usage of resources are monitored. The portal is aimed for both expert users and people who are not familiar with a specific domain in life sciences. The new web interface provides, in particular, visual guidance for newcomers to ExPASy.

A plant-derived flavonoid inhibits entry of all HCV genotypes into human hepatocytes.

BACKGROUND & AIMS: Interferon-based therapies for hepatitis C virus (HCV) infection are limited by side effects and incomplete response rates, particularly among transplant recipients. We screened a library of plant-derived small molecules to identify HCV inhibitors with novel mechanisms. METHODS: We isolated phenolic compounds from Marrubium peregrinum L (Lamiaceae). Replication of HCV RNA, virus production, and cell entry were monitored using replicons and infectious HCV. Inhibition of HCV was measured in hepatoma cells and primary human hepatocytes using luciferase reporter gene assays, core enzyme-linked immunosorbent assays, or infectivity titration. We tested the bioavailability of the compound in mice. RESULTS: We identified a flavonoid, ladanein (BJ486K), with unreported antiviral activity and established its oral bioavailability in mice. Natural and synthetic BJ486K inhibited a post-attachment entry step, but not RNA replication or assembly; its inhibitory concentration 50% was 2.5 µm. BJ486K was effective against all major HCV genotypes, including a variant that is resistant to an entry inhibitor; it prevented infection of primary human hepatocytes. Combined administration of BJ486K and cyclosporine A had a synergistic effect in inhibition of HCV infection. CONCLUSIONS: BJ486K has oral bioavailability and interferes with entry of HCV into cultured human hepatocytes. It synergizes with cyclosporine A to inhibit HCV infection. Its inhibitory effects are independent of HCV genotype, including a variant that is resistant to an entry inhibitor against scavenger receptor class B type I. Flavonoid derivatives therefore might be developed as components of combination therapies because they are potent, broadly active inhibitors of HCV entry that could prevent graft reinfection after liver transplantation.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use.

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

Hepatitis C viral protein NS5A induces EMT and participates in oncogenic transformation of primary hepatocyte precursors.

BACKGROUND & AIMS: Apicobasal polarity, which is essential for epithelial structure and function, is targeted by several tumour-related pathogens and is generally perturbed in the course of carcinogenesis. Hepatitis C virus (HCV) infection is associated with a strong risk of hepatocellular carcinoma, typically preceded by dysplastic alterations of cell morphology. We investigated the molecular mechanisms and the functional consequences of HCV-driven perturbations of epithelial polarity. METHODS: We used biochemical, genetic, and cell biology approaches to assess the impact of hepatitis C viral protein NS5A on the polarity and function of hepatocytes and hepatic progenitors. Transgenic animals and xenograft models served for in vivo validation of the results obtained in cell culture. RESULTS: We found that expression of HCV-NS5A in primary hepatic precursors and in immortalized hepatocyte cell lines gave rise to profound modifications of cell polarity, leading to epithelial to mesenchymal transition (EMT). NS5A, either alone or in the context of the full complement of viral proteins in the course of infection, acted through activating Twist2, a transcriptional regulator of EMT. The effects of NS5A were additive to those of TGF-ß, a cytokine abundant in diseased liver and highly relevant to HCV-related pathology. Moreover, NS5A cooperates with oncogenic Ras, giving rise to transformed, invasive cells that are highly tumorigenic in vivo. CONCLUSIONS: Our data suggest that in the context of HCV infection, NS5A favors formation of preneoplastic lesions by disrupting cell polarity and additional oncogenic events cooperate with the viral protein to give rise to motile and invasive tumour cells.

Hepatitis C virus (HCV) evades NKG2D-dependent NK cell responses through NS5A-mediated imbalance of inflammatory cytokines.

Understanding how hepatitis C virus (HCV) induces and circumvents the host's natural killer (NK) cell-mediated immunity is of critical importance in efforts to design effective therapeutics. We report here the decreased expression of the NKG2D activating receptor as a novel strategy adopted by HCV to evade NK-cell mediated responses. We show that chronic HCV infection is associated with expression of ligands for NKG2D, the MHC class I-related Chain (MIC) molecules, on hepatocytes. However, NKG2D expression is downmodulated on circulating NK cells, and consequently NK cell-mediated cytotoxic capacity and interferon-γ production are impaired. Using an endotoxin-free recombinant NS5A protein, we show that NS5A stimulation of monocytes through Toll-like Receptor 4 (TLR4) promotes p38- and PI3 kinase-dependent IL-10 production, while inhibiting IL-12 production. In turn, IL-10 triggers secretion of TGFß which downmodulates NKG2D expression on NK cells, leading to their impaired effector functions. Moreover, culture supernatants of HCV JFH1 replicating Huh-7.5.1 cells reproduce the effect of recombinant NS5A on NKG2D downmodulation. Exogenous IL-15 can antagonize the TGFß effect and restore normal NKG2D expression on NK cells. We conclude that NKG2D-dependent NK cell functions are modulated during chronic HCV infection, and demonstrate that this alteration can be prevented by exogenous IL-15, which could represent a meaningful adjuvant for therapeutic intervention.

The purinergic P2Y14 receptor links hepatocyte death to hepatic stellate cell activation and fibrogenesis in the liver.

Fibrosis contributes to ~45% of deaths in western countries. In chronic liver disease, fibrosis is a major factor determining outcomes, but efficient antifibrotic therapies are lacking. Although platelet-derived growth factor and transforming growth factor-ß constitute key fibrogenic mediators, they do not account for the well-established link between cell death and fibrosis in the liver. Here, we hypothesized that damage-associated molecular patterns (DAMPs) may link epithelial cell death to fibrogenesis in the injured liver. DAMP receptor screening identified purinergic receptor P2Y14 among several candidates as highly enriched in hepatic stellate cells (HSCs), the main fibrogenic cell type of the liver. Conversely, P2Y14 ligands uridine 5'-diphosphate (UDP)-glucose and UDP-galactose were enriched in hepatocytes and were released upon different modes of cell death. Accordingly, ligand-receptor interaction analysis that combined proteomic and single-cell RNA sequencing data revealed P2Y14 ligands and P2Y14 receptor as a link between dying cells and HSCs, respectively. Treatment with P2Y14 ligands or coculture with dying hepatocytes promoted HSC activation in a P2Y14-dependent manner. P2Y14 ligands activated extracellular signal-regulated kinase (ERK) and Yes-associated protein (YAP) signaling in HSCs, resulting in ERK-dependent HSC activation. Global and HSC-selective P2Y14 deficiency attenuated liver fibrosis in multiple mouse models of liver injury. Functional expression of P2Y14 was confirmed in healthy and diseased human liver and human HSCs. In conclusion, P2Y14 ligands and their receptor constitute a profibrogenic DAMP pathway that directly links cell death to fibrogenesis.

Addressing trypsin bias in large scale (phospho)proteome analysis by size exclusion chromatography and secondary digestion of large post-trypsin peptides.

In the vast majority of bottom-up proteomics studies, protein digestion is performed using only mammalian trypsin. Although it is clearly the best enzyme available, the sole use of trypsin rarely leads to complete sequence coverage, even for abundant proteins. It is commonly assumed that this is because many tryptic peptides are either too short or too long to be identified by RPLC-MS/MS. We show through in silico analysis that 20-30% of the total sequence of three proteomes (Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Homo sapiens) is expected to be covered by Large post-Trypsin Peptides (LpTPs) with M(r) above 3000 Da. We then established size exclusion chromatography to fractionate complex yeast tryptic digests into pools of peptides based on size. We found that secondary digestion of LpTPs followed by LC-MS/MS analysis leads to a significant increase in identified proteins and a 32-50% relative increase in average sequence coverage compared to trypsin digestion alone. Application of the developed strategy to analyze the phosphoproteomes of S. pombe and of a human cell line identified a significant fraction of novel phosphosites. Overall our data indicate that specific targeting of LpTPs can complement standard bottom-up workflows to reveal a largely neglected portion of the proteome.

Production of infectious hepatitis C virus in primary cultures of human adult hepatocytes.

BACKGROUND & AIMS: Although hepatitis C virus (HCV) can be grown in the hepatocarcinoma-derived cell line Huh-7, a cell-culture model is needed that supports its complete, productive infection cycle in normal, quiescent, highly differentiated human hepatocytes. We sought to develop such a system. METHODS: Primary cultures of human adult hepatocytes were inoculated with HCV derived from Huh-7 cell culture (HCVcc) and monitored for expression of hepatocyte differentiation markers and replication of HCV. Culture supernatants were assayed for HCV RNA, core antigen, and infectivity titer. The buoyant densities of input and progeny virus were compared in iodixanol gradients. RESULTS: While retaining expression of differentiation markers, primary hepatocytes supported the complete infectious cycle of HCV, including production of significant titers of new infectious progeny virus, which was called primary-culture-derived virus (HCVpc). Compared with HCVcc, HCVpc had lower average buoyant density and higher specific infectivity; this was similar to the characteristics of virus particles associated with the very-low-density lipoproteins that are produced during in vivo infection. These properties were lost after re-culture of HCVpc in poorly differentiated Huh-7 cells, suggesting that authentic virions can be produced only by normal hepatocytes that secrete authentic very-low-density lipoproteins. CONCLUSIONS: We have established a cell-culture-based system that allows production of infectious HCV in physiologically relevant human hepatocytes. This provides a useful tool for the study of HCV interactions with its natural host cell and for the development of antiviral therapies.

Benchmarking joint multi-omics dimensionality reduction approaches for the study of cancer.

High-dimensional multi-omics data are now standard in biology. They can greatly enhance our understanding of biological systems when effectively integrated. To achieve proper integration, joint Dimensionality Reduction (jDR) methods are among the most efficient approaches. However, several jDR methods are available, urging the need for a comprehensive benchmark with practical guidelines. We perform a systematic evaluation of nine representative jDR methods using three complementary benchmarks. First, we evaluate their performances in retrieving ground-truth sample clustering from simulated multi-omics datasets. Second, we use TCGA cancer data to assess their strengths in predicting survival, clinical annotations and known pathways/biological processes. Finally, we assess their classification of multi-omics single-cell data. From these in-depth comparisons, we observe that intNMF performs best in clustering, while MCIA offers an effective behavior across many contexts. The code developed for this benchmark study is implemented in a Jupyter notebook-multi-omics mix (momix)-to foster reproducibility, and support users and future developers.