Transcriptome assembly for a colour-polymorphic grasshopper (Gomphocerus sibiricus) with a very large genome size

Characterizing Hox genes in mayflies (Ephemeroptera), with Hexagenia limbata as a new mayfly model

BACKGROUND

Hox genes are key regulators of appendage development in the insect body plan. The body plan of mayfly (Ephemeroptera) nymphs differs due to the presence of abdominal appendages called gills. Despite mayflies' phylogenetic position in Paleoptera and novel morphology amongst insects, little is known of their developmental genetics, such as the appendage-regulating Hox genes. To address this issue we present an annotated, early instar transcriptome and embryonic expression profiles for Antennapedia, Ultrabithorax, and Abdominal A proteins in the mayfly Hexagenia limbata, identify putative Hox protein sequences in the mayflies H. limbata, Cloeon dipterum, and Ephemera danica, and describe the genomic organization of the Hox gene cluster in E. danica.

RESULTS

Transcriptomic sequencing of early instar H. limbata nymphs yielded a high-quality assembly of 83,795 contigs, of which 22,975 were annotated against Folsomia candida, Nilaparvata lugens, Zootermopsis nevadensis and UniRef90 protein databases. Homeodomain protein phylogeny and peptide annotations identified coding sequences for eight of the ten canonical Hox genes (excluding zerknüllt/Hox3 and fushi tarazu) in H. limbata and C. dipterum, and all ten in E. danica. Mayfly Hox protein sequences and embryonic expression patterns of Antp, Ubx, and Abd-A appear highly conserved with those seen in other non-holometabolan insects. Similarly, the genomic organization of the Hox cluster in E. danica resembles that seen in most insects.

CONCLUSIONS

We present evidence that mayfly Hox peptide sequences and the embryonic expression patterns for Antp, Ubx, and Abd-A are extensively conserved with other insects, as is organization of the mayfly Hox gene cluster. The protein data suggest mayfly Antp, Ubx, and Abd-A play appendage promoting and repressing roles during embryogenesis in the thorax and abdomen, respectively, as in other insects. The identified expression of eight Hox genes, including Ubx and abd-A, in early instar nymphs further indicates a post-embryonic role, possibly in gill development. These data provide a basis for H. limbata as a complementary Ephemeridae model to the growing repertoire of mayfly model species and molecular techniques.

Systematic analysis of the lysine malonylome in Sanghuangporus sanghuang

Background

Sanghuangporus sanghuang is a well-known traditional medicinal mushroom associated with mulberry. Despite the properties of this mushroom being known for many years, the regulatory mechanisms of bioactive compound biosynthesis in this medicinal mushroom are still unclear. Lysine malonylation is a posttranslational modification that has many critical functions in various aspects of cell metabolism. However, at present we do not know its role in S. sanghuang. In this study, a global investigation of the lysine malonylome in S. sanghuang was therefore carried out.

Results

In total, 714 malonyl modification sites were matched to 255 different proteins. The analysis indicated that malonyl modifications were involved in a wide range of cellular functions and displayed a distinct subcellular localization. Bioinformatics analysis indicated that malonylated proteins were engaged in different metabolic pathways, including glyoxylate and dicarboxylate metabolism, glycolysis/gluconeogenesis, and the tricarboxylic acid (TCA) cycle. Notably, a total of 26 enzymes related to triterpene and polysaccharide biosynthesis were found to be malonylated, indicating an indispensable role of lysine malonylation in bioactive compound biosynthesis in S. sanghuang.

Conclusions

These findings suggest that malonylation is associated with many metabolic pathways, particularly the metabolism of the bioactive compounds triterpene and polysaccharide. This paper represents the first comprehensive survey of malonylation in S. sanghuang and provides important data for further study on the physiological function of lysine malonylation in S. sanghuang and other medicinal mushrooms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08120-0.

Genome Size Estimation and Full-Length Transcriptome of Sphingonotus tsinlingensis: Genetic Background of a Drought-Adapted Grasshopper

Sphingonotus Fieber, 1852 (Orthoptera: Acrididae), is a grasshopper genus comprising approximately 170 species, all of which prefer dry environments such as deserts, steppes, and stony benchlands. In this study, we aimed to examine the adaptation of grasshopper species to arid environments. The genome size of Sphingonotus tsinlingensis was estimated using flow cytometry, and the first high-quality full-length transcriptome of this species was produced. The genome size of S. tsinlingensis is approximately 12.8 Gb. Based on 146.98 Gb of PacBio sequencing data, 221.47 Mb full-length transcripts were assembled. Among these, 88,693 non-redundant isoforms were identified with an N50 value of 2,726 bp, which was markedly longer than previous grasshopper transcriptome assemblies. In total, 48,502 protein-coding sequences were identified, and 37,569 were annotated using public gene function databases. Moreover, 36,488 simple tandem repeats, 12,765 long non-coding RNAs, and 414 transcription factors were identified. According to gene functions, 61 cytochrome P450 (CYP450) and 66 heat shock protein (HSP) genes, which may be associated with drought adaptation of S. tsinlingensis, were identified. We compared the transcriptomes of S. tsinlingensis and two other grasshopper species which were less tolerant to drought, namely Mongolotettix japonicus and Gomphocerus licenti. We observed the expression of CYP450 and HSP genes in S. tsinlingensis were higher. We produced the first full-length transcriptome of a Sphingonotus species that has an ultra-large genome. The assembly characteristics were better than those of all known grasshopper transcriptomes. This full-length transcriptome may thus be used to understand the genetic background and evolution of grasshoppers.

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.

The green-brown polymorphism of the club-legged grasshopper Gomphocerus sibiricus is heritable and appears genetically simple

Background

Local coexistence of distinct, genetically determined color morphs can be unstable and transitional. Stable, long-term coexistence requires some form of balancing selection to protect morphs from getting lost by directional selection or genetic drift. However, not all phenotypic polymorphism need to have a genetic basis. We here report on the genetic basis of two color polymorphisms in the club-legged grasshopper Gomphocerus sibiricus: a green-brown polymorphism that is phylogenetically and geographically widespread among orthopteran insects and a pied-brown pattern polymorphism that is shared among many gomphocerine grasshoppers.

Results

We found a remarkably clear outcome of matings within and between morph that suggest not only that the green-brown polymorphism is heritable in this species, but that results can be most parsimoniously explained by a single autosomal locus with two alleles in which the green allele is dominant over the brown allele. A few individuals did not match this pattern and suggest the existence of genetic modifiers and/or developmental phenocopies. We also show that the pied-brown polymorphism is highly heritable, although the evidence for the involvement of one or more loci is less clear-cut.

Conclusions

Overall, our data demonstrate that the two polymorphisms are heritable in the club-legged grasshopper and appear genetically simple, at least with respect to green morphs. The results are consistent with the idea that the synthesis or transport of a pigment involved in the production of green coloration (likely biliverdin) is lost by homozygosity for loss-of-function alleles in brown individuals. The apparently simple genetic architecture of the green-brown polymorphism offer potential for studying balancing selection in the field and for genetic mapping in this species.