Is developmental plasticity triggered by DNA methylation changes in the invasive cane toad (Rhinella marina)?

Many organisms can adjust their development according to environmental conditions, including the presence of conspecifics. Although this developmental plasticity is common in amphibians, its underlying molecular mechanisms remain largely unknown. Exposure during development to either 'cannibal cues' from older conspecifics, or 'alarm cues' from injured conspecifics, causes reduced growth and survival in cane toad (Rhinella marina) tadpoles. Epigenetic modifications, such as changes in DNA methylation patterns, are a plausible mechanism underlying these developmental plastic responses. Here we tested this hypothesis, and asked whether cannibal cues and alarm cues trigger the same DNA methylation changes in developing cane toads. We found that exposure to both cannibal cues and alarm cues was associated with local changes in DNA methylation patterns. These DNA methylation changes affected genes putatively involved in developmental processes, but in different genomic regions for different conspecific-derived cues. Genetic background explains most of the epigenetic variation among individuals. Overall, the molecular mechanisms triggered by exposure to cannibal cues seem to differ from those triggered by alarm cues. Studies linking epigenetic modifications to transcriptional activity are needed to clarify the proximate mechanisms that regulate developmental plasticity in cane toads.

Effects of Psychotropic Drugs on Ribosomal Genes and Protein Synthesis.

Altered protein synthesis has been implicated in the pathophysiology of several neuropsychiatric disorders, particularly schizophrenia. Ribosomes are the machinery responsible for protein synthesis. However, there remains little information on whether current psychotropic drugs affect ribosomes and contribute to their therapeutic effects. We treated human neuronal-like (NT2-N) cells with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), valproate (0.5 mM) or vehicle control for 24 h. Transcriptomic and gene set enrichment analysis (GSEA) identified that the ribosomal pathway was altered by these drugs. We found that three of the eight drugs tested significantly decreased ribosomal gene expression, whilst one increased it. Most changes were observed in the components of cytosolic ribosomes and not mitochondrial ribosomes. Protein synthesis assays revealed that aripiprazole, clozapine and lithium all decreased protein synthesis. Several currently prescribed psychotropic drugs seem to impact ribosomal gene expression and protein synthesis. This suggests the possibility of using protein synthesis inhibitors as novel therapeutic agents for neuropsychiatric disorders.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina).

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

Integrative Analyses of Transcriptomes to Explore Common Molecular Effects of Antipsychotic Drugs.

There is little understanding of the underlying molecular mechanism(s) involved in the clinical efficacy of antipsychotics for schizophrenia. This study integrated schizophrenia-associated transcriptional perturbations with antipsychotic-induced gene expression profiles to detect potentially relevant therapeutic targets shared by multiple antipsychotics. Human neuronal-like cells (NT2-N) were treated for 24 h with one of the following antipsychotic drugs: amisulpride, aripiprazole, clozapine, risperidone, or vehicle controls. Drug-induced gene expression patterns were compared to schizophrenia-associated transcriptional data in post-mortem brain tissues. Genes regulated by each of four antipsychotic drugs in the reverse direction to schizophrenia were identified as potential therapeutic-relevant genes. A total of 886 genes were reversely expressed between at least one drug treatment (versus vehicle) and schizophrenia (versus healthy control), in which 218 genes were commonly regulated by all four antipsychotic drugs. The most enriched biological pathways include Wnt signaling and action potential regulation. The protein-protein interaction (PPI) networks found two main clusters having schizophrenia expression quantitative trait loci (eQTL) genes such as PDCD10, ANK2, and AKT3, suggesting further investigation on these genes as potential novel treatment targets.

Epithelial de-differentiation triggered by co-ordinate epigenetic inactivation of the EHF and CDX1 transcription factors drives colorectal cancer progression.

Colorectal cancers (CRCs) often display histological features indicative of aberrant differentiation but the molecular underpinnings of this trait and whether it directly drives disease progression is unclear. Here, we identify co-ordinate epigenetic inactivation of two epithelial-specific transcription factors, EHF and CDX1, as a mechanism driving differentiation loss in CRCs. Re-expression of EHF and CDX1 in poorly-differentiated CRC cells induced extensive chromatin remodelling, transcriptional re-programming, and differentiation along the enterocytic lineage, leading to reduced growth and metastasis. Strikingly, EHF and CDX1 were also able to reprogramme non-colonic epithelial cells to express colonic differentiation markers. By contrast, inactivation of EHF and CDX1 in well-differentiated CRC cells triggered tumour de-differentiation. Mechanistically, we demonstrate that EHF physically interacts with CDX1 via its PNT domain, and that these transcription factors co-operatively drive transcription of the colonic differentiation marker, VIL1. Compound genetic deletion of Ehf and Cdx1 in the mouse colon disrupted normal colonic differentiation and significantly enhanced colorectal tumour progression. These findings thus reveal a novel mechanism driving epithelial de-differentiation and tumour progression in CRC.

Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia.

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk-taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads' dispersal ability. Common-garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal-related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA-seq to gain a better understanding of the molecular mechanisms underlying dispersal-related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

Biological Mechanism(s) Underpinning the Association between Antipsychotic Drugs and Weight Gain.

Weight gain and consequent metabolic alterations are common side-effects of many antipsychotic drugs. Interestingly, several studies have suggested that improvement in symptoms and adverse metabolic effects are correlated. We used next generation sequencing data from NT-2 (human neuronal) cells treated with aripiprazole, amisulpride, risperidone, quetiapine, clozapine, or vehicle control, and compared with the Pillinger P-score (ranked from 0 to 1, indicating greater increase in weight gain and related metabolic parameters) to identify the genes most associated with the drugs' propensity to cause weight gain. The top 500 genes ranked for their correlation with the drugs' propensity to cause weight gain were subjected to pathway analysis using DAVID (NIH). We further investigated transcription factors (TFs) that are more likely to regulate the genes involved in these processes using the prediction tool of key TFs from TRRUST. The results suggest an enrichment for genes involved in lipid biosynthesis and metabolism, which are of interest for mechanisms underpinning weight-gain. The list of genes involved in the lipid pathways that correlated with weight gain was enriched for genes transcriptionally regulated by SREBF1 and SREBF2. Furthermore, quetiapine significantly increased the expression of SREBF1 and SREBF2 in NT-2 cells. Our results suggest that the effects of these antipsychotic drugs on lipid metabolism may be mediated, at least in part, via regulation of SREBF1/SREBF2 expression, with evidence of a direct effect of quetiapine on the expression of SREBF1/2. The effects of antipsychotic drugs on lipid metabolism may influence white matter structure (therapeutic effect) and the risk of weight gain, lipid disturbances, and, consequently, metabolic syndrome (adverse effects). Understanding the different molecular effects of these drugs could inform a personalized medicine approach in treating patients with schizophrenia.

Discovery of Novel Viruses Associated With the Invasive Cane Toad (Rhinella marina) in Its Native and Introduced Ranges.

Cane toads (Rhinella marina) are notoriously successful invaders: from 101 individuals brought to Australia in 1935, poisonous toads now cover an area >1.2 million km2 with adverse effects on native fauna. Despite extensive research on the role of macroparasites in cane toad invasion, viral research is lagging. We compared viral prevalence and diversity between toads in their native range (French Guiana, n=25) and two introduced ranges: Australia (n=151) and Hawai'i (n=10) with a metatranscriptomic and metagenomic approach combined with PCR screening. Australian toads almost exclusively harbor one of seven viruses detected globally. Rhimavirus-A (Picornaviridae) exhibited low genetic diversity and likely actively infected 9% of sampled Australian toads extending across ~2,000km of Northern Australia and up to the current invasion front. In native range cane toads, we identified multiple phylogenetically distinct viruses (Iridoviridae, Picornaviridae, Papillomaviridae, and Nackedna-like virus). None of the same viruses was detected in both ranges, suggesting that Australian cane toads have largely escaped the viral infection experienced by their native range counterparts. The novel native range viruses described here are potential biocontrol agents, as Australian toads likely lack prior immunological exposure to these viruses. Overall, our evidence suggests that there may be differences between viruses infecting cane toads in their native vs. introduced ranges, which lays the groundwork for further studies on how these viruses have influenced the toads' invasion history.

Transcriptional Modulation of the Hippo Signaling Pathway by Drugs Used to Treat Bipolar Disorder and Schizophrenia.

Recent reports suggest a link between positive regulation of the Hippo pathway with bipolar disorder (BD), and the Hippo pathway is known to interact with multiple other signaling pathways previously associated with BD and other psychiatric disorders. In this study, neuronal-like NT2 cells were treated with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), valproate (0.5 mM), or vehicle control for 24 h. Genome-wide mRNA expression was quantified and analyzed using gene set enrichment analysis (GSEA), with genes belonging to Hippo, Wnt, Notch, TGF- ß, and Hedgehog retrieved from the KEGG database. Five of the eight drugs downregulated the genes of the Hippo pathway and modulated several genes involved in the interacting pathways. We speculate that the regulation of these genes, especially by aripiprazole, clozapine, and quetiapine, results in a reduction of MAPK and NFκB pro-inflammatory signaling through modulation of Hippo, Wnt, and TGF-ß pathways. We also employed connectivity map analysis to identify compounds that act on these pathways in a similar manner to the known psychiatric drugs. Thirty-six compounds were identified. The presence of antidepressants and antipsychotics validates our approach and reveals possible new targets for drug repurposing.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species.

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

Effects of psychoactive drugs on cellular bioenergetic pathways.

OBJECTIVES: To investigate the actions of lithium, valproate, lamotrigine and quetiapine on bioenergetic pathways in cultured NT2-N neuronal-like cells and C8-B4 microglial cells. METHODS: NT2-N and C8-B4 cells were cultured and treated with lithium (2.5 mM), valproate (0.5 mM), quetiapine (0.05 mM) or lamotrigine (0.05 mM) for 24 hours. Gene expression and the mitochondrial bioenergetic profile were measured in both cell lines. RESULTS: In NT2-N cells, valproate increased oxidative phosphorylation (OXPHOS) gene expression, mitochondrial uncoupling and maximal respiratory capacity, while quetiapine decreased OXPHOS gene expression and respiration linked to ATP turnover, as well as decreasing the expression of genes in the citric acid cycle. Lamotrigine decreased OXPHOS gene expression but had no effect on respiration, while lithium reduced the expression of genes in the citric acid cycle. In C8-B4 cells, valproate and lithium increased OXPHOS gene expression, and valproate increased basal respiratory rate and maximal and spare respiratory capacities. In contrast, quetiapine significantly reduced basal respiratory rate and maximal and spare respiratory capacities. CONCLUSIONS: Overall our data suggest that some drugs used to treat neuropsychiatric and affective disorders have actions on a range of cellular bioenergetic processes, which could impact their effects in patients.

Transcriptional Effects of Psychoactive Drugs on Genes Involved in Neurogenesis.

Although neurogenesis is affected in several psychiatric diseases, the effects and mechanisms of action of psychoactive drugs on neurogenesis remain unknown and/or controversial. This study aims to evaluate the effects of psychoactive drugs on the expression of genes involved in neurogenesis. Neuronal-like cells (NT2-N) were treated with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), or valproate (0.5 mM) for 24 h. Genome wide mRNA expression was quantified and analysed using gene set enrichment analysis, with the neurogenesis gene set retrieved from the Gene Ontology database and the Mammalian Adult Neurogenesis Gene Ontology (MANGO) database. Transcription factors that are more likely to regulate these genes were investigated to better understand the biological processes driving neurogenesis. Targeted metabolomics were performed using gas chromatography-mass spectrometry. Six of the eight drugs decreased the expression of genes involved in neurogenesis in both databases. This suggests that acute treatment with these psychoactive drugs negatively regulates the expression of genes involved in neurogenesis in vitro. SOX2 and three of its target genes (CCND1, BMP4, and DKK1) were also decreased after treatment with quetiapine. This can, at least in part, explain the mechanisms by which these drugs decrease neurogenesis at a transcriptional level in vitro. These results were supported by the finding of increased metabolite markers of mature neurons following treatment with most of the drugs tested, suggesting increased proportions of mature relative to immature neurons consistent with reduced neurogenesis.

Increased Adaptive Variation Despite Reduced Overall Genetic Diversity in a Rapidly Adapting Invader.

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this "genetic paradox". The invasive cane toad (Rhinella marina) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations.

A field guide for the compositional analysis of any-omics data.

BACKGROUND: Next-generation sequencing (NGS) has made it possible to determine the sequence and relative abundance of all nucleotides in a biological or environmental sample. A cornerstone of NGS is the quantification of RNA or DNA presence as counts. However, these counts are not counts per se: their magnitude is determined arbitrarily by the sequencing depth, not by the input material. Consequently, counts must undergo normalization prior to use. Conventional normalization methods require a set of assumptions: they assume that the majority of features are unchanged and that all environments under study have the same carrying capacity for nucleotide synthesis. These assumptions are often untestable and may not hold when heterogeneous samples are compared. RESULTS: Methods developed within the field of compositional data analysis offer a general solution that is assumption-free and valid for all data. Herein, we synthesize the extant literature to provide a concise guide on how to apply compositional data analysis to NGS count data. CONCLUSIONS: In highlighting the limitations of total library size, effective library size, and spike-in normalizations, we propose the log-ratio transformation as a general solution to answer the question, "Relative to some important activity of the cell, what is changing?"

Chromosome-Level Alpaca Reference Genome VicPac3.1 Improves Genomic Insight Into the Biology of New World Camelids.

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex (MHC), the region involved in the Minute Chromosome Syndrome (MCS) and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

Immune and environment-driven gene expression during invasion: An eco-immunological application of RNA-Seq.

Host-pathogen associations change rapidly during a biological invasion and are predicted to impose strong selection on immune function. It has been proposed that the invader may experience an abrupt reduction in pathogen-mediated selection ("enemy release"), thereby favoring decreased investment into "costly" immune responses. Across plants and animals, there is mixed support for this prediction. Pathogens are not the only form of selection imposed on invaders; differences in abiotic environmental conditions between native and introduced ranges are also expected to drive rapid evolution. Here, we use RNA-Seq to assess the expression patterns of immune and environmentally associated genes in the cane toad (Rhinella marina) across its invasive Australian range. Transcripts encoding mediators of costly immune responses (inflammation, cytotoxicity) showed a curvilinear relationship with invasion history, with highest expression in toads from oldest and newest colonized areas. This pattern is surprising given theoretical expectations of density dynamics in invasive species and may be because density influences both intraspecific competition and parasite transmission, generating conflicting effects on the strength of immune responses. Alternatively, this expression pattern may be the result of other evolutionary forces, such as spatial sorting and genetic drift, working simultaneously with natural selection. Our findings do not support predictions about immune function based on the enemy release hypothesis and suggest instead that the effects of enemy release are difficult to isolate in wild populations, especially in the absence of information regarding parasite and pathogen infection. Additionally, expression patterns of genes underlying putatively environmentally associated traits are consistent with previous genetic studies, providing further support that Australian cane toads have adapted to novel abiotic challenges.

Draft genome assembly of the invasive cane toad, Rhinella marina.

Background: The cane toad (Rhinella marina formerly Bufo marinus) is a species native to Central and South America that has spread across many regions of the globe. Cane toads are known for their rapid adaptation and deleterious impacts on native fauna in invaded regions. However, despite an iconic status, there are major gaps in our understanding of cane toad genetics. The availability of a genome would help to close these gaps and accelerate cane toad research. Findings: We report a draft genome assembly for R. marina, the first of its kind for the Bufonidae family. We used a combination of long-read Pacific Biosciences RS II and short-read Illumina HiSeq X sequencing to generate 359.5 Gb of raw sequence data. The final hybrid assembly of 31,392 scaffolds was 2.55 Gb in length with a scaffold N50 of 168 kb. BUSCO analysis revealed that the assembly included full length or partial fragments of 90.6% of tetrapod universal single-copy orthologs (n = 3950), illustrating that the gene-containing regions have been well assembled. Annotation predicted 25,846 protein coding genes with similarity to known proteins in Swiss-Prot. Repeat sequences were estimated to account for 63.9% of the assembly. Conclusions: The R. marina draft genome assembly will be an invaluable resource that can be used to further probe the biology of this invasive species. Future analysis of the genome will provide insights into cane toad evolution and enrich our understanding of their interplay with the ecosystem at large.

Benchmarking differential expression analysis tools for RNA-Seq: normalization-based vs. log-ratio transformation-based methods.

BACKGROUND: Count data generated by next-generation sequencing assays do not measure absolute transcript abundances. Instead, the data are constrained to an arbitrary "library size" by the sequencing depth of the assay, and typically must be normalized prior to statistical analysis. The constrained nature of these data means one could alternatively use a log-ratio transformation in lieu of normalization, as often done when testing for differential abundance (DA) of operational taxonomic units (OTUs) in 16S rRNA data. Therefore, we benchmark how well the ALDEx2 package, a transformation-based DA tool, detects differential expression in high-throughput RNA-sequencing data (RNA-Seq), compared to conventional RNA-Seq methods such as edgeR and DESeq2. RESULTS: To evaluate the performance of log-ratio transformation-based tools, we apply the ALDEx2 package to two simulated, and two real, RNA-Seq data sets. One of the latter was previously used to benchmark dozens of conventional RNA-Seq differential expression methods, enabling us to directly compare transformation-based approaches. We show that ALDEx2, widely used in meta-genomics research, identifies differentially expressed genes (and transcripts) from RNA-Seq data with high precision and, given sufficient sample sizes, high recall too (regardless of the alignment and quantification procedure used). Although we show that the choice in log-ratio transformation can affect performance, ALDEx2 has high precision (i.e., few false positives) across all transformations. Finally, we present a novel, iterative log-ratio transformation (now implemented in ALDEx2) that further improves performance in simulations. CONCLUSIONS: Our results suggest that log-ratio transformation-based methods can work to measure differential expression from RNA-Seq data, provided that certain assumptions are met. Moreover, these methods have very high precision (i.e., few false positives) in simulations and perform well on real data too. With previously demonstrated applicability to 16S rRNA data, ALDEx2 can thus serve as a single tool for data from multiple sequencing modalities.

Understanding sequencing data as compositions: an outlook and review.

Motivation: Although seldom acknowledged explicitly, count data generated by sequencing platforms exist as compositions for which the abundance of each component (e.g. gene or transcript) is only coherently interpretable relative to other components within that sample. This property arises from the assay technology itself, whereby the number of counts recorded for each sample is constrained by an arbitrary total sum (i.e. library size). Consequently, sequencing data, as compositional data, exist in a non-Euclidean space that, without normalization or transformation, renders invalid many conventional analyses, including distance measures, correlation coefficients and multivariate statistical models. Results: The purpose of this review is to summarize the principles of compositional data analysis (CoDA), provide evidence for why sequencing data are compositional, discuss compositionally valid methods available for analyzing sequencing data, and highlight future directions with regard to this field of study. Supplementary information: Supplementary data are available at Bioinformatics online.

Improving amphibian genomic resources: a multitissue reference transcriptome of an iconic invader.

Background: Cane toads (Rhinella marina) are an iconic invasive species introduced to 4 continents and well utilized for studies of rapid evolution in introduced environments. Despite the long introduction history of this species, its profound ecological impacts, and its utility for demonstrating evolutionary principles, genetic information is sparse. Here we produce a de novo transcriptome spanning multiple tissues and life stages to enable investigation of the genetic basis of previously identified rapid phenotypic change over the introduced range. Findings: Using approximately 1.9 billion reads from developing tadpoles and 6 adult tissue-specific cDNA libraries, as well as a transcriptome assembly pipeline encompassing 100 separate de novo assemblies, we constructed 62 202 transcripts, of which we functionally annotated ∼50%. Our transcriptome assembly exhibits 90% full-length completeness of the Benchmarking Universal Single-Copy Orthologs data set. Robust assembly metrics and comparisons with several available anuran transcriptomes and genomes indicate that our cane toad assembly is one of the most complete anuran genomic resources available. Conclusions: This comprehensive anuran transcriptome will provide a valuable resource for investigation of genes under selection during invasion in cane toads, but will also greatly expand our general knowledge of anuran genomes, which are underrepresented in the literature. The data set is publically available in NCBI and GigaDB to serve as a resource for other researchers.