A practical guide to build de-novo assemblies for single tissues of non-model organisms: the example of a Neotropical frog

De novo transcriptome assembly and identification of G-Protein-Coupled-Receptors (GPCRs) in two species of monogenean parasites of fish

Genomic resources for Platyhelminthes of the class Monogenea are scarce, despite the diversity of these parasites, some species of which are highly pathogenic to their fish hosts. This work aimed to generate de novo-assembled transcriptomes of two monogenean species, Scutogyrus longicornis (Dactylogyridae) and Rhabdosynochus viridisi (Diplectanidae), providing a protocol for cDNA library preparation with low input samples used in single cell transcriptomics. This allowed us to work with sub-microgram amounts of total RNA with success. These transcriptomes consist of 25,696 and 47,187 putative proteins, respectively, which were further annotated according to the Swiss-Prot, Pfam, GO, KEGG, and COG databases. The completeness values of these transcriptomes evaluated with BUSCO against Metazoa databases were 54.1% and 73%, respectively, which is in the range of other monogenean species. Among the annotations, a large number of terms related to G-protein-coupled receptors (GPCRs) were found. We identified 109 GPCR-like sequences in R. viridisi, and 102 in S. longicornis, including family members specific for Platyhelminthes. Rhodopsin was the largest family according to GRAFS classification. Two putative melatonin receptors found in S. longicornis represent the first record of this group of proteins in parasitic Platyhelminthes. Forty GPCRs of R. viridisi and 32 of S. longicornis that were absent in Vertebrata might be potential drug targets. The present study provides the first publicly available transcriptomes for monogeneans of the subclass Monopisthocotylea, which can serve as useful genomic datasets for functional genomic research of this important group of parasites.

De novo assembly and functional annotation of blood transcriptome of loggerhead turtle, and in silico characterization of peroxiredoxins and thioredoxins

The aim of this study was to generate and analyze the atlas of the loggerhead turtle blood transcriptome by RNA-seq, as well as identify and characterize thioredoxin (Tnxs) and peroxiredoxin (Prdxs) antioxidant enzymes of the greatest interest in the control of peroxide levels and other biological functions. The transcriptome of loggerhead turtle was sequenced using the Illumina Hiseq 2000 platform and de novo assembly was performed using the Trinity pipeline. The assembly comprised 515,597 contigs with an N50 of 2,631 bp. Contigs were analyzed with CD-Hit obtaining 374,545 unigenes, of which 165,676 had ORFs encoding putative proteins longer than 100 amino acids. A total of 52,147 (31.5%) of these transcripts had significant homology matches in at least one of the five databases used. From the enrichment of GO terms, 180 proteins with antioxidant activity were identified, among these 28 Prdxs and 50 putative Tnxs. The putative proteins of loggerhead turtles encoded by the genes Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip were predicted and characterized in silico. When comparing Prdxs and Txns of loggerhead turtle with homologous human proteins, they showed 18 (9%), 52 (18%) 94 (43%), 36 (16%), 35 (33%) and 74 (19%) amino acid mutations respectively. However, they showed high conservation in active sites and structural motifs (98%), with few specific modifications. Of these, Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip presented 0, 25, 18, three, six and two deleterious changes. This study provides a high quality blood transcriptome and functional annotation of loggerhead sea turtles.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

A tree frog (Boana pugnax) dataset of skin transcriptome for the identification of biomolecules with potential antimicrobial activities

Increases in the prevalence of multiply resistant microbes have necessitated the search for new molecules with antimicrobial properties. One noteworthy avenue in this search is inspired by the presence of native antimicrobial peptides in the skin of amphibians. Having the second highest diversity of frogs worldwide, Colombian anurans represent an extensive natural reservoir that could be tapped in this search. Among this diversity, species such as Boana pugnax (the Chirique-Flusse Treefrog) are particularly notable, in that they thrive in a diversity of marginal habitats, utilize both aquatic and arboreal habitats, and are members of one of few genera that are known to mount a robust immunological response against the fungus Batrachochytrium dendrobatidis, which has decimated the population of frogs worldwide. To search for molecules with potential antimicrobial activity, we have assembled and annotated a reference transcriptome from the skin of four wild captured B. pugnax from Antioquia, Colombia. Analysis of potential antimicrobial and immunological components was performed using ontology analyses, we identified several antimicrobial chemokines with particularly strong potential for exhibiting broadscale antimicrobial activities, as well as several genes related to rapid alteration of transcriptional (KRAB zinc finger protein) and phosphorylation (MAPK) responses to exogenous stressors. We also found eight families of transmembrane transport proteins, including sodium, potassium and voltage-dependent calcium channels, which will be invaluable in future studies aimed at more precisely defining the diversity and function of cationic antimicrobial peptides with alpha-helical structures. These data highlight the utility of frogs such as Boana pugnax in the search of new antimicrobial molecules. Moreover, the molecular datasets presented here allow us to expand our knowledge of this species and illustrate the importance of preserving the vast potential of Colombian biodiversity for the identification of useful biomolecules.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Methodological Insight Into Mosquito Microbiome Studies

Symbiotic bacteria affect competence for pathogen transmission in insect vectors, including mosquitoes. However, knowledge on mosquito-microbiome-pathogen interactions remains limited, largely due to methodological reasons. The current, cost-effective practice of sample pooling used in mosquito surveillance and epidemiology prevents correlation of individual traits (i.e., microbiome profile) and infection status. Moreover, many mosquito studies employ laboratory-reared colonies that do not necessarily reflect the natural microbiome composition and variation in wild populations. As a consequence, epidemiological and microbiome studies in mosquitoes are to some extent uncoupled, and the interactions among pathogens, microbiomes, and natural mosquito populations remain poorly understood. This study focuses on the effect the pooling practice poses on mosquito microbiome profiles, and tests different approaches to find an optimized low-cost methodology for extensive sampling while allowing for accurate, individual-level microbiome studies. We tested the effect of pooling by comparing wild-caught, individually processed mosquitoes with pooled samples. With individual mosquitoes, we also tested two methodological aspects that directly affect the cost and feasibility of broad-scale molecular studies: sample preservation and tissue dissection. Pooling affected both alpha- and beta-diversity measures of the microbiome, highlighting the importance of using individual samples when possible. Both RNA and DNA yields were higher when using inexpensive reagents such as NAP (nucleic acid preservation) buffer or absolute ethanol, without freezing for short-term storage. Microbiome alpha- and beta-diversity did not show overall significant differences between the tested treatments compared to the controls (freshly extracted samples or dissected guts). However, the use of standardized protocols is highly recommended to avoid methodological bias in the data.

Next generation sequencing and RNA-seq characterization of adipose tissue in the Nile crocodile (Crocodylus niloticus) in South Africa: Possible mechanism(s) of pathogenesis and pathophysiology of pansteatitis

BACKGROUND

Concerted efforts to identify the pathogenesis and mechanism(s) involved in pansteatitis, (a generalized inflammation of the adipose tissue), that was attributed to the recent crocodile die off in the Olifants River and Loskop Dam in Kruger National Park, Mpumalanga, South Africa have been in the forefront of research in recent time. As part of the efforts, molecular characterization of healthy and pansteatitis adipose tissue was carried out by RNA sequencing (RNA-Seq) using Next Generation Sequencing (NGS) and de novo assembly of the adipose transcriptome, followed by differential gene expression analysis.

METHODOLOGY

Healthy adipose tissue consisting of fifty samples was collected from the subcutaneous, visceral, intermuscular adipose tissues and the abdominal fat body of ten 4 years old juvenile crocodiles from a local crocodile farm in Pretoria, South Africa. Ten pansteatitis samples were collected from visceral and intermuscular adipose tissues of five crocodiles that were dying of pansteatitis.

RESULTS

Forty-two thousand, two hundred and one (42,201) transcripts were assembled, out of which 37, 835 had previously been characterized. The de novo assembled transcriptome had an N50 (average sequence) of 436 bp, percentage GC content of 43.92, which compared well with previously assembled transcripts in the saltwater crocodile. Seventy genes were differentially expressed and upregulated in pansteatitis. These included genes coding for extracellular matrix (ECM) signaling ligands, inflammatory cytokines and tumour necrosis factor alpha (TNFα) receptors, fatty acid synthase and fatty acid binding proteins, peroxisome proliferator-activated receptor gamma (PPARγ), nuclear factor and apoptosis signaling ligands, and mitogen activated protein kinase enzymes among others. Majority (88.6%) of the upregulated genes were found to be involved in hypoxia inducible pathways for activation of NFkβ and inflammation, apoptosis, Toll-like receptor pathway and PPARγ. Bicaudal homologous 2 Drosophila gene (BICD2) associated with spinal and lower extremity muscle atrophy was also upregulated in pansteatitis while Sphingosine -1-phosphate phosphatase 2 (SGPP2) involved in Sphingosine -1- phosphate metabolism was downregulated. Futhermore, Doublesex-mab-related transcription factor 1 (DMRT1) responsible for sex gonad development and germ cell differentiation was also downregulated.

CONCLUSION

Thus, from the present study, based on differentially expressed genes in pansteatitis, affected Nile crocodiles might have died partly due to their inability to utilize stored triglycerides as a result of inflammation induced insulin resistance, leading to starvation in the midst of plenty. Affected animals may have also suffered muscular atrophy of the lower extremities and poor fertility.