Next-generation sequencing of transcriptomes: a guide to RNA isolation in nonmodel animals

High-throughput molecular gut content analysis of aphids identifies plants relevant for potato virus Y epidemiology

Aphids are phloem-feeding insects that reduce crop productivity due to feeding and transmission of plant viruses. When aphids disperse across the landscape to colonize new host plants, they will often probe on a wide variety of nonhost plants before settling on a host suitable for feeding and reproduction. There is limited understanding of the diversity of plants that aphids probe on within a landscape, and characterizing this diversity can help us better understand host use patterns of aphids. Here, we used gut content analysis (GCA) to identify plant genera that were probed by aphid vectors of potato virus Y (PVY). Aphids were trapped weekly near potato fields during the growing seasons of 2020 and 2021 in San Luis Valley in Colorado. High-throughput sequencing of plant barcoding genes, trnF and ITS2, from 200 individual alate (i.e., winged) aphids representing nine vector species of PVY was performed using the PacBio sequencing platform, and sequences were identified to genus using NCBI BLASTn. We found that 34.7% of aphids probed upon presumed PVY host plants and that two of the most frequently detected plant genera, Solanum and Brassica, represent important crops and weeds within the study region. We found that 75% of aphids frequently probed upon PVY nonhosts including many species that are outside of their reported host ranges. Additionally, 19% of aphids probed upon more than one plant species. This study provides the first evidence from high-throughput molecular GCA of aphids and reveals host use patterns that are relevant for PVY epidemiology.

Comparison of human dental tissue RNA extraction methods for RNA sequencing

OBJECTIVE

The purpose of this study was to identify an efficient RNA extraction method for periodontal ligament (PDL) and dental pulp (DP) tissues to be used in RNA sequencing studies, given the increased use of these techniques in dental research and the lack of standard protocols.

DESIGN

PDL and DP were harvested from extracted third molars. Total RNA was extracted with four RNA extraction kits. RNA concentration, purity and integrity were assessed by means of NanoDrop and Bioanalyzer and statistically compared.

RESULTS

RNA from PDL was more likely to be degraded than that of DP. The TRIzol method yielded the highest RNA concentration from both tissues. All methods harvested RNA with A260/A280 close to 2.0 and with A260/A230 above 1.5, except for the A260/A230 from PDL obtained with the RNeasy Mini kit. For RNA integrity, the RNeasy Fibrous Tissue Mini kit yielded the highest RIN values and 28 S/18 S from PDL, while the RNeasy Mini kit obtained relatively high RIN values with an appropriate 28 S/18 S for DP.

CONCLUSION

Significantly different results were obtained for PDL and DP when using the RNeasy Mini kit. The RNeasy Mini kit provided the highest RNA yields and quality for DP, while the RNeasy Fibrous Tissue Mini kit obtained the highest quality RNA from PDL.

Influence of RNA-Seq library construction, sampling methods, and tissue harvesting time on gene expression estimation

RNA sequencing (RNA-Seq) is popular for measuring gene expression in non-model organisms, including wild populations. While RNA-Seq can detect gene expression variation among wild-caught individuals and yield important insights into biological function, sampling methods can also affect gene expression estimates. We examined the influence of multiple technical variables on estimated gene expression in a non-model fish, the westslope cutthroat trout (Oncorhynchus clarkii lewisi), using two RNA-Seq library types: 3' RNA-Seq (QuantSeq) and whole mRNA-Seq (NEB). We evaluated effects of dip netting versus electrofishing, and of harvesting tissue immediately versus 5 min after euthanasia on estimated gene expression in blood, gill, and muscle. We found no significant differences in gene expression between sampling methods or tissue collection times with either library type. When library types were compared using the same blood samples, 58% of genes detected by both NEB and QuantSeq showed significantly different expression between library types, and NEB detected 31% more genes than QuantSeq. Although the two library types recovered different numbers of genes and expression levels, results with NEB and QuantSeq were consistent in that neither library type showed differences in gene expression between sampling methods and tissue harvesting times. Our study suggests that researchers can safely rely on different fish sampling strategies in the field. In addition, while QuantSeq is more cost effective, NEB detects more expressed genes. Therefore, when it is crucial to detect as many genes as possible (especially low expressed genes), when alternative splicing is of interest, or when working with an organism lacking good genomic resources, whole mRNA-Seq is more powerful.

Purification methods to reduce interference by dextran sodium sulfate with quantification of gene expression in intestinal tissue samples from a piglet model of colitis

Dextran sodium sulfate (DSS) is commonly used to induce intestinal (i.e., colonic) inflammation in a variety of animal models. However, DSS is known to cause interference when using quantitative-real time polymerase chain reaction (qRT-PCR) methods, thereby invalidating accurate and precise measurement of tissue gene expression. Therefore, the goal of this study was to determine whether different mRNA purification methods would reduce DSS-interference. Colonic tissue samples were collected at postnatal days (PND) 27 or 28 from pigs that had not been administered DSS (Control), and two independent groups of pigs that received 1.25 g of DSS/kg of BW/d (DSS-1 and DSS-2) from PND 14 to 18. Tissue samples collected were subsequently stratified into three purification methods (i.e., 9 total treatment × method combinations), including: 1) no purification, 2) purification with lithium chloride (LiCl), or 3) purification using spin column filtration. All data were analyzed using a one-way ANOVA in the Mixed procedure of SAS. The average RNA concentrations across all treatments were between 1,300 and 1,800 μg/μL for all three in vivo groups. Although there were statistical differences among purification methods, the 260/280 and 260/230 ratios fell between acceptable limits of 2.0 to 2.1 and 2.0 to 2.2, respectively, for all treatment groups. This confirms the RNA quality was adequate and not influenced by purification method in addition to suggesting the absence of phenol, salts, and carbohydrate contamination. For pigs in the Control group that did not receive DSS, qRT-PCR Ct values of four cytokines were achieved, though these values were not altered by purification method. For pigs that had undergone DSS dosing, those tissues subjected to either no purification or purification using LiCl did not generate applicable Ct values. However, when tissues derive from DSS-treated pigs underwent spin column purification, half of the samples from DSS-1 and DSS-2 groups generated appropriate Ct estimates. Therefore, spin column purification appeared to be more effective than LiCl purification, but no method was 100% effective, so caution should be exercised when interpreting gene expression results from studies where animals are exposed to DSS-induced colitis.

Lineage-specific rearrangement of chromatin loops and epigenomic features during adipocytes and osteoblasts commitment

Human mesenchymal stem cells (hMSCs) can be differentiated into adipocytes and osteoblasts. The processes are driven by the rewiring of chromatin architectures and transcriptomic/epigenomic changes. Here, we induced hMSCs to adipogenic and osteogenic differentiation, and performed 2 kb resolution Hi-C experiments for chromatin loops detection. We also generated matched RNA-seq, ChIP-seq and ATAC-seq data for integrative analysis. After comprehensively comparing adipogenesis and osteogenesis, we quantitatively identified lineage-specific loops and screened out lineage-specific enhancers and open chromatin. We reveal that lineage-specific loops can activate gene expression and facilitate cell commitment through combining enhancers and accessible chromatin in a lineage-specific manner. We finally proposed loop-mediated regulatory networks and identified the controlling factors for adipocytes and osteoblasts determination. Functional experiments validated the lineage-specific regulation networks towards IRS2 and RUNX2 that are associated with adipogenesis and osteogenesis, respectively. These results are expected to help better understand the chromatin conformation determinants of hMSCs fate commitment.

Single-chromosome fission yeast models reveal the configuration robustness of a functional genome

In eukaryotic organisms, genetic information is usually carried on multiple chromosomes. Whether and how the number and configuration of chromosomes affect organismal fitness and speciation remain unclear. Here, we have successfully established several single-chromosome fission yeast Schizosaccharomyces pombe strains, in which the three natural chromosomes have been fused into one giant chromosome in different orders. Chromosome fusions accompanied by the deletions of telomeres and centromeres result in the loss of chromosomal interactions and a drastic change of global chromosome organization, but alter gene expression marginally. The single-chromosome strains display little defects in cell morphology, mitosis, genotoxin sensitivity, and meiosis. Crosses between a wild-type strain and a single-chromosome strain or between two single-chromosome strains with different fusion orders suffer defective meiosis and poor spore viability. We conclude that eukaryotic genomes are robust against dramatic chromosomal reconfiguration, and stochastic changes in chromosome number and genome organization during evolution underlie reproductive isolation and speciation.

De Novo Transcriptome Assembly and Analysis of the Flat Oyster Pathogenic Protozoa Bonamia Ostreae

The flat oyster Ostrea edulis is an oyster species native to Europe. It has declined to functional extinction in many areas of the NE Atlantic for several decades. Factors explaining this decline include over-exploitation of natural populations and diseases like bonamiosis, regulated across both the EU and the wider world and caused by the intracellular protozoan parasite Bonamia ostreae. To date, very limited sequence data are available for this Haplosporidian species. We present here the first transcriptome of B. ostreae. As this protozoan is not yet culturable, it remains extremely challenging to obtain high-quality -omic data. Thanks to a specific parasite isolation protocol and a dedicated bioinformatic pipeline, we were able to obtain a high-quality transcriptome for an intracellular marine micro-eukaryote, which will be very helpful to better understand its biology and to consider the development of new relevant diagnostic tools.

Assessment of RNA extraction protocols from cladocerans

The usage of cladocerans as non-model organisms in ecotoxicological and risk assessment studies has intensified in recent years due to their ecological importance in aquatic ecosystems. The molecular assessment such as gene expression analysis has been introduced in ecotoxicological and risk assessment to link the expression of specific genes to a biological process in the cladocerans. The validity and accuracy of gene expression analysis depends on the quantity, quality and integrity of extracted ribonucleic acid (RNA) of the sample. However, the standard methods of RNA extraction from the cladocerans are still lacking. This study evaluates the extraction of RNA from tropical freshwater cladocerans Moina micrura using two methods: the phenol-chloroform extraction method (QIAzol) and a column-based kit (Qiagen Micro Kit). Glycogen was introduced in both approaches to enhance the recovery of extracted RNA and the extracted RNA was characterised using spectrophotometric analysis (NanoDrop), capillary electrophoresis (Bioanalyzer). Then, the extracted RNA was analysed with reverse transcription polymerase chain reaction (RT-PCR) to validate the RNA extraction method towards downstream gene expression analysis. The results indicate that the column-based kit is most suitable for the extraction of RNA from M. micrura, with the quantity (RNA concentration = 26.90 ± 6.89 ng/μl), quality (A260:230 = 1.95 ± 0.15, A280:230 = 1.85 ± 0.09) and integrity (RNA integrity number, RIN = 7.20 ± 0.16). The RT-PCR analysis shows that the method successfully amplified both alpha tubulin and actin gene at 33-35 cycles (i.e. Ct = 32.64 to 33.48). The results demonstrate that the addition of glycogen is only suitable for the phenol-chloroform extraction method. RNA extraction with high and comprehensive quality control assessment will increase the accuracy and reliability of downstream gene expression, thus providing more ecotoxicological data at the molecular biological level on other freshwater zooplankton species.

Vitellogenin gene expression in marine mussels exposed to ethinylestradiol: No induction at the transcriptional level

Vitellogenin (Vtg), a large multidomain protein precursor of egg-yolk proteins, is used as an endocrine disruption biomarker in fish, and in the last decades, its use has been extended to invertebrates like mollusks. However, it remains unclear whether invertebrate endocrine system produces Vtg in response to estrogens, like it occurs in oviparous vertebrates. In a previous study, no evidence of induction of Vtg expression at protein level was found in gonads of the marine mussel Mytilus galloprovincialis after exposure to the estrogenic chemical 17α-ethinylestradiol (EE2). In the present follow-up study, it was investigated whether there is any effect of EE2 on Vtg abundance at transcriptional level in M. galloprovincialis gonads. To this aim, RT-qPCR analysis targeting three different domains of Vtg transcript was performed on gonads of mussels that were exposed either 4 or 24 days to 100 ng/L EE2. In addition, several reference genes were analysed and a selection of these for potential use in further RT-qPCR analyses on mussel male and female gonads is provided. Results showed higher expression in females than in males for the three analysed Vtg domains, and no evidence of Vtg mRNA induction due to EE2 either in females or males. The present results, together with those obtained from previous analysis at protein level, support that Vtg is not an adequate biomarker for xenoestrogenicity in marine mussels. Additionally, nucleotide sequences of Vtg transcripts of three closely-related species from Mytilus edulis complex (M. galloprovincialis, M. edulis and M. trossulus) are provided and compared with Vtg sequences from other mollusk species to assess the level of conservation and evolutionary relationships among species.

Bias in RNA-seq Library Preparation: Current Challenges and Solutions

Although RNA sequencing (RNA-seq) has become the most advanced technology for transcriptome analysis, it also confronts various challenges. As we all know, the workflow of RNA-seq is extremely complicated and it is easy to produce bias. This may damage the quality of RNA-seq dataset and lead to an incorrect interpretation for sequencing result. Thus, our detailed understanding of the source and nature of these biases is essential for the interpretation of RNA-seq data, finding methods to improve the quality of RNA-seq experimental, or development bioinformatics tools to compensate for these biases. Here, we discuss the sources of experimental bias in RNA-seq. And for each type of bias, we discussed the method for improvement, in order to provide some useful suggestions for researcher in RNA-seq experimental.

A toolkit for studying Varroa genomics and transcriptomics: preservation, extraction, and sequencing library preparation

BACKGROUND

The honey bee parasite, Varroa destructor, is a leading cause of honey bee population declines. In addition to being an obligate ectoparasitic mite, Varroa carries several viruses that infect honey bees and act as the proximal cause of colony collapses. Nevertheless, until recently, studies of Varroa have been limited by the paucity of genomic tools. Lab- and field-based methods exploiting such methods are still nascent. This study developed a set of methods for preserving Varroa DNA and RNA from the field to the lab and processing them into sequencing libraries. We performed preservation experiments in which Varroa mites were immersed in TRIzol, RNAlater, and absolute ethanol for preservation periods up to 21 days post-treatment to assess DNA and RNA integrity.

RESULTS

For both DNA and RNA, mites preserved in TRIzol and RNAlater at room temperature degraded within 10 days post-treatment. Mites preserved in ethanol at room temperature and 4 °C remained intact through 21 days. Varroa mite DNA and RNA libraries were created and sequenced for ethanol preserved samples, 15 and 21 days post-treatment. All DNA sequences mapped to the V. destructor genome at above 95% on average, while RNA sequences mapped to V. destructor, but also sometimes to high levels of the deformed-wing virus and to various organisms.

CONCLUSIONS

Ethanolic preservation of field-collected mites is inexpensive and simple, and allows them to be shipped and processed successfully in the lab for a wide variety of sequencing applications. It appears to preserve RNA from both Varroa and at least some of the viruses it vectors.

Testing the Stability of Plasma Protein and Whole Blood RNA in Archived Blood of Loggerhead Sea Turtles, Caretta caretta

Sample storage conditions can affect accuracy and reproducibility of biological measurements. Storing samples rapidly at the lowest available temperatures is considered ideal but is not always feasible when sampling in remote and logistically challenging field conditions, as is often the case with sea turtles. The objective of this study was to examine the stability of plasma proteins and quality of whole blood RNA from loggerhead sea turtle samples collected as part of an eighteen-year-long curated specimen collection. These biological variables are often used to assess sea turtle health; therefore, it is necessary to maintain the integrity of these components during storage. Protein electrophoresis was conducted on heparinized plasma from individual turtles collected in 2018 (n = 3), 2008 (n = 3), and 2001 (n = 3). Plasma was also pooled from four turtles sampled in 2018 and subjected to various storage temperatures. Whole blood was collected in blood collection tubes containing sodium heparin or PAXgene tubes with an RNA preservative. These were subjected to different storage treatments that can possibly occur during logistically difficult field sampling. Following various treatments, plasma proteins showed minor differences across collection years and no differences among storage treatments were observed, even when exposed to 38°C for three hours. RNA quality was assessed from whole blood using an RNA integrity number (RIN). RINs were poor from sodium heparin tubes that were frozen and from PAXgene tubes after an extended thaw. High-quality RNA was obtained from sodium heparin tubes that were never frozen and from PAXgene tubes with freezing delayed by up to 11 days. Overall, these results indicate that plasma proteins remain stable over time and when exposed to undesirable storage conditions, and RNA degrades rapidly in sea turtle blood after freezing and when not properly preserved. These aspects are important to consider when planning sampling protocols and logistics for optimal long-term sample preservation.

Identification of Genic SSRs Provide a Perspective for Studying Environmental Adaptation in the Endemic Shrub Tetraena mongolica

Tetraena mongolica is a xerophytic shrub endemic to desert regions in Inner Mongolia. This species has evolved distinct survival strategies that allow it to adapt to hyper-drought and heterogeneous habitats. Simple sequence repeats (SSRs) may provide a molecular basis in plants for fast adaptation to environmental change. Thus, identifying SSRs and their possible effects on gene behavior has the potential to provide valuable information for studies of adaptation. In this study, we sequenced six individual transcriptomes of T. mongolica from heterogeneous habitats, focused on SSRs located in genes, and identified 811 polymorphic SSRs. Of the identified SSRs, 172, 470, and 76 were located in 5' UTRs, CDSs, and 3' UTRs in 591 transcripts; and AG/CT, AAC/GTT, and AT/AT were the most abundant repeats in each gene region. Functional annotation showed that many of the identified polymorphic SSRs were in genes that were enriched in several GO terms and KEGG pathways, suggesting the functional significance of these genes in the environmental adaptation process. The identification of polymorphic genic SSRs in our study lays a foundation for future studies investigating the contribution of SSRs to regulation of genes in natural populations of T. mongolica and their importance for adaptive evolution of this species.

Fatty acid-binding protein 4 downregulation drives calcification in the development of kidney stone disease

Nephrolithiasis is a significant source of morbidity, and its incidence has increased significantly over the last decades. This rise has been attributed to concurrent increasing rates of obesity, associated with a 3-time risk of developing NL. To date, the mechanism by which obesity is linked to stone formation has not been elucidated. We aimed to utilize a transcriptomics approach to discover the missing link between these two epidemic diseases. We investigated gene expression profiling of nephrolithiasis patients by two RNA-sequencing approaches: comparison between renal papilla tissue with and without the presence of calcified Randall's plaques (RP), and comparison between the papilla, medulla, and cortex regions from within a single recurrent stone forming kidney. Results were overlaid between differently expressed genes found in the patient cohort and in the severely lithogenic kidney to identify common genes. Overlay of these two RNA-sequencing datasets demonstrated there is impairment of lipid metabolism in renal papilla tissue containing RP linked to downregulation of fatty acid binding protein (FABP) 4. Immunohistochemistry of human kidney specimens and microarray analysis of renal tissue from a nephrolithiasis mouse model confirmed that FABP4 downregulation is associated with renal stone formation. In a FABP4 knockout mouse model, FABP4 deficiency resulted in development of both renal and urinary crystals. Our study revealed that FABP4 plays an important, previously unrecognized role in kidney stone formation, providing a feasible mechanism to explain the link between nephrolithiasis and metabolic syndrome.

Comparison of procedures for RNA-extraction from peripheral blood mononuclear cells

RNA quality and quantity are important factors for ensuring the accuracy of gene expression analysis and other RNA-based downstream applications. Thus far, only a limited number of methodological studies have compared sample storage and RNA extraction procedures for human cells. We compared three commercially available RNA extraction kits, i.e., (NucliSENS) easyMAG, RNeasy (Mini Kit) and RiboPure (RNA Purification Kit–blood). In addition, additional conditions, such as storage medium and storage temperature of human peripheral blood mononuclear cells were evaluated, i.e., 4 °C for RNAlater or -80 °C for QIAzol and for the respective cognate lysis buffers; easyMAG, RNeasy or RiboPure. RNA was extracted from aliquots that had been stored for one day (Run 1) or 83 days (Run 2). After DNase treatment, quantity and quality of RNA were assessed by means of a NanoDrop spectrophotometer, 2100 Bioanalyzer and RT-qPCR for the ACTB reference gene. We observed that high-quality RNA can be obtained using RNeasy and RiboPure, regardless of the storage medium, whereas samples stored in RNAlater resulted in the least amount of RNA extracted. In addition, RiboPure combined with storage of samples in its cognate lysis buffer yielded twice as much RNA as all other procedures. These results were supported by RT-qPCR and by the reproducibility observed for two independent extraction runs.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

The utility of reptile blood transcriptomes in molecular ecology

Reptiles and other nonmammalian vertebrates have transcriptionally active nucleated red blood cells. If blood transcriptomes can provide quantitative data to address questions relevant to molecular ecology, this could circumvent the need to euthanize animals to assay tissues. This would allow longitudinal sampling of animals' responses to treatments, as well as sampling of protected taxa. We developed and annotated blood transcriptomes from six reptile species and found on average 25,000 proteins are being transcribed in the blood, and there is a CORE group of 9,282 orthogroups that are found in at least four of six species. In comparison to liver transcriptomes from the same taxa, approximately two-thirds of the orthogroups were found in both blood and liver; and a similar percentage of ecologically relevant gene groups (insulin and insulin-like signalling, electron transport chain, oxidative stress, glucocorticoid receptors) were found transcribed in both blood and liver. As a resource, we provide a user-friendly database of gene ids identified in each blood transcriptome. Although on average 37% of reads mapped to haemoglobin, importantly, the majority of nonhaemoglobin transcripts had sufficient depth (e.g., 97% at ≥10 reads) to be included in differential gene expression analysis. Thus, we demonstrate that RNAseq blood transcriptomes from a very small blood sample (<10 μl) is a minimally invasive option in nonmammalian vertebrates for quantifying expression of a large number of ecologically relevant genes that would allow longitudinal sampling and sampling of protected populations.

Whole transcriptome analysis of an estuarine amphipod exposed to highway road dust

Urban road dust can potentially have adverse effects on aquatic and benthic ecosystems if discharged into receiving waters; however, little is known about the mode of action of road dust toxicity within aquatic organisms. With an aim to reveal the biological effects of road dust on benthic crustacean species, we performed a de novo transcriptome analysis of the estuarine amphipod Grandidierella japonica exposed to road dust collected from highways around Tokyo. A transcriptome analysis by Illumina HiSeq 2500 identified differentially expressed genes related to the gamma-aminobutyric acid (GABA) signaling pathway, oxidative damage, and cuticle metabolism. Among these, a GABA_B receptor subunit showed down-regulation in the road dust treatment, but a constant expression in the treatment of road dust with a carbonaceous resin XAD-4, which can reduce the acute toxicity of road dust to G. japonica. These results and the time course expressions of the related genes were partially confirmed by quantitative PCR (qPCR) experiments. Although the linkage between acute lethal toxicity and the molecular initiating events induced by road dust was still unclear, our findings provide lines of evidence to identify the causative toxicants in urban road dust.

A Phylogenomic Supertree of Birds

It has long been appreciated that analyses of genomic data (e.g., whole genome sequencing or sequence capture) have the potential to reveal the tree of life, but it remains challenging to move from sequence data to a clear understanding of evolutionary history, in part due to the computational challenges of phylogenetic estimation using genome-scale data. Supertree methods solve that challenge because they facilitate a divide-and-conquer approach for large-scale phylogeny inference by integrating smaller subtrees in a computationally efficient manner. Here, we combined information from sequence capture and whole-genome phylogenies using supertree methods. However, the available phylogenomic trees had limited overlap so we used taxon-rich (but not phylogenomic) megaphylogenies to weave them together. This allowed us to construct a phylogenomic supertree, with support values, that included 707 bird species (~7% of avian species diversity). We estimated branch lengths using mitochondrial sequence data and we used these branch lengths to estimate divergence times. Our time-calibrated supertree supports radiation of all three major avian clades (Palaeognathae, Galloanseres, and Neoaves) near the Cretaceous-Paleogene (K-Pg) boundary. The approach we used will permit the continued addition of taxa to this supertree as new phylogenomic data are published, and it could be applied to other taxa as well.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

Comparative transcriptome analysis reveals potential evolutionary differences in adaptation of temperature and body shape among four Percidae species

Considering the divergent temperature habitats and morphological traits of four Percidae species: yellow perch (Perca flavescens), Eurasian perch (Perca fluviatilis), pike perch (Sander lucioperca), and ruffe (Gymnocephalus cernua), we stepped into the transcriptome level to discover genes and mechanisms that drive adaptation to different temperature environments and evolution in body shape. Based on 93,566 to 181,246 annotated unigenes of the four species, we identified 1,117 one-to-one orthologous genes and subsequently constructed the phylogenetic trees that are consistent with previous studies. Together with the tree, the ratios of nonsynonymous to synonymous substitutions presented decreased evolutionary rates from the D. rerio branch to the sub-branch clustered by P. flavescens and P. fluviatilis. The specific 93 fast-evolving genes and 57 positively selected genes in P. flavescens, compared with 22 shared fast-evolving genes among P. fluviatilis, G. cernua, and S. lucioperca, showed an intrinsic foundation that ensure its adaptation to the warmer Great Lakes and farther south, especially in functional terms like “Cul4-RING E3 ubiquitin ligase complex.” Meanwhile, the specific 78 fast-evolving genes and 41 positively selected genes in S. lucioperca drew a clear picture of how it evolved to a large and elongated body with camera-type eyes and muscle strength so that it could occupy the highest position in the food web. Overall, our results uncover genetic basis that support evolutionary adaptation of temperature and body shape in four Percid species, and could furthermore assist studies on environmental adaptation in fishes.

RNA-seq data from mature male gonads of marine mussels Mytilus edulis and M. galloprovincialis

The mussels Mytilus edulis and Mytilus galloprovincialis are marine organisms with external fertilization able to hybridize where their distributions overlap allowing the study of reproductive isolation mechanisms in nature. We provide raw data of a transcriptomic analysis of mature male gonads from these two Mytilus spp. using NGS (Illumina) technology and a preliminary list of transcript that were functionally annotated showing species-specific differential expression. A shortlist including some of these genes and their corresponding proteins have been thoroughly analysed and discussed in Romero et al. (2018, Submitted for publication).

De novo transcriptome assembly and differential gene expression analysis of the calanoid copepod Acartia tonsa exposed to nickel nanoparticles

The calanoid copepod Acartia tonsa is a reference species in standardized ecotoxicology bioassay. Despite this interest, there is a lack of knowledge on molecular responses of A. tonsa to contaminants. We generated a de novo assembled transcriptome of A. tonsa exposed 4 days to 8.5 and 17 mg/L nickel nanoparticles (NiNPs), which have been shown to reduce egg hatching success and larval survival but had no effects on the adults. Aims of our study were to 1) improve the knowledge on the molecular responses of A. tonsa copepod and 2) increase the genomic resources of this copepod for further identification of potential biomarkers of NP exposure. The de novo assembled transcriptome of A. tonsa consisted of 53,619 unigenes, which were further annotated to nr, GO, KOG and KEGG databases. In particular, most unigenes were assigned to Metabolic and Cellular processes (34-45%) GO terms, and to Human disease (28%) and Organismal systems (23%) KEGG categories. Comparison among treatments showed that 373 unigenes were differentially expressed in A. tonsa exposed to NiNPs at 8.5 and 17 mg/L, with respect to control. Most of these genes were downregulated and took part in ribosome biogenesis, translation and protein turnover, thus suggesting that NiNPs could affect the copepod ribosome synthesis machinery and functioning. Overall, our study highlights the potential of toxicogenomic approach in gaining more mechanistic and functional information about the mode of action of emerging compounds on marine organisms, for biomarker discovering in crustaceans.

Molecular and cellular characterization of apoptosis in flat oyster a key mechanisms at the heart of host-parasite interactions

Bonamia ostreae has been associated with the decline of flat oyster Ostrea edulis populations in some European countries. This obligatory intracellular parasite persists and multiplies into hemocytes. Previous in vitro experiments showed that apoptosis is activated in hemocytes between 1 h and 4 h of contact with the parasite. The flat oyster uses the apoptosis pathway to defend against B. ostreae. However, the parasite might be also able to modulate this response in order to survive in its host. In order to investigate this hypothesis the apoptotic response of the host was evaluated using flow cytometry, transmission electron microscopy and by measuring the response of genes involved in the apoptotic pathway after 4 h. In parallel, the parasite response was investigated by measuring the expression of B. ostreae genes involved in different biological functions including cell cycle and cell death. Obtained results allow describing molecular apoptotic pathways in O. edulis and confirm that apoptosis is early activated in hemocytes after a contact with B. ostreae. Interestingly, at cellular and molecular levels this process appeared downregulated after 44 h of contact. Concurrently, parasite gene expression appeared reduced suggesting that the parasite could inhibit its own metabolism to escape the immune response.

Creating a functional single-chromosome yeast

Eukaryotic genomes are generally organized in multiple chromosomes. Here we have created a functional single-chromosome yeast from a Saccharomyces cerevisiae haploid cell containing sixteen linear chromosomes, by successive end-to-end chromosome fusions and centromere deletions. The fusion of sixteen native linear chromosomes into a single chromosome results in marked changes to the global three-dimensional structure of the chromosome due to the loss of all centromere-associated inter-chromosomal interactions, most telomere-associated inter-chromosomal interactions and 67.4% of intra-chromosomal interactions. However, the single-chromosome and wild-type yeast cells have nearly identical transcriptome and similar phenome profiles. The giant single chromosome can support cell life, although this strain shows reduced growth across environments, competitiveness, gamete production and viability. This synthetic biology study demonstrates an approach to exploration of eukaryote evolution with respect to chromosome structure and function.

Phylogenomics of palearctic Formica species suggests a single origin of temporary parasitism and gives insights to the evolutionary pathway toward slave-making behaviour

Background

The ants of the Formica genus are classical model species in evolutionary biology. In particular, Darwin used Formica as model species to better understand the evolution of slave-making, a parasitic behaviour where workers of another species are stolen to exploit their workforce. In his book “On the Origin of Species” (1859), Darwin first hypothesized that slave-making behaviour in Formica evolved in incremental steps from a free-living ancestor.

Methods

The absence of a well-resolved phylogenetic tree of the genus prevent an assessment of whether relationships among Formica subgenera are compatible with this scenario. In this study, we resolve the relationships among the 4 palearctic Formica subgenera (Formica str. s., Coptoformica, Raptiformica and Serviformica) using a phylogenomic dataset of 945 genes for 16 species.

Results

We provide a reference tree resolving the relationships among the main Formica subgenera with high bootstrap supports.

Discussion

The branching order of our tree suggests that the free-living lifestyle is ancestral in the Formica genus and that parasitic colony founding could have evolved a single time, probably acting as a pre-adaptation to slave-making behaviour.

Conclusion

This phylogenetic tree provides a solid backbone for future evolutionary studies in the Formica genus and slave-making behaviour.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1159-4) contains supplementary material, which is available to authorized users.

Tidal and diel orchestration of behaviour and gene expression in an intertidal mollusc

Intertidal inhabitants are exposed to the 24-hour solar day, and the 12.4 hour rising and falling of the tides. One or both of these cycles govern intertidal organisms' behaviour and physiology, yet little is known about the molecular clockworks of tidal rhythmicity. Here, we show that the limpet Cellana rota exhibits robust tidally rhythmic behaviour and gene expression. We assembled a de-novo transcriptome, identifying novel tidal, along with known circadian clock genes. Surprisingly, most of the putative circadian clock genes, lack a typical rhythmicity. We identified numerous tidally rhythmic genes and pathways commonly associated with the circadian clock. We show that not only is the behaviour of an intertidal organism in tune with the tides, but so too are many of its genes and pathways. These findings highlight the plasticity of biological timekeeping in nature, strengthening the growing notion that the role of 'canonical' circadian clock genes may be more fluid than previously thought, as exhibited in an organism which has evolved in an environment where tidal oscillations are the dominant driving force.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Optimizing RNA Extraction of Renal Papilla Biopsy Tissue in Kidney Stone Formers: A New Methodology for Genomic Study

INTRODUCTION

Endoscopic tools have provided versatile examination and treatment for kidney stone procedures. Despite endourologists researching urinary stone disease using endoscopes to collect tissue, this tissue collection method is limited. Endoscopically removed tissues are small in size, restricting the types of genome-based examination possible. We investigated a new method of renal papilla biopsy and RNA extraction to establish a genomic research methodology for kidney stone disease.

PATIENTS AND METHODS

We conducted a prospective multi-institutional study and collected renal papilla specimens from consecutive percutaneous nephrolithotomy and ureteroscopy (URS) cases performed for removal of upper urinary tract stones. Renal papilla tissue was extracted using ureteroscopic biopsy forceps after stone removal. RNA was extracted using two different extraction kits, and their quantity and quality were examined. Additionally, the impact of biopsy on surgical complications was compared between cases performed with and without biopsy by matched case-control analysis adjusted for age, gender, body mass index, bilaterality, and stone burden.

RESULTS

A total of 90 biopsies from 49 patients were performed, and the median duration between specimen collection and RNA extraction was 61 days. Both univariate and multivariate analyses showed BIGopsy^® forceps usage significantly increased the total yield (p = 0.004) and quality (p = 0.001 for A260/280, p = 0.004 for A260/A230) of extracted RNA. Extraction using the RNeasy Micro Kit^® also improved A260/A230, whereas reduced RNA integrity number of extracted RNA by univariate and multivariate analyses (p = 0.002 and p < 0.001, respectively). Moreover, matched case-control study demonstrated that endoscopic renal papilla biopsy caused no significant surgical complications, including bleeding, decreased stone clearance and hematocrit, and renal dysfunction. Biopsies during URS imparted an average of 20 minutes of procedure time over nonbiopsy cases.

CONCLUSIONS

We demonstrate a safe methodology for optimal RNA extraction of renal papilla tissue. This technique will accelerate advanced genomic studies for kidney stone formers by facilitating larger tissue yields.

Application of cDNA-AFLP to biomarker exploration in a non-model species Grandidierella japonica

Biomarkers of exposure can be used to identify specific contaminants that are adversely affecting aquatic organisms. However, it remains prohibitively costly to investigate multiple novel biomarkers of exposure in a non-model species, despite the development of next-generation sequencing technology. In this study, we focused on the use of cDNA-amplified fragment length polymorphism (AFLP) as a cost-effective biomarker discovery tool to test whether it could identify biomarkers of exposure in the non-model amphipod species Grandidierella japonica. Loci were identified that were differentially expressed in amphipods exposed to reference chemicals (Cu, Zn, and nicotine) and to an environmental sample (road dust) at sublethal concentrations. Eight loci were shown to respond consistently to nicotine at different concentrations, but not to Cu or Zn. Some of the loci also responded to an environmental road dust sample containing nicotine. These findings suggest that loci identified using cDNA-AFLP could be used as biomarkers of nicotine exposure in environmental samples with complex matrices. Further studies with other organisms and toxicants are needed, but we have demonstrated that the use of cDNA-AFLP to identify biomarkers for ecotoxicological studies of non-model species is at least feasible.

Patterns of cross-contamination in a multispecies population genomic project: detection, quantification, impact, and solutions

Background

Contamination is a well-known but often neglected problem in molecular biology. Here, we investigated the prevalence of cross-contamination among 446 samples from 116 distinct species of animals, which were processed in the same laboratory and subjected to subcontracted transcriptome sequencing.

Results

Using cytochrome oxidase 1 as a barcode, we identified a minimum of 782 events of between-species contamination, with approximately 80% of our samples being affected. An analysis of laboratory metadata revealed a strong effect of the sequencing center: nearly all the detected events of between-species contamination involved species that were sent the same day to the same company. We introduce new methods to address the amount of within-species, between-individual contamination, and to correct for this problem when calling genotypes from base read counts.

Conclusions

We report evidence for pervasive within-species contamination in this data set, and show that classical population genomic statistics, such as synonymous diversity, the ratio of non-synonymous to synonymous diversity, inbreeding coefficient F_IT, and Tajima’s D, are sensitive to this problem to various extents. Control analyses suggest that our published results are probably robust to the problem of contamination. Recommendations on how to prevent or avoid contamination in large-scale population genomics/molecular ecology are provided based on this analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0366-6) contains supplementary material, which is available to authorized users.

High-quality RNA extraction from the sea urchin Paracentrotus lividus embryos

The sea urchin Paracentrotus lividus (Lamarck, 1816) is a keystone herbivore in the Mediterranean Sea due to its ability to transform macroalgal-dominated communities into barren areas characterized by increased cover of bare substrates and encrusting coralline algae, reduced biodiversity and altered ecosystem functions. P. lividus is also an excellent animal model for toxicology, physiology and biology investigations having been used for more than a century as a model for embryological studies with synchronously developing embryos which are easy to manipulate and analyze for morphological aberrations. Despite its importance for the scientific community, the complete genome is still not fully annotated. To date, only a few molecular tools are available and a few Next Generation Sequencing (NGS) studies have been performed. Here we aimed at setting-up an RNA extraction method to obtain high quality and sufficient quantity of RNA for NGS from P. lividus embryos at the pluteus stage. We compared five different RNA extraction protocols from four different pools of plutei (500, 1000, 2500 and 5000 embryos): TRIzol®, and four widely-used Silica Membrane kits, GenElute™ Mammalian Total RNA Miniprep Kit, RNAqueous® Micro Kit, RNeasy® Micro Kit and Aurum™ Total RNA Mini Kit. The quantity of RNA isolated was evaluated using NanoDrop. The quality, considering the purity, was measured as A260/A280 and A260/230 ratios. The integrity was measured by RNA Integrity Number (RIN). Our results demonstrated that the most efficient procedures were GenElute, RNeasy and Aurum, producing a sufficient quantity of RNA for NGS. The Bioanalyzer profiles and RIN values revealed that the most efficient methods guaranteeing for RNA integrity were RNeasy and Aurum combined with an initial preservation in RNAlater. This research represents the first attempt to standardize a method for high-quality RNA extraction from sea urchin embryos at the pluteus stage, providing a new resource for this established model marine organism.

Positive selection on sperm ion channels in a brooding brittle star: consequence of life-history traits evolution

Closely related species are key models to investigate mechanisms leading to reproductive isolation and early stages of diversification, also at the genomic level. The brittle star cryptic species complex Ophioderma longicauda encompasses the sympatric broadcast-spawning species C3 and the internal brooding species C5. Here, we used de novo transcriptome sequencing and assembly in two closely related species displaying contrasting reproductive modes to compare their genetic diversity and to investigate the role of natural selection in reproductive isolation. We reconstructed 20 146 and 22 123 genes for C3 and C5, respectively, and characterized a set of 12 229 orthologs. Genetic diversity was 1.5-2 times higher in C3 compared to C5, confirming that species with low parental investment display higher levels of genetic diversity. Forty-eight genes were the targets of positive diversifying selection during the evolution of the two species. Notably, two genes (NHE and TetraKCNG) are sperm-specific ion channels involved in sperm motility. Ancestral sequence reconstructions show that natural selection targeted the two genes in the brooding species. This may result from an adaptation to the novel environmental conditions surrounding sperm in the brooding species, either directly affecting sperm or via an increase in male/female conflict. This phenomenon could have promoted prezygotic reproductive isolation between C3 and C5. Finally, the sperm receptors to egg chemoattractants differed between C3 and C5 in the ligand-binding region. We propose that mechanisms of species-specific gamete recognition in brittle stars occur during sperm chemotaxis (sperm attraction towards the eggs), contrary to other marine invertebrates where prezygotic barriers to interspecific hybridization typically occur before sperm-egg fusion.

Convergent evolution of social hybridogenesis in Messor harvester ants

Sexual reproduction generally requires no more than two partners. Here, we show convergent evolution of social hybridogenesis, a reproductive system requiring three reproductive partners in harvester ants. In this unorthodox reproductive system, two distinct genetic lineages live in sympatry and queens have to mate with males of their own lineage to produce queens along with males of the alternative lineage to produce workers. Using a large transcriptomic data set of nine species, we show that social hybridogenesis evolved at least three times independently in the genus Messor. Moreover, a study of 13 populations of Messor barbarus revealed that this mode of reproduction is fixed in the whole range of this ecologically dominant species. Finally, we show that workers can produce males carrying genes of the two genetic lineages, raising the possibility of rare gene flow between lineages contributing to the long-term maintenance of pairs of interdependent lineages. These results emphasize the evolutionary importance of social hybridogenesis, a major transition possibly linked to the peculiar ecology of harvester ants.

Hemizygosity Enhances Purifying Selection: Lack of Fast-Z Evolution in Two Satyrine Butterflies

The fixation probability of a recessive beneficial mutation is increased on the X or Z chromosome, relative to autosomes, because recessive alleles carried by X or Z are exposed to selection in the heterogametic sex. This leads to an increased dN/dS ratio on sex chromosomes relative to autosomes, a pattern called the “fast-X” or “fast-Z” effect. Besides positive selection, the strength of genetic drift and the efficacy of purifying selection, which affect the rate of molecular evolution, might differ between sex chromosomes and autosomes. Disentangling the complex effects of these distinct forces requires the genome-wide analysis of polymorphism, divergence and gene expression data in a variety of taxa. Here we study the influence of hemizygosity of the Z chromosome in Maniola jurtina and Pyronia tithonus, two species of butterflies (Lepidoptera, Nymphalidae, Satyrinae). Using transcriptome data, we compare the strength of positive and negative selection between Z and autosomes accounting for sex-specific gene expression. We show that M. jurtina and P. tithonus do not experience a faster, but rather a slightly slower evolutionary rate on the Z than on autosomes. Our analysis failed to detect a significant difference in adaptive evolutionary rate between Z and autosomes, but comparison of male-biased, unbiased and female-biased Z-linked genes revealed an increased efficacy of purifying selection against recessive deleterious mutations in female-biased Z-linked genes. This probably contributes to the lack of fast-Z evolution of satyrines. We suggest that the effect of hemizygosity on the fate of recessive deleterious mutations should be taken into account when interpreting patterns of molecular evolution in sex chromosomes vs. autosomes.

Discovery of parvovirus-related sequences in an unexpected broad range of animals

Our knowledge of the genetic diversity and host ranges of viruses is fragmentary. This is particularly true for the Parvoviridae family. Genetic diversity studies of single stranded DNA viruses within this family have been largely focused on arthropod- and vertebrate-infecting species that cause diseases of humans and our domesticated animals: a focus that has biased our perception of parvovirus diversity. While metagenomics approaches could help rectify this bias, so too could transcriptomics studies. Large amounts of transcriptomic data are available for a diverse array of animal species and whenever this data has inadvertently been gathered from virus-infected individuals, it could contain detectable viral transcripts. We therefore performed a systematic search for parvovirus-related sequences (PRSs) within publicly available transcript, genome and protein databases and eleven new transcriptome datasets. This revealed 463 PRSs in the transcript databases of 118 animals. At least 41 of these PRSs are likely integrated within animal genomes in that they were also found within genomic sequence databases. Besides illuminating the ubiquity of parvoviruses, the number of parvoviral sequences discovered within public databases revealed numerous previously unknown parvovirus-host combinations; particularly in invertebrates. Our findings suggest that the host-ranges of extant parvoviruses might span the entire animal kingdom.

A replicated climate change field experiment reveals rapid evolutionary response in an ecologically important soil invertebrate

Whether species can respond evolutionarily to current climate change is crucial for the persistence of many species. Yet, very few studies have examined genetic responses to climate change in manipulated experiments carried out in natural field conditions. We examined the evolutionary response to climate change in a common annelid worm using a controlled replicated experiment where climatic conditions were manipulated in a natural setting. Analyzing the transcribed genome of 15 local populations, we found that about 12% of the genetic polymorphisms exhibit differences in allele frequencies associated to changes in soil temperature and soil moisture. This shows an evolutionary response to realistic climate change happening over short-time scale, and calls for incorporating evolution into models predicting future response of species to climate change. It also shows that designed climate change experiments coupled with genome sequencing offer great potential to test for the occurrence (or lack) of an evolutionary response.

Adaptive Protein Evolution in Animals and the Effective Population Size Hypothesis

The rate at which genomes adapt to environmental changes and the prevalence of adaptive processes in molecular evolution are two controversial issues in current evolutionary genetics. Previous attempts to quantify the genome-wide rate of adaptation through amino-acid substitution have revealed a surprising diversity of patterns, with some species (e.g. Drosophila) experiencing a very high adaptive rate, while other (e.g. humans) are dominated by nearly-neutral processes. It has been suggested that this discrepancy reflects between-species differences in effective population size. Published studies, however, were mainly focused on model organisms, and relied on disparate data sets and methodologies, so that an overview of the prevalence of adaptive protein evolution in nature is currently lacking. Here we extend existing estimators of the amino-acid adaptive rate by explicitly modelling the effect of favourable mutations on non-synonymous polymorphism patterns, and we apply these methods to a newly-built, homogeneous data set of 44 non-model animal species pairs. Data analysis uncovers a major contribution of adaptive evolution to the amino-acid substitution process across all major metazoan phyla—with the notable exception of humans and primates. The proportion of adaptive amino-acid substitution is found to be positively correlated to species effective population size. This relationship, however, appears to be primarily driven by a decreased rate of nearly-neutral amino-acid substitution because of more efficient purifying selection in large populations. Our results reveal that adaptive processes dominate the evolution of proteins in most animal species, but do not corroborate the hypothesis that adaptive substitutions accumulate at a faster rate in large populations. Implications regarding the factors influencing the rate of adaptive evolution and positive selection detection in humans vs. other organisms are discussed.

The rate at which species adapt to environmental changes is a controversial topic. The theory predicts that adaptation is easier in large than in small populations, and the genomic studies of model organisms have revealed a much higher adaptive rate in large population-sized flies than in small population-sized humans and apes. Here we build and analyse a large data set of protein-coding sequences made of thousands of genes in 44 pairs of species from various groups of animals including insects, molluscs, annelids, echinoderms, reptiles, birds, and mammals. Extending and improving existing data analysis methods, we show that adaptation is a major process in protein evolution across all phyla of animals: the proportion of amino-acid substitutions that occurred adaptively is above 50% in a majority of species, and reaches up to 90%. Our analysis does not confirm that population size, here approached through species genetic diversity and ecological traits, does influence the rate of adaptive molecular evolution, but points to human and apes as a special case, compared to other animals, in terms of adaptive genomic processes.

Molecular biomarkers for chronological age in animal ecology

The chronological age of an individual animal predicts many of its biological characteristics, and these in turn influence population-level ecological processes. Animal age information can therefore be valuable in ecological research, but many species have no external features that allow age to be reliably determined. Molecular age biomarkers provide a potential solution to this problem. Research in this area of molecular ecology has so far focused on a limited range of age biomarkers. The most commonly tested molecular age biomarker is change in average telomere length, which predicts age well in a small number of species and tissues, but performs poorly in many other situations. Epigenetic regulation of gene expression has recently been shown to cause age-related modifications to DNA and to cause changes in abundance of several RNA types throughout animal lifespans. Age biomarkers based on these epigenetic changes, and other new DNA-based assays, have already been applied to model organisms, humans and a limited number of wild animals. There is clear potential to apply these marker types more widely in ecological studies. For many species, these new approaches will produce age estimates where this was previously impractical. They will also enable age information to be gathered in cross-sectional studies and expand the range of demographic characteristics that can be quantified with molecular methods. We describe the range of molecular age biomarkers that have been investigated to date and suggest approaches for developing the newer marker types as age assays in nonmodel animal species.

Comparative Analysis of Transcriptomes from Secondary Reproductives of Three Reticulitermes Termite Species

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.

Evaluation of the impact of RNA preservation methods of spiders for de novo transcriptome assembly

With advances in high-throughput sequencing technologies, de novo transcriptome sequencing and assembly has become a cost-effective method to obtain comprehensive genetic information of a species of interest, especially in nonmodel species with large genomes such as spiders. However, high-quality RNA is essential for successful sequencing, and sample preservation conditions require careful consideration for the effective storage of field-collected samples. To this end, we report a streamlined feasibility study of various storage conditions and their effects on de novo transcriptome assembly results. The storage parameters considered include temperatures ranging from room temperature to -80°C; preservatives, including ethanol, RNAlater, TRIzol and RNAlater-ICE; and sample submersion states. As a result, intact RNA was extracted and assembly was successful when samples were preserved at low temperatures regardless of the type of preservative used. The assemblies as well as the gene expression profiles were shown to be robust to RNA degradation, when 30 million 150-bp paired-end reads are obtained. The parameters for sample storage, RNA extraction, library preparation, sequencing and in silico assembly considered in this work provide a guideline for the study of field-collected samples of spiders.