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

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Metatranscriptomic holobiont analysis of carbohydrate-active enzymes in the millipede Telodeinopus aoutii (Diplopoda, Spirostreptida)

As important decomposers of soil organic matter, millipedes contribute to lignocellulose decomposition and nutrient cycling. The degradation of lignocellulose requires the action of several carbohydrate-active enzymes (CAZymes) and, in most invertebrates, depends on the activity of mutualistic gut microorganisms. To address the question of the importance of the microbiota and endogenous (host) enzymes in digestive processes in millipedes, we analyzed metatranscriptomic data from the tropical millipede Telodeinopus aoutii at the holobiont level. Functional annotation included identification of expressed CAZymes (CAZy families and EC terms) in the host and its intestinal microbiota, foregut, midgut, and hindgut, compared to non-intestinal tissues. Most of the 175 CAZy families were expressed exclusively in the gut microbiota and more than 50% of these microbial families were expressed exclusively in the hindgut. The greatest diversity of expressed endogenous CAZymes from all gut sections was found in the midgut (77 families). Bacteria were the major microbial producers of CAZymes, Proteobacteria dominating in the midgut and Bacteriodetes with Firmicutes in the hindgut. The contribution of the eukaryotic microbiota to CAZymes production was negligible. Functional classification of expressed CAZy families confirmed a broad functional spectrum of CAZymes potentially expressed in the holobiont. Degradation of lignocellulose in the digestive tract of the millipede T. aoutii depends largely on bacterial enzymes expressed in the hindgut. Endogenous cellulases were not detected, except for the potentially cellulolytic family AA15, but an expression of cellulolytic enzymes of this family was not confirmed at the EC-number level. The midgut had the greatest diversity of expressed endogenous CAZymes, mainly amylases, indicating the importance of digesting α-glucosidases for the millipede. In contrast, bacterial lignocellulolytic enzymes are sparsely expressed here. The hindgut was the hotspot of microbial degradation of cellulose and hemicellulases. The gain of the millipede from the microbial lignocellulose degradation in the gut, and consequently the mutualistic status of the relationship between the millipede and its cellulolytic gut bacteria, depends on the ability of the millipede to take up microbial metabolites as nutrients through the hindgut wall. Enzymes expressed in the intestine can degrade all components of lignocellulose except lignin. Assuming that soil microbiota is partially degraded lignin in the millipede diet, T. aoutii can be considered a decomposer of soil organic matter relying primarily on its gut bacteria. The deposition of millipede fecal pellets containing an organic matter modified by the hindgut bacterial community could be of ecological significance.

Differential Selection on Caste-Associated Genes in a Subterranean Termite

Analyzing the information-rich content of RNA can help uncover genetic events associated with social insect castes or other social polymorphisms. Here, we exploit a series of cDNA libraries previously derived from whole-body tissue of different castes as well as from three behaviourally distinct populations of the Eastern subterranean termite Reticulitermes flavipes. We found that the number (~0.5 M) of single nucleotide variants (SNVs) was roughly equal between nymph, worker and soldier caste libraries, but dN/dS (ratio of nonsynonymous to synonymous substitutions) analysis suggested that some of these variants confer a caste-specific advantage. Specifically, the dN/dS ratio was high (~4.3) for genes expressed in the defensively specialized soldier caste, relative to genes expressed by other castes (~1.7−1.8) and regardless of the North American population (Toronto, Raleigh, Boston) from which the castes were sampled. The populations, meanwhile, did show a large difference in SNV count but not in the manner expected from known demographic and behavioural differences; the highly invasive unicolonial population from Toronto was not the least diverse and did not show any other unique substitution patterns, suggesting any past bottleneck associated with invasion or with current unicoloniality has become obscured at the RNA level. Our study raises two important hypotheses relevant to termite sociobiology. First, the positive selection (dN/dS > 1) inferred for soldier-biased genes is presumably indirect and of the type mediated through kin selection, and second, the behavioural changes that accompany some social insect urban invasions (i.e., ‘unicoloniality’) may be detached from the loss-of-diversity expected from invasion bottlenecks.

A Century of Synergy in Termite Symbiosis Research: Linking the Past with New Genomic Insights

Termites have long been studied for their symbiotic associations with gut microbes. In the late nineteenth century, this relationship was poorly understood and captured the interest of parasitologists such as Joseph Leidy; this research led to that of twentieth-century biologists and entomologists including Cleveland, Hungate, Trager, and Lüscher. Early insights came via microscopy, organismal, and defaunation studies, which led to descriptions of microbes present, descriptions of the roles of symbionts in lignocellulose digestion, and early insights into energy gas utilization by the host termite. Focus then progressed to culture-dependent microbiology and biochemical studies of host-symbiont complementarity, which revealed specific microhabitat requirements for symbionts and noncellulosic mechanisms of symbiosis (e.g., N₂ fixation). Today, knowledge on termite symbiosis has accrued exponentially thanks to omic technologies that reveal symbiont identities, functions, and interdependence, as well as intricacies of host-symbiont complementarity. Moving forward, the merging of classical twentieth-century approaches with evolving omic tools should provide even deeper insights into host-symbiont interplay.

The de novo transcriptome of workers head of the higher group termite Globitermes sulphureus Haviland (Blattodea: Termitidae)

The subterranean termite Globitermus sulphureus is an important Southeast Asian pest with limited genomic resources that causes damages to agriculture crops and building structures. Therefore, the main goal of this study was to survey the G. sulphureus transcriptome composition. Here, we performed de novo transcriptome for G. sulphureus workers’ heads using Illumina HiSeq paired-end sequencing technology. A total of 88, 639, 408 clean reads were collected and assembled into 243, 057 transcripts and 193, 344 putative genes. The transcripts were annotated with the Trinotate pipeline. In total, 27, 061 transcripts were successfully annotated using BLASTX against the SwissProt database and 17, 816 genes were assigned to 47, 598 GO terms. We classified 14, 223 transcripts into COG classification, resulting in 25 groups of functional annotations. Next, a total of 12, 194 genes were matched in the KEGG pathway and 392 metabolic pathways were predicted based on the annotation. Moreover, we detected two endogenous cellulases in the sequences. The RT-qPCR analysis showed that there were significant differences in the expression levels of two genes β-glucosidase and endo-β-1,4-glucanase between worker and soldier heads of G. sulphureus. This is the first study to characterize the complete head transcriptome of a higher termite G. sulphureus using a high-throughput sequencing. Our study may provide an overview and comprehensive molecular resource for comparative studies of the transcriptomics and genomics of termites.

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.

Comparative transcriptome and potential antiviral signaling pathways analysis of the gills in the red swamp crayfish, Procambarus clarkii infected with White Spot Syndrome Virus (WSSV)

Red swamp crayfish is an important model organism for research of the invertebrate innate immunity mechanism. Its excellent disease resistance against bacteria, fungi, and viruses is well-known. However, the antiviral mechanisms of crayfish remain unclear. In this study, we obtained high-quality sequence reads from normal and white spot syndrome virus (WSSV)-challenged crayfish gills. For group normal (GN), 39,390,280 high-quality clean reads were randomly assembled to produce 172,591 contigs; whereas, 34,011,488 high-quality clean reads were randomly assembled to produce 182,176 contigs for group WSSV-challenged (GW). After GO annotations analysis, a total of 35,539 (90.01%), 14,931 (37.82%), 28,221 (71.48%), 25,290 (64.05%), 15,595 (39.50%), and 13,848 (35.07%) unigenes had significant matches with sequences in the Nr, Nt, Swiss-Prot, KEGG, COG and GO databases, respectively. Through the comparative analysis between GN and GW, 12,868 genes were identified as differentially up-regulated DEGs, and 9,194 genes were identified as differentially down-regulated DEGs. Ultimately, these DEGs were mapped into different signaling pathways, including three important signaling pathways related to innate immunity responses. These results could provide new insights into crayfish antiviral immunity mechanism.

In-depth comparative transcriptome analysis of intestines of red swamp crayfish, Procambarus clarkii, infected with WSSV

Crayfish has become one of the most important farmed aquatic species in China due to its excellent disease resistance against bacteria and viruses. However, the antiviral mechanism of crayfish is still not very clear. In the present study, many high-quality sequence reads from crayfish intestine were obtained using Illumina-based transcriptome sequencing. For the normal group (GN), 44,600,142 high-quality clean reads were randomly assembled to produce 125,394 contigs. For the WSSV-challenged group (GW), 47,790,746 high-quality clean reads were randomly assembled to produce 148,983 contigs. After GO annotation, 39,482 unigenes were annotated into three ontologies: biological processes, cellular components, and molecular functions. In addition, 15,959 unigenes were mapped to 25 different COG categories. Moreover, 7,000 DEGs were screened out after a comparative analysis between the GN and GW samples, which were mapped into 250 KEGG pathways. Among these pathways, 36 were obviously changed (P-values < 0.05) and 28 pathways were extremely significantly changed (P-values < 0.01). Finally, five key DEGs involved in the JAK-STAT signaling pathway were chosen for qRT-PCR. The results showed that these five DEGs were obviously up-regulated at 36 h post WSSV infection in crayfish intestine. These results provide new insight into crayfish antiviral immunity mechanisms.