Differential Gene Expression Analysis of the Epacromius coerulipes (Orthoptera: Acrididae) Transcriptome

Transcriptome of the pygmy grasshopper Formosatettix qinlingensis (Orthoptera: Tetrigidae)

Formosatettix qinlingensis (Zheng, 1982) is a tiny grasshopper endemic to Qinling in China. For further study of its transcriptomic features, we obtained RNA-Seq data by Illumina HiSeq X Ten sequencing platform. Firstly, transcriptomic analysis showed that transcriptome read numbers of two female and one male samples were 25,043,314, 24,429,905, and 25,034,457, respectively. We assembled 65,977 unigenes, their average length was 1,072.09 bp, and the length of N50 was 2,031 bp. The average lengths of F. qinlingensis female and male unigenes were 911.30 bp, and 941.82 bp, and the N50 lengths were 1,745 bp and 1,735 bp, respectively. Eight databases were used to annotate the functions of unigenes, and 23,268 functional unigenes were obtained. Besides, we also studied the body color, immunity and insecticide resistance of F. qinlingensis. Thirty-nine pigment-related genes were annotated. Some immunity genes and signaling pathways were found, such as JAK-STAT and Toll-LIKE receptor signaling pathways. There are also some insecticide resistance genes and signal pathways, like nAChR, GST and DDT. Further, some of these genes were differentially expressed in female and male samples, including pigment, immunity and insecticide resistance. The transcriptomic study of F. qinlingensis will provide data reference for gene prediction and molecular expression study of other Tetrigidae species in the future. Differential genetic screening of males and females provides a basis for studying sex and immune balance in insects.

The susceptibility of Lymantria dispar larvae to Beauveria bassiana under Cd stress: A multi-omics study

Insect susceptibility to entomopathogenic microorganisms under heavy metal stress, as well as its regulatory mechanism is still poorly understood. This study aims to investigate the susceptibility of gypsy moth larvae to Beauveria bassiana under cadmium (Cd) stress (at 3.248 or 44.473 mg Cd/kg fresh food), and reveal the potential molecular mechanisms underlying the Cd effect on the larval susceptibility to B. bassiana via combined transcriptome and proteome analyses. Our results showed that pre-exposure to Cd increased the susceptibility of gypsy moth larvae to B. bassiana, and there was an additive effect between Cd exposure and B. bassiana infection on the larval mortality. Under the Cd stress at low and high concentrations, 138 and 899 differentially expressed genes (DEGs), as well as 514 and 840 differentially expressed proteins (DEPs) were identified, respectively. Immunotoxic effects induced by Cd exposure at the transcription level increased in a negative dose-response manner, with no immunity-related DEGs obtained at the low Cd concentration and a high number of immunity-related DEGs down-regulated at the high Cd concentration. In contrast, a potentially suppressed or stimulated trend in the Toll and Imd signaling pathway at protein level was revealed under low or high concentration of Cd treatment. Analysis of xenobiotics biodegradation-related pathways at both transcription and translation levels revealed that the gypsy moth larvae possessed an efficient homeostasis regulatory mechanism to the low-level Cd exposure, but exhibited a reduced xenobiotics biodegradation capability to the Cd stress at high levels. Together, these findings demonstrate Cd contamination promote the microbial-based biocontrol efficacy, and unravel the molecular regulatory network of heavy metal exposures that affects susceptibility of insects to pathogenic diseases.

Characterization and analysis of full-length transcriptomes from two grasshoppers, Gomphocerus licenti and Mongolotettix japonicus

Acrididae are diverse in size, body shape, behavior, ecology and life history; widely distributed; easy to collect; and important to agriculture. They represent promising model candidates for functional genomics, but their extremely large genomes have hindered this research; establishing a reference transcriptome for a species is the primary means of obtaining genetic information. Here, two Acrididae species, Gomphocerus licenti and Mongolotettix japonicus, were selected for full-length (FL) PacBio transcriptome sequencing. For G. licenti and M. japonicus, respectively, 590,112 and 566,165 circular consensus sequences (CCS) were generated, which identified 458,131 and 428,979 full-length nonchimeric (FLNC) reads. After isoform-level clustering, next-generation sequencing (NGS) short sequences were used for error correction, and remove redundant sequences with CD-HIT, 17,970 and 16,766 unigenes were generated for G. licenti and M. japonicus. In addition, we obtained 17,495 and 16,373 coding sequences, 1,082 and 813 transcription factors, 11,840 and 10,814 simple sequence repeats, and 905 and 706 long noncoding RNAs by analyzing the transcriptomes of G. licenti and M. japonicus, respectively, and 15,803 and 14,846 unigenes were annotated in eight functional databases. This is the first study to sequence FL transcriptomes of G. licenti and M. japonicus, providing valuable genetic resources for further functional genomics research.

Role of long non-coding RNA in DEET- and fipronil-mediated alteration of transcripts associated with Phase I and Phase II xenobiotic metabolism in human primary hepatocytes

Human exposure to environmental chemicals both individually and in combination occurs frequently world-wide most often with unknown consequences. Use of molecular approaches to aide in the assessment of risk involved in chemical exposure is a growing field in toxicology. In this study, we examined the impact of two environmental chemicals used in and around homes, the insect repellent DEET (N,N-diethyl-m-toluamide) and the phenylpyrazole insecticide fipronil (fluocyanobenpyrazole) on transcript levels of enzymes potentially involved in xenobiotic metabolism and on long non-coding RNAs (lncRNAs). Primary human hepatocytes were treated with these two chemicals both individually and in combination. Using RNA-Seq, we found that 10 major enzyme categories involved in phase 1 and phase 2 xenobiotic metabolism were significantly (α = 0.05) up- and down-regulated (i.e., 100 μM DEET-19 transcripts, 89% up and 11% down; 10 μM fipronil-52 transcripts, 53% up and 47% down; and 100 μM DEET +10 μM fipronil-69 transcripts, 43% up and 57% down). The altered genes were then mapped to the human genome and their proximity (within 1,000,000 bp) to lncRNAs examined. Unique proximities were discovered between altered lncRNA and altered P450s (CYP) and other enzymes (DEET, 2 CYP; Fipronil, 6 CYP and 15 other; and DEET + fipronil, 7 CYP and 21 other). Many of the altered P450 transcripts were in multiple clusters in the genome with proximal altered lncRNAs, suggesting a regulator function for the lncRNA. At the gene level there was high percent identity for lncRNAs near P450 clusters, but this relationship was not found at the transcript level. The role of these altered lncRNAs associated with xenobiotic induction, human diseases and chemical mixtures is discussed.

De Novo Assembly and Characterization of the Xenocatantops brachycerus Transcriptome

Grasshoppers are common pests but also have high nutritional and commercial potential. Xenocatantops brachycerus Willemse (Orthoptera: Acrididae) is an economically important grasshopper species that is reared in China. Using the IlluminaHiSeqTM 4000 platform, three transcriptomes of the adult male, adult female, and nymph of X. brachycerus were sequenced. A total of 133,194,848 clean reads were obtained and de novo assembled into 43,187 unigenes with an average length of 964 bp (N50 of 1799 bp); of these, 24,717 (57.23%) unigenes matched known proteins. Based on these annotations, many putative transcripts related to X. brachycerus growth, development, environmental adaptability, and metabolism of nutritional components and bioactive components were identified. In addition, the expression profiles of all three transcriptome datasets were analyzed, and many differentially expressed genes were detected using RSEM and PossionDis. Unigenes. Unigenes with functions associated with growth and development exhibited higher transcript levels at the nymph stage, and unigenes associated with environmental adaptability showed increased transcription in the adults. These comprehensive X. brachycerus transcriptomic data will provide a useful molecular resource for gene prediction, molecular marker development, and studies on signaling pathways in this species and will serve as a reference for the efficient use of other grasshoppers.

Comparative transcriptomic analysis of the evolution and development of flower size in Saltugilia (Polemoniaceae)

BACKGROUND

Flower size varies dramatically across angiosperms, representing innovations over the course of >130 million years of evolution and contributing substantially to relationships with pollinators. However, the genetic underpinning of flower size is not well understood. Saltugilia (Polemoniaceae) provides an excellent non-model system for extending the genetic study of flower size to interspecific differences that coincide with variation in pollinators.

RESULTS

Using targeted gene capture methods, we infer phylogenetic relationships among all members of Saltugilia to provide a framework for investigating the genetic control of flower size differences via RNA-Seq de novo assembly. Nuclear concatenation and species tree inference methods provide congruent topologies. The inferred evolutionary trajectory of flower size is from small flowers to larger flowers. We identified 4 to 10,368 transcripts that are differentially expressed during flower development, with many unigenes associated with cell wall modification and components of the auxin and gibberellin pathways.

CONCLUSIONS

Saltugilia is an excellent model for investigating covarying floral and pollinator evolution. Four candidate genes from model systems (BIG BROTHER, BIG PETAL, GASA, and LONGIFOLIA) show differential expression during development of flowers in Saltugilia, and four other genes (FLOWERING-PROMOTING FACTOR 1, PECTINESTERASE, POLYGALACTURONASE, and SUCROSE SYNTHASE) fit into hypothesized organ size pathways. Together, these gene sets provide a strong foundation for future functional studies to determine their roles in specifying interspecific differences in flower size.