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.

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.

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.

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.

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.

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.

Complex selection on 5' splice sites in intron-rich organisms.

In contrast to the typically streamlined genomes of prokaryotes, many eukaryotic genomes are riddled with long intergenic regions, spliceosomal introns, and repetitive elements. What explains the persistence of these and other seemingly suboptimal structures? There are three general hypotheses: (1) the structures in question are not actually suboptimal but optimal, being favored by selection, for unknown reasons; (2) the structures are not suboptimal, but of (essentially) equal fitness to "optimal" ones; or (3) the structures are truly suboptimal, but selection is too weak to systematically eliminate them. The 5' splice sites of introns offer a rare opportunity to directly test these hypotheses. Intron-poor species show a clear consensus splice site; most introns begin with the same six nucleotide sequence (typically GTAAGT or GTATGT), indicating efficient selection for this consensus sequence. In contrast, intron-rich species have much less pronounced boundary consensus sequences, and only small minorities of introns in intron-rich species share the same boundary sequence. We studied rates of evolutionary change of 5' splice sites in three groups of closely related intron-rich species--three primates, five Drosophila species, and four Cryptococcus fungi. Surprisingly, the results indicate that changes from consensus-to-variant nucleotides are generally disfavored by selection, but that changes from variant to consensus are neither favored nor disfavored. This evolutionary pattern is consistent with selective differences across introns, for instance, due to compensatory changes at other sites within the gene, which compensate for the otherwise suboptimal consensus-to-variant changes in splice boundaries.

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.

A proteomics approach to decipher the molecular nature of planarian stem cells.

BACKGROUND: In recent years, planaria have emerged as an important model system for research into stem cells and regeneration. Attention is focused on their unique stem cells, the neoblasts, which can differentiate into any cell type present in the adult organism. Sequencing of the Schmidtea mediterranea genome and some expressed sequence tag projects have generated extensive data on the genetic profile of these cells. However, little information is available on their protein dynamics. RESULTS: We developed a proteomic strategy to identify neoblast-specific proteins. Here we describe the method and discuss the results in comparison to the genomic high-throughput analyses carried out in planaria and to proteomic studies using other stem cell systems. We also show functional data for some of the candidate genes selected in our proteomic approach. CONCLUSIONS: We have developed an accurate and reliable mass-spectra-based proteomics approach to complement previous genomic studies and to further achieve a more accurate understanding and description of the molecular and cellular processes related to the neoblasts.

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.

Smed454 dataset: unravelling the transcriptome of Schmidtea mediterranea.

BACKGROUND: Freshwater planarians are an attractive model for regeneration and stem cell research and have become a promising tool in the field of regenerative medicine. With the availability of a sequenced planarian genome, the recent application of modern genetic and high-throughput tools has resulted in revitalized interest in these animals, long known for their amazing regenerative capabilities, which enable them to regrow even a new head after decapitation. However, a detailed description of the planarian transcriptome is essential for future investigation into regenerative processes using planarians as a model system. RESULTS: In order to complement and improve existing gene annotations, we used a 454 pyrosequencing approach to analyze the transcriptome of the planarian species Schmidtea mediterranea Altogether, 598,435 454-sequencing reads, with an average length of 327 bp, were assembled together with the ~10,000 sequences of the S. mediterranea UniGene set using different similarity cutoffs. The assembly was then mapped onto the current genome data. Remarkably, our Smed454 dataset contains more than 3 million novel transcribed nucleotides sequenced for the first time. A descriptive analysis of planarian splice sites was conducted on those Smed454 contigs that mapped univocally to the current genome assembly. Sequence analysis allowed us to identify genes encoding putative proteins with defined structural properties, such as transmembrane domains. Moreover, we annotated the Smed454 dataset using Gene Ontology, and identified putative homologues of several gene families that may play a key role during regeneration, such as neurotransmitter and hormone receptors, homeobox-containing genes, and genes related to eye function. CONCLUSIONS: We report the first planarian transcript dataset, Smed454, as an open resource tool that can be accessed via a web interface. Smed454 contains significant novel sequence information about most expressed genes of S. mediterranea. Analysis of the annotated data promises to contribute to identification of gene families poorly characterized at a functional level. The Smed454 transcriptome data will assist in the molecular characterization of S. mediterranea as a model organism, which will be useful to a broad scientific community.

The diploid genome sequence of an individual human.

Presented here is a genome sequence of an individual human. It was produced from approximately 32 million random DNA fragments, sequenced by Sanger dideoxy technology and assembled into 4,528 scaffolds, comprising 2,810 million bases (Mb) of contiguous sequence with approximately 7.5-fold coverage for any given region. We developed a modified version of the Celera assembler to facilitate the identification and comparison of alternate alleles within this individual diploid genome. Comparison of this genome and the National Center for Biotechnology Information human reference assembly revealed more than 4.1 million DNA variants, encompassing 12.3 Mb. These variants (of which 1,288,319 were novel) included 3,213,401 single nucleotide polymorphisms (SNPs), 53,823 block substitutions (2-206 bp), 292,102 heterozygous insertion/deletion events (indels)(1-571 bp), 559,473 homozygous indels (1-82,711 bp), 90 inversions, as well as numerous segmental duplications and copy number variation regions. Non-SNP DNA variation accounts for 22% of all events identified in the donor, however they involve 74% of all variant bases. This suggests an important role for non-SNP genetic alterations in defining the diploid genome structure. Moreover, 44% of genes were heterozygous for one or more variants. Using a novel haplotype assembly strategy, we were able to span 1.5 Gb of genome sequence in segments >200 kb, providing further precision to the diploid nature of the genome. These data depict a definitive molecular portrait of a diploid human genome that provides a starting point for future genome comparisons and enables an era of individualized genomic information.

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).

Computational gene annotation in new genome assemblies using GeneID.

The sequence of many eukaryotic genomes is nowadays available from a personal computer to any researcher in the world-wide scientific community. However, the sequences are worthless without the adequate annotation of the biological meaningful elements. The annotation of the genes, in particular, is a challenging task that can not be tackled without the aid of specific bioinformatics tools. We present in this chapter a simple protocol mainly based on the combination of the program GeneID and other computational tools to annotate the location of a gene, which was previously annotated in D. melanogaster, in the recently assembled genome of D. yakuba.