The Oyster River Protocol: a multi-assembler and kmer approach for de novo transcriptome assembly

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

Characterization of the transcriptional effects of the plastic additive dibutyl phthalate alone and in combination with microplastic on the green-lipped mussel Perna canaliculus

The presence and persistence of microplastics (MPs) in diverse aquatic environments are of global concern. Microplastics can impact marine organisms via direct physical interaction and the release of potentially harmful chemical additives incorporated into the plastic. These chemicals are physically bound to the plastic matrix and can leach out. The hazards associated with chemical additives to exposed organisms is not well characterized. We investigated the hazards of plastic additives leaching from plastic. We used the common plasticizer dibutyl phthalate (DBP) as a chemical additive proxy and the New Zealand green-lipped mussel (Perna canaliculus) as a model. We used early-adult P. canaliculus exposed to combinations of virgin and DBP-spiked polyvinyl chloride (PVC), MPs, and DBP alone for 7 days. Whole transcriptome sequencing (RNA-seq) was conducted to assess whether leaching of DBP from MPs poses a hazard. The differences between groups were evaluated using pairwise permutational multivariate analysis of variance (PERMANOVA), and all treatments were significantly different from controls. In addition, a significant difference was seen between DBP and PVC MP treatment. Transcriptome analysis revealed that mussels exposed to DBP alone had the most differentially expressed genes (914), followed by PVC MP + DBP (448), and PVC MP (250). Gene ontology functional analysis revealed that the most enriched pathway types were in cellular metabolism, immune response, and endocrine disruption. Microplastic treatments enriched numerous pathways related to cellular metabolism and immune response. The combined exposure of PVC MP + DBP appears to cause combined effects, suggesting that DBP is bioavailable to the exposed mussels in the PVC MP + DBP treatment. Our results support the hypothesis that chemical additives are potentially an important driver of MP toxicity. Environ Toxicol Chem 2024;43:1604-1614. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Egg-laying hormone expression in identified neurons across developmental stages and reproductive states of the nudibranch Berghia stephanieae

Neuropeptides play essential roles in coordinating reproduction. Egg-laying hormone (ELH) is conserved in genetic sequence and behavioral function across molluscs, where neuronal clusters secrete ELH to modulate and induce egg-laying. Here we investigated ELH in the nudibranch mollusc, Berghia stephanieae. ELH preprohormone gene orthologs, which showed clade-specific differences at the C-terminus of the predicted bioactive peptide, were identified in brain transcriptomes across several nudipleuran species, including B. stephanieae. ELH shares deep homology with the corticotropin-releasing hormone gene family, which has roles broadly in stress response. Injection of synthesized B. stephanieae ELH peptide into mature individuals induced egg-laying. ELH gene expression in the brain and body was mapped using in-situ hybridization chain reaction. Across the adult brain, 300-400 neurons expressed ELH. Twenty-one different cell types were identified in adults, three of which were located unilaterally on the right side, which corresponds to the location of the reproductive organs. Ten cell types were present in pre-reproductive juvenile stages. An asymmetric cluster of approximately 100 small neurons appeared in the right pedal ganglion of late-stage juveniles. Additional neurons in the pleural and pedal ganglia expressed ELH only in adults that were actively laying eggs and sub-adults that were on the verge of doing so, implicating their direct role in reproduction. Outside the brain, ELH was expressed on sensory appendages, including in presumptive sensory neurons. Its widespread expression in the nudibranch B. stephanieae suggests that ELH plays a role beyond reproduction in gastropod molluscs.

Cellular-resolution gene expression mapping reveals organization in the head ganglia of the gastropod, Berghia stephanieae

Gastropod molluscs such as Aplysia, Lymnaea, and Tritonia have been important for determining fundamental rules of motor control, learning, and memory because of their large, individually identifiable neurons. Yet only a small number of gastropod neurons have known molecular markers, limiting the ability to establish brain-wide structure-function relations. Here we combine high-throughput, single-cell RNA sequencing with in situ hybridization chain reaction in the nudibranch Berghia stephanieae to identify and visualize the expression of markers for cell types. Broad neuronal classes were characterized by genes associated with neurotransmitters, like acetylcholine, glutamate, serotonin, and GABA, as well as neuropeptides. These classes were subdivided by other genes including transcriptional regulators and unannotated genes. Marker genes expressed by neurons and glia formed discrete, previously unrecognized regions within and between ganglia. This study provides the foundation for understanding the fundamental cellular organization of gastropod nervous systems.

Regulatory and evolutionary impact of DNA methylation in two songbird species and their naturally occurring F1 hybrids

BACKGROUND

Regulation of transcription by DNA methylation in 5'-CpG-3' context is a widespread mechanism allowing differential expression of genetically identical cells to persist throughout development. Consequently, differences in DNA methylation can reinforce variation in gene expression among cells, tissues, populations, and species. Despite a surge in studies on DNA methylation, we know little about the importance of DNA methylation in population differentiation and speciation. Here we investigate the regulatory and evolutionary impact of DNA methylation in five tissues of two Ficedula flycatcher species and their naturally occurring F1 hybrids.

RESULTS

We show that the density of CpG in the promoters of genes determines the strength of the association between DNA methylation and gene expression. The impact of DNA methylation on gene expression varies among tissues with the brain showing unique patterns. Differentially expressed genes between parental species are predicted by genetic and methylation differentiation in CpG-rich promoters. However, both these factors fail to predict hybrid misexpression suggesting that promoter mismethylation is not a main determinant of hybrid misexpression in Ficedula flycatchers. Using allele-specific methylation estimates in hybrids, we also determine the genome-wide contribution of cis- and trans effects in DNA methylation differentiation. These distinct mechanisms are roughly balanced in all tissues except the brain, where trans differences predominate.

CONCLUSIONS

Overall, this study provides insight on the regulatory and evolutionary impact of DNA methylation in songbirds.

Transcriptional Profiling of Early Defense Response to White Pine Blister Rust Infection in Pinus albicaulis (Whitebark Pine)

Pathogen perception generates the activation of signal transduction cascades to host defense. White pine blister rust (WPBR) is caused by Cronartium ribicola J.C. Fisch and affects a number of species of Pinus. One of the most severely affected species is Pinus albicaulis Engelm (whitebark pine). WPBR resistance in the species is a polygenic and complex trait that requires an optimized immune response. We identified early responses in 2-year-old seedlings after four days of fungal inoculation and compared the underlying transcriptomic response with that of healthy non-inoculated individuals. A de novo transcriptome assembly was constructed with 56,796 high quality-annotations derived from the needles of susceptible and resistant individuals in a resistant half-sib family. Differential expression analysis identified 599 differentially expressed transcripts, from which 375 were upregulated and 224 were downregulated in the inoculated seedlings. These included components of the initial phase of active responses to abiotic factors and stress regulators, such as those involved in the first steps of flavonoid biosynthesis. Four days after the inoculation, infected individuals showed an overexpression of chitinases, reactive oxygen species (ROS) regulation signaling, and flavonoid intermediates. Our research sheds light on the first stage of infection and emergence of disease symptoms among whitebark pine seedlings. RNA sequencing (RNA-seq) data encoding hypersensitive response, cell wall modification, oxidative regulation signaling, programmed cell death, and plant innate immunity were differentially expressed during the defense response against C. ribicola.

A genome sequence for the threatened whitebark pine

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

Phototactic preference and its genetic basis in the planulae of the colonial Hydrozoan Hydractinia symbiolongicarpus

Background

Marine organisms with sessile adults commonly possess motile larval stages that make settlement decisions based on integrating environmental sensory cues. Phototaxis, the movement toward or away from light, is a common behavioral characteristic of aquatic and marine metazoan larvae, and of algae, protists, and fungi. In cnidarians, behavioral genomic investigations of motile planulae larvae have been conducted in anthozoans (corals and sea anemones) and scyphozoans (true jellyfish), but such studies are presently lacking in hydrozoans. Here, we examined the behavioral genomics of phototaxis in planulae of the hydrozoan Hydractinia symbiolongicarpus.

Results

A behavioral phototaxis study of day 3 planulae indicated preferential phototaxis to green (523 nm) and blue (470 nm) wavelengths of light, but not red (625 nm) wavelengths. A developmental transcriptome study where planula larvae were collected from four developmental time points for RNA-seq revealed that many genes critical to the physiology and development of ciliary photosensory systems are dynamically expressed in planula development and correspond to the expression of phototactic behavior. Microscopical investigations using immunohistochemistry and in situ hybridization demonstrated that several transcripts with predicted function in photoreceptors, including cnidops class opsin, CNG ion channel, and CRX-like transcription factor, localize to ciliated bipolar sensory neurons of the aboral sensory neural plexus, which is associated with the direction of phototaxis and the site of settlement.

Conclusions

The phototactic preference displayed by planulae is consistent with the shallow sandy marine habitats they experience in nature. Our genomic investigations add further evidence of similarities between cnidops-mediated photoreceptors of hydrozoans and other cnidarians and ciliary photoreceptors as found in the eyes of humans and other bilaterians, suggesting aspects of their shared evolutionary history.

Fatty acid production and associated gene pathways are altered by increased salinity and dimethyl sulfoxide treatments during cryopreservation of Symbiodinium pilosum (Symbiodiniaceae)

The Symbiodinium genus is ancestral among other Symbiodiniaceae lineages with species that are both symbiotic and free living. Changes in marine ecosystems threaten their existence and crucial ecological roles. Cryopreservation offers an avenue for their long-term storage for future habitat restoration after coral bleaching. In our previous study we demonstrated that high salinity treatments of Symbiodiniaceae isolates led to changes in their fatty acid (FA) profiles and higher cell viabilities after cryopreservation. In this study, we investigated the role of increased salinity on FA production and the genes involved in FA biosynthesis and degradation pathways during the cryopreservation of Symbiodinium pilosum. Overall, there was a twofold increase in mass of FAs produced by S. pilosum after being cultured in medium with increased salinity (54 parts per thousand; ppt). Dimethyl sulfoxide (Me2SO) led to a ninefold increase of FAs in standard salinity (SS) treatment, compared to a fivefold increase in increased salinity (IS) treatments. The mass of the FA classes returned to baseline during recovery. Transcriptomic analyses showed an acyl carrier protein gene was significantly upregulated after Me2SO treatment in the SS cultures. Cytochrome P450 reductase genes were significantly down regulated after Me2SO addition in SS treatment preventing FA degradation. These changes in the expression of FA biosynthesis and degradation genes contributed to more FAs in SS treated isolates. Understanding how increased salinity changes FA production and the roles of specific genes in regulating FA pathways will help improve current freezing protocols for Symbiodiniaceae and other marine microalgae.

Differential responses of selectively bred mussels (Perna canaliculus) to heat stress-survival, immunology, gene expression and microbiome diversity

New Zealand's green-lipped mussel (Perna canaliculus) is an ecologically and economically important species. Marine heatwaves are increasing in frequency around NZ's coastline, and these events are correlated with increased stress and mortality of some aquaculture species. This study aimed to identify general biomarkers of heat stress in P. canaliculus and to assess whether responses differed between genetically distinct selectively bred mussels. We exposed three families of selectively bred mussels (families A, B and C) to three seawater temperature regimes in the laboratory: 1) a "control" treatment (ambient 12°C), 2) a 26°C heat challenge with a subsequent recovery period, and 3) a sustained 26°C heat challenge with no recovery. We investigated whether the survival, immune response (hemocyte concentration and viability, oxidative stress and total antioxidant capacity), hemocyte gene expression and gill microbiome differed between the families during the temperature challenges. In the sustained heat-stress treatment, family A had the highest survival rate (42% compared with 25% and 5% for families C and B, respectively). Gene expression levels significantly shifted during thermal stress and differed between families, with family A more dissimilar than families B and C. Family C had substantially more genes impacted by temperature treatment and timepoint than the other families, while family B had very little genes/pathways that responded to thermal stress. Genes related to heat shock proteins and immune responses (e.g., AIF1, CTSC, TOLL8, CASP9, FNTA, AHCY, CRYAB, PPIF) were upregulated in all families during heat stress. Microbiome species-richness differed between families before and during heat-stress, with family A having a distinctly different microbiome flora than the other families. Microbial diversity changed similarly in all families exposed to prolonged heat-stress, with species of Vibrio and Campylobacter increasing in these mussels. Our study highlights the use of non-lethal sampling of hemocytes as a diagnostic tool to explore the immune response and gene expression of selectively bred mussels, to predict their response to ocean warming. This approach can identify potential thermotolerant candidates for further selective breeding, which may increase the resilience of the mussel aquaculture industry in a warming ocean.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

A Genome Sequence for the Threatened Whitebark Pine

Whitebark pine (WBP, Pinus albicaulis ) is a white pine of subalpine regions in western contiguous US and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola ) and additional threats from mountain pine beetle ( Dendroctonus ponderosae ), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short-reads of haploid megametophyte tissue and Oxford Nanopore long-reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gbp of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gbp). Approximately 87.2% (24.0 Gbp) of total sequence was placed on the twelve WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich-repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the three subclasses of NLRs (TNL, CNL, RNL). Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo assembled transcriptomes.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality

BACKGROUND

Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected.

RESULTS

Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification.

CONCLUSIONS

Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Genome assembly of the foot-flagging frog, Staurois parvus: a resource for understanding mechanisms of behavior

Elaborate and skilled movements of the body have been selected in a variety of species as courtship and rivalry signals. One roadblock in studying these behaviors has been a lack of resources for understanding how they evolved at the genetic level. The Bornean rock frog (Staurois parvus) is an ideal species in which to address this issue. Males wave their hindlimbs in a "foot-flagging" display when competing for mates. The evolution of foot flagging in S. parvus and other species is accompanied by increases in the expression of the androgen receptor gene within its neuromuscular system, but it remains unclear what genetic or transcriptional changes are associated with this behavioral phenotype. We have now assembled the genome of S. parvus, resulting in 3.98 Gbp of 22,402 contigs with an N50 of 611,229 bp. The genome will be a resource for finding genes related to the physiology underlying foot flagging and to adaptations of the neuromuscular system. As a first application of the genome, we also began work in comparative genomics and differential gene expression analysis. We show that the androgen receptor is diverged from other anuran species, and we identify unique expression patterns of genes in the spinal cord and leg muscle that are important for axial patterning, cell specification and morphology, or muscle contraction. This genome will continue to be an important tool for future -omics studies to understand the evolution of elaborate signaling behaviors in this and potentially related species.

The evolution of antimicrobial peptides in Chiroptera

High viral tolerance coupled with an extraordinary regulation of the immune response makes bats a great model to study host-pathogen evolution. Although many immune-related gene gains and losses have been previously reported in bats, important gene families such as antimicrobial peptides (AMPs) remain understudied. We built an exhaustive bioinformatic pipeline targeting the major gene families of defensins and cathelicidins to explore AMP diversity and analyze their evolution and distribution across six bat families. A combination of manual and automated procedures identified 29 AMP families across queried species, with α-, β-defensins, and cathelicidins representing around 10% of AMP diversity. Gene duplications were inferred in both α-defensins, which were absent in five species, and three β-defensin gene subfamilies, but cathelicidins did not show significant shifts in gene family size and were absent in Anoura caudifer and the pteropodids. Based on lineage-specific gains and losses, we propose diet and diet-related microbiome evolution may determine the evolution of α- and β-defensins gene families and subfamilies. These results highlight the importance of building species-specific libraries for genome annotation in non-model organisms and shed light on possible drivers responsible for the rapid evolution of AMPs. By focusing on these understudied defenses, we provide a robust framework for explaining bat responses to pathogens.

Increasing temperature-driven changes in life history traits and gene expression of an Antarctic tardigrade species

The Antarctic region has been experiencing some of the planet's strongest climatic changes, including an expected increase of the land temperature. The potential effects of this warming trend will lead ecosystems to a risk of losing biodiversity. Antarctic mosses and lichens host different microbial groups, micro-arthropods and meiofaunal organisms (e.g., tardigrades, rotifers). The eutardigrade Acutuncus antarcticus is considered a model animal to study the effect of increasing temperature due to global warming on Antarctic terrestrial communities. In this study, life history traits and fitness of this species are analyzed by rearing specimens at two different and increasing temperatures (5°C vs. 15°C). Moreover, the first transcriptome analysis on A. antarcticus is performed, exposing adult animals to a gradual increase of temperature (5°C, 10°C, 15°C, and 20°C) to find differentially expressed genes under short- (1 day) and long-term (15 days) heat stress. Acutuncus antarcticus specimens reared at 5°C live longer (maximum life span: 686 days), reach sexual maturity later, lay more eggs (which hatch in longer time and in lower percentage) compared with animals reared at 15°C. The fitness decreases in animals belonging to the second generation at both rearing temperatures. The short-term heat exposure leads to significant changes at transcriptomic level, with 67 differentially expressed genes. Of these, 23 upregulated genes suggest alterations of mitochondrial activity and oxido-reductive processes, and two intrinsically disordered protein genes confirm their role to cope with heat stress. The long-term exposure induces alterations limited to 14 genes, and only one annotated gene is upregulated in response to both heat stresses. The decline in transcriptomic response after a long-term exposure indicates that the changes observed in the short-term are likely due to an acclimation response. Therefore, A. antarcticus could be able to cope with increasing temperature over time, including the future conditions imposed by global climate change.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Reverse engineering environmental metatranscriptomes clarifies best practices for eukaryotic assembly

BACKGROUND

Diverse communities of microbial eukaryotes in the global ocean provide a variety of essential ecosystem services, from primary production and carbon flow through trophic transfer to cooperation via symbioses. Increasingly, these communities are being understood through the lens of omics tools, which enable high-throughput processing of diverse communities. Metatranscriptomics offers an understanding of near real-time gene expression in microbial eukaryotic communities, providing a window into community metabolic activity.

RESULTS

Here we present a workflow for eukaryotic metatranscriptome assembly, and validate the ability of the pipeline to recapitulate real and manufactured eukaryotic community-level expression data. We also include an open-source tool for simulating environmental metatranscriptomes for testing and validation purposes. We reanalyze previously published metatranscriptomic datasets using our metatranscriptome analysis approach.

CONCLUSION

We determined that a multi-assembler approach improves eukaryotic metatranscriptome assembly based on recapitulated taxonomic and functional annotations from an in-silico mock community. The systematic validation of metatranscriptome assembly and annotation methods provided here is a necessary step to assess the fidelity of our community composition measurements and functional content assignments from eukaryotic metatranscriptomes.

Latitudinal variation and plasticity in response to temperature in Geukensia demissa

As global temperatures warm, species must adapt to a changing climate or transition to a different location suitable for their survival. Understanding the extent to which species are able to do so, particularly keystone species, is imperative to ensuring the survival of key ecosystems. The ribbed mussel Geukensia demissa is an integral part of salt marshes along the Atlantic coast of North America. Spatial patterns of genomic and phenotypic divergence have been previously documented, although their link with coastal environmental variation is unknown. Here, we study how populations of G. demissa in the northern (Massachusetts) and southern (Georgia) portions of the species range respond to changes in temperature. We combine assays of variation in oxygen consumption and RNA transcriptomic data with genomic divergence analyses to identify how separate populations of G. demissa may vary in distinct thermal environments. Our results show differences in constitutive oxygen consumption between mussels from Georgia and Massachusetts, as well as shared and disparate patterns of gene expression across temperature profiles. We also find that metabolic genes seem to be a strong component of divergence between these two populations. Our analysis highlights the importance of studying integrative patterns of genomic and phenotypic variation in species that are key for particular ecosystems, and how they might respond to further changes in climate.

Killer Knots: Molecular Evolution of Inhibitor Cystine Knot Toxins in Wandering Spiders (Araneae: Ctenidae)

Venom expressed by the nearly 50,000 species of spiders on Earth largely remains an untapped reservoir of a diverse array of biomolecules with potential for pharmacological and agricultural applications. A large fraction of the noxious components of spider venoms are a functionally diverse family of structurally related polypeptides with an inhibitor cystine knot (ICK) motif. The cysteine-rich nature of these toxins makes structural elucidation difficult, and most studies have focused on venom components from the small handful of medically relevant spider species such as the highly aggressive Brazilian wandering spider Phoneutria nigriventer. To alleviate difficulties associated with the study of ICK toxins in spiders, we devised a comprehensive approach to explore the evolutionary patterns that have shaped ICK functional diversification using venom gland transcriptomes and proteomes from phylogenetically distinct lineages of wandering spiders and their close relatives. We identified 626 unique ICK toxins belonging to seven topological elaborations. Phylogenetic tests of episodic diversification revealed distinct regions between cysteine residues that demonstrated differential evidence of positive or negative selection, which may have structural implications towards the specificity and efficacy of these toxins. Increased taxon sampling and whole genome sequencing will provide invaluable insights to further understand the evolutionary processes that have given rise to this diverse class of toxins.

Ecological constraints on highly evolvable olfactory receptor genes and morphology in neotropical bats

Although evolvability of genes and traits may promote specialization during species diversification, how ecology subsequently restricts such variation remains unclear. Chemosensation requires animals to decipher a complex chemical background to locate fitness-related resources, and thus the underlying genomic architecture and morphology must cope with constant exposure to a changing odorant landscape; detecting adaptation amidst extensive chemosensory diversity is an open challenge. In phyllostomid bats, an ecologically diverse clade that evolved plant visiting from a presumed insectivorous ancestor, the evolution of novel food detection mechanisms is suggested to be a key innovation, as plant-visiting species rely strongly on olfaction, supplementarily using echolocation. If this is true, exceptional variation in underlying olfactory genes and phenotypes may have preceded dietary diversification. We compared olfactory receptor (OR) genes sequenced from olfactory epithelium transcriptomes and olfactory epithelium surface area of bats with differing diets. Surprisingly, although OR evolution rates were quite variable and generally high, they are largely independent of diet. Olfactory epithelial surface area, however, is relatively larger in plant-visiting bats and there is an inverse relationship between OR evolution rates and surface area. Relatively larger surface areas suggest greater reliance on olfactory detection and stronger constraint on maintaining an already diverse OR repertoire. Instead of the typical case in which specialization and elaboration are coupled with rapid diversification of associated genes, here the relevant genes are already evolving so quickly that increased reliance on smell has led to stabilizing selection, presumably to maintain the ability to consistently discriminate among specific odorants-a potential ecological constraint on sensory evolution.

De novo assembly and annotation of the transcriptome of the endangered seagrass Zostera capensis: Insights from differential gene expression under thermal stress

Seagrasses are important marine ecosystem engineers but anthropogenic impacts and climate change have led to numerous population declines globally. In South Africa, Zostera capensis is endangered due to fragmented populations and heavy anthropogenic pressures on estuarine ecosystems that house the core of the populations. Addressing questions of how pressures such as climate change affect foundational species, including Z. capensis are crucial to supporting their conservation and underpin restoration efforts. Here we use ecological transcriptomics to study key functional responses of Z. capensis through quantification of gene expression after thermal stress and present the first reference transcriptome of Z. capensis. Four de novo reference assemblies (Trinity, IDBA-tran, RNAspades, SOAPdenovo) filtered through the EvidentialGene pipeline resulted in 153,755 transcripts with a BUSCO score of 66.1% for completeness. Differential expression analysis between heat stressed (32 °C for three days) and pre-warming plants identified genes involved in photosynthesis, oxidative stress, translation, metabolic and biosynthetic processes in the Z. capensis thermal stress response. This reference transcriptome is a significant contribution to the limited available genomic resources for Z. capensis and represents a vital tool for addressing questions around the species restoration and potential functional responses to warming marine environments.

Characterizing Hox genes in mayflies (Ephemeroptera), with Hexagenia limbata as a new mayfly model

BACKGROUND

Hox genes are key regulators of appendage development in the insect body plan. The body plan of mayfly (Ephemeroptera) nymphs differs due to the presence of abdominal appendages called gills. Despite mayflies' phylogenetic position in Paleoptera and novel morphology amongst insects, little is known of their developmental genetics, such as the appendage-regulating Hox genes. To address this issue we present an annotated, early instar transcriptome and embryonic expression profiles for Antennapedia, Ultrabithorax, and Abdominal A proteins in the mayfly Hexagenia limbata, identify putative Hox protein sequences in the mayflies H. limbata, Cloeon dipterum, and Ephemera danica, and describe the genomic organization of the Hox gene cluster in E. danica.

RESULTS

Transcriptomic sequencing of early instar H. limbata nymphs yielded a high-quality assembly of 83,795 contigs, of which 22,975 were annotated against Folsomia candida, Nilaparvata lugens, Zootermopsis nevadensis and UniRef90 protein databases. Homeodomain protein phylogeny and peptide annotations identified coding sequences for eight of the ten canonical Hox genes (excluding zerknüllt/Hox3 and fushi tarazu) in H. limbata and C. dipterum, and all ten in E. danica. Mayfly Hox protein sequences and embryonic expression patterns of Antp, Ubx, and Abd-A appear highly conserved with those seen in other non-holometabolan insects. Similarly, the genomic organization of the Hox cluster in E. danica resembles that seen in most insects.

CONCLUSIONS

We present evidence that mayfly Hox peptide sequences and the embryonic expression patterns for Antp, Ubx, and Abd-A are extensively conserved with other insects, as is organization of the mayfly Hox gene cluster. The protein data suggest mayfly Antp, Ubx, and Abd-A play appendage promoting and repressing roles during embryogenesis in the thorax and abdomen, respectively, as in other insects. The identified expression of eight Hox genes, including Ubx and abd-A, in early instar nymphs further indicates a post-embryonic role, possibly in gill development. These data provide a basis for H. limbata as a complementary Ephemeridae model to the growing repertoire of mayfly model species and molecular techniques.

Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components

Modern venomics is increasing its focus on hymenopterans such as honeybees, bumblebees, parasitoid wasps, ants and true wasps. However solitary bees remain understudied in comparison and the few available venom studies focus on short melittin-like sequences and antimicrobial peptides. Herein we describe the first comprehensive venom profile of a solitary bee, the violet carpenter bee Xylocopa violacea, by using proteo-transcriptomics. We reveal a diverse and complex venom profile with 43 different protein families identified from dissected venom gland extracts of which 32 are also detected in the defensively injected venom. Melittin and apamin are the most highly secreted components, followed by Phospholipase A2, Icarapin, Secapin and three novel components. Other components, including eight novel protein families, are rather lowly expressed. We further identify multiple forms of apamin-like peptides. The melittin-like sequences of solitary bees separate into two clades, one comprised most sequences from solitary bees including xylopin (the variant in Xylocopa), while sequences from Lasioglossa appear closer related to melittin-like peptides from Bombus (Bombolittins). Our study suggests that more proteo-transcriptomic data from other solitary bees should be complemented with corresponding genome data to fully understand the evolution and complexity of bee venom proteins, and is of a particular need to disentangle the ambiguous phylogenetic relations of short peptides.

•

Venoms of solitary bees are so far surprisingly little studied.

•

We analyse here the venom profile from the solitary carpenter bee X. violacea.

•

The venom of X. violacea is similar complex as in honey bees.

•

Besides new melittin variants eight novel venom components are identified.

•

Our results contribute to understand better the evoluiton of bee venoms.

No evidence of DUI in the Mediterranean alien species Brachidontes pharaonis (P. Fisher, 1870) despite mitochondrial heteroplasmy

Two genetically different mitochondrial haplogroups of Brachidontes pharaonis (p-distance 6.8%) have been identified in the Mediterranean Sea. This hinted at a possible presence of doubly uniparental inheritance in this species. To ascertain this possibility, we sequenced two complete mitogenomes of Brachidontes pharaonis mussels and performed a qPCR analysis to measure the relative mitogenome copy numbers of both mtDNAs. Despite the presence of two very similar regions composed entirely of repetitive sequences in the two haplogroups, no recombination between mitogenomes was detected. In heteroplasmic individuals, both mitogenomes were present in the generative tissues of both sexes, which argues against the presence of doubly uniparental inheritance in this species.

Modern venomics-Current insights, novel methods, and future perspectives in biological and applied animal venom research

Venoms have evolved >100 times in all major animal groups, and their components, known as toxins, have been fine-tuned over millions of years into highly effective biochemical weapons. There are many outstanding questions on the evolution of toxin arsenals, such as how venom genes originate, how venom contributes to the fitness of venomous species, and which modifications at the genomic, transcriptomic, and protein level drive their evolution. These questions have received particularly little attention outside of snakes, cone snails, spiders, and scorpions. Venom compounds have further become a source of inspiration for translational research using their diverse bioactivities for various applications. We highlight here recent advances and new strategies in modern venomics and discuss how recent technological innovations and multi-omic methods dramatically improve research on venomous animals. The study of genomes and their modifications through CRISPR and knockdown technologies will increase our understanding of how toxins evolve and which functions they have in the different ontogenetic stages during the development of venomous animals. Mass spectrometry imaging combined with spatial transcriptomics, in situ hybridization techniques, and modern computer tomography gives us further insights into the spatial distribution of toxins in the venom system and the function of the venom apparatus. All these evolutionary and biological insights contribute to more efficiently identify venom compounds, which can then be synthesized or produced in adapted expression systems to test their bioactivity. Finally, we critically discuss recent agrochemical, pharmaceutical, therapeutic, and diagnostic (so-called translational) aspects of venoms from which humans benefit.

A Survey on the Distribution of Ovothiol and ovoA Gene Expression in Different Tissues and Cells: A Comparative Analysis in Sea Urchins and Mussels

Ovothiols are histidine-derived thiols produced by a variety of marine invertebrates, protists and bacteria. These compounds, which are among the strongest natural antioxidants, are involved in controlling the cellular redox balance due to their redox exchange with glutathione. Although ovothiols were initially reported as protective agents against environmental stressors, new evidence suggests that they can also act as pheromones and participate in fundamental biological processes such as embryogenesis. To get further insight into the biological roles of ovothiols, we compared ovothiol biosynthesis in the sea urchin Paracentrotus lividus and in the mussel Mytilus galloprovincialis, the two species that represent the richest sources of these compounds among marine invertebrates. Ovothiol content was measured in different tissues and in the immune cells from both species and the expression levels of ovoA, the gene responsible for ovothiol biosynthesis, was inferred from publicly available transcriptomes. A comparative analysis of ovothiol biosynthesis in the two species allowed the identification of the tissues and cells synthesizing the metabolite and highlighted analogies and differences between sea urchins and mussels. By improving our knowledge on the biological roles of ovothiols and pointing out the existence of sustainable natural sources for their isolation, this study provides the basis for future biotechnological investigations on these valuable compounds.

TransPi - a comprehensive TRanscriptome ANalysiS PIpeline for de novo transcriptome assembly

The use of RNA-Seq data and the generation of de novo transcriptome assemblies have been pivotal for studies in ecology and evolution. This is distinctly true for non-model organisms, where no genome information is available. In such organisms, studies of differential gene expression, DNA enrichment baits design, and phylogenetics can all be accomplished with de novo transcriptome assemblies. Multiple tools are available for transcriptome assembly, however, no single tool can provide the best assembly for all datasets. Therefore, a multi assembler approach, followed by a reduction step, is often sought to generate an improved representation of the assembly. To reduce errors in these complex analyses while at the same time attaining reproducibility and scalability, automated workflows have been essential in the analysis of RNA-Seq data. However, most of these tools are designed for species where genome data is used as reference for the assembly process, limiting their use in non-model organisms. We present TransPi, a comprehensive pipeline for de novo transcriptome assembly, with minimum user input but without losing the ability of a thorough analysis. A combination of different model organisms, k-mer sets, read lengths, and read quantities were used for assessing the tool. Furthermore, a total of 49 non-model organisms, spanning different phyla, were also analysed. Compared to approaches using single assemblers only, TransPi produces higher BUSCO completeness percentages, and a concurrent significant reduction in duplication rates. TransPi is easy to configure and can be deployed seamlessly using Conda, Docker and Singularity.

A simple guide to de novo transcriptome assembly and annotation

A transcriptome constructed from short-read RNA sequencing (RNA-seq) is an easily attainable proxy catalog of protein-coding genes when genome assembly is unnecessary, expensive or difficult. In the absence of a sequenced genome to guide the reconstruction process, the transcriptome must be assembled de novo using only the information available in the RNA-seq reads. Subsequently, the sequences must be annotated in order to identify sequence-intrinsic and evolutionary features in them (for example, protein-coding regions). Although straightforward at first glance, de novo transcriptome assembly and annotation can quickly prove to be challenging undertakings. In addition to familiarizing themselves with the conceptual and technical intricacies of the tasks at hand and the numerous pre- and post-processing steps involved, those interested must also grapple with an overwhelmingly large choice of tools. The lack of standardized workflows, fast pace of development of new tools and techniques and paucity of authoritative literature have served to exacerbate the difficulty of the task even further. Here, we present a comprehensive overview of de novo transcriptome assembly and annotation. We discuss the procedures involved, including pre- and post-processing steps, and present a compendium of corresponding tools.

Salivary and Intestinal Transcriptomes Reveal Differential Gene Expression in Starving, Fed and Trypanosoma cruzi-Infected Rhodnius neglectus

Rhodnius neglectus is a potential vector of Trypanosoma cruzi (Tc), the causative agent of Chagas disease. The salivary glands (SGs) and intestine (INT) are actively required during blood feeding. The saliva from SGs is injected into the vertebrate host, modulating immune responses and favoring feeding for INT digestion. Tc infection significantly alters the physiology of these tissues; however, studies that assess this are still scarce. This study aimed to gain a better understanding of the global transcriptional expression of genes in SGs and INT during fasting (FA), fed (FE), and fed in the presence of Tc (FE + Tc) conditions. In FA, the expression of transcripts related to homeostasis maintenance proteins during periods of stress was predominant. Therefore, the transcript levels of Tret1-like and Hsp70Ba proteins were increased. Blood appeared to be responsible for alterations found in the FE group, as most of the expressed transcripts, such as proteases and cathepsin D, were related to digestion. In FE + Tc group, there was a decreased expression of blood processing genes for insect metabolism (e.g., Antigen-5 precursor, Pr13a, and Obp), detoxification (Sult1) in INT and acid phosphatases in SG. We also found decreased transcriptional expression of lipocalins and nitrophorins in SG and two new proteins, pacifastin and diptericin, in INT. Several transcripts of unknown proteins with investigative potential were found in both tissues. Our results also show that the presence of Tc can change the expression in both tissues for a long or short period of time. While SG homeostasis seems to be re-established on day 9, changes in INT are still evident. The findings of this study may be used for future research on parasite-vector interactions and contribute to the understanding of food physiology and post-meal/infection in triatomines.

Production of Agrocinopine A by Ipomoea batatas Agrocinopine Synthase in Transgenic Tobacco and Its Effect on the Rhizosphere Microbial Community

Agrobacterium tumefaciens is a bacterial pathogen that causes crown gall disease on a wide range of eudicot plants by genetic transformation. Besides T-DNA integrated by natural transformation of plant vegetative tissues by pathogenic Agrobacterium spp., previous reports have indicated that T-DNA sequences originating from an ancestral Agrobacterium sp. are present in the genomes of all cultivated sweet potato (Ipomoea batatas) varieties analyzed. Expression of an Agrobacterium-derived agrocinopine synthase (ACS) gene was detected in leaf and root tissues of sweet potato, suggesting that the plant can produce agrocinopine, a sugar-phosphodiester opine considered to be utilized by some strains of Agrobacterium spp. in crown gall. To validate the product synthesized by Ipomoea batatas ACS (IbACS), we introduced IbACS into tobacco under a constitutive promoter. High-voltage paper electrophoresis followed by alkaline silver nitrate staining detected the production of an agrocinopine-like substance in IbACS1-expressing tobacco, and further mass spectrometry and nuclear magnetic resonance analyses of the product confirmed that IbACS can produce agrocinopine A from natural plant substrates. The partially purified compound was biologically active in an agrocinopine A bioassay. A 16S ribosomal RNA amplicon sequencing and meta-transcriptome analysis revealed that the rhizosphere microbial community of tobacco was affected by the expression of IbACS. A new species of Leifsonia (actinobacteria) was isolated as an enriched bacterium in the rhizosphere of IbACS1-expressing tobacco. This Leifsonia sp. can catabolize agrocinopine A produced in tobacco, indicating that the production of agrocinopine A attracts rhizosphere bacteria that can utilize this sugar-phosphodiester. These results suggest a potential role of IbACS conserved among sweet potato cultivars in manipulating their microbial community.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Relaxed selection on male mitochondrial genes in DUI bivalves eases the need for mitonuclear coevolution

Mitonuclear coevolution is an important prerequisite for efficient energy production in eukaryotes. However, many bivalve taxa experience doubly uniparental inheritance (DUI) and have sex-specific mitochondrial (mt) genomes, providing a challenge for mitonuclear coevolution. We examined possible mechanisms to reconcile mitonuclear coevolution with DUI. No nuclear-encoded, sex-specific OXPHOS paralogs were found in the DUI clam Ruditapes philippinarum, refuting OXPHOS paralogy as a solution in this species. It is also unlikely that mt changes causing disruption of nuclear interactions are strongly selected against because sex-specific mt-residues or those under positive selection in M mt genes were not depleted for contacting nuclear-encoded residues. However, M genomes showed consistently higher d_N /d_S ratios compared to putatively ancestral F genomes in all mt OXPHOS genes and across all DUI species. Further analyses indicated that this was consistently due to relaxed, not positive selection on M vs. F mt OXPHOS genes. Similarly, selection was relaxed on the F genome of DUI species compared to species with strict maternal inheritance. Coupled with recent physiological and molecular evolution studies, we suggest that relaxed selection on M mt function limits the need to maintain mitonuclear interactions in M genomes compared to F genomes. We discuss our findings with regard to OXPHOS function and the origin of DUI.

De novo Assembly, Annotation, and Analysis of Transcriptome Data of the Ladakh Ground Skink Provide Genetic Information on High-Altitude Adaptation

The Himalayan Arc is recognized as a global biodiversity hotspot. Among its numerous cryptic and undiscovered organisms, this composite high-mountain ecosystem harbors many taxa with adaptations to life in high elevations. However, evolutionary patterns and genomic features have been relatively rarely studied in Himalayan vertebrates. Here, we provide the first well-annotated transcriptome of a Greater Himalayan reptile species, the Ladakh Ground skink Asymblepharus ladacensis (Squamata: Scincidae). Based on tissues from the brain, an embryonic disc, and pooled organ material, using pair-end Illumina NextSeq 500 RNAseq, we assembled ~77,000 transcripts, which were annotated using seven functional databases. We tested ~1600 genes, known to be under positive selection in anurans and reptiles adapted to high elevations, and potentially detected positive selection for 114 of these genes in Asymblepharus. Even though the strength of these results is limited due to the single-animal approach, our transcriptome resource may be valuable data for further studies on squamate reptile evolution in the Himalayas as a hotspot of biodiversity.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which multiple species share the same conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which vivid colour and pattern are produced in a Müllerian mimicry complex of poison frogs. We did this by first assembling a high-quality de novo genome assembly for the mimic poison frog Ranitomeya imitator. This assembled genome is 6.8 Gbp in size, with a contig N50 of 300 Kbp R. imitator and two colour morphs from both Ranitomeya fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes that are differentially expressed throughout development, many of them related to melanocyte development, melanin synthesis, iridophore development and guanine synthesis. Polytypic differences within species may be the result of differences in expression and/or timing of expression, whereas convergence for colour pattern between species does not appear to be due to the same changes in gene expression. In addition, we identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species. Finally, we hypothesize that genes in the keratin family are important for producing different structural colours within these frogs.

Novel doublesex duplication associated with sexually dimorphic development of dogface butterfly wings

Sexually dimorphic development is responsible for some of the most remarkable phenotypic variation found in nature. Alternative splicing of the transcription factor gene doublesex (dsx) is a highly conserved developmental switch controlling the expression of sex-specific pathways. Here, we leverage sex-specific differences in butterfly wing color pattern to characterize the genetic basis of sexually dimorphic development. We use RNA-seq, immunolocalization, and motif binding site analysis to test specific predictions about the role of dsx in the development of structurally based ultraviolet (UV) wing patterns in Zerene cesonia (Southern Dogface). Unexpectedly, we discover a novel duplication of dsx that shows a sex-specific burst of expression associated with the sexually dimorphic UV coloration. The derived copy consists of a single exon that encodes a DNA binding but no protein-binding domain and has experienced rapid amino-acid divergence. We propose the novel dsx paralog may suppress UV scale differentiation in females, which is supported by an excess of Dsx-binding sites at cytoskeletal and chitin-related genes with sex-biased expression. These findings illustrate the molecular flexibility of the dsx gene in mediating the differentiation of secondary sexual characteristics.

Transcriptomic Analysis of Skin Color in Anole Lizards

Color and color pattern are critical for animal camouflage, reproduction, and defense. Few studies, however, have attempted to identify candidate genes for color and color pattern in squamate reptiles, a colorful group with over 10,000 species. We used comparative transcriptomic analyses between white, orange, and yellow skin in a color-polymorphic species of anole lizard to 1) identify candidate color and color-pattern genes in squamates and 2) assess if squamates share an underlying genetic basis for color and color pattern variation with other vertebrates. Squamates have three types of chromatophores that determine color pattern: guanine-filled iridophores, carotenoid- or pteridine-filled xanthophores/erythrophores, and melanin-filled melanophores. We identified 13 best candidate squamate color and color-pattern genes shared with other vertebrates: six genes linked to pigment synthesis pathways, and seven genes linked to chromatophore development and maintenance. In comparisons of expression profiles between pigment-rich and white skin, pigment-rich skin upregulated the pteridine pathway as well as xanthophore/erythrophore development and maintenance genes; in comparisons between orange and yellow skin, orange skin upregulated the pteridine and carotenoid pathways as well as melanophore maintenance genes. Our results corroborate the predictions that squamates can produce similar colors using distinct color-reflecting molecules, and that both color and color-pattern genes are likely conserved across vertebrates. Furthermore, this study provides a concise list of candidate genes for future functional verification, representing a first step in determining the genetic basis of color and color pattern in anoles.

Pincho: A Modular Approach to High Quality De Novo Transcriptomics

Transcriptomic reconstructions without reference (i.e., de novo) are common for data samples derived from non-model biological systems. These assemblies involve massive parallel short read sequence reconstructions from experiments, but they usually employ ad-hoc bioinformatic workflows that exhibit limited standardization and customization. The increasing number of transcriptome assembly software continues to provide little room for standardization which is exacerbated by the lack of studies on modularity that compare the effects of assembler synergy. We developed a customizable management workflow for de novo transcriptomics that includes modular units for short read cleaning, assembly, validation, annotation, and expression analysis by connecting twenty-five individual bioinformatic tools. With our software tool, we were able to compare the assessment scores based on 129 distinct single-, bi- and tri-assembler combinations with diverse k-mer size selections. Our results demonstrate a drastic increase in the quality of transcriptome assemblies with bi- and tri- assembler combinations. We aim for our software to improve de novo transcriptome reconstructions for the ever-growing landscape of RNA-seq data derived from non-model systems. We offer guidance to ensure the most complete transcriptomic reconstructions via the inclusion of modular multi-assembly software controlled from a single master console.

Mediterranean Pine Vole, Microtus duodecimcostatus: A Paradigm of an Opportunistic Breeder

Simple Summary

In temperate zones of the Earth, some mammalian species reproduce seasonally whereas others do it continuously. Other species are summer breeders in the north and winter breeders in the south. Thus, the reproductive pattern seems not to be a species-specific but a population-specific trait. We investigated the reproduction pattern of the Mediterranean pine vole, Microtus duodecimcostatus, in the area around the city of Granada in Southern Spain, and found that individuals living in wastelands reproduce seasonally whereas those living in close poplar plantations (just 8 km apart) reproduce throughout the year, as did voles captured in wastelands and kept in captivity. These animals represent thus a paradigm of an opportunistic breeder as particular individuals stop breeding or not, depending on the environmental conditions they face at any moment. Sexually inactive male voles undergo complete testis inactivation and their sperm production is halted. The immune system in active testes is depressed, a phenomenon known as “immune privilege” that protect germ cells from autoimmune attack. We studied gene activity in active and inactive testes and our results indicate that such an immune privilege is lost in inactive testes, suggesting an important role for this process during testis regression.

Abstract

Most mammalian species of the temperate zones of the Earth reproduce seasonally, existing a non-breeding period in which the gonads of both sexes undergo functional regression. It is widely accepted that photoperiod is the principal environmental cue controlling these seasonal changes, although several exceptions have been described in other mammalian species in which breeding depends on cues such as food or water availability. We studied the circannual reproductive cycle in males of the Mediterranean pine vole, Microtus duodecimcostatus, in the Southeastern Iberian Peninsula. Morphological, hormonal, functional, molecular and transcriptomic analyses were performed. As reported for populations of other species from the same geographic area, male voles captured in wastelands underwent seasonal testis regression in summer whereas, surprisingly, those living either in close poplar plantations or in our animal house reproduced throughout the year, showing that it is the microenvironment of a particular vole subpopulation what determines its reproductive status and that these animals are pure opportunistic, photoperiod-independent breeders. In addition, we show that several molecular pathways, including MAPK, are deregulated and that the testicular “immune privilege” is lost in the inactive testes, providing novel mechanisms linking seasonal testosterone reduction and testis regression.

The integrin-mediated adhesive complex in the ancestor of animals, fungi, and amoebae

Integrins are transmembrane receptors that activate signal transduction pathways upon extracellular matrix binding. The integrin-mediated adhesive complex (IMAC) mediates various cell physiological processes. Although the IMAC was thought to be specific to animals, in the past ten years these complexes were discovered in other lineages of Obazoa, the group containing animals, fungi, and several microbial eukaryotes. Very recently, many genomes and transcriptomes from Amoebozoa (the eukaryotic supergroup sister to Obazoa), other obazoans, orphan protist lineages, and the eukaryotes' closest prokaryotic relatives, have become available. To increase the resolution of where and when IMAC proteins exist and have emerged, we surveyed these newly available genomes and transcriptomes for the presence of IMAC proteins. Our results highlight that many of these proteins appear to have evolved earlier in eukaryote evolution than previously thought and that co-option of this apparently ancient protein complex was key to the emergence of animal-type multicellularity. The role of the IMACs in amoebozoans is unknown, but they play critical adhesive roles in at least some unicellular organisms.

Combined Proteotranscriptomic-Based Strategy to Discover Novel Antimicrobial Peptides from Cone Snails

Despite their impressive diversity and already broad therapeutic applications, cone snail venoms have received less attention as a natural source in the investigation of antimicrobial peptides than other venomous animals such as scorpions, spiders, or snakes. Cone snails are among the largest genera (Conus sp.) of marine invertebrates, with more than seven hundred species described to date. These predatory mollusks use their sophisticated venom apparatus to capture prey or defend themselves. In-depth studies of these venoms have unraveled many biologically active peptides with pharmacological properties of interest in the field of pain management, the treatment of epilepsy, neurodegenerative diseases, and cardiac ischemia. Considering sequencing efficiency and affordability, cone snail venom gland transcriptome analyses could allow the discovery of new, promising antimicrobial peptides. We first present here the need for novel compounds like antimicrobial peptides as a viable alternative to conventional antibiotics. Secondly, we review the current knowledge on cone snails as a source of antimicrobial peptides. Then, we present the current state of the art in analytical methods applied to crude or milked venom followed by how antibacterial activity assay can be implemented for fostering cone snail antimicrobial peptides studies. We also propose a new innovative profile Hidden Markov model-based approach to annotate full venom gland transcriptomes and speed up the discovery of potentially active peptides from cone snails.

Signal, bias, and the role of transcriptome assembly quality in phylogenomic inference

BACKGROUND

Phylogenomic approaches have great power to reconstruct evolutionary histories, however they rely on multi-step processes in which each stage has the potential to affect the accuracy of the final result. Many studies have empirically tested and established methodology for resolving robust phylogenies, including selecting appropriate evolutionary models, identifying orthologs, or isolating partitions with strong phylogenetic signal. However, few have investigated errors that may be initiated at earlier stages of the analysis. Biases introduced during the generation of the phylogenomic dataset itself could produce downstream effects on analyses of evolutionary history. Transcriptomes are widely used in phylogenomics studies, though there is little understanding of how a poor-quality assembly of these datasets could impact the accuracy of phylogenomic hypotheses. Here we examined how transcriptome assembly quality affects phylogenomic inferences by creating independent datasets from the same input data representing high-quality and low-quality transcriptome assembly outcomes.

RESULTS

By studying the performance of phylogenomic datasets derived from alternative high- and low-quality assembly inputs in a controlled experiment, we show that high-quality transcriptomes produce richer phylogenomic datasets with a greater number of unique partitions than low-quality assemblies. High-quality assemblies also give rise to partitions that have lower alignment ambiguity and less compositional bias. In addition, high-quality partitions hold stronger phylogenetic signal than their low-quality transcriptome assembly counterparts in both concatenation- and coalescent-based analyses.

CONCLUSIONS

Our findings demonstrate the importance of transcriptome assembly quality in phylogenomic analyses and suggest that a portion of the uncertainty observed in such studies could be alleviated at the assembly stage.

Gene family amplification facilitates adaptation in freshwater unionid bivalve Megalonaias nervosa

Freshwater unionid bivalves currently face severe anthropogenic challenges. Over 70% of species in the United States are threatened, endangered or extinct due to pollution, damming of waterways and overfishing. These species are notable for their unusual life history strategy, parasite-host co-evolution and biparental mitochondrial inheritance. Among this clade, the washboard mussel Megalonaias nervosa is one species that remains prevalent across the Southeastern United States, with robust population sizes. We have created a reference genome for M. nervosa to determine how genome content has evolved in the face of these widespread environmental challenges. We observe dynamic changes in genome content, with a burst of recent transposable element proliferation causing a 382 Mb expansion in genome content. Birth-death models suggest rapid expansions among gene families, with a mutation rate of 1.16 × 10^-8 duplications per gene per generation. Cytochrome P450 gene families have experienced exceptional recent amplification beyond expectations based on genome-wide birth-death processes. These genes are associated with increased rates of amino acid changes, a signature of selection driving evolution of detox genes. Fitting evolutionary models of adaptation from standing genetic variation, we can compare adaptive potential across species and mutation types. The large population size in M. nervosa suggests a 4.7-fold advantage in the ability to adapt from standing genetic variation compared with a low diversity endemic E. hopetonensis. Estimates suggest that gene family evolution may offer an exceptional substrate of genetic variation in M. nervosa, with P_sgv  = 0.185 compared with P_sgv  = 0.067 for single nucleotide changes. Hence, we suggest that gene family evolution is a source of 'hopeful monsters' within the genome that may facilitate adaptation when selective pressures shift. These results suggest that gene family expansion is a key driver of adaptive evolution in this key species of freshwater Unionidae that is currently facing widespread environmental challenges. This work has clear implications for conservation genomics on freshwater bivalves as well as evolutionary theory. This genome represents a first step to facilitate reverse ecological genomics in Unionidae and identify the genetic underpinnings of phenotypic diversity.

Comparative transcriptomics of spotted seatrout (Cynoscion nebulosus) populations to cold and heat stress

Resilience to climate change depends on a species' adaptive potential and phenotypic plasticity. The latter can enhance survival of individual organisms during short periods of extreme environmental perturbations, allowing genetic adaptation to take place over generations. Along the U.S. East Coast, estuarine-dependent spotted seatrout (Cynoscion nebulosus) populations span a steep temperature gradient that provides an ideal opportunity to explore the molecular basis of phenotypic plasticity. Genetically distinct spotted seatrout sampled from a northern and a southern population were exposed to acute cold and heat stress (5 biological replicates in each treatment and control group), and their transcriptomic responses were compared using RNA-sequencing (RNA-seq). The southern population showed a larger transcriptomic response to acute cold stress, whereas the northern population showed a larger transcriptomic response to acute heat stress compared with their respective population controls. Shared transcripts showing significant differences in expression levels were predominantly enriched in pathways that included metabolism, transcriptional regulation, and immune response. In response to heat stress, only the northern population significantly upregulated genes in the apoptosis pathway, which could suggest greater vulnerability to future heat waves in this population as compared to the southern population. Genes showing population-specific patterns of expression, including hpt, acot, hspa5, and hsc71, are candidates for future studies aiming to monitor intraspecific differences in temperature stress responses in spotted seatrout. Our findings contribute to the current understanding of phenotypic plasticity and provide a basis for predicting the response of a eurythermal fish species to future extreme temperatures.

Giant lungfish genome elucidates the conquest of land by vertebrates

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Isolating the Role of Corticosterone in the Hypothalamic-Pituitary-Gonadal Transcriptomic Stress Response

Investigation of the negative impacts of stress on reproduction has largely centered around the effects of the adrenal steroid hormone, corticosterone (CORT), and its influence on a system of tissues vital for reproduction-the hypothalamus of the brain, the pituitary gland, and the gonads (the HPG axis). Research on the action of CORT on the HPG axis has predominated the stress and reproductive biology literature, potentially overshadowing other influential mediators. To gain a more complete understanding of how elevated CORT affects transcriptomic activity of the HPG axis, we experimentally examined its role in male and female rock doves (Columba livia). We exogenously administrated CORT to mimic circulating levels during the stress response, specifically 30 min of restraint stress, an experimental paradigm known to increase circulating CORT in vertebrates. We examined all changes in transcription within each level of the HPG axis as compared to both restraint-stressed birds and vehicle-injected controls. We also investigated the differential transcriptomic response to CORT and restraint-stress in each sex. We report causal and sex-specific effects of CORT on the HPG transcriptomic stress response. Restraint stress caused 1567 genes to uniquely differentially express while elevated circulating CORT was responsible for the differential expression of 304 genes. Only 108 genes in females and 8 in males differentially expressed in subjects that underwent restraint stress and those who were given exogenous CORT. In response to elevated CORT and restraint-stress, both sexes shared the differential expression of 5 genes, KCNJ5, CISH, PTGER3, CEBPD, and ZBTB16, all located in the pituitary. The known functions of these genes suggest potential influence of elevated CORT on immune function and prolactin synthesis. Gene expression unique to each sex indicated that elevated CORT affected more gene transcription in females than males (78 genes versus 3 genes, respectively). To our knowledge, this is the first study to isolate the role of CORT in HPG genomic transcription during a stress response. We present an extensive and openly accessible view of the role corticosterone in the HPG transcriptomic stress response. Because the HPG system is well conserved across vertebrates, these data have the potential to inspire new therapeutic strategies for reproductive dysregulation in multiple vertebrate systems, including our own.

Novel Allergen Discovery through Comprehensive De Novo Transcriptomic Analyses of Five Shrimp Species

Shellfish allergy affects 2% of the world's population and persists for life in most patients. The diagnosis of shellfish allergy, in particular shrimp, is challenging due to the similarity of allergenic proteins from other invertebrates. Despite the clinical importance of immunological cross-reactivity among shellfish species and between allergenic invertebrates such as dust mites, the underlying molecular basis is not well understood. Here we mine the complete transcriptome of five frequently consumed shrimp species to identify and compare allergens with all known allergen sources. The transcriptomes were assembled de novo, using Trinity, from raw RNA-Seq data of the whiteleg shrimp (Litopenaeus vannamei), black tiger shrimp (Penaeus monodon), banana shrimp (Fenneropenaeus merguiensis), king shrimp (Melicertus latisulcatus), and endeavour shrimp (Metapenaeus endeavouri). BLAST searching using the two major allergen databases, WHO/IUIS Allergen Nomenclature and AllergenOnline, successfully identified all seven known crustacean allergens. The analyses revealed up to 39 unreported allergens in the different shrimp species, including heat shock protein (HSP), alpha-tubulin, chymotrypsin, cyclophilin, beta-enolase, aldolase A, and glyceraldehyde-3-phosphate dehydrogenase (G3PD). Multiple sequence alignment (Clustal Omega) demonstrated high homology with allergens from other invertebrates including mites and cockroaches. This first transcriptomic analyses of allergens in a major food source provides a valuable resource for investigating shellfish allergens, comparing invertebrate allergens and future development of improved diagnostics for food allergy.

RNAflow: An Effective and Simple RNA-Seq Differential Gene Expression Pipeline Using Nextflow

RNA-Seq enables the identification and quantification of RNA molecules, often with the aim of detecting differentially expressed genes (DEGs). Although RNA-Seq evolved into a standard technique, there is no universal gold standard for these data's computational analysis. On top of that, previous studies proved the irreproducibility of RNA-Seq studies. Here, we present a portable, scalable, and parallelizable Nextflow RNA-Seq pipeline to detect DEGs, which assures a high level of reproducibility. The pipeline automatically takes care of common pitfalls, such as ribosomal RNA removal and low abundance gene filtering. Apart from various visualizations for the DEG results, we incorporated downstream pathway analysis for common species as Homo sapiens and Mus musculus. We evaluated the DEG detection functionality while using qRT-PCR data serving as a reference and observed a very high correlation of the logarithmized gene expression fold changes.

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

Semimytilus algosus: first known hermaphroditic mussel with doubly uniparental inheritance of mitochondrial DNA

Doubly uniparental inheritance (DUI) of mitochondrial DNA is a rare phenomenon occurring in some freshwater and marine bivalves and is usually characterized by the mitochondrial heteroplasmy of male individuals. Previous research on freshwater Unionida mussels showed that hermaphroditic species do not have DUI even if their closest gonochoristic counterparts do. No records showing DUI in a hermaphrodite have ever been reported. Here we show for the first time that the hermaphroditic mussel Semimytilus algosus (Mytilida), very likely has DUI, based on the complete sequences of both mitochondrial DNAs and the distribution of mtDNA types between male and female gonads. The two mitogenomes show considerable divergence (34.7%). The presumably paternal M type mitogenome dominated the male gonads of most studied mussels, while remaining at very low or undetectable levels in the female gonads of the same individuals. If indeed DUI can function in the context of simultaneous hermaphroditism, a change of paradigm regarding its involvement in sex determination is needed. It is apparently associated with gonadal differentiation rather than with sex determination in bivalves.