Transcriptome assembly and identification of genes and SNPs associated with growth traits in largemouth bass (Micropterus salmoides)

Transcriptome analysis of large yellow croaker (Larimichthys crocea) at different growth rates

The unsynchronized growth of the large yellow croaker (Larimichthys crocea), which impacts growth efficiency, poses a challenge for aquaculture practitioners. In our study, juvenile stocks of large yellow croaker were sorted by size after being cultured in offshore cages for 4 months. Subsequently, individuals from both the fast-growing (FG) and slow-growing (SG) groups were sampled for analysis. High-throughput RNA-Seq was employed to identify genes and pathways that are differentially expressed during varying growth rates, which could suggest potential physiological mechanisms that influence growth rate. Our transcriptome analysis identified 382 differentially expressed genes (DEGs), comprising 145 upregulated and 237 downregulated genes in comparison to the SG group. GO and KEGG enrichment analyses indicated that these DEGs are predominantly involved in signal transduction and biochemical metabolic pathways. Quantitative PCR (qPCR) results demonstrated that cat, fasn, idh1, pgd, fgf19, igf2, and fads2 exhibited higher expression levels, whereas gadd45b and gadd45g showed lower expression compared to the slow-growing group. In conclusion, the differential growth rates of large yellow croaker are intricately associated with cellular proliferation, metabolic rates of the organism, and immune regulation. These findings offer novel insights into the molecular mechanisms and regulatory aspects of growth in large yellow croaker and enhance our understanding of growth-related genes.

Comparative transcriptomic analyses of brain-liver-muscle in channel catfish (Ictalurus punctatus) with differential growth rate

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

Muscle Fiber Characteristics and Transcriptome Analysis in Slow- and Fast-Growing Megalobrama amblycephala

Growth is an important trait in aquaculture that is influenced by various factors, among which genetic regulation plays a crucial role. Megalobrama amblycephala, one of the most important freshwater species in China, exhibits wide variations in body mass among individuals of the same age within the same pool. But the molecular mechanisms underlying wide variation in body mass remain unclear. Here, we performed muscle histological and transcriptome analysis of muscle tissues from Fast-Growing (FG) and Slow-Growing (SG) M. amblycephala at the age of 4 months old (4 mo) and 10 months old (10 mo) to elucidate its muscle development and growth mechanism. The muscle histological analysis showed smaller diameter and higher total number of muscle fibers in FG compared to SG at 4 mo, while larger diameter and total number of muscle fibers were detected in FG at 10 mo. The transcriptome analysis of muscle tissue detected 1171 differentially expressed genes (DEGs) between FG and SG at 4 mo, and 718 DEGs between FG and SG at 10 mo. Furthermore, 44 DEGs were consistently up-regulated in FG at both 4 mo and 10 mo. Up-regulated DEGs in FG at 4 mo were mainly enriched in the pathways related to cell proliferation, while down-regulated DEGs were significantly enriched in cell fusion and muscle contraction. Up-regulated DEGs in FG at 10 mo were mainly enriched in the pathways related to cell proliferation and protein synthesis. Therefore, these results provide novel insights into the molecular mechanism of M. amblycephala muscle growth at different stages, and will be of great guiding significance to promote the fast growth of M. amblycephala.

Single nucleotide polymorphisms in the MRFs gene family associated with growth in Nile tilapia

BACKGROUND

Muscle occupies most of the fish body, promoting the proliferation of fish muscle fibers can facilitate rapid growth and increase the body weight of fish. Some studiesSeveral previous suggest that Myogenic regulatory factors (MRFs) play an important role in the growth of fish.

OBJECTIVE

To investigate the association between the polymorphism of MRFs gene family and growth traits in Nile tilapia (Oreochromis niloticus), get more molecular markers for growth.

METHODS

Amplified the Nile tilapia MRFs family gene, including Myogenic determination 1 (Myod1), Myogenic determination 2 (Myod2), Myogenin (Myog), Myogenic factor 5 (Myf5), and Myogenic factor 6 (Myf6), single nucleotide polymorphism (SNP) were screened by Sanger sequencing.

RESULTS

A total of 16 SNP loci were screened, including six for Myf5, six for Myf6, one for Myog, one for Myod1 and two for Myod2. The growth traits were analyzed in relation to these 16 SNP loci, and the results indicated significant associations between all 16 SNP loci and the growth traits (P < 0.05). The linkage disequilibrium analysis revealed that D1 and D2 diplotypes of Myf5 gene, E1, E2, E3 and E4 of Myf6 gene, and F1 diplotype of Myod2 gene were significantly associated with superior growth traits.

CONCLUSION

There were 6, 6, 1, 1 and 2 growth-related molecular markers in Myf5, Myf6, Myog, Myod1 and Myod2 genes, respectively, which could be applied to the breeding of Nile tilapia.

Single nucleotide polymorphism SNP19140160 A > C is a potential breeding locus for fast-growth largemouth bass (Micropterus salmoides)

BACKGROUND

Largemouth bass (Micropterus salmoides) has significant economic value as a high-yielding fish species in China's freshwater aquaculture industry. Determining the major genes related to growth traits and identifying molecular markers associated with these traits serve as the foundation for breeding strategies involving gene pyramiding. In this study, we screened restriction-site associated DNA sequencing (RAD-seq) data to identify single nucleotide polymorphism (SNP) loci potentially associated with extreme growth differences between fast-growth and slow-growth groups in the F1 generation of a largemouth bass population.

RESULTS

We subsequently identified associations between these loci and specific candidate genes related to four key growth traits (body weight, body length, body height, and body thickness) based on SNP genotyping. In total, 4,196,486 high-quality SNPs were distributed across 23 chromosomes. Using a population-specific genotype frequency threshold of 0.7, we identified 30 potential SNPs associated with growth traits. Among the 30 SNPs, SNP19140160, SNP9639603, SNP9639605, and SNP23355498 showed significant associations; three of them (SNP9639603, SNP9639605, and SNP23355498) were significantly associated with one trait, body length, in the F1 generation, and one (SNP19140160) was significantly linked with four traits (body weight, height, length, and thickness) in the F1 generation. The markers SNP19140160 and SNP23355498 were located near two growth candidate genes, fam174b and ppip5k1b, respectively, and these candidate genes were closely linked with growth, development, and feeding. The average body weight of the group with four dominant genotypes at these SNP loci in the F1 generation population (703.86 g) was 19.63% higher than that of the group without dominant genotypes at these loci (588.36 g).

CONCLUSIONS

Thus, these four markers could be used to construct a population with dominant genotypes at loci related to fast growth. These findings demonstrate how markers can be used to identify genes related to fast growth, and will be useful for molecular marker-assisted selection in the breeding of high-quality largemouth bass.

Muscle transcriptome analysis provides new insights into the growth gap between fast- and slow-growing Sinocyclocheilus grahami

Sinocyclocheilus grahami is an economically valuable and famous fish in Yunnan Province, China. However, given its slow growth (40 g/2 years) and large growth differences among individuals, its growth performance needs to be improved for sustainable future use, in which molecular breeding technology can play an important role. In the current study, we conducted muscle transcriptomic analysis to investigate the growth gaps among individuals and the mechanism underlying growth within 14 fast- and 14 slow-growth S. grahami. In total, 1,647 differentially expressed genes (DEGs) were obtained, including 947 up-regulated and 700 down-regulated DEGs in fast-growth group. Most DEGs were significantly enriched in ECM-receptor interaction, starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, amino acids biosynthesis and metabolism, peroxisome, and PPAR signaling pathway. Some genes related to glycogen degradation, glucose transport, and glycolysis (e.g., adipoq, prkag1, slc2a1, agl, pygm, pgm1, pfkm, gapdh, aldoa, pgk1, pgam2, bpgm, and eno3) were up-regulated, while some genes related to fatty acid degradation and transport (e.g., acox1, acaa1, fabp1b.1, slc27a1, and slc27a2) and amino acid metabolism (e.g., agxt, shmt1, glula, and cth) were down-regulated in the fast-growth group. Weighted gene co-expression network analysis identified col1a1, col1a2, col5a1, col6a2, col10a1, col26a1, bglap, and krt15 as crucial genes for S. grahami growth. Several genes related to bone and muscle growth (e.g., bmp2, bmp3, tgfb1, tgfb2, gdf10, and myog) were also up-regulated in the fast-growth group. These results suggest that fast-growth fish may uptake adequate energy (e.g., glucose, fatty acid, and amino acids) from fodder, with excess energy substances used to synthesize collagen to accelerate bone and muscle growth after normal life activities are maintained. Moreover, energy uptake may be the root cause, while collagen synthesis may be the direct reason for the growth gap between fast- and slow-growth fish. Hence, improving food intake and collagen synthesis may be crucial for accelerating S. grahami growth, and further research is required to fully understand and confirm these associations.

Identification of Shewanella putrefaciens as a novel pathogen of the largemouth bass (Micropterus salmoides) and histopathological analysis of diseased fish

The largemouth bass (Micropterus salmoides) is an economically important aquaculture species in China, and its production has increased rapidly in recent years. Although Shewanella putrefaciens is known to infect several fish species, its role in infecting M. salmoides is relatively unknown. Here, we isolated a gram-negative bacterial strain (termed XX2021) from farmed largemouth bass. Based on the results of 16S rRNA sequencing and phylogenetic analyses, the isolate was identified as S. putrefaciens. The virulence of XX2021 was dependent on water temperature, such as the LD50 values were 4.21×104, 7.26×105, and 2.47×106 CFU/g fish weight at 10°C, 18°C, and 25°C, respectively. Four virulent genes—including dksA, hem, lonR, and fur—were screened through a PCR assay. The results of an antibiotic resistance test showed that XX2021 was sensitive to kanamycin, cefotaxime, doxycycline, sulfamethoxazole, florfenicol, tetracycline, and gentamicin; showed intermediate susceptibility to streptomycin, ampicillin, and norfloxacin; and was resistant to nalidixic acid and penicillin. XX2021-infected fish showed clinical symptoms typical of S. putrefaciens infection. In addition, we re-isolated XX2021 from infected fish and confirmed its identity using 16S rRNA sequencing. Histopathological changes were observed in the intestine, head kidney, spleen, and liver of diseased fish. This study presents the first report of the pathogenic effects of S. putrefaciens in farmed largemouth bass. Our findings may help develop effective disease control strategies for aquaculture fish and prevent disease outbreaks under low water temperatures.

ATP synthase subunit e is a shrimp growth-associated breeding marker

Penaeus vannamei is one of the most popular aquaculture species in the world. This species is featured with its fast-growing and delicious taste, which drives people develop various strains. During this process identification of trait-associated markers could effectively increase breeding efficiency. Driven by this, we tried to screen fast-growing key regulators via a FACS-based high throughput method, in which 2-NBDG was applied as a fluorescent indicator for direct glucose uptake measurement. Totally six candidate genes were screened out followed by in vitro validation in 293T cells. After that, the correlation between these genes and shrimp growing was further verified in a hybrid lineage. The expression level of two genes including ATP synthase subunit e and inhibitor of apoptosis protein showed some correlation with shrimp growth speed. Furthermore, we tested these two candidate markers in various lineages and confirmed that ATP synthase subunit e could be a shrimp growth-associated breeding marker.

Development of EST-Molecular Markers from RNA Sequencing for Genetic Management and Identification of Growth Traits in Potato Grouper (Epinephelus tukula)

Simple Summary

The potato grouper is a novel aquaculture species in Taiwan. Due to the lack of genetic information concerning this species, we have developed molecular markers based on transcriptome sequencing and further characterized their association with gene diversity and growth traits of this species. Ultimately, these markers could be utilized as accurate and efficient tools for genetic management and marker-assisted selection of potato grouper with distinct growth traits.

Abstract

The accuracy and efficiency of marker-assisted selection (MAS) has been proven for economically critical aquaculture species. The potato grouper (Epinephelus tukula), a novel cultured grouper species in Taiwan, shows large potential in aquaculture because of its fast growth rate among other groupers. Because of the lack of genetic information for the potato grouper, the first transcriptome and expressed sequence tag (EST)-derived simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were developed. Initially, the transcriptome was obtained from seven cDNA libraries by using the Illumina platform. De novo transcriptome of the potato grouper yielded 51.34 Gb and 111,490 unigenes. The EST-derived SSR and SNP markers were applied in genetic management, in parentage analysis, and to discover the functional markers of economic traits. The F₁ juveniles were identified as siblings from one pair of parents (80 broodstocks). Fast- and slow-growth individuals were analyzed using functional molecular markers and through their association with growth performance. The results revealed that two SNPs were correlated with growth traits. The transcriptome database obtained in this study and its derived SSR and SNP markers may be applied not only for MAS but also to maintain functional gene diversity in the novel cultured grouper.

Verticillium dahliae VdBre1 is required for cotton infection by modulating lipid metabolism and secondary metabolites

The soil-borne ascomycete Verticillium dahliae causes wilt disease in more than two hundred dicotyledonous plants including the economically important crop cotton, and results in a severe reduction in cotton fiber yield and quality. During infection, V. dahliae secretes numerous secondary metabolites, which act as toxic factors to promote the infection process. However, the mechanism underlying how V. dahliae secondary metabolites regulate cotton infection remains largely unexplored. In this study, we report that VdBre1, an ubiquitin ligase (E3) enzyme to modify H2B, regulates radial growth and conidia production of V. dahliae. The VdBre1 deletion strains show nonpathogenic symptoms on cotton, and microscopic inspection and penetration assay indicated that penetration ability of the ∆VdBre1 strain was dramatically reduced. RNA-seq revealed that a total of 1643 differentially expressed genes between the ∆VdBre1 strain and the wild type strain V592, among which genes related to lipid metabolism were significantly overrepresented. Remarkably, the volume of lipid droplets in the ∆VdBre1 conidia was shown to be smaller than that of wild-type strains. Further metabolomics analysis revealed that the pathways of lipid metabolism and secondary metabolites, such as steroid biosynthesis and metabolism of terpenoids and polyketides, have dramatically changed in the ∆VdBre1 metabolome. Taken together, these results indicate that VdBre1 plays crucial roles in cotton infection and pathogenecity, by globally regulating lipid metabolism and secondary metabolism of V. dahliae.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

Identification of Candidate Growth-Related SNPs and Genes Using GWAS in Brown-Marbled Grouper (Epinephelus fuscoguttatus)

Brown-marbled grouper, Epinephelus fuscoguttatus, is not only an important commercial fish species, but also an important crossbreeding parent in grouper industry. Improvement of growth traits of this species contributes to the development of grouper breeding. Currently, the development of molecular marker associated with growth of brown-marbled grouper is rare. Thus, we performed the first genome-wide association study (GWAS) for five growth traits in 172 brown-marbled groupers with 43,688 SNPs detected by ddRAD-seq. We identified a total of 5 significant and 18 suggestive QTLs located in multiple chromosomes associated with growth traits. In the 20 kb window of the significant SNPs and suggestive SNPs, 5 and 14 potential candidate genes affecting growth were detected, respectively. Five potential candidate genes near the significantly associated SNPs were selected for expression analysis. Among of which, bmp2k, wasf1, and acyp2 involved in bone development, maintenance of mitochondrion structure, and metabolism were differentially expressed. Interestingly, the SNP 23:29601315 located in the intron of bmp2k was significantly associated with body weight, body length, body height, and body thickness and suggestively associated with total length. We verified the locus using another new group including 123 individuals. The results showed that individuals with CC genotype have better growth traits comparing other individuals. Our findings not only contribute to understanding the molecular mechanism of growth regulation, but also promote the advance of marker-assisted selection in brown-marbled grouper.

A novel transcriptome-derived SNPs array for tench (Tinca tinca L.)

Tench (Tinca tinca L.) has great economic potential due to its high rate of fecundity and long-life span. Population genetic studies based on allozymes, microsatellites, PCR-RFLP and sequence analysis of genes and DNA fragments have revealed the presence of Eastern and Western phylogroups. However, the lack of genomic resources for this species has complicated the development of genetic markers. In this study, the tench transcriptome and genome were sequenced by high-throughput sequencing. A total of 60,414 putative SNPs were identified in the tench transcriptome using a computational pipeline. A set of 96 SNPs was selected for validation and a total of 92 SNPs was validated, resulting in the highest conversion and validation rate for a non-model species obtained to date (95.83%). The validated SNPs were used to genotype 140 individuals belonging to two tench breeds (Tabor and Hungarian), showing low (FST = 0.0450) but significant (<0.0001) genetic differentiation between the two tench breeds. This implies that set of validated SNPs array can be used to distinguish the tench breeds and that it might be useful for studying a range of associations between DNA sequence and traits of importance. These genomic resources created for the tench will provide insight into population genetics, conservation fish stock management, and aquaculture.

The transcriptome sequencing and functional analysis of eyestalk ganglions in Chinese mitten crab (Eriocheir sinensis) treated with different photoperiods

Photoperiod plays an important role in individual growth, development, and metabolism in crustaceans. The growth and reproduction of crabs are closely related to the photoperiod. However, as of yet, there are still no transcriptomic reports of eyestalk ganglions treated under different photoperiods in the Chinese mitten crab (Eriocheir sinensis), which is a benthonic crab with high commercial value in Asia. In this study, we collected the eyestalk ganglions of crabs that were reared under different photoperiods, including a control group (L: D = 12 h: 12 h, named CC), a constant light group (L: D = 24 h: 0 h, named LL) and a constant darkness group (L: D = 0 h: 24 h, named DD). RNA sequencing was performed on these tissues in order to examine the effects of different photoperiods. The total numbers of clean reads from the CC, LL and DD groups were 48,772,584 bp, 53,943,281 bp and 53,815,178 bp, respectively. After de novo assembly, 161,380 unigenes were obtained and were matched with different databases. The DEGs were significantly enriched in phototransduction and energy metabolism pathways. Results from RT-qPCR showed that TRP channel protein (TRP) in the phototransduction pathway had a significantly higher level of expression in LL and DD groups than in the CC group. We found that the downregulation of the pyruvate dehydrogenase complex (PDC) gene and the upregulation phosphoenolpyruvate carboxykinase (PPC) gene were involved in energy metabolism processes in LL or DD. In addition, we also found that the upregulation of the expression level of the genes Gαq, pyruvate kinase (PK), NADH peroxidase (NADH) and ATPase is involved in phototransduction and energy metabolism. These results may shed some light on the molecular mechanism underlying the effect of photoperiod in physiological activity of E. sinensis.

Analysis of the transcriptome data in Litopenaeus vannamei reveals the immune basis and predicts the hub regulation-genes in response to high-pH stress

Soil salinization erodes the farmlands and poses a serious threat to human life, reuse of the saline-alkali lands as cultivated resources becomes increasingly prominent. Pacific white shrimp (Litopenaeus vannamei) is an important farmed aquatic species for the development and utilization of the saline-alkali areas. However, little is known about the adaptation mechanism of this species in terms of high-pH stress. In the present study, a transcriptome analysis on the gill tissues of L. vannamei in response to high-pH stress (pH 9.3 ± 0.1) was conducted. After analyzing, the cyclic nucleotide gated channel-Ca2+ (CNGC-Ca2+) and patched 1 (Ptc1) were detected as the majority annotated components in the cAMP signaling pathway (KO04024), indicating that the CNGC-Ca2+ and Ptc1 might be the candidate components for transducing and maintaining the high-pH stress signals, respectively. The immunoglobulin superfamily (IgSF), heat shock protein (HSP), glutathione s-transferase (GST), prophenoloxidase/phenoloxidase (proPO/PO), superoxide dismutase (SOD), anti-lipopolysaccharide factor (ALF) and lipoprotein were discovered as the major transcribed immune factors in response to high-pH stress. To further detect hub regulation-genes, protein-protein interaction (PPI) networks were constructed; the genes/proteins "Polymerase (RNA) II (DNA directed) polypeptide A" (POLR2A), "Histone acetyltransferase p300" (EP300) and "Heat shock 70kDa protein 8" (HSPA8) were suggested as the top three hub regulation-genes in response to acute high-pH stress; the genes/proteins "Heat shock 70kDa protein 4" (HSPA4), "FBJ murine osteosarcoma viral oncogene homolog" (FOS) and "Nucleoporin 54kDa" (NUP54) were proposed as the top three hub regulation-genes involved in adapting endurance high-pH stress; the protein-interactions of "EP300-HSPA8" and "HSPA4-NUP54" were detected as the most important biological interactions in response to the high-pH stress; and the HSP70 family genes might play essential roles in the adaptation of the high-pH stress environment in L. vannamei. These findings provide the first insight into the molecular and immune basis of L. vannamei in terms of high-pH environments, and the construction of a PPI network might improve our understanding in revealing the hub regulation-genes in response to abiotic stress in shrimp species and might be beneficial for further studies.

Gene co-expression network analysis reveals coordinated regulation of three characteristic secondary biosynthetic pathways in tea plant (Camellia sinensis)

BACKGROUND

The leaves of tea plants (Camellia sinensis) are used to produce tea, which is one of the most popular beverages consumed worldwide. The nutritional value and health benefits of tea are mainly related to three abundant characteristic metabolites; catechins, theanine and caffeine. Weighted gene co-expression network analysis (WGCNA) is a powerful system for investigating correlations between genes, identifying modules among highly correlated genes, and relating modules to phenotypic traits based on gene expression profiling. Currently, relatively little is known about the regulatory mechanisms and correlations between these three secondary metabolic pathways at the omics level in tea.

RESULTS

In this study, levels of the three secondary metabolites in ten different tissues of tea plants were determined, 87,319 high-quality unigenes were assembled, and 55,607 differentially expressed genes (DEGs) were identified by pairwise comparison. The resultant co-expression network included 35 co-expression modules, of which 20 modules were significantly associated with the biosynthesis of catechins, theanine and caffeine. Furthermore, we identified several hub genes related to these three metabolic pathways, and analysed their regulatory relationships using RNA-Seq data. The results showed that these hub genes are regulated by genes involved in all three metabolic pathways, and they regulate the biosynthesis of all three metabolites. It is notable that light was identified as an important regulator for the biosynthesis of catechins.

CONCLUSION

Our integrated omics-level WGCNA analysis provides novel insights into the potential regulatory mechanisms of catechins, theanine and caffeine metabolism, and the identified hub genes provide an important reference for further research on the molecular biology of tea plants.

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.