Search

Skeletal muscle regeneration after injury is a complex process involving inflammatory signaling and myoblast activation. Pro-inflammatory cytokines like tumor necrosis factor-alpha (TNF-α) are key mediators, but their effects on gene expression in proliferating myoblasts are unclear. We performed the RNA sequencing of TNF-α treated C2C12 myoblasts to elucidate the signaling pathways and gene networks regulated by TNF-α during myoblast proliferation. The TNF-α (10 ng/mL) treatment of C2C12 cells led to 958 differentially expressed genes compared to the controls. Pathway analysis revealed significant regulation of TNF-α signaling, along with the chemokine and IL-17 pathways. Key upregulated genes included cytokines (e.g., IL-6), chemokines (e.g., CCL7), and matrix metalloproteinases (MMPs). TNF-α increased myogenic factor 5 (Myf5) but decreased MyoD protein levels and stimulated the release of MMP-9, MMP-10, and MMP-13. TNF-α also upregulates versican and myostatin mRNA. Overall, our study demonstrates the TNF-α modulation of distinct gene expression patterns and signaling pathways that likely contribute to enhanced myoblast proliferation while suppressing premature differentiation after muscle injury. Elucidating the mechanisms involved in skeletal muscle regeneration can aid in the development of regeneration-enhancing therapeutics.

Although non-front fanged snakes account for almost two-thirds of snake diversity, most studies on venom composition and evolution focus exclusively on front-fanged species, which comprise most of the clinically relevant accidents. Comprehensive reports on venom composition of non-front fanged snakes are still scarce for several groups. In this study, we address such shortage of knowledge by providing new insights about the venom composition among species of Phalotris, a poorly studied Neotropical dipsadid genus. Phalotris are known for their specialized venom delivery system and toxic venoms, which can cause life-threatening accidents in humans. We evaluate the venom-gland transcriptome of Phalotris, comparing the following three South American species: P. reticulatus for the Araucaria Pine forests, P. lemniscatus for the Pampa grasslands, and P. mertensi for the Brazilian Cerrado. Our results indicate similar venom profiles, in which they share a high expression level of Kunitz-type inhibitors (KUNZ). On the other hand, comparative analyses revealed substantial differences in the expression levels of C-type lectins (CTL) and snake venom metalloproteinases (SVMP). The diverse set of SVMP and CTL isoforms shows signals of positive selection, and we also identified truncated forms of type III SVMPs, which resemble type II and type I SVMPs of viperids. Additionally, we identified a CNP precursor hosting a proline-rich region containing a BPP motif resembling those commonly detected in viperid venoms with hypotensive activity. Altogether, our results suggest an evolutionary history favoring high expression levels of few KUNZ isoforms in Phalotris venoms, contrasting with a highly diverse set of SVMP and CTL isoforms. Such diversity can be comparable with the venom variability observed in some viperids. Our findings highlight the extreme phenotypic diversity of non-front fanged snakes and the importance to allocate greater effort to study neglected groups of Colubroidea., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests:Pedro Nachtigall reports financial support was provided by State of Sao Paulo Research Foundation. Omar M. Entiauspe-Neto reports financial support was provided by State of Sao Paulo Research Foundation. Inacio Loiola Meirelles Junqueira de Azevedo reports financial support was provided by National Council for Scientific and Technological Development. Felipe Gobbi Grazziotin reports financial support was provided by National Council for Scientific and Technological Development. Felipe Gobbi Grazziotin reports financial support was provided by State of Sao Paulo Research Foundation. Inacio Loiola Meirelles Junqueira de Azevedo reports financial support was provided by State of Sao Paulo Research Foundation. Omar M. Entiauspe-Neto reports financial support was provided by Coordination of Higher Education Personnel Improvement. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Snake venoms have evolved in several families of Caenophidae, and their toxins have been assumed to be biochemical weapons with a role as a trophic adaptation. However, it remains unclear how venom contributes to the success of venomous species for adaptation to different environments. Here we compared the venoms from Bothrocophias hyoprora , Bothrops taeniatus , Bothrops bilineatus smaragdinus , Bothrops brazili , and Bothrops atrox collected in the Amazon Rainforest, aiming to understand the ecological and toxinological consequences of venom composition. Transcriptomic and proteomic analyses indicated that the venoms presented the same toxin groups characteristic from bothropoids, but with distinct isoforms with variable qualitative and quantitative abundances, contributing to distinct enzymatic and toxic effects. Despite the particularities of each venom, commercial Bothrops antivenom recognized the venom components and neutralized the lethality of all species. No clear features could be observed between venoms from arboreal and terrestrial habitats, nor in the dispersion of the species throughout the Amazon habitats, supporting the notion that venom composition may not shape the ecological or toxinological characteristics of these snake species and that other factors influence their foraging or dispersal in different ecological niches.

Background: The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process., Results: Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models., Conclusions: ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome., (© The Author(s) 2024. Published by Oxford University Press GigaScience.)

Oral and other cephalic glands have been surveyed by several studies with distinct purposes. Despite the wide diversity and medical relevance of the New World coral snakes, studies focusing on understanding the biological roles of the glands within this group are still scarce. Specifically, the venom glands of some coral snakes were previously investigated but all other cephalic glands remain uncharacterized. In this sense, performing morphological and molecular analysis of these glands may help better understand their biological role. Here, we studied the morphology of the venom, infralabial, rictal, and harderian glands of thirteen species of Micrurus and Micruroides euryxanthus. We also performed a molecular characterization of these glands from selected species of Micrurus using transcriptomic and proteomic approaches. We described substantial morphological variation in the cephalic glands of New World coral snakes and structural evidence for protein-secreting cells in the inferior rictal glands. Our molecular analysis revealed that the venom glands, as expected, are majorly devoted to toxin production, however, the infralabial and inferior rictal glands also expressed some toxin genes at low to medium levels, despite the marked morphological differences. On the other hand, the harderian glands were dominated by the expression of lipocalins, but do not produce toxins. Our integrative analysis, including the prediction of biological processes and pathways, helped decipher some important traits of cephalic glands and better understand their biology., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Snake venoms harbor a wide and diverse array of enzymatic and nonenzymatic toxic components, allowing them to exert myriad effects on their prey. However, they appear to trend toward a few optimal compositional scaffolds, dominated by four major toxin classes: SVMPs, SVSPs, 3FTxs, and PLA2s. Nevertheless, the latter appears to be restricted to vipers and elapids, as it has never been reported as a major venom component in rear-fanged species. Here, by investigating the original transcriptomes from 19 species distributed in eight genera from the Pseudoboini tribe (Dipsadidae: Xenodontinae) and screening among seven additional tribes of Dipsadidae and three additional families of advanced snakes, we discovered that a novel type of venom PLA2, resembling a PLA2-IIE, has been recruited to the venom of some species of the Pseudoboini tribe, where it is a major component. Proteomic and functional analyses of these venoms further indicate that these PLA2s play a relevant role in the venoms from this tribe. Moreover, we reconstructed the phylogeny of PLA2s across different snake groups and show that different types of these toxins have been recruited in at least five independent events in caenophidian snakes. Additionally, we present the first compositional profiling of Pseudoboini venoms. Our results demonstrate how relevant phenotypic traits are convergently recruited by different means and from homologous and nonhomologous genes in phylogenetically and ecologically divergent snake groups, possibly optimizing venom composition to overcome diverse adaptative landscapes., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Snakes of the Philodryadini tribe are included in the Dipsadidae family, which is a diverse group of rear-fanged snakes widespread in different ecological conditions, including habitats and diet. However, little is known about the composition and effects of their venoms despite their relevance for understanding the evolution of these snakes or even their impact on the occasional cases of human envenoming. In this study, we integrated venom gland transcriptomics, venom proteomics and functional assays to characterize the venoms from eight species of the Philodryadini tribe, which includes the genus Philodryas , Chlorosoma and Xenoxybelis . The most abundant components identified in the venoms were snake venom metalloproteinases (SVMPs), cysteine-rich secretory proteins (CRISPs), C-type lectins (CTLs), snake endogenous matrix metalloproteinases type 9 (seMMP-9) and snake venom serinoproteinases (SVSPs). These protein families showed a variable expression profile in each genus. SVMPs were the most abundant components in Philodryas , while seMMP-9 and CRISPs were the most expressed in Chlorosoma and Xenoxybelis , respectively. Lineage-specific differences in venom composition were also observed among Philodryas species, whereas P. olfersii presented the highest amount of SVSPs and P. agassizii was the only species to express significant amounts of 3FTx. The variability observed in venom composition was confirmed by the venom functional assays. Philodryas species presented the highest SVMP activity, whereas Chlorosoma species showed higher levels of gelatin activity, which may correlate to the seMMP-9 enzymes. The variability observed in the composition of these venoms may be related to the tribe phylogeny and influenced by their diets. In the presented study, we expanded the set of venomics studies of the Philodryadini tribe, which paves new roads for further studies on the evolution and ecology of Dipsadidae snakes.

Biological specialization reduces the size of niche space while increasing efficiency in the use of available resources. Specialization often leads to phenotypic changes via natural selection aligning with niche space constraints. Commonly observed changes are in size, shape, behavior, and traits associated with feeding. One often selected trait for dietary specialization is venom, which, in snakes, often shows variation dependent on diet across and within species. The Neotropical Blunt-headed Treesnake (Imantodes cenchoa) is a highly specialized, rear-fanged, arboreal, lizard hunter that displays a long thin body, enlarged eyes, and a large Duvernoy's gland. However, toxin characterization of I. cenchoa has never been completed. Here, we use RNA-seq and mass spectrometry to assemble, annotate, and analyze the venom gland transcriptomes of four I. cenchoa from across their range. We find a lack of significant venom variation at the sequence and expression levels, suggesting venom conservation across the species. We propose this conservation provides evidence of a specialized venom repertoire, adapted to maximize efficiency of capturing and processing lizards. Importantly, this study provides the most complete venom gland transcriptomes of I. cenchoa and evidence of venom specialization in a rear-fanged snake, giving insight into selective pressures of venom across all snake species., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Interspecific differences in snake venom compositions can result from distinct regulatory mechanisms acting in each species. However, comparative analyses focusing on identifying regulatory elements and patterns that led to distinct venom composition are still scarce. Among venomous snakes, Bothrops cotiara and Bothrops fonsecai represent ideal models to complement our understanding of the regulatory mechanisms of venom production. These recently diverged species share a similar specialized diet, habitat, and natural history, but each presents a distinct venom phenotype. Here, we integrated data from the venom gland transcriptome and miRNome and the venom proteome of B. fonsecai and B. cotiara to better understand the regulatory mechanisms that may be acting to produce differing venom compositions. We detected not only the presence of similar toxin isoforms in both species but also distinct expression profiles of phospholipases A2 (PLA2) and some snake venom metalloproteinases (SVMPs) and snake venom serine proteinases (SVSPs) isoforms. We found evidence of modular expression regulation of several toxin isoforms implicated in venom divergence and observed correlated expression of several transcription factors. We did not find strong evidence for miRNAs shaping interspecific divergence of the venom phenotypes, but we identified a subset of toxin isoforms whose final expression may be fine-tuned by specific miRNAs. Sequence analysis on orthologous toxins showed a high rate of substitutions between PLA2s, which indicates that these toxins may be under strong positive selection or represent paralogous toxins in these species. Our results support other recent studies in suggesting that gene regulation is a principal mode of venom evolution across recent timescales, especially among species with conserved ecotypes.

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using μCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of 'harmless' non-front-fanged snakes., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

DNA methylation is one of the epigenetic modifications that configures gene transcription programs. This study describes the DNA methylation profile of HIV-infected individuals with distinct characteristics related to natural and artificial viremia control. Sheared DNA from circulating mononuclear cells was subjected to target enrichment bisulfite sequencing designed to cover CpG-rich genomic regions. Gene expression was assessed through RNA-seq. Hypermethylation in virologic responders was highly distributed closer to Transcription Start Sites (p-value = 0.03). Hyper and hypomethylation levels within TSS adjacencies varied according to disease progression status (Kruskal-Wallis, p < 0.001), and specific differentially methylated regions associated genes were identified for each group. The lower the promoter methylation, the higher the gene expression in subjects undergoing virologic failure (R = - 0.82, p = 0.00068). Among the inversely correlated genes, those supporting glycolysis and its related pathways were hypomethylated and up-regulated in virologic failures. Disease progression heterogeneity was associated with distinct DNA methylation patterns in terms of rates and distribution. Methylation was associated with the expression of genes sustaining intracellular glucose metabolism in subjects undergoing antiretroviral virologic failure. Our findings highlight that DNA methylation is associated with latency, disease progression, and fundamental cellular processes., (© 2021. The Author(s).)

Motivation: Next-generation sequencing has become exceedingly common and has transformed our ability to explore nonmodel systems. In particular, transcriptomics has facilitated the study of venom and evolution of toxins in venomous lineages; however, many challenges remain. Primarily, annotation of toxins in the transcriptome is a laborious and time-consuming task. Current annotation software often fails to predict the correct coding sequence and overestimates the number of toxins present in the transcriptome. Here, we present ToxCodAn, a python script designed to perform precise annotation of snake venom gland transcriptomes. We test ToxCodAn with a set of previously curated transcriptomes and compare the results to other annotators. In addition, we provide a guide for venom gland transcriptomics to facilitate future research and use Bothrops alternatus as a case study for ToxCodAn and our guide., Results: Our analysis reveals that ToxCodAn provides precise annotation of toxins present in the transcriptome of venom glands of snakes. Comparison with other annotators demonstrates that ToxCodAn has better performance with regard to run time ($>20x$ faster), coding sequence prediction ($>3x$ more accurate) and the number of toxins predicted (generating $>4x$ less false positives). In this sense, ToxCodAn is a valuable resource for toxin annotation. The ToxCodAn framework can be expanded in the future to work with other venomous lineages and detect novel toxins., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Motivation: Over the past decade, the field of next-generation sequencing (NGS) has seen dramatic advances in methods and a decrease in costs. Consequently, a large expansion of data has been generated by NGS, most of which have originated from RNA-sequencing (RNA-seq) experiments. Because mitochondrial genes are expressed in most eukaryotic cells, mitochondrial mRNA sequences are usually co-sequenced within the target transcriptome, generating data that are commonly underused or discarded. Here, we present MITGARD, an automated pipeline that reliably recovers the mitochondrial genome from RNA-seq data from various sources. The pipeline identifies mitochondrial sequence reads based on a phylogenetically related reference, assembles them into contigs, and extracts a complete mtDNA for the target species., Results: We demonstrate that MITGARD can reconstruct the mitochondrial genomes of several species throughout the tree of life. We noticed that MITGARD can recover the mitogenomes in different sequencing schemes and even in a scenario of low-sequencing depth. Moreover, we showed that the use of references from congeneric species diverging up to 30 million years ago (MYA) from the target species is sufficient to recover the entire mitogenome, whereas the use of species diverging between 30 and 60 MYA allows the recovery of most mitochondrial genes. Additionally, we provide a case study with original data in which we estimate a phylogenetic tree of snakes from the genus Bothrops, further demonstrating that MITGARD is suitable for use on biodiversity projects. MITGARD is then a valuable tool to obtain high-quality information for studies focusing on the phylogenetic and evolutionary aspects of eukaryotes and provides data for easily identifying a sample using barcoding, and to check for cross-contamination using third-party tools., (© The Author(s) 2021. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca , a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A 2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option., Competing Interests: The authors declare no competing interest.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Ontogenetic changes in venom composition have been described in Bothrops snakes, but only a few studies have attempted to identify the targeted paralogues or the molecular mechanisms involved in modifications of gene expression during ontogeny. In this study, we decoded B. jararacussu venom gland transcripts from six specimens of varying sizes and analyzed the variability in the composition of independent venom proteomes from 19 individuals. We identified 125 distinct putative toxin transcripts, and of these, 73 were detected in venom proteomes and only 10 were involved in the ontogenetic changes. Ontogenetic variability was linearly related to snake size and did not correspond to the maturation of the reproductive stage. Changes in the transcriptome were highly predictive of changes in the venom proteome. The basic myotoxic phospholipases A 2 (PLA 2 s) were the most abundant components in larger snakes, while in venoms from smaller snakes, PIII-class SVMPs were the major components. The snake venom metalloproteinases (SVMPs) identified corresponded to novel sequences and conferred higher pro-coagulant and hemorrhagic functions to the venom of small snakes. The mechanisms modulating venom variability are predominantly related to transcriptional events and may consist of an advantage of higher hematotoxicity and more efficient predatory function in the venom from small snakes.

Pararamosis is a disease that occurs due to contact with the hairs of the larval stage of the Brazilian moth Premolis semirufa . Envenomation induces osteoarticular alterations with cartilage impairment that resembles joint synovitis. Thus, the toxic venom present in the caterpillar hairs interferes with the phenotype of the cells present in the joints, resulting in inflammation and promoting tissue injury. Therefore, to address the inflammatory mechanisms triggered by envenomation, we studied the effects of P . semirufa hair extract on human chondrocytes. We have selected for the investigation, cytokines, chemokines, matrix metalloproteinases (MMPs), complement components, eicosanoids, and extracellular matrix (ECM) components related to OA and RA. In addition, for measuring protein-coding mRNAs of some molecules associated with osteoarthritis (OA) and rheumatoid arthritis (RA), reverse transcription (RT) was performed followed by quantitative real-time PCR (RT-qPCR) and we performed the RNA-sequencing (RNA-seq) analysis of the chondrocytes transcriptome. In the supernatant of cell cultures treated with the extract, we observed increased IL-6, IL-8, MCP-1, prostaglandin E2, metalloproteinases (MMP-1, MMP-2, MMP-3 and MMP-13), and complement system components (C3, C4, and C5). We noticed a significant decrease in both aggrecan and type II collagen and an increase in HMGB1 protein in chondrocytes after extract treatment. RNA-seq analysis of the chondrocyte transcriptome allowed us to identify important pathways related to the inflammatory process of the disease, such as the inflammatory response, chemotaxis of immune cells and extracellular matrix (ECM) remodeling. Thus, these results suggest that components of Premolis semirufa hair have strong inflammatory potential and are able to induce cartilage degradation and ECM remodeling, promoting a disease with an osteoarthritis signature. Modulation of the signaling pathways that were identified as being involved in this pathology may be a promising approach to develop new therapeutic strategies for the control of pararamosis and other inflammatory joint diseases., (Copyright © 2020 Villas-Boas, Pidde, Lichtenstein, Ching, Junqueira-de-Azevedo, DeOcesano-Pereira, Madureira Trufen, Chudzinski-Tavassi, Morais and Tambourgi.)

Haemocyanins (Hcs) are copper-containing, respiratory proteins that occur in the haemolymph of many arthropod species. Here, we report the presence of Hcs in the chilopode Myriapoda, demonstrating that these proteins are more widespread among the Arthropoda than previously thought. The analysis of transcriptome of S. subspinipes subpinipes reveals the presence of two distinct subunits of Hc, where the signal peptide is present, and six of prophenoloxidase (PPO), where the signal peptide is absent, in the 75 kDa range. Size exclusion chromatography profiles indicate different quaternary organization for Hc of both species, which was corroborated by TEM analysis: S. viridicornis Hc is a 6 × 6-mer and S. subspinipes Hc is a 3 × 6-mer, which resembles the half-structure of the 6 × 6-mer but also includes the presence of phenoloxidases, since the 1 × 6-mer quaternary organization is commonly associated with hexamers of PPO. Studies with Chelicerata showed that PPO activity are exclusively associated with the Hcs. This study indicates that Scolopendra may have different proteins playing oxygen transport (Hc) and PO function, both following the hexameric oligomerization observed in Hcs.

Coffea arabica L. is an important agricultural commodity, accounting for 60% of traded coffee worldwide. Nitrogen (N) is a macronutrient that is usually limiting to plant yield; however, molecular mechanisms of plant acclimation to N limitation remain largely unknown in tropical woody crops. In this study, we investigated the transcriptome of coffee roots under N starvation, analyzing poly-A+ libraries and small RNAs. We also evaluated the concentration of selected amino acids and N-source preferences in roots. Ammonium was preferentially taken up over nitrate, and asparagine and glutamate were the most abundant amino acids observed in coffee roots. We obtained 34,654 assembled contigs by mRNA sequencing, and validated the transcriptional profile of 12 genes by RT-qPCR. Illumina small RNA sequencing yielded 8,524,332 non-redundant reads, resulting in the identification of 86 microRNA families targeting 253 genes. The transcriptional pattern of eight miRNA families was also validated. To our knowledge, this is the first catalog of differentially regulated amino acids, N sources, mRNAs, and sRNAs in Arabica coffee roots.

Elucidating the molecular mechanisms underlying snake venom variability provides important clues for understanding how the biological functions of this powerful toxic arsenal evolve. We analyzed in detail individual transcripts and venom protein isoforms produced by five specimens of a venomous snake (Bothrops atrox) from two nearby but genetically distinct populations from the Brazilian Amazon rainforest which show functional similarities in venom properties. Individual variation was observed among the venoms of these specimens, but the overall abundance of each general toxin family was conserved both in transcript and in venom protein levels. However, when expression of independent paralogues was analyzed, remarkable differences were observed within and among each toxin group, both between individuals and between populations. Transcripts for functionally essential venom proteins ("core function" proteins) were highly expressed in all specimens and showed similar transcription/translation rates. In contrast, other paralogues ("adaptive" proteins) showed lower expression levels and the toxins they coded for varied among different individuals. These results provide support for the inferences that (a) expression and translational differences play a greater role in defining adaptive variation in venom phenotypes than does sequence variation in protein coding genes and (b) convergent adaptive venom phenotypes can be generated through different molecular mechanisms., Significance: Analysis of individual transcripts and venom protein isoforms produced by specimens of a venomous snake (Bothrops atrox), from the Brazilian Amazon rainforest, revealed that transcriptional and translational mechanisms contribute to venom phenotypic variation. Our finding of evidence for high expression of toxin proteins with conserved function supports the hypothesis that the venom phenotype consists of two kinds of proteins: conserved "core function" proteins that provide essential functional activities with broader relevance and less conserved "adaptive" proteins that vary in expression and may permit customization of protein function. These observations allowed us to suggest that genetic mechanisms controlling venom variability are not restricted to selection of gene copies or mutations in structural genes but also to selection of the mechanisms controlling gene expression, contributing to the plasticity of this important phenotype for venomous snakes., (Copyright © 2018. Published by Elsevier B.V.)

Background: Except for the northern region, where the Amazonian black scorpion, T. obscurus, represents the predominant and most medically relevant scorpion species, Tityus serrulatus, the Brazilian yellow scorpion, is widely distributed throughout Brazil, causing most envenoming and fatalities due to scorpion sting. In order to evaluate and compare the diversity of venom components of Tityus obscurus and T. serrulatus, we performed a transcriptomic investigation of the telsons (venom glands) corroborated by a shotgun proteomic analysis of the venom from the two species., Results: The putative venom components represented 11.4% and 16.7% of the total gene expression for T. obscurus and T. serrulatus, respectively. Transcriptome and proteome data revealed high abundance of metalloproteinases sequences followed by sodium and potassium channel toxins, making the toxin core of the venom. The phylogenetic analysis of metalloproteinases from T. obscurus and T. serrulatus suggested an intraspecific gene expansion, as we previously observed for T. bahiensis, indicating that this enzyme may be under evolutionary pressure for diversification. We also identified several putative venom components such as anionic peptides, antimicrobial peptides, bradykinin-potentiating peptide, cysteine rich protein, serine proteinases, cathepsins, angiotensin-converting enzyme, endothelin-converting enzyme and chymotrypsin like protein, proteinases inhibitors, phospholipases and hyaluronidases., Conclusion: The present work shows that the venom composition of these two allopatric species of Tityus are considerably similar in terms of the major classes of proteins produced and secreted, although their individual toxin sequences are considerably divergent. These differences at amino acid level may reflect in different epitopes for the same protein classes in each species, explaining the basis for the poor recognition of T. obscurus venom by the antiserum raised against other species.

Ornithodoros mimon is an argasid tick that parasitizes bats, birds and opossums and is also harmful to humans. Knowledge of the transcripts present in the tick gut helps in understanding the role of vital molecules in the digestion process and parasite-host relationship, while also providing information about the evolution of arthropod hematophagy. Thus, the present study aimed to know and ascertain the main molecules expressed in the gut of argasid after their blood meal, through analysis on the gut transcriptome of engorged females of O. mimon using 454-based RNA sequencing. The gut transcriptome analysis reveals several transcripts associated with hemoglobin digestion, such as serine, cysteine, aspartic proteases and metalloenzymes. The phylogenetic analysis on the peptidases confirmed that most of them are clustered with other tick genes. We recorded the presence a cathepsin O peptidase-coding transcript in ticks. The topology of the phylogenetic inferences, based on transcripts of inferred families of homologues, was similar to that of previous reports based on mitochondrial genome and nuclear rRNA sequences. We deposited 2,213 sequence of O. mimon to the public databases. Our findings may help towards better understanding of important argasid metabolic processes, such as digestion, nutrition and immunity.

The complete mitochondrial genome, containing 17,526 bp, was determined from the pitviper Bothrops jararaca . It is the first mitogenome for the most medically important genus of snake in Latin America. This mitogenome has common snake mitochondrial features such as a duplicated control region that has nearly identical sequences at two different locations of the mitogenome and a translocation of tRNA-Leu (UUR). Besides, we found a translocation of the tRNA-Pro compared to Colubridae snakes. Finally, an unusual possible duplication containing a tRNA-Phe was observed for the first time and may represent a marker of the genus., Competing Interests: The authors have declared that no conflict of interest exists., (© 2016 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)