CAPTVRED: an automated pipeline for viral tracking and discovery from capture-based metagenomics samples.

Summary: Target Enrichment Sequencing or Capture-based metagenomics has emerged as an approach of interest for viral metagenomics in complex samples. However, these datasets are usually analyzed with standard downstream Bioinformatics analyses. CAPTVRED (Capture-based metagenomics Analysis Pipeline for tracking ViRal species from Environmental Datasets), has been designed to assess the virome present in complex samples, specially focused on those obtained by Target Enrichment Sequencing approach. This work aims to provide a user-friendly tool that complements this sequencing approach for the total or partial virome description, especially from environmental matrices. It includes a setup module which allows preparation and adjustment of the pipeline to any capture panel directed to a set of species of interest. The tool also aims to reduce time and computational cost, as well as to provide comprehensive, reproducible, and accessible results while being easy to costume, set up, and install. Availability and implementation: Source code and test datasets are freely available at github repository: https://github.com/CompGenLabUB/CAPTVRED.git.

Gene regulatory dynamics during the development of a paleopteran insect, the mayfly Cloeon dipterum.

The evolution of insects has been marked by the appearance of key body plan innovations that promoted the outstanding ability of this lineage to adapt to new habitats, boosting the most successful radiation in animals. To understand the evolution of these new structures, it is essential to investigate which genes and gene regulatory networks participate during the embryonic development of insects. Great efforts have been made to fully understand gene expression and gene regulation during the development of holometabolous insects, in particular Drosophila melanogaster. Conversely, functional genomics resources and databases in other insect lineages are scarce. To provide a new platform to study gene regulation in insects, we generated ATAC-seq for the first time during the development of the mayfly Cloeon dipterum, which belongs to Paleoptera, the sister group to all other winged insects. With these comprehensive datasets along six developmental stages, we characterized pronounced changes in accessible chromatin between early and late embryogenesis. The application of ATAC-seq in mayflies provides a fundamental resource to understand the evolution of gene regulation in insects.

Target enrichment metaviromics enables comprehensive surveillance of coronaviruses in environmental and animal samples.

The COVID-19 pandemic has underscored the importance of understanding the role of animals in the transmission of coronaviruses (CoVs) and their impact on human health. A One Health approach, integrating human, animal, and environmental health, is essential for effective CoVs control. Next-generation sequencing has played a pivotal role in identifying and monitoring the evolution of novel CoVs strains, like SARS-CoV-2. However, viral occurrence and diversity studies in environmental and animal samples are challenging because of the complexity of viral communities and low abundance of viruses in these samples. Target enrichment sequencing (TES) has emerged as a valuable tool for investigating viral families in challenging samples. This approach involves the specific capture and enrichment of viral genomes using sequence-specific probes, thereby enhancing the efficiency of detection and characterization. In this study, we aimed to develop and validate a TES panel to study CoVs in various complex environmental and animal derived samples. The results demonstrated the panel's effectiveness in capturing and sequencing a wide diversity of CoVs providing valuable insights into their abundance and host diversity in urban wastewater, farm animal corpses lixiviates and bat guano samples. In sewage samples, CoVs were detected solely when TES was employed while in guano samples, sequencing of CoVs species was achieved in 2 out of 4 samples showing an almost three-logarithmic increase in the number of reads obtained in comparison with the untargeted approach. For animal lixiviates, only the TES application enabled the acquisition of CoVs reads. The information obtained can significantly contribute to early detection, surveillance, and control measures for CoVs, including viral discovery and potential spillover events. Additionally, this sequencing panel shows potential for studying other significant viral families and monitoring viral diversity in different animal populations.

Wnt/ß-catenin signalling is required for pole-specific chromatin remodeling during planarian regeneration.

For successful regeneration, the identity of the missing tissue must be specified according to the pre-existing tissue. Planarians are ideal for the study of the mechanisms underlying this process; the same field of cells can regrow a head or a tail according to the missing body part. After amputation, the differential activation of the Wnt/ß-catenin signal specifies anterior versus posterior identity. Initially, both wnt1 and notum (Wnt inhibitor) are expressed in all wounds, but 48 hours later they are restricted to posterior or anterior facing wounds, respectively, by an unknown mechanism. Here we show that 12 hours after amputation, the chromatin accessibility of cells in the wound region changes according to the polarity of the pre-existing tissue in a Wnt/ß-catenin-dependent manner. Genomic analyses suggest that homeobox transcription factors and chromatin-remodeling proteins are direct Wnt/ß-catenin targets, which trigger the expression of posterior effectors. Finally, we identify FoxG as a wnt1 up-stream regulator, probably via binding to its first intron enhancer region.

Oxidative Stress Is Associated with Overgrowth in Drosophila l(3)mbt Mutant Imaginal Discs.

The loss-of-function conditions for an l(3)malignant brain tumour (l(3)mbt) in larvae reared at 29 °C results in malignant brain tumours and hyperplastic imaginal discs. Unlike the former that have been extensively characterised, little is known about the latter. Here we report the results of a study of the hyperplastic l(3)mbt mutant wing imaginal discs. We identify the l(3)mbt wing disc tumour transcriptome and find it to include genes involved in reactive oxygen species (ROS) metabolism. Furthermore, we show the presence of oxidative stress in l(3)mbt hyperplastic discs, even in apoptosis-blocked conditions, but not in l(3)mbt brain tumours. We also find that chemically blocking oxidative stress in l(3)mbt wing discs reduces the incidence of wing disc overgrowths. Our results reveal the involvement of oxidative stress in l(3)mbt wing discs hyperplastic growth.

Gene Network of Susceptibility to Atypical Femoral Fractures Related to Bisphosphonate Treatment.

Atypical femoral fractures (AFF) are rare fragility fractures in the subtrocantheric or diaphysis femoral region associated with long-term bisphosphonate (BP) treatment. The etiology of AFF is still unclear even though a genetic basis is suggested. We performed whole exome sequencing (WES) analysis of 12 patients receiving BPs for at least 5 years who sustained AFFs and 4 controls, also long-term treated with BPs but without any fracture. After filtration and prioritization of rare variants predicted to be damaging and present in genes shared among at least two patients, a total of 272 variants in 132 genes were identified. Twelve of these genes were known to be involved in bone metabolism and/or AFF, highlighting DAAM2 and LRP5, both involved in the Wnt pathway, as the most representative. Afterwards, we intersected all mutated genes with a list of 34 genes obtained from a previous study of three sisters with BP-related AFF, identifying nine genes. One of these (MEX3D) harbored damaging variants in two AFF patients from the present study and one shared among the three sisters. Gene interaction analysis using the AFFNET web suggested a complex network among bone-related genes as well as with other mutated genes. BinGO biological function analysis highlighted cytoskeleton and cilium organization. In conclusion, several genes and their interactions could provide genetic susceptibility to AFF, that along with BPs treatment and in some cases with glucocorticoids may trigger this so feared complication.

Exploring the diversity of coronavirus in sewage during COVID-19 pandemic: Don't miss the forest for the trees.

In the wake of the COVID-19 pandemic, the use of next generation sequencing (NGS) has proved to be an important tool for the genetic characterization of SARS-CoV-2 from clinical samples. The use of different available NGS tools applied to wastewater samples could be the key for an in-depth study of the excreted virome, not only focusing on SARS-CoV-2 circulation and typing, but also to detect other potentially pandemic viruses within the same family. With this aim, 24-hours composite wastewater samples from March and July 2020 were sequenced by applying specific viral NGS as well as target enrichment NGS. The full virome of the analyzed samples was obtained, with human Coronaviridae members (CoV) present in one of those samples after applying the enrichment. One contig was identified as HCoV-OC43 and 8 contigs as SARS-CoV-2. CoVs from other animal hosts were also detected when applying this technique. These contigs were compared with those obtained from contemporary clinical specimens by applying the same target enrichment approach. The results showed that there is a co-circulation in urban areas of human and animal coronaviruses infecting domestic animals and rodents. NGS enrichment-based protocols might be crucial to describe the occurrence and genetic characteristics of SARS-CoV-2 and other Coronaviridae family members within the excreted virome present in wastewater.

The Deubiquitinating Enzyme Ataxin-3 Regulates Ciliogenesis and Phagocytosis in the Retina.

Expansion of a CAG repeat in ATXN3 causes the dominant polyglutamine disease spinocerebellar ataxia type 3 (SCA3), yet the physiological role of ATXN3 remains unclear. Here, we focus on unveiling the function of Ataxin-3 (ATXN3) in the retina, a neurological organ amenable to morphological and physiological studies. Depletion of Atxn3 in zebrafish and mice causes morphological and functional retinal alterations and, more precisely, photoreceptor cilium and outer segment elongation, cone opsin mislocalization, and cone hyperexcitation. ATXN3 localizes at the basal body and axoneme of the cilium, supporting its role in regulating ciliary length. Abrogation of Atxn3 expression causes decreased levels of the regulatory protein KEAP1 in the retina and delayed phagosome maturation in the retinal pigment epithelium. We propose that ATXN3 regulates two relevant biological processes in the retina, namely, ciliogenesis and phagocytosis, by modulating microtubule polymerization and microtubule-dependent retrograde transport, thus positing ATXN3 as a causative or modifier gene in retinal/macular dystrophies.

Characterisation of the sewage virome: comparison of NGS tools and occurrence of significant pathogens.

NGS techniques are excellent tools to monitor and identify viral pathogens circulating among the population with some limitations that need to be overcome, especially in complex matrices. Sewage contains a high amount of other microorganisms that could interfere when trying to sequence viruses for which random PCR amplifications are needed before NGS. The selection of appropriate NGS tools is important for reliable identification of viral diversity among the population. We have compared different NGS methodologies (Untargeted Viral Metagenomics, Target Enrichment Sequencing and Amplicon Deep Sequencing) for the detection and characterisation of viruses in urban sewage, focusing on three important human pathogens: papillomaviruses, adenoviruses and enteroviruses. A full picture of excreted viruses was obtained by applying Untargeted Viral Metagenomics, which detected members of four different vertebrate viral families in addition to bacteriophages, plant viruses and viruses infecting other hosts. Target Enrichment Sequencing, using specific vertebrate viral probes, allowed the detection of up to eight families containing human viruses, with high variety of types within the families and with a high genome coverage. By applying Amplicon Deep Sequencing, the diversity of enteroviruses, adenoviruses and papillomaviruses observed was higher than when applying the other two strategies and this technique allowed the subtyping of an enterovirus A71 C1 strain related to a brainstem encephalitis outbreak occurring at the same time in the sampling area. From the data obtained, we concluded that the different strategies studied provided different levels of analysis: TES is the best strategy to obtain a broad picture of human viruses present in complex samples such as sewage. Other NGS strategies are useful for studying the virome of complex samples when also targeting viruses infecting plants, bacteria, invertebrates or fungi (Untargeted Viral Metagenomics) or when observing the variety within a sole viral family is the objective of the study (Amplicon Deep Sequencing).

RPGeNet v2.0: expanding the universe of retinal disease gene interactions network.

RPGeNet offers researchers a user-friendly queriable tool to visualize the interactome network of visual disorder genes, thus enabling the identification of new potential causative genes and the assignment of novel candidates to specific retinal or cellular pathways. This can be highly relevant for clinical applications as retinal dystrophies affect 1:3000 people worldwide, and the causative genes are still unknown for 30% of the patients. RPGeNet is a refined interaction network interface that limits its skeleton network to the shortest paths between each and every known causative gene of inherited syndromic and non-syndromic retinal dystrophies. RPGeNet integrates interaction information from STRING, BioGRID and PPaxe, along with retina-specific expression data and associated genetic variants, over a Cytoscape.js web interface. For the new version, RPGeNet v2.0, the database engine was migrated to Neo4j graph database manager, which speeds up the initial queries and can handle whole interactome data for new ways to query the network. Further, user facilities have been introduced as the capability of saving and restoring a researcher customized network layout or as novel features to facilitate navigation and data projection on the network explorer interface. Responsiveness has been further improved by transferring some functionality to the client side.

Ask1 and Akt act synergistically to promote ROS-dependent regeneration in Drosophila.

How cells communicate to initiate a regenerative response after damage has captivated scientists during the last few decades. It is known that one of the main signals emanating from injured cells is the Reactive Oxygen Species (ROS), which propagate to the surrounding tissue to trigger the replacement of the missing cells. However, the link between ROS production and the activation of regenerative signaling pathways is not yet fully understood. We describe here the non-autonomous ROS sensing mechanism by which living cells launch their regenerative program. To this aim, we used Drosophila imaginal discs as a model system due to its well-characterized regenerative ability after injury or cell death. We genetically-induced cell death and found that the Apoptosis signal-regulating kinase 1 (Ask1) is essential for regenerative growth. Ask1 senses ROS both in dying and living cells, but its activation is selectively attenuated in living cells by Akt1, the core kinase component of the insulin/insulin-like growth factor pathway. Akt1 phosphorylates Ask1 in a secondary site outside the kinase domain, which attenuates its activity. This modulation of Ask1 activity results in moderate levels of JNK signaling in the living tissue, as well as in activation of p38 signaling, both pathways required to turn on the regenerative response. Our findings demonstrate a non-autonomous activation of a ROS sensing mechanism by Ask1 and Akt1 to replace the missing tissue after damage. Collectively, these results provide the basis for understanding the molecular mechanism of communication between dying and living cells that triggers regeneration.

SiNoPsis: Single Nucleotide Polymorphisms selection and promoter profiling.

MOTIVATION: The selection of a single nucleotide polymorphism (SNP) using bibliographic methods can be a very time-consuming task. Moreover, a SNP selected in this way may not be easily visualized in its genomic context by a standard user hoping to correlate it with other valuable information. Here we propose a web form built on top of Circos that can assist SNP-centered screening, based on their location in the genome and the regulatory modules they can disrupt. Its use may allow researchers to prioritize SNPs in genotyping and disease studies. RESULTS: SiNoPsis is bundled as a web portal. It focuses on the different structures involved in the genomic expression of a gene, especially those found in the core promoter upstream region. These structures include transcription factor binding sites (for promoter and enhancer signals), histones and promoter flanking regions. Additionally, the tool provides eQTL and linkage disequilibrium (LD) properties for a given SNP query, yielding further clues about other indirectly associated SNPs. Possible disruptions of the aforementioned structures affecting gene transcription are reported using multiple resource databases. SiNoPsis has a simple user-friendly interface, which allows single queries by gene symbol, genomic coordinates, Ensembl gene identifiers, RefSeq transcript identifiers and SNPs. It is the only portal providing useful SNP selection based on regulatory modules and LD with functional variants in both textual and graphic modes (by properly defining the arguments and parameters needed to run Circos). AVAILABILITY AND IMPLEMENTATION: SiNoPsis is freely available at https://compgen.bio.ub.edu/SiNoPsis/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evaluation of Methods for the Concentration and Extraction of Viruses from Sewage in the Context of Metagenomic Sequencing.

Viral sewage metagenomics is a novel field of study used for surveillance, epidemiological studies, and evaluation of waste water treatment efficiency. In raw sewage human waste is mixed with household, industrial and drainage water, and virus particles are, therefore, only found in low concentrations. This necessitates a step of sample concentration to allow for sensitive virus detection. Additionally, viruses harbor a large diversity of both surface and genome structures, which makes universal viral genomic extraction difficult. Current studies have tackled these challenges in many different ways employing a wide range of viral concentration and extraction procedures. However, there is limited knowledge of the efficacy and inherent biases associated with these methods in respect to viral sewage metagenomics, hampering the development of this field. By the use of next generation sequencing this study aimed to evaluate the efficiency of four commonly applied viral concentrations techniques (precipitation with polyethylene glycol, organic flocculation with skim milk, monolithic adsorption filtration and glass wool filtration) and extraction methods (Nucleospin RNA XS, QIAamp Viral RNA Mini Kit, NucliSENS® miniMAG®, or PowerViral® Environmental RNA/DNA Isolation Kit) to determine the viriome in a sewage sample. We found a significant influence of concentration and extraction protocols on the detected viriome. The viral richness was largest in samples extracted with QIAamp Viral RNA Mini Kit or PowerViral® Environmental RNA/DNA Isolation Kit. Highest viral specificity were found in samples concentrated by precipitation with polyethylene glycol or extracted with Nucleospin RNA XS. Detection of viral pathogens depended on the method used. These results contribute to the understanding of method associated biases, within the field of viral sewage metagenomics, making evaluation of the current literature easier and helping with the design of future studies.

The coelomic epithelium transcriptome from a clonal sea star, Coscinasterias muricata.

Coscinasterias is a cosmopolitan genus of large asteroid sea stars with the ability of somatic fission as a clonal reproductive strategy. During fission, the animals tear themselves apart across their central disc, where the lost body parts are regenerated afterwards. Here, we have sequenced and subsequently analysed the transcriptome of the coelomic epithelium of a clonal Coscinasterias muricata specimen from New Zealand. Out of the total 389,768 raw reads, 11,344 contigs were assembled and grouped into functions. Raw read and assembled contig sequences are available at NCBI (BioSample: SAMN03371637), while the annotated assembly can be accessed through the project transcriptome browser (compgen.bio.ub.edu/gbrowse/starfish_transcriptome/). Our data is valuable for future detailed exploration of the coelomic epithelium functions as well as for a better understanding of sea star physiology.

Distilling a Visual Network of Retinitis Pigmentosa Gene-Protein Interactions to Uncover New Disease Candidates.

BACKGROUND: Retinitis pigmentosa (RP) is a highly heterogeneous genetic visual disorder with more than 70 known causative genes, some of them shared with other non-syndromic retinal dystrophies (e.g. Leber congenital amaurosis, LCA). The identification of RP genes has increased steadily during the last decade, and the 30% of the cases that still remain unassigned will soon decrease after the advent of exome/genome sequencing. A considerable amount of genetic and functional data on single RD genes and mutations has been gathered, but a comprehensive view of the RP genes and their interacting partners is still very fragmentary. This is the main gap that needs to be filled in order to understand how mutations relate to progressive blinding disorders and devise effective therapies. METHODOLOGY: We have built an RP-specific network (RPGeNet) by merging data from different sources: high-throughput data from BioGRID and STRING databases, manually curated data for interactions retrieved from iHOP, as well as interactions filtered out by syntactical parsing from up-to-date abstracts and full-text papers related to the RP research field. The paths emerging when known RP genes were used as baits over the whole interactome have been analysed, and the minimal number of connections among the RP genes and their close neighbors were distilled in order to simplify the search space. CONCLUSIONS: In contrast to the analysis of single isolated genes, finding the networks linking disease genes renders powerful etiopathological insights. We here provide an interactive interface, RPGeNet, for the molecular biologist to explore the network centered on the non-syndromic and syndromic RP and LCA causative genes. By integrating tissue-specific expression levels and phenotypic data on top of that network, a more comprehensive biological view will highlight key molecular players of retinal degeneration and unveil new RP disease candidates.

Digital gene expression approach over multiple RNA-Seq data sets to detect neoblast transcriptional changes in Schmidtea mediterranea.

BACKGROUND: The freshwater planarian Schmidtea mediterranea is recognised as a valuable model for research into adult stem cells and regeneration. With the advent of the high-throughput sequencing technologies, it has become feasible to undertake detailed transcriptional analysis of its unique stem cell population, the neoblasts. Nonetheless, a reliable reference for this type of studies is still lacking. RESULTS: Taking advantage of digital gene expression (DGE) sequencing technology we compare all the available transcriptomes for S. mediterranea and improve their annotation. These results are accessible via web for the community of researchers. Using the quantitative nature of DGE, we describe the transcriptional profile of neoblasts and present 42 new neoblast genes, including several cancer-related genes and transcription factors. Furthermore, we describe in detail the Smed-meis-like gene and the three Nuclear Factor Y subunits Smed-nf-YA, Smed-nf-YB-2 and Smed-nf-YC. CONCLUSIONS: DGE is a valuable tool for gene discovery, quantification and annotation. The application of DGE in S. mediterranea confirms the planarian stem cells or neoblasts as a complex population of pluripotent and multipotent cells regulated by a mixture of transcription factors and cancer-related genes.

The nervous system of Xenacoelomorpha: a genomic perspective.

Xenacoelomorpha is, most probably, a monophyletic group that includes three clades: Acoela, Nemertodermatida and Xenoturbellida. The group still has contentious phylogenetic affinities; though most authors place it as the sister group of the remaining bilaterians, some would include it as a fourth phylum within the Deuterostomia. Over the past few years, our group, along with others, has undertaken a systematic study of the microscopic anatomy of these worms; our main aim is to understand the structure and development of the nervous system. This research plan has been aided by the use of molecular/developmental tools, the most important of which has been the sequencing of the complete genomes and transcriptomes of different members of the three clades. The data obtained has been used to analyse the evolutionary history of gene families and to study their expression patterns during development, in both space and time. A major focus of our research is the origin of 'cephalized' (centralized) nervous systems. How complex brains are assembled from simpler neuronal arrays has been a matter of intense debate for at least 100 years. We are now tackling this issue using Xenacoelomorpha models. These represent an ideal system for this work because the members of the three clades have nervous systems with different degrees of cephalization; from the relatively simple sub-epithelial net of Xenoturbella to the compact brain of acoels. How this process of 'progressive' cephalization is reflected in the genomes or transcriptomes of these three groups of animals is the subject of this paper.

Assessment of transcript reconstruction methods for RNA-seq.

We evaluated 25 protocol variants of 14 independent computational methods for exon identification, transcript reconstruction and expression-level quantification from RNA-seq data. Our results show that most algorithms are able to identify discrete transcript components with high success rates but that assembly of complete isoform structures poses a major challenge even when all constituent elements are identified. Expression-level estimates also varied widely across methods, even when based on similar transcript models. Consequently, the complexity of higher eukaryotic genomes imposes severe limitations on transcript recall and splice product discrimination that are likely to remain limiting factors for the analysis of current-generation RNA-seq data.