De novo transcriptome sequencing of Cryptotermes domesticus and comparative analysis of gene expression in response to different wood species

Exploration into Natural Variation Genes Associated with Determinate and Capitulum-like Inflorescence in Brassica napus

Brassica napus is a globally important vegetable and oil crop. The research is meaningful for the yield and plant architecture of B. napus. In this study, one natural mutant line with determinate and capitulum-like inflorescence was chosen for further study. Genetic analysis indicated that the segregation patterns of inflorescences in the F2 populations supported a digenic inheritance model, which was further approved via the BSA-Seq technique. The BSA-Seq method detected two QTL regions on C02 (14.27-18.41 Mb) and C06 (32.98-33.68 Mb) for the genetic control of determinate inflorescences in MT plants. In addition, the expression profile in MT compared with WT was analyzed, and a total of 133 candidate genes for regulating the flower development (75 genes, 56.4%), shoot meristem development (29 genes, 21.8%), and inflorescence meristem development (13 genes, 9.8%) were identified. Then one joint analysis combing BSA-Seq and RNA-Seq identified two candidate genes of BnaTFL1 and BnaAP1 for regulating the MT phenotype. Furthermore, the potential utilization of the MT plants was also discussed.

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.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Exploring the basis of 2-propenyl and 3-butenyl glucosinolate synthesis by QTL mapping and RNA-sequencing in Brassica juncea

Brassica juncea is used as a condiment, as vegetables and as an oilseed crop, especially in semiarid areas. In the present study, we constructed a genetic map using one recombinant inbred line (RIL) of B. juncea. A total of 304 ILP (intron length polymorphism) markers were mapped to 18 linkage groups designated LG01-LG18 in B. juncea. The constructed map covered a total genetic length of 1671.13 cM with an average marker interval of 5.50 cM. The QTLs for 2-propenyl glucosinolates (GSLs) colocalized with the QTLs for 3-butenyl GSLs between At1g26180 and BnapPIP1580 on LG08 in the field experiments of 2016 and 2017. These QTLs accounted for an average of 42.3% and 42.6% phenotypic variation for 2-propenyl and 3-butenyl GSLs, respectively. Furthermore, the Illumina RNA-sequencing technique was used to excavate the genes responsible for the synthesis of GSLs in the siliques of the parental lines of the RIL mapping population, because the bulk of the seed GSLs might originate from the siliques. Comparative analysis and annotation by gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) revealed that 324 genes were involved in GSL metabolism, among which only 24 transcripts were differentially expressed genes (DEGs). Among those DEGs, 15 genes were involved in the biosynthesis and transport of aliphatic GSLs, and their expression patterns were further validated by qRT-PCR analysis. Joint QTL mapping and RNA-sequencing analyses reveal one candidate gene of IIL1 (LOC106416451) for GSL metabolism in B. juncea. These results will be helpful for further fine mapping, gene cloning and genetic mechanisms of 2-propenyl and 3-butenyl GSLs in B. juncea.

Phylogeny and evolution of the cholesterol transporter NPC1 in insects

Sterols are essential nutrients for eukaryotes. Insects are obligate sterol auxotrophs and must acquire this key nutrient from their diets. The digestive tract is the organ for absorbing nutrients as well as sterols from food. In mice, the Niemann-Pick type C1 Like 1 (NPC1L1) gene is highly expressed in the intestine and is critical for cholesterol absorption. In contrast, the molecular mechanisms for the absorption of dietary sterols in insects have not been well studied. We annotated NPC1 genes in 39 insects from 10 orders using available genomic and transcriptomic information and inferred phylogenetic relationships. Insect NPC1 genes were grouped into two sister-clades, NPC1a and NPC1b, suggesting a likely duplication in the ancestor of insects. The former exhibited weaker gut-biased expression or a complete lack of tissue-biased expression, depending on the species, while the latter was highly enriched in the gut of three lepidopteran species. This result is similar to previous findings in Drosophila melanogaster. In insects, NPC1a accumulated non-synonymous substitutions at a lower rate than NPC1b. This pattern was consistent across orders, indicating that NPC1a evolved under stronger molecular constraint than NPC1b.

Optimization of a metatranscriptomic approach to study the lignocellulolytic potential of the higher termite gut microbiome

Background

Thanks to specific adaptations developed over millions of years, the efficiency of lignin, cellulose and hemicellulose decomposition of higher termite symbiotic system exceeds that of many other lignocellulose utilizing environments. Especially, the examination of its symbiotic microbes should reveal interesting carbohydrate-active enzymes, which are of primary interest for the industry. Previous metatranscriptomic reports (high-throughput mRNA sequencing) highlight the high representation and overexpression of cellulose and hemicelluloses degrading genes in the termite hindgut digestomes, indicating the potential of this technology in search for new enzymes. Nevertheless, several factors associated with the material sampling and library preparation steps make the metatranscriptomic studies of termite gut prokaryotic symbionts challenging.

Methods

In this study, we first examined the influence of the sampling strategy, including the whole termite gut and luminal fluid, on the diversity and the metatranscriptomic profiles of the higher termite gut symbiotic bacteria. Secondly, we evaluated different commercially available kits combined in two library preparative pipelines for the best bacterial mRNA enrichment strategy.

Results

We showed that the sampling strategy did not significantly impact the generated results, both in terms of the representation of the microbes and their transcriptomic profiles. Nevertheless collecting luminal fluid reduces the co-amplification of unwanted RNA species of host origin. Furthermore, for the four studied higher termite species, the library preparative pipeline employing Ribo-Zero Gold rRNA Removal Kit “Epidemiology” in combination with Poly(A) Purist MAG kit resulted in a more efficient rRNA and poly-A-mRNAdepletion (up to 98.44% rRNA removed) than the pipeline utilizing MICROBExpress and MICROBEnrich kits. High correlation of both Ribo-Zero and MICROBExpresse depleted gene expression profiles with total non-depleted RNA-seq data has been shown for all studied samples, indicating no systematic skewing of the studied pipelines.

Conclusions

We have extensively evaluated the impact of the sampling strategy and library preparation steps on the metatranscriptomic profiles of the higher termite gut symbiotic bacteria. The presented methodological approach has great potential to enhance metatranscriptomic studies of the higher termite intestinal flora and to unravel novel carbohydrate-active enzymes.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4076-9) contains supplementary material, which is available to authorized users.