Search

Author(s): The FANTOM Consortium and the RIKEN Genome Exploration Research Group Phase I & II Team*; FANTOM Consortium:; Y. Okazaki [1, 2]; M. Furuno [1]; T. Kasukawa [1, 3]; J. [...]

In order to investigate the relationships between Italian wild boar and major pig breeds, we studied the genetic variability of four wild boar populations in Italy (Arezzo, Pisa, Parma, Bergamo) using a 533-bp fragment of the mitochondrial control region. Sixty-nine wild boar samples were analysed, allowing the identification of 10 distinct haplotypes, which involve a total of 15 single nucleotide polymorphisms. Phylogenetic and network analyses were performed also considering several sequences of wild and domesticated forms available in the databases. The Bayesian phylogenetic tree and the Median-Joining network analyses show three main groups: the Italian (IT), European (EU) and Asian (AS) clades. The IT clade corresponds to the Maremma endemic wild boar population and also includes Sardinian individuals, while the EU and AS groups include wild boars as well as domestic pig breeds. Only two individuals from Pisa cluster in the IT group, whereas two haplotypes from Bergamo cluster in the AS group and all other samples cluster in the EU clade. These findings suggest that in Italy wild boar populations have a mixed origin, both EU and AS, and that an interbreeding between wild and domesticated strains has probably occurred. Eight of the 10 wild boars coming from the Migliarino-San Rossore-Massaciuccoli Regional Park (Pisa) belong to H2 and H3 haplotypes, and cluster into the EU clade, suggesting that this regional park is not anymore exclusive of the endemic Maremma wild boar. [ABSTRACT FROM AUTHOR]

The base composition of 25 complete mammalian mitochondrial (mt) genomes has been analyzed taking into account all three codon positions (P1230 and fourfold degenerate sites (P4FD) of H-strand genes. In the nontranscribed L strand, G is the less represented base and A is the most represented one in all cases, while C and T differ among species. H-strand protein-coding genes show an asymmetric distribution of the four bases between the two strands. The asymmetry indexes AT and GC skews on P4FD are much higher than those on P123, suggesting the existence of asymmetrical directional mutation pressure. Relationships between the compositional features and transcription of replication processes have been investigated in order to find a possible mechanism that could explain the origin of this asymmetry. AT and GC skews, the base composition in fourfold degenerate sites, and the number of variable sites for each gene are significantly correlated with the duration of single-stranded state of the H-stranded genes during replication. We tested different replication-related hypotheses, such as the existence of biased dNTP pools, gamma DNA polymerase mispairing, and the asymmetric replication itself. Most of them failed to explain the observed results, hydrolytic deaminations being the only one in agreement with our data. Thus, we hypothesize that one of the crucial processes for the origin of asymmetric and biased base composition of mammalian mitochondrial genomes is the spontaneous deamination of C and A in the H strand during replication. [ABSTRACT FROM PUBLISHER]

Abstract Background Recent discoveries have highlighted the fact that alternative splicing and alternative transcripts are the rule, rather than the exception, in metazoan genes. Since multiple transcript and protein variants expressed by the same gene are, by definition, structurally distinct and need not to be functionally equivalent, the concept of gene orthology should be extended to the transcript level in order to describe evolutionary relationships between structurally similar transcript variants. In other words, the identification of true orthology relationships between gene products now should progress beyond primary sequence and "splicing orthology", consisting in ancestrally shared exon-intron structures, is required to define orthologous isoforms at transcript level. Results As a starting step in this direction, in this work we performed a large scale human- mouse gene comparison with a twofold goal: first, to assess if and to which extent traditional gene annotations such as RefSeq capture genuine splicing orthology; second, to provide a more detailed annotation and quantification of true human-mouse orthologous transcripts defined as transcripts of orthologous genes exhibiting the same splicing patterns. Conclusions We observed an identical exon/intron structure for 32% of human and mouse orthologous genes. This figure increases to 87% using less stringent criteria for gene structure similarity, thus implying that for about 13% of the human RefSeq annotated genes (and about 25% of the corresponding transcripts) we could not identify any mouse transcript showing sufficient similarity to be confidently assigned as a splicing ortholog. Our data suggest that current gene and transcript data may still be rather incomplete - with several splicing variants still unknown. The observation that alternative splicing produces large numbers of alternative transcripts and proteins, some of them conserved across species and others truly species-specific, suggests that, still maintaining the conventional definition of gene orthology, a new concept of "splicing orthology" can be defined at transcript level.

Abstract Background Within Chordata, the subphyla Vertebrata and Cephalochordata (lancelets) are characterized by a remarkable stability of the mitochondrial (mt) genome, with constancy of gene content and almost invariant gene order, whereas the limited mitochondrial data on the subphylum Tunicata suggest frequent and extensive gene rearrangements, observed also within ascidians of the same genus. Results To confirm this evolutionary trend and to better understand the evolutionary dynamics of the mitochondrial genome in Tunicata Ascidiacea, we have sequenced and characterized the complete mt genome of two congeneric ascidian species, Phallusia mammillata and Phallusia fumigata (Phlebobranchiata, Ascidiidae). The two mtDNAs are surprisingly rearranged, both with respect to one another and relative to those of other tunicates and chordates, with gene rearrangements affecting both protein-coding and tRNA genes. The new data highlight the extraordinary variability of ascidian mt genome in base composition, tRNA secondary structure, tRNA gene content, and non-coding regions (number, size, sequence and location). Indeed, both Phallusia genomes lack the trnD gene, show loss/acquisition of DHU-arm in two tRNAs, and have a G+C content two-fold higher than other ascidians. Moreover, the mt genome of P. fumigata presents two identical copies of trnI, an extra tRNA gene with uncertain amino acid specificity, and four almost identical sequence regions. In addition, a truncated cytochrome b, lacking a C-terminal tail that commonly protrudes into the mt matrix, has been identified as a new mt feature probably shared by all tunicates. Conclusion The frequent occurrence of major gene order rearrangements in ascidians both at high taxonomic level and within the same genus makes this taxon an excellent model to study the mechanisms of gene rearrangement, and renders the mt genome an invaluable phylogenetic marker to investigate molecular biodiversity and speciation events in this largely unexplored group of basal chordates.

Abstract Background Phylogenetic relationships between Lagomorpha, Rodentia and Primates and their allies (Euarchontoglires) have long been debated. While it is now generally agreed that Rodentia constitutes a monophyletic sister-group of Lagomorpha and that this clade (Glires) is sister to Primates and Dermoptera, higher-level relationships within Rodentia remain contentious. Results We have sequenced and performed extensive evolutionary analyses on the mitochondrial genome of the scaly-tailed flying squirrel Anomalurus sp., an enigmatic rodent whose phylogenetic affinities have been obscure and extensively debated. Our phylogenetic analyses of the coding regions of available complete mitochondrial genome sequences from Euarchontoglires suggest that Anomalurus is a sister taxon to the Hystricognathi, and that this clade represents the most basal divergence among sampled Rodentia. Bayesian dating methods incorporating a relaxed molecular clock provide divergence-time estimates which are consistently in agreement with the fossil record and which indicate a rapid radiation within Glires around 60 million years ago. Conclusion Taken together, the data presented provide a working hypothesis as to the phylogenetic placement of Anomalurus, underline the utility of mitochondrial sequences in the resolution of even relatively deep divergences and go some way to explaining the difficulty of conclusively resolving higher-level relationships within Glires with available data and methodologies.

Abstract Background To gain insight into the evolutionary features of the huntingtin (htt) gene in Chordata, we have sequenced and characterized the full-length htt mRNA in the ascidian Ciona intestinalis, a basal chordate emerging as new invertebrate model organism. Moreover, taking advantage of the availability of genomic and EST sequences, the htt gene structure of a number of chordate species, including the cogeneric ascidian Ciona savignyi, and the vertebrates Xenopus and Gallus was reconstructed. Results The C. intestinalis htt transcript exhibits some peculiar features, such as spliced leader trans-splicing in the 98 nt-long 5' untranslated region (UTR), an alternative splicing in the coding region, eight alternative polyadenylation sites, and no similarities of both 5' and 3'UTRs compared to homologs of the cogeneric C. savignyi. The predicted protein is 2946 amino acids long, shorter than its vertebrate homologs, and lacks the polyQ and the polyP stretches found in the the N-terminal regions of mammalian homologs. The exon-intron organization of the htt gene is almost identical among vertebrates, and significantly conserved between Ciona and vertebrates, allowing us to hypothesize an ancestral chordate gene consisting of at least 40 coding exons. Conclusion During chordate diversification, events of gain/loss, sliding, phase changes, and expansion of introns occurred in both vertebrate and ascidian lineages predominantly in the 5'-half of the htt gene, where there is also evidence of lineage-specific evolutionary dynamics in vertebrates. On the contrary, the 3'-half of the gene is highly conserved in all chordates at the level of both gene structure and protein sequence. Between the two Ciona species, a fast evolutionary rate and/or an early divergence time is suggested by the absence of significant similarity between UTRs, protein divergence comparable to that observed between mammals and fishes, and different distribution of repetitive elements.

Cyclophilin (CyP) D is a regulator of the mitochondrial F-ATP synthase. Here we report the discovery of a form of CyPD lacking the first 10 (mouse) or 13 (human) N-terminal residues (ΔN-CyPD), a protein region with species-specific features. NMR studies on recombinant human full-length CyPD (FL-CyPD) and ΔN-CyPD form revealed that the N-terminus is highly flexible, in contrast with the rigid globular part. We have studied the interactions of FL and ΔN-CyPD with F-ATP synthase at the OSCP subunit, a site where CyPD binding inhibits catalysis and favors the transition of the enzyme complex to the permeability transition pore. At variance from FL-CyPD, ΔN-CyPD binds OSCP in saline media, indicating that the N-terminus substantially decreases the binding affinity for OSCP. We also provide evidence that calpain 1 is responsible for generation of ΔN-CyPD in cells. Altogether, our work suggests the existence of a novel mechanism of modulation of CyPD through cleavage of its N-terminus that may have significant pathophysiological implications., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

A genomic database of all Earth's eukaryotic species could contribute to many scientific discoveries; however, only a tiny fraction of species have genomic information available. In 2018, scientists across the world united under the Earth BioGenome Project (EBP), aiming to produce a database of high-quality reference genomes containing all ~1.5 million recognized eukaryotic species. As the European node of the EBP, the European Reference Genome Atlas (ERGA) sought to implement a new decentralised, equitable and inclusive model for producing reference genomes. For this, ERGA launched a Pilot Project establishing the first distributed reference genome production infrastructure and testing it on 98 eukaryotic species from 33 European countries. Here we outline the infrastructure and explore its effectiveness for scaling high-quality reference genome production, whilst considering equity and inclusion. The outcomes and lessons learned provide a solid foundation for ERGA while offering key learnings to other transnational, national genomic resource projects and the EBP., (© 2024. The Author(s).)

Sequencing the mitochondrial genome of the tunicate Oikopleura dioica is a challenging task due to the presence of long poly-A/T homopolymer stretches, which impair sequencing and assembly. Here, we report on the sequencing and annotation of the majority of the mitochondrial genome of O. dioica by means of combining several DNA and amplicon reads obtained by Illumina and MinIon Oxford Nanopore Technologies with public RNA sequences. We document extensive RNA editing, since all homopolymer stretches present in the mitochondrial DNA correspond to 6U-regions in the mitochondrial RNA. Out of the 13 canonical protein-coding genes, we were able to detect eight, plus an unassigned open reading frame that lacked sequence similarity to canonical mitochondrial protein-coding genes. We show that the nad3 gene has been transferred to the nucleus and acquired a mitochondria-targeting signal. In addition to two very short rRNAs, we could only identify a single tRNA (tRNA-Met), suggesting multiple losses of tRNA genes, supported by a corresponding loss of mitochondrial aminoacyl-tRNA synthetases in the nuclear genome. Based on the eight canonical protein-coding genes identified, we reconstructed maximum likelihood and Bayesian phylogenetic trees and inferred an extreme evolutionary rate of this mitochondrial genome. The phylogenetic position of appendicularians among tunicates, however, could not be accurately determined., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Ribonucleotides represent the most common non-canonical nucleotides found in eukaryotic genomes. The sources of chromosome-embedded ribonucleotides and the mechanisms by which unrepaired rNMPs trigger genome instability and human pathologies are not fully understood. The available sequencing technologies only allow to indirectly deduce the genomic location of rNMPs. Oxford Nanopore Technologies (ONT) may overcome such limitation, revealing the sites of rNMPs incorporation in genomic DNA directly from raw sequencing signals. We synthesized two types of DNA molecules containing rNMPs at known or random positions and we developed data analysis pipelines for DNA-embedded ribonucleotides detection by ONT. We report that ONT can identify all four ribonucleotides incorporated in DNA by capturing rNMPs-specific alterations in nucleotide alignment features, current intensity, and dwell time. We propose that ONT may be successfully employed to directly map rNMPs in genomic DNA and we suggest a strategy to build an ad hoc basecaller to analyse native genomes., (© 2024. The Author(s).)

The 16S rRNA amplicon-based sequencing approach represents the most common and cost-effective strategy with great potential for microbiome profiling. The use of second-generation sequencing (NGS) technologies has led to protocols based on the amplification of one or a few hypervariable regions, impacting the outcome of the analysis. Nowadays, comparative studies are necessary to assess different amplicon-based approaches, including the full-locus sequencing currently feasible thanks to third-generation sequencing (TGS) technologies. This study compared three different methods to achieve the deepest microbiome taxonomic characterization: (a) the single-region approach, (b) the multiplex approach, covering several regions of the target gene/region, both based on NGS short reads, and (c) the full-length approach, which analyzes the whole length of the target gene thanks to TGS long reads. Analyses carried out on benchmark microbiome samples, with a known taxonomic composition, highlighted a different classification performance, strongly associated with the type of hypervariable regions and the coverage of the target gene. Indeed, the full-length approach showed the greatest discriminating power, up to species level, also on complex real samples. This study supports the transition from NGS to TGS for the study of the microbiome, even if experimental and bioinformatic improvements are still necessary.

Accurate and timely monitoring of the evolution of SARS-CoV-2 is crucial for identifying and tracking potentially more transmissible/virulent viral variants, and implement mitigation strategies to limit their spread. Here we introduce HaploCoV, a novel software framework that enables the exploration of SARS-CoV-2 genomic diversity through space and time, to identify novel emerging viral variants and prioritize variants of potential epidemiological interest in a rapid and unsupervised manner. HaploCoV can integrate with any classification/nomenclature and incorporates an effective scoring system for the prioritization of SARS-CoV-2 variants. By performing retrospective analyses of more than 11.5 M genome sequences we show that HaploCoV demonstrates high levels of accuracy and reproducibility and identifies the large majority of epidemiologically relevant viral variants - as flagged by international health authorities - automatically and with rapid turn-around times.Our results highlight the importance of the application of strategies based on the systematic analysis and integration of regional data for rapid identification of novel, emerging variants of SARS-CoV-2. We believe that the approach outlined in this study will contribute to relevant advances to current and future genomic surveillance methods., (© 2023. The Author(s).)

Rotator cuff tendon (RCT) disease results from multifactorial mechanisms, in which inflammation plays a key role. Pro-inflammatory cytokines and tendon stem cell/progenitor cells (TSPCs) have been shown to participate in the inflammatory response. However, the underlying molecular mechanism is still not clear. In this study, flow cytometry analyses of different subpopulations of RCT-derived TSPCs demonstrate that after three days of administration, TNFα alone or in combination with IFNγ significantly decreases the percentage of CD146+CD49d+ and CD146+CD49f+ but not CD146+CD109+ TSPCs populations. In parallel, the same pro-inflammatory cytokines upregulate the expression of CD200 in the CD146+ TSPCs population. Additionally, the TNFα/IFNγ combination modulates the protein expression of STAT1, STAT3, and MMP9, but not fibromodulin. At the gene level, IRF1 , CAAT (CAAT/EBPbeta) , and DOK2 but not NF-κb , TGRF2 (TGFBR2) , and RAS-GAP are modulated. In conclusion, although our study has several important limitations, the results highlight a new potential role of CD200 in regulating inflammation during tendon injuries. In addition, the genes analyzed here might be new potential players in the inflammatory response of TSPCs.

In vertebrates, the four transcription factors Sox2 , c-Myc , Pou5f1 and Klf4 are involved in the differentiation of several tissues during vertebrate embryogenesis; moreover, they are normally co-expressed in embryonic stem cells and play roles in pluripotency, self-renewal, and maintenance of the undifferentiated state in adult cells. The in vitro forced co-expression of these factors, named Yamanaka factors (YFs), induces pluripotency in human or mouse fibroblasts. Botryllus schlosseri is a colonial tunicate undergoing continuous stem cell-mediated asexual development, providing a valuable model system for the study of pluripotency in the closest living relatives of vertebrates. In this study, we identified B. schlosseri orthologs of human Sox2 and c-Myc genes , as well as the closest homologs of the vertebrate-specific Pou5f1 gene, through an in-depth evolutionary analysis of the YF gene families in tunicates and other deuterostomes. Then, we studied the expression of these genes during the asexual cycle of B. schlosseri using in situ hybridization in order to investigate their possible involvement in tissue differentiation and in pluripotency maintenance. Our results show a shared spatio-temporal expression pattern consistent with the reported functions of these genes in invertebrate and vertebrate embryogenesis. Moreover, Myc , SoxB1 and Pou3 were expressed in candidate stem cells residing in their niches, while Pou2 was found expressed exclusively in the immature previtellogenic oocytes, both in gonads and circulating in the colonial vascular system. Our data suggest that Myc , SoxB1 and Pou3 may be individually involved in the differentiation of the same territories seen in other chordates, and that, together, they may play a role in stemness even in this colonial ascidian., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Vanni, Salonna, Gasparini, Martini, Anselmi, Gissi and Manni.)

Aquaporins (AQPs) are a family of membrane channels facilitating diffusion of water and small solutes into and out of cells. Despite their biological relevance in osmoregulation and ubiquitous distribution throughout metazoans, the presence of AQPs in annelids has been poorly investigated. Here, we searched and annotated Aqp sequences in public genomes and transcriptomes of annelids, inferred their evolutionary relationships through phylogenetic analyses and discussed their putative physiological relevance. We identified a total of 401 Aqp sequences in 27 annelid species, including 367 sequences previously unrecognized as Aqps . Similar to vertebrates, phylogenetic tree reconstructions clustered these annelid Aqps in four clades: AQP1-like, AQP3-like, AQP8-like and AQP11-like. We found no clear indication of the existence of paralogs exclusive to annelids; however, several gene duplications seem to have occurred in the ancestors of some Sedentaria annelid families, mainly in the AQP1-like clade. Three of the six Aqps annotated in Alitta succinea , an estuarine annelid showing high salinity tolerance, were validated by RT-PCR sequencing, and their similarity to human AQPs was investigated at the level of "key" conserved residues and predicted three-dimensional structure. Our results suggest a diversification of the structures and functions of AQPs in Annelida comparable to that observed in other taxa.

Seahorses are considered a flagship species for conservation efforts and due to their conservation status, improving knowledge on their dietary composition while applying a non-invasive approach, could be useful. Using Hippocampus guttulatus as a case study, the present study represents pioneering research into investigating the diet of seahorses by NGS-based DNA metabarcoding of fecal samples. The study developed and tested the protocol for fecal DNA metabarcoding during the feeding trials where captive seahorses were fed on a diet of known composition; the process was subsequently applied on fecal samples collected from wild individuals. The analysis of samples collected during the feeding trials indicated the reliability of the applied molecular approach by allowing the characterization of the effectively ingested prey. In the field study, among detected prey species, results revealed that the majority of the seahorse samples contained taxa such as Amphipoda, Decapoda, Isopoda, and Calanoida, while less common prey taxa were Gastropoda and Polyplacophora. As only a small amount of starting fecal material is needed and the sampling procedure is neither invasive nor lethal. The present study indicates DNA metabarcoding as useful for investigating seahorse diet and could help define management and conservation actions.

Effective systems for the analysis of molecular data are fundamental for monitoring the spread of infectious diseases and studying pathogen evolution. The rapid identification of emerging viral strains, and/or genetic variants potentially associated with novel phenotypic features is one of the most important objectives of genomic surveillance of human pathogens and represents one of the first lines of defense for the control of their spread. During the COVID 19 pandemic, several taxonomic frameworks have been proposed for the classification of SARS-Cov-2 isolates. These systems, which are typically based on phylogenetic approaches, represent essential tools for epidemiological studies as well as contributing to the study of the origin of the outbreak. Here, we propose an alternative, reproducible, and transparent phenetic method to study changes in SARS-CoV-2 genomic diversity over time. We suggest that our approach can complement other systems and facilitate the identification of biologically relevant variants in the viral genome. To demonstrate the validity of our approach, we present comparative genomic analyses of more than 175,000 genomes. Our method delineates 22 distinct SARS-CoV-2 haplogroups, which, based on the distribution of high-frequency genetic variants, fall into four major macrohaplogroups. We highlight biased spatiotemporal distributions of SARS-CoV-2 genetic profiles and show that seven of the 22 haplogroups (and of all of the four haplogroup clusters) showed a broad geographic distribution within China by the time the outbreak was widely recognized-suggesting early emergence and widespread cryptic circulation of the virus well before its isolation in January 2020. General patterns of genomic variability are remarkably similar within all major SARS-CoV-2 haplogroups, with UTRs consistently exhibiting the greatest variability, with s2m, a conserved secondary structure element of unknown function in the 3'-UTR of the viral genome showing evidence of a functional shift. Although several polymorphic sites that are specific to one or more haplogroups were predicted to be under positive or negative selection, overall our analyses suggest that the emergence of novel types is unlikely to be driven by convergent evolution and independent fixation of advantageous substitutions, or by selection of recombined strains. In the absence of extensive clinical metadata for most available genome sequences, and in the context of extensive geographic and temporal biases in the sampling, many questions regarding the evolution and clinical characteristics of SARS-CoV-2 isolates remain open. However, our data indicate that the approach outlined here can be usefully employed in the identification of candidate SARS-CoV-2 genetic variants of clinical and epidemiological importance., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Botryllids are colonial ascidians widely studied for their potential invasiveness and as model organisms, however the morphological description and discrimination of these species is very problematic, leading to frequent specimen misidentifications. To facilitate species discrimination and detection of cryptic/new species, we developed new barcoding primers for the amplification of a COI fragment of about 860 bp (860-COI), which is an extension of the common Folmer's barcode region. Our 860-COI was successfully amplified in 177 worldwide-sampled botryllid colonies. Combined with morphological analyses, 860-COI allowed not only discriminating known species, but also identifying undescribed and cryptic species, resurrecting old species currently in synonymy, and proposing the assignment of clade D of the model organism Botryllus schlosseri to Botryllus renierii. Importantly, within clade A of B. schlosseri, 860-COI recognized at least two candidate species against only one recognized by the Folmer's fragment, underlining the need of further genetic investigations on this clade. This result also suggests that the 860-COI could have a greater ability to diagnose cryptic/new species than the Folmer's fragment at very short evolutionary distances, such as those observed within clade A. Finally, our new primers simplify the amplification of 860-COI even in non-botryllid ascidians, suggesting their wider usefulness in ascidians.

Rotator cuff tears (RCTs) and rotator cuff disease (RCD) are important causes of disability in middle-aged individuals affected by nontraumatic shoulder dysfunctions. Our previous studies have demonstrated that four different hyaluronic acid preparations (HAPs), including Artrosulfur ® hyaluronic acid (HA) (Alfakjn S.r.l., Garlasco, Italy), may exert a protective effect in human RCT-derived tendon cells undergoing oxidative stress damage. Recently, methylsulfonylmethane (MSM) (Barentz, Paderno Dugnano, Italy) has proven to have anti-inflammatory properties and to cause pain relief in patients affected by tendinopathies. This study aims at evaluating three preparations (Artrosulfur ® HA, MSM, and Artrosulfur ® MSM + HA) in the recovery from hydrogen peroxide-induced oxidative stress damage in human tenocyte. Cell proliferation, Lactate Dehydrogenase (LDH) release, and inducible nitric oxide synthases (iNOS) and prostaglandin E2 (PGE2) modulation were investigated. In parallel, expression of metalloproteinases 2 (MMP2) and 14 (MMP14) and collagen types I and III were also examined. Results demonstrate that Artrosulfur ® MSM + HA improves cell escape from oxidative stress by decreasing cytotoxicity and by reducing iNOS and PGE2 secretion. Furthermore, it differentially modulates MMP2 and MMP14 levels and enhances collagen III expression after 24 h, proteins globally related to rapid acceleration of the extracellular matrix (ECM) remodelling and thus tendon healing. By improving the anti-cytotoxic effect of HA, the supplementation of MSM may represent a feasible strategy to ameliorate cuff tendinopathies.

Conventional central chondrosarcoma (CCC) is a malignant bone tumor that is characterized by the production of chondroid tissue. Since radiation therapy and chemotherapy have limited effects on CCC, treatment of most patients depends on surgical resection. This study aimed to identify the expression profiles of microRNAs (miRNAs) and isomiRs in CCC tissues to highlight their possible participation to the regulation of pathways critical for the formation and growth of this type of tumor. Our study analyzed miRNAs and isomiRs from Grade I (GI), Grade II (GII), and Grade III (GIII) histologically validated CCC tissue samples. While the different histological grades shared a similar expression profile for the top abundant miRNAs, we found several microRNAs and isomiRs showing a strong different modulation in GII + GIII vs GI grade samples and their involvement in tumor biology could be consistently hypothesized. We then in silico validated these differently expressed miRNAs in a larger chondrosarcoma public dataset and confirmed the expression trend for 17 out of 34 miRNAs. Our results clearly suggests that the contribution of miRNA deregulation, and their targeted pathways, to the progression of CCC could be relevant and strongly indicates that when studying miRNA deregulation in tumors, not only the canonical miRNAs, but the whole set of corresponding isomiRs should be taken in account. Improving understanding of the precise roles of miRNAs and isomiRs over the course of central chondrosarcoma progression could help identifying possible targets for precision medicine therapeutic intervention., Competing Interests: Conflict of interestThe authors declare that they have no conflict of interest., (© The Author(s) 2020.)

Mesenchymal stromal/stem cells (MSCs) are increasingly employed for tissue regeneration, largely mediated through paracrine actions. Currently, extracellular vesicles (EVs) released by MSCs are major mediators of these paracrine effects. We evaluated whether rat-bone-marrow-MSC-derived EVs (rBMSCs-EVs) can ameliorate tendon injury in an in vivo rat model. Pro-collagen1A2 and MMP14 protein are expressed in rBMSC-EVs, and are important factors for extracellular-matrix tendon-remodeling. In addition, we found pro-collagen1A2 in rBMSC-EV surface-membranes by dot blot. In vitro on cells isolated from Achilles tendons, utilized as rBMSC -EVs recipient cells, EVs at both low and high doses induce migration of tenocytes; at higher concentration, they induce proliferation and increase expression of Collagen type I in tenocytes. Pretreatment with trypsin abrogate the effect of EVs on cell proliferation and migration, and the expression of collagen I. When either low- or high-dose rBMSCs-EVs were injected into a rat-Achilles tendon injury-model (immediately after damage), at 30 days, rBMSC-EVs were found to have accelerated the remodeling stage of tendon repair in a dose-dependent manner. At histology and histomorphology evaluation, high doses of rBMSCs-EVs produced better restoration of tendon architecture, with optimal tendon-fiber alignment and lower vascularity. Higher EV-concentrations demonstrated greater expression of collagen type I and lower expression of collagen type III. BMSC-EVs hold promise as a novel cell-free modality for the management of tendon injuries., Competing Interests: The authors have declared that no competing interests exist.

Non-traumatic rotator cuff tears (RCTs) are a frequent and potentially disabling injury. There is growing evidence that hyaluronic acid (HA) is effective for pain relief and to counteract inflammation in RCTs, however, its effective role in tendinopathies remains poorly studied. This study aims to disclose a possible molecular mechanism underlying the cytoprotective effects of four different HA preparations (Artrosulfur HA ® , Synolis VA ® , Hyalgan ® and Hyalubrix ® ) under H 2 O 2 -induced oxidative stress. Expression-levels of Lactate dehydrogenase (LDH) released were quantified in cell supernatants, CD44 expression levels were analysed by fluorescence microscopy, the mitochondrial membrane depolarisation (TMRE assay) was measured by flow cytometry and the role of the transcription factor Nrf2 was investigated as a potential therapeutic target for RCT treatment. The modulation of extracellular matrix- (ECM) related protein expression (Integrin β1, pro-collagen 1A2 and collagen 1A1) and autophagy occurrence (Erk 1/2 and phosphoErk 1/2 and LC3B), were all investigated by Western Blot. Results demonstrate that Artrosulfur HA, Hyalubrix and Hyalgan improve cell escape from H 2 O 2 -induced oxidative stress, decreasing cytotoxicity, reducing Nrf2 expression and enhancing catalase recovery. This study lays the grounds for further investigations insight novel pharmaceutical strategies targeting key effectors involved in the molecular cascade triggered by HA.

In order to accelerate the isolation and characterization of structurally new or novel secondary metabolites, it is crucial to develop efficient strategies that prioritize samples with greatest promise early in the workflow so that resources can be utilized in a more efficient and cost-effective manner. We have developed a metrics-based prioritization approach using exact LC-HRMS, which uses data for 24 618 marine natural products held in the PharmaSea database. Each sample was evaluated and allocated a metric score by a software algorithm based on the ratio of new masses over the total (sample novelty), ratio of known masses over the total (chemical novelty), number of peaks above a defined peak area threshold (sample complexity), and peak area (sample diversity). Samples were then ranked and prioritized based on these metric scores. To validate the approach, eight marine sponges and six tunicate samples collected from the Fiji Islands were analyzed, metric scores calculated, and samples targeted for isolation and characterization of new compounds. Structures of new compounds were elucidated by spectroscopic techniques, including 1D and 2D NMR, MS, and MS/MS. Structures were confirmed by computer-assisted structure elucidation methods (CASE) using the ACD/Structure Elucidator Suite.

Background: This pilot study aimed to ascertain whether the local application of ascorbic acid (AA), of T 3 , and of rat (r) bone marrow mesenchymal stem cells (BMSCs), alone or in all possible combinations, promoted healing after an Achilles tendon injury in a rat model., Methods: An Achilles tendon defect was produced in 24 6-8-week-old male inbred Lewis rats. The animals were then randomly divided into eight groups of three rats each. The tendon defect was filled with 50 μL of phosphate-buffered saline (PBS) containing (1) 50 μg/mL AA (AA group), (2) 10 -7  M T 3 (T 3 group), (3) 4 × 10 6 rBMSCs (rBMSC group), (4) 50 μg/mL AA + 10 -7  M T 3 (AA + T 3 group), (5) 4 × 10 6 rBMSCs + 50 μg/mL AA (rBMSC + AA group), (6) 4 × 10 6 rBMSCs + 10 -7  M T 3 (rBMSC + T 3 group), (7) 4 × 10 6 rBMSCS + 50 μg/mL AA + 10 -7  M T 3 (rBMSC + AA + T 3 group), and (8) PBS only (control group: CTRL). All treatments were administered by local injection immediately after the tendons had been damaged; additionally, AA was injected also on the second and fourth day from the first injection (for groups 1, 4, 5, and 7), and T 3 was injected again every day for 4 days (for groups 2, 4, 6, and 7). At 30 days from initial treatment, tendon samples were harvested, and the quality of tendon repair was evaluated using histological and histomorphological analysis. The structure and morphology of the injured Achilles tendons were evaluated using the modified Svensson, Soslowsky, and Cook score, and the collagen type I and III ratio was calculated., Results: The group treated with AA combined with T 3 displayed the lowest Svensson, Soslowsky, and Cook total score value of all tissue sections at histopathological examination, with fiber structure close to regular orientation, normal-like tendon vasculature, and no cartilage formation. AA + T 3 also showed the highest collagen I and the lowest collagen III values compared to all other treatments including the CTRL., Conclusion: There are potential benefits using a combination of AA and T 3 to accelerate tendon healing.

Background: Complete mitochondrial (mt) genomes have been sequenced for thousands of animals and represent a molecule of choice for many evolutionary studies. Nevertheless, some animal groups have remained under-sampled. Ctenophora (comb jellies) is one such example, with only two complete mt sequences determined hitherto for this phylum, which encompasses ca. 150-200 described species. This lack of data derives from the extremely fast mt evolutionary rate in this lineage, complicating primer design and DNA amplification. Indeed, in the two ctenophore mt genomes sequenced to date, i.e. those of Mnemiopsis leidyi (order Lobata) and Pleurobrachia bachei (order Cydippida), both rRNA and protein coding genes exhibit an extraordinary size reduction and have highly derived sequences. Additionally, all tRNAs, and the atp6 and atp8 genes are absent. In order to determine whether these characteristics are shared by other ctenophores, we obtained the complete mt genomes of three benthic ctenophores belonging to the so far unsampled order of Platyctenida: Coeloplana loyai, Coeloplana yulianicorum and Vallicula multiformis., Results: The mt genomes of benthic ctenophores reveal the same peculiarities found in Mnemiopsis and Pleurobrachia, demonstrating that the fast evolutionary rate is a general trait of the ctenophore mt genomes. Our results also indicate that this high evolutionary rate not only affects the nucleotide substitution but also gene rearrangements. Indeed, gene order was highly rearranged among representatives of the different taxonomic orders in which it was close to random, but also quite variable within Platyctenida, in which the genera Coeloplana and Vallicula share only four conserved synteny blocks. However, the two congeneric Coeloplana species display exactly the same gene order. Because of the extreme evolutionary rate, our phylogenetic analyses were unable to resolve the phylogenetic position of ctenophores within metazoans or the relationships among the different Ctenophora orders. Comparative sequence-analyses allowed us to correct the annotation of the Pleurobrachia mt genome, confirming the absence of tRNAs, the presence of both rRNA genes, and the existence of a reassignment of codon TGA from tryptophan to serine for this species., Conclusions: Since Platyctenida is an early diverging lineage among Ctenophora, our findings suggest that the mt traits described above are ancestral characteristics of this phylum.

Botryllus schlosseri is a widespread colonial ascidian commonly considered cosmopolitan and amply used as model for researches ranging from developmental biology to immunobiology. Recently, molecular data lead to hypothesize that the species named B. schlosseri may consist of more than a single taxon. Indeed, five highly divergent clades, named A-E, have been genetically identified and are referred as cryptic species. In this context, and lacking both a type and a detailed morphological description, we believe that it is necessary, as a taxonomic reference point, to designate a neotype and re-describe the species. Therefore, a sample from the Lagoon of Venice (Adriatic Sea, Italy) was deposited as neotype in the Natural History Museum of Venice (Italy), preserved both in formalin and in 90% ethanol. Here we provide a morphological description of the suggested neotype of B. schlosseri that takes into account several developmental stages (oozooid, zooid of first blastogenic generations, and mature zooid) and is carefully compared with the previous descriptions of samples coming from other European and non-European localities. Finally, we associate our morphological description to a "DNA barcode", consisting in a long fragment of the mitochondrial COI gene. Our description is associated to clade A, although at now we cannot guarantee that this association is univocal.

Background: Hyaluronic Acid (HA) has been already approved by Food and Drug Administration (FDA) for osteoarthritis (OA), while its use in the treatment of tendinopathy is still debated. The aim of this study was to evaluate the effects of two different HA on human rotator cuff tendon derived cells in terms of cell viability, proliferation and apoptosis., Methods: An in vitro model was developed on human tendon derived cells from rotator cuff tears to study the effects of two different HA preparations: Sinovial HL® (High-Low molecular weight) (MW: 80-100 kDa) and KDa Sinovial Forte SF (MW: 800-1200), at various concentrations. Tendon derived cells morphology was evaluated after 0, 7 and 14 d of culture. Viability and proliferation were analyzed after 0, 24, and 48 h of culture and apoptosis occurrence was assessed after 24 h of culture., Results: All the HAPs tested here increased viability and proliferation, in a dose-dependent manner and they reduced apoptosis at early stages (24 h) compared to control cells (without HAPs)., Conclusions: HAPs enhanced viability and proliferation and counteracted apoptosis in tendon derived cells., Competing Interests: Conflict of interest The Authors declare no conflict of interest.

Liver cancer (LC) is one of the most common cancers and represents the third highest cause of cancer-related deaths worldwide. Extracellular vesicle (EVs) cargoes, which are selectively enriched in RNA, offer great promise for the diagnosis, prognosis and treatment of LC. Our study analyzed the RNA cargoes of EVs derived from 4 liver-cancer cell lines: HuH7, Hep3B, HepG2 (hepato-cellular carcinoma) and HuH6 (hepatoblastoma), generating two different sets of sequencing libraries for each. One library was size-selected for small RNAs and the other targeted the whole transcriptome. Here are reported genome wide data of the expression level of coding and non-coding transcripts, microRNAs, isomiRs and snoRNAs providing the first comprehensive overview of the extracellular-vesicle RNA cargo released from LC cell lines. The EV-RNA expression profiles of the four liver cancer cell lines share a similar background, but cell-specific features clearly emerge showing the marked heterogeneity of the EV-cargo among the individual cell lines, evident both for the coding and non-coding RNA species., Competing Interests: CONFLICTS OF INTEREST The authors declare that there are no conflicts of interest.

The order Aplousobranchia (Chordata, Ascidiacea) contains approximately 1500 species distributed worldwide. Their phylogeny, however, remains unclear, with unresolved family relationships. While most Aplousobranchia are colonial, debates exist concerning the phylogenetic position of families such as the Diazonidae and Cionidae, which exhibit a solitary lifestyle and share morphological characteristics with both Aplousobranchia and Phlebobranchia orders. To clarify the phylogenetic position of the Diazonidae and Cionidae, we determined the complete mitochondrial sequence of the solitary diazonid Rhopalaea idoneta. The phylogenetic reconstruction based on the 13 mitochondrial protein coding genes strongly supports a positioning of Diazonidae well-nested within the Aplousobranchia rather than a positioning as a sister clade of the Aplousobranchia. In addition, we examined the regenerative ability of R. idoneta. Similar to colonial Aplousobranchia, R. idoneta was found to be able to completely regenerate its thorax. Ciona, also known to possess high regenerative abilities, is the Aplousobranchia sister clade rather than a member of the Phlebobranchia. Our results thus indicate that the colonial lifestyle was acquired in the Aplousobranchia, starting from a Ciona-like solitary ancestor and secondarily lost in Diazonidae representatives such as Rhopalaea. The solitary lifestyle of Rhopalaea is thus a derived characteristic rather than an ancestral trait., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Ascidians belong to the tunicates, the sister group of vertebrates and are recognized model organisms in the field of embryonic development, regeneration and stem cells. ANISEED is the main information system in the field of ascidian developmental biology. This article reports the development of the system since its initial publication in 2010. Over the past five years, we refactored the system from an initial custom schema to an extended version of the Chado schema and redesigned all user and back end interfaces. This new architecture was used to improve and enrich the description of Ciona intestinalis embryonic development, based on an improved genome assembly and gene model set, refined functional gene annotation, and anatomical ontologies, and a new collection of full ORF cDNAs. The genomes of nine ascidian species have been sequenced since the release of the C. intestinalis genome. In ANISEED 2015, all nine new ascidian species can be explored via dedicated genome browsers, and searched by Blast. In addition, ANISEED provides full functional gene annotation, anatomical ontologies and some gene expression data for the six species with highest quality genomes. ANISEED is publicly available at: http://www.aniseed.cnrs.fr., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The cosmopolitan ascidian Ciona intestinalis is the most common model species of Tunicata, the sister-group of Vertebrata, and widely used in developmental biology, genomics and evolutionary studies. Recently, molecular studies suggested the presence of cryptic species hidden within the C. intestinalis species, namely C. intestinalis type A and type B. So far, no substantial morphological differences have been identified between individuals belonging to the two types. Here we present morphometric, immunohistochemical, and histological analyses, as well as 3-D reconstructions, of late larvae obtained by cross-fertilization experiments of molecularly determined type A and type B adults, sampled in different seasons and in four different localities. Our data point to quantitative and qualitative differences in the trunk shape of larvae belonging to the two types. In particular, type B larvae exhibit a longer pre-oral lobe, longer and relatively narrower total body length, and a shorter ocellus-tail distance than type A larvae. All these differences were found to be statistically significant in a Discriminant Analysis. Depending on the number of analyzed parameters, the obtained discriminant function was able to correctly classify > 93% of the larvae, with the remaining misclassified larvae attributable to the existence of intra-type seasonal variability. No larval differences were observed at the level of histology and immunohistochemical localization of peripheral sensory neurons. We conclude that type A and type B are two distinct species that can be distinguished on the basis of larval morphology and molecular data. Since the identified larval differences appear to be valid diagnostic characters, we suggest to raise both types to the rank of species and to assign them distinct names.

Eukaryotic cells contain a population of mitochondria, variable in number and shape, which in turn contain multiple copies of a tiny compact genome (mtDNA) whose expression and function is strictly coordinated with the nuclear one. mtDNA copy number varies between different cell or tissues types, both in response to overall metabolic and bioenergetics demands and as a consequence or cause of specific pathological conditions. Here we present a novel and reliable methodology to assess the effective mtDNA copy number per diploid genome by investigating off-target reads obtained by whole-exome sequencing (WES) experiments. We also investigate whether and how mtDNA copy number correlates with mitochondrial mass, respiratory activity and expression levels. Analyzing six different tissues from three age- and sex-matched human individuals, we found a highly significant linear correlation between mtDNA copy number estimated by qPCR and the frequency of mtDNA off target WES reads. Furthermore, mtDNA copy number showed highly significant correlation with mitochondrial gene expression levels as measured by RNA-Seq as well as with mitochondrial mass and respiratory activity. Our methodology makes thus feasible, at a large scale, the investigation of mtDNA copy number in diverse cell-types, tissues and pathological conditions or in response to specific treatments., (Copyright © 2014 © Elsevier B.V. and Mitochondria Research Society. Published by Elsevier B.V. All rights reserved.)

Ascidians are a fascinating group of filter-feeding marine chordates characterized by rapid evolution of both sequences and structure of their nuclear and mitochondrial genomes. Moreover, they include several model organisms used to investigate complex biological processes in chordates. To study the evolutionary dynamics of ascidians at short phylogenetic distances, we sequenced 13 new mitogenomes and analyzed them, together with 15 other available mitogenomes, using a novel approach involving detailed whole-mitogenome comparisons of conspecific and congeneric pairs. The evolutionary rate was quite homogeneous at both intraspecific and congeneric level, and the lowest congeneric rates were found in cryptic (morphologically undistinguishable) and in morphologically very similar species pairs. Moreover, congeneric nonsynonymous rates (dN) were up to two orders of magnitude higher than in intraspecies pairs. Overall, a clear-cut gap sets apart conspecific from congeneric pairs. These evolutionary peculiarities allowed easily identifying an extraordinary intraspecific variability in the model ascidian Botryllus schlosseri, where most pairs show a dN value between that observed at intraspecies and congeneric level, yet consistently lower than that of the Ciona intestinalis cryptic species pair. These data suggest ongoing speciation events producing genetically distinct B. schlosseri entities. Remarkably, these ongoing speciation events were undetectable by the cox1 barcode fragment, demonstrating that, at low phylogenetic distances, the whole mitogenome has a higher resolving power than cox1. Our study shows that whole-mitogenome comparative analyses, performed on a suitable sample of congeneric and intraspecies pairs, may allow detecting not only cryptic species but also ongoing speciation events.

Botryllus schlosseri is a colonial urochordate that follows the chordate plan of development following sexual reproduction, but invokes a stem cell-mediated budding program during subsequent rounds of asexual reproduction. As urochordates are considered to be the closest living invertebrate relatives of vertebrates, they are ideal subjects for whole genome sequence analyses. Using a novel method for high-throughput sequencing of eukaryotic genomes, we sequenced and assembled 580 Mbp of the B. schlosseri genome. The genome assembly is comprised of nearly 14,000 intron-containing predicted genes, and 13,500 intron-less predicted genes, 40% of which could be confidently parceled into 13 (of 16 haploid) chromosomes. A comparison of homologous genes between B. schlosseri and other diverse taxonomic groups revealed genomic events underlying the evolution of vertebrates and lymphoid-mediated immunity. The B. schlosseri genome is a community resource for studying alternative modes of reproduction, natural transplantation reactions, and stem cell-mediated regeneration. DOI:http://dx.doi.org/10.7554/eLife.00569.001.

The β(1,3)-glucanosyltransferases of the GH72 family are redundant enzymes that are essential for the formation and dynamic remodeling of the fungal wall during different stages of the life cycle. Four putative genes encoding glycosylphosphatidylinositol (GPI)-anchored β(1,3)-glucanosyltransferases, designated TmelGEL1, TmelGEL2, TmelGEL4 and TmelGAS4, have been annotated in the genome of Tuber melanosporum, an ectomycorrhizal fungus that also produces a hypogeous fruiting body (FB) of great commercial value (black truffle). This work focuses on the characterization and expression of this multigene family by taking advantage of a laser microdissection (LMD) technology that has been used to separate two distinct compartments in the FB, the hyphae and the asci containing the ascospores. Of the four genes, TmelGEL1 was the most up-regulated in the FB compared to the free-living mycelium. Inside the FB, the expression of TmelGEL1 was restricted to the hyphal compartment. A phylogenetic analysis of the Gel/Gas protein family of T. melanosporum was also carried out. A total of 237 GH72 proteins from 51 Ascomycotina and 3 Basidiomycota (outgroup) species were analyzed. The resulting tree provides insight into the evolution of the T. melanosporum proteins and identifies new GH72 paralogs/subfamilies. Moreover, it represents a starting point to formulate new hypotheses on the significance of the striking GH72 gene redundancy in fungal biology., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Ascidians or sea squirts form a diverse group within chordates, which includes a few thousand members of marine sessile filter-feeding animals. Their mitochondrial genomes are characterized by particularly high evolutionary rates and rampant gene rearrangements. This extreme variability complicates standard polymerase chain reaction (PCR) based techniques for molecular characterization studies, and consequently only a few complete Ascidian mitochondrial genome sequences are available. Using the standard PCR and Sanger sequencing approach, we produced the mitochondrial genome of Ascidiella aspersa only after a great effort. In contrast, we produced five additional mitogenomes (Botrylloides aff. leachii, Halocynthia spinosa, Polycarpa mytiligera, Pyura gangelion, and Rhodosoma turcicum) with a novel strategy, consisting in sequencing the pooled total DNA samples of these five species using one Illumina HiSeq 2000 flow cell lane. Each mitogenome was efficiently assembled in a single contig using de novo transcriptome assembly, as de novo genome assembly generally performed poorly for this task. Each of the new six mitogenomes presents a different and novel gene order, showing that no syntenic block has been conserved at the ordinal level (in Stolidobranchia and in Phlebobranchia). Phylogenetic analyses support the paraphyly of both Ascidiacea and Phlebobranchia, with Thaliacea nested inside Phlebobranchia, although the deepest nodes of the Phlebobranchia-Thaliacea clade are not well resolved. The strategy described here thus provides a cost-effective approach to obtain complete mitogenomes characterized by a highly plastic gene order and a fast nucleotide/amino acid substitution rate.