Genome Survey Sequencing of the Mole Cricket Gryllotalpa orientalis.

The mole cricket Gryllotalpa orientalis is an evolutionarily, medicinal, and agriculturally significant insect that inhabits underground environments and is distributed globally. This study measured genome size by flow cytometry and k-mer based on low-coverage sequencing, and nuclear repetitive elements were also identified. The haploid genome size estimate is 3.14 Gb by flow cytometry, 3.17 Gb, and 3.77 Gb-based two k-mer methods, respectively, which is well within the range previously reported for other species of the suborder Ensifera. 56% of repetitive elements were found in G. orientalis, similar to 56.83% in Locusta migratoria. However, the great size of repetitive sequences could not be annotated to specific repeat element families. For the repetitive elements that were annotated, Class I-LINE retrotransposon elements were the most common families and more abundant than satellite and Class I-LTR. These results based on the newly developed genome survey could be used in the taxonomic study and whole genome sequencing to improve the understanding of the biology of G. orientalis.

De Novo Assembly and Characterization of the Transcriptome of an Omnivorous Camel Cricket (Tachycines meditationis).

Tachycines meditationis (Orthoptera: Rhaphidophoridae: Tachycines) is a widely distributed insect in eastern Asia. This species is common in urban environments, and its unique omnivorous diet may contribute to its success in various habitats. However, molecular studies on the species are scarce. Here, we obtained the first transcriptome sequence of T. meditationis and performed preliminary analyses to test whether the evolution of coding sequences fits the expectations based on the species' ecology. We retrieved 476,495 effective transcripts and annotated 46,593 coding sequences (CDS). We analysed the codon usage and found that directional mutation pressure was the leading cause of codon usage bias in this species. This genome-wide relaxed codon usage pattern in T. meditationis is surprising, given the potentially large population size of this species. Moreover, despite the omnivorous diet, the chemosensory genes of this species do not exhibit codon usage deviating significantly from the genome-level pattern. They also do not seem to experience more gene family expansion than other cave cricket species do. A thorough search for rapidly evolved genes using the dN/dS value showed that genes associated with substance synthesis and metabolic pathways, such as retinol metabolism, aminoacyl-tRNA biosynthesis, and fatty acid metabolism, underwent species-specific positive selection. While some results seem to contradict the species ecology, our transcriptome assembly provides a valuable molecular resource for future studies on camel cricket evolution and molecular genetics for feeding ecology in insects, in general.

The Chromosome-Level Genome of Hestina assimilis (Lepidoptera: Nymphalidae) Reveals the Evolution of Saprophagy-Related Genes in Brush-Footed Butterflies.

Most butterflies feed on nectar, while some saprophagous butterflies forage on various non-nectar foods. To date, little is known about the genomic and molecular shifts associated with the evolution of the saprophagous feeding strategy. Here, we assembled the high-quality chromosome-level genome of Hestina assimilis to explore its saprophagous molecular and genetic mechanisms. This chromosome-level genome of H. assimilis is 412.82 Mb, with a scaffold N50 of 15.70 Mb. In total, 98.11% of contigs were anchored to 30 chromosomes. Compared with H. assimilis and other Nymphalidae butterflies, the genes of metabolism and detoxification experienced expansions. We annotated 80 cytochrome P450 (CYP) genes in the H. assimilis genome, among which genes belonging to the CYP4 subfamily were significantly expanded (p < 0.01). These P450 genes were unevenly distributed and mainly concentrated on chromosomes 6-9. We identified 33 olfactory receptor (OR), 20 odorant-binding protein (OBP), and six gustatory receptor (GR) genes in the H. assimilis genome, which were fewer than in the nectarivorous Danaus plexippus. A decreased number of OBP, OR, and GR genes implied that H. assimilis should resort less to olfaction and gustation than their nectarivorous counterparts, which need highly specialized olfactory and gustatory functions. Moreover, we found one site under positive selection occurred in residue 996 (phenylalanine) of GR genes exclusive to H. assimilis, which is conservative in most lineages. Our study provides support for the adaptive evolution of feeding habits in butterflies.

Evolution of Gene Arrangements in the Mitogenomes of Ensifera and Characterization of the Complete Mitogenome of Schizodactylus jimo.

Gene arrangement (relative location of genes) is another evolutionary marker of the mitogenome that can provide extensive information on the evolutionary mechanism. To explore the evolution of gene arrangements in the mitogenome of diversified Ensifera, we sequenced the mitogenome of the unique dune cricket species found in China and used it for phylogenetic analysis, in combination with 84 known Ensiferan mitogenomes. The mitogenome of Schizodactylus jimo is a 16,428-bp circular molecule that contains 37 genes. We identified eight types of gene arrangement in the 85 ensiferan mitogenomes. The gene location changes (i.e., gene translocation and duplication) were in three gene blocks: I-Q-M-ND2, rrnl-rns-V, and ND3-A-R-N-S-E-F. From the phylogenetic tree, we found that Schizodactylus jimo and most other species share a typical and ancient gene arrangement type (Type I), while Grylloidea has two types (Types II and III), and the other five types are rare and scattered in the phylogenetic tree. We deduced that the tandem replication-random loss model is the evolutionary mechanism of gene arrangements in Ensifera. Selection pressure analysis revealed that purifying selection dominated the evolution of the ensiferan mitochondrial genome. This study suggests that most gene rearrangements in the ensiferan mitogenome are rare accidental events.

Bone marrow mesenchymal stem cells alleviate the formation of pathological scars in rats.

Introduction: Although bone marrow-derived mesenchymal stem cells (BMSCs) have attracted increasing attention because of their pivotal functions in the process of wound healing and fibrosis alleviation, the underlying molecular mechanisms have been poorly understood. Moreover, transforming growth factor beta 1 (TGF-ß1) is positively correlated with scar formation, whereas TGF-ß3 inhibits the pathological scar formation process. However, the relation of TGF-ß1, TGF-ß3, and the TGF-ß/Smad signaling pathway with BMSCs is unknown and requires further investigation. Methods: A cell co-culture platform was used to examine the relationship between BMSCs and dermal fibroblasts (DFs). EdU labelling and cell cycle detection were carried out to examine the viability of DF cells. Transwell and wound healing assays were used to test the cell migration of DFs. The expression of TGF-ß pathway components and collagens were determined by RT-qPCR and western blotting. A damaged skin rat model was applied to test the effects of BMSC treatment on skin wound healing. Results: The results showed that BMSC secretion could inhibit the viability and migration of DFs. Moreover, we observed that the TGF-ß-induced expression of TGF-ß1, Smad2, Smad3, COLI and COLIII was attenuated upon BMSC treatment in DFs, while the decrease in TGF-ß3 expression was enhanced by BMSCs. Furthermore, BMSC treatment accelerated wound healing and attenuated skin collagen deposition in a damaged skin rat model, leading to the mitigation of cell proliferation and enhancement of cell apoptosis. In addition, the expression of alpha-smooth muscle actin (α-SMA), COLI, and COLII was alleviated by BMSC treatment. Conclusions: Our results indicate that BMSCs can promote wound healing and inhibit skin collagen deposition, which is associated with the TGF-ß/Smad signaling pathway.

The Comparison of Gut Bacteria Communities and the Functions Among the Sympatric Grasshopper Species From the Loess Plateau.

Gut bacteria exert effects on the health and fitness of their insect hosts. Grasshoppers are an important part of the grassland ecosystem and provide important ecosystem services. As the most valuable feature in grassland ecosystem, the compositions and potential influences of gut bacterial in herbivorous grasshoppers in the same ecological environment are essential but undetermined. To facilitate such studies, we collected nine species of grasshoppers (n = 110) from a rebuild grassland on the Loess Plateau in northern Shaanxi, China, which is a representative area of ecosystem restoration model. We characterized the composition and function of the gut bacteria. We found that 326 OTUs were exhibited in all grasshoppers in which Enterobacter, Pantoea, Bacillus, and Spiroplsma are dominant. Among them, 18 OTUs were shared across all nine species of grasshoppers. The predicted function showed that the majority function of those OTUs were involved in survival dependent processes including membrane transport, carbohydrate metabolism, amino acid metabolism, and DNA replication and repair. The composition of gut bacteria is specific to each grasshopper species, and the bacteria community is most various in Trilophidia annulata. These results highlight the gut bacterial community diversity in different grasshopper species. Our findings are necessary for better understanding the relationships between this important herbivorous insect and their microbiomes and have the potential contribution of evaluating the revegetation and ecosystem management in this area.

Description of three new species of Gomphomastacinae (Orthoptera: Eumastacoidea) from the east of Qinghai-Tibet plateau, China.

The grasshopper Subfamily Gomphomastacinae (Orthoptera: Acridoidea: Eumastacidae) is widely distributed in the Qinghai-Tibetan Plateau and Central Asia. This paper describes three new species: Ptygomastax nihilsulcus Ge, sp. nov. Phytomastax pentaspinula Ge, sp. nov. and Pentaspinula unispinula Ge, sp. nov. of the subfamily from the eastern Qinghai-Tibetan Plateau, along with an updated identification key of species of the related genera. Detail figures of the body and genitalia of the new species are also provided.

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.

Description of Heteracris lecoqi Yetchom-Fondjo amp; Kekeunou sp. nov. and Heteracris hannai Wandji amp; Kekeunou sp. nov. and redescription of Heteracris guineensis (Krauss, 1890) (Orthoptera: Acrididae) with comments on its ecology in the Southern part of Cameroon.

The specimens studied were collected with sweep net and pitfall in the forests, agro-forests, herbaceous fallows, and crop fields of 14 localities in the southern part of Cameroon, from August 2015 to February 2018. The results show that Heteracris hannai sp. nov. and Heteracris lecoqi sp. nov. differ from previously known species of the same genus by the pattern of coloration and the details or shape of the phallic complex. H. lecoqi sp. nov. is characterized by light brownish body; male cercus with flattened, downcurved and obtuse apex; lophus strongly curved; interlophal space with V shape; apodeme of cingulum bow, convergent, with V-shape; valve of cingulum in lateral view longer than apical valve of penis; ramus in ventral view joined. H. hannai sp. nov. is distinct to other species by brown to grey body; male cercus with apex rounded, curved inside; lophus curved; interlophal space with U-shape; apodeme of cingulum thick, slightly parallel, with U shape; valve of cingulum in lateral view hardly longer than apical valve of penis; ramus in ventral view opened. The characteristics of H. guineensis are closer to H. hannai sp. nov. than to H. lecoqi sp. nov. H. lecoqi sp. nov. was collected only in the fallows while H. hannai sp. nov. was collected in the forests, agro-forests, fallows, and crop fields. Compared to H. guineensis, both new species are scarce in the natural vegetation and their distribution area is limited to two and four localities respectively for H. lecoqi sp. nov. and H. hannai sp. nov. All these three grasshopper's species were recorded as accidental species in all types of vegetation.

A synoptic review of Gomphomastacinae (Orthoptera: Eumastacoidea) from the Qinghai-Tibet plateau with a description of new species.

With unique climate, topography, and vegetation, the Qinghai-Tibet plateau is a special biogeographic region with richness of endemic species. However, the taxonomy and distribution of many insect groups on this plateau are still poorly known. Here, we synthesized a species checklist of subfamily Gomphomastacinae (Orthoptera: Acridoidea: Eumastacidae) for this region, while describing a new species: Myrmeleomastax wideis Qiu, sp. nov.. Images of the new species and the distribution map of all the known species are provided. Type specimens for the new species are deposited in the Zoological and Botanical Museum, Shaanxi Normal University, Xi'an, China (SNNU).

Comparison of the Complete Eragrostis pilosa Chloroplast Genome with Its Relatives in Eragrostideae (Chloridoideae; Poaceae).

Eragrostis of the tribe Eragrostideae is a taxonomically complex genus, because of its polyploid nature and the presence of similar morphological characters among its species. However, the relationship between these morphologically indistinguishable species at the genomic level has not yet been investigated. Here, we report the complete chloroplast genome of E. pilosa and compare its genome structures, gene contents, simple sequence repeats (SSRs), sequence divergence, codon usage bias, and Kimura 2-parameter (K2P) interspecific genetic distances with those of other Eragrostideae species. The E. pilosa chloroplast genome was 134,815 bp in length and contained 132 genes and four regions, including a large single-copy region (80,100 bp), a small single-copy region (12,661 bp), and a pair of inverted repeats (21,027 bp). The average nucleotide diversity between E. pilosa and E. tef was estimated to be 0.011, and 0.01689 among all species. The minimum and maximum K2P interspecific genetic distance values were identified in psaA (0.007) and matK (0.029), respectively. Of 45 SSRs, eight were shared with E. tef, all of which were in the LSC region. Phylogenetic analysis resolved the monophyly of the sampled Eragrostis species and confirmed the close relationship between E. pilosa and E. tef. This study provides useful chlorophyll genomic information for further species identification and phylogenetic reconstruction of Eragrostis species.

Cadmium-induced genome-wide DNA methylation changes in growth and oxidative metabolism in Drosophila melanogaster.

BACKGROUND: Cadmium (Cd)-containing chemicals can cause serious damage to biological systems. In animals and plants, Cd exposure can lead to metabolic disorders or death. However, for the most part the effects of Cd on specific biological processes are not known. DNA methylation is an important mechanism for the regulation of gene expression. In this study we examined the effects of Cd exposure on global DNA methylation in a living organism by whole-genome bisulfite sequencing (WGBS) using Drosophila melanogaster as model. RESULTS: A total of 71 differentially methylated regions and 63 differentially methylated genes (DMGs) were identified by WGBS. A total of 39 genes were demethylated in the Cd treatment group but not in the control group, whereas 24 showed increased methylation in the former relative to the latter. In most cases, demethylation activated gene expression: genes such as Cdc42 and Mekk1 were upregulated as a result of demethylation. There were 37 DMGs that overlapped with differentially expressed genes from the digital expression library including baz, Act5C, and ss, which are associated with development, reproduction, and energy metabolism. CONCLUSIONS: DNA methylation actively regulates the physiological response to heavy metal stress in Drosophila in part via activation of apoptosis.

Transcriptomics and differential gene expression in Whitmania pigra (Annelida: Clitellata: Hirudinida: Hirudinidae): Contrasting feeding and fasting modes.

The medicinal utility of leeches has been demonstrated through decades of use in modern hospital settings, mainly as relievers of venous congestion following flap or digit replantation surgery. In the present study, we sequence and annotate (through BLAST- and Gene Ontology-based approaches) the salivary transcriptome of the nonblood feeding hirudinid Whitmania pigra and assess the differential gene expression of anticoagulation factors (through both quantitative real-time PCR [qRT-PCR] and in silico-based methods) during feeding and fasting conditions. This was done in order to evince the diversity of putative anticoagulation factors, as well as estimate the levels of upregulation of genes immediately after feeding. In total, we found sequences with demonstrated orthology (via both phylogenetic analyses and BLAST-based approaches) to seven different proteins that have previously been linked to anticoagulatory capabilities-eglin C, bdellin, granulin, guamerin, hyaluronidase, destabilase I, and lipocalin. All of these were recovered from leeches both in the fasting and in the feeding conditions, but all show signs of upregulation in the feeding leeches. Interestingly, our RNA-seq effort, coupled with a hypergeometric test, indicated that the differentially expressed genes were disproportionately involved in three main immunological pathways (endocytosis, peroxisome regulation, and lysosome regulation). The results and implications of the finding of anticoagulants in this nonblood feeding leech and the putative upregulation of anticoagulation factors after feeding are briefly discussed in an evolutionary context.

Different mitogenomic codon usage patterns between damselflies and dragonflies and nine complete mitogenomes for odonates.

Damselflies and dragonflies, of the order Odonata, have distinct body plans and predatory abilities. Knowledge of their various evolutionary histories will allow for an understanding of the genetic and phenotypic evolution of insects. Mitogenomes are suitable materials to elucidate this, but the mitogenome of only a few odonates have been annotated. Herein, we report the complete mitogenome of nine odonates, including seven dragonflies and two damselflies, and a comprehensive analysis of the codon usage in 31 Odonata mitogenomes with the aim to estimate their evolutionary characteristics. Overall, a weak codon bias exists among odonate mitogenomes, although this favours AT-ending codons. Damselflies have a weaker codon usage bias than dragonflies, and 37 codons have significantly different usages. Both directional mutation and purifying selection shape damselfly and dragonfly mitogenomes. Although inevitable, directional mutation bias plays a minor role, whereas purifying selection pressure is the dominant evolutionary force. A higher selection pressure is observed in dragonflies than in damselflies, but it mainly acts on codon usage patterns rather than amino acid translation. Our findings suggest that dragonflies might have more efficient mitochondrial gene expression levels than damselflies, producing more proteins that support their locomotion and predatory abilities.

Analysis of synonymous codon usage pattern of genes in unique non-blood-sucking leech Whitmania pigra.

Whitmania pigra is a unique, fluid-sucking ectoparasite and an anticoagulant medical leech. The codon usage bias (CUB) is the nonuniform usage of synonymous codons in which some codons are more preferred than others. Here, we performed a comprehensive analysis of CUB of genes in W. pigra, analyzing 140 780 transcripts, 59 553 unigenes, and 20 304 qualified coding sequences (CDSs) from the transcriptomic data of W. pigra. The effective number of codons values suggested that the CUB was low in these genes. We recognized profoundly favored codons in W. pigra that have a G/C-ending. Parity rule two-bias plots suggested that both mutation pressure and natural selection might have influenced the CUB. However, neutrality plots revealed that natural selection might have played a major role while mutation pressure might have played a minor role in shaping the CUB. We applied principal component analysis to relative synonymous codon usage values for divided CDSs based on GC content and codon-ending bases. Codon usage in W. pigra had a general inclination toward C-ending codons and natural selection rather than mutation pressure is the dominant force in the genetic evolution of W. pigra. To our knowledge, this is the first study to describe a complete codon usage analysis of W. pigra; this will increase the understanding of CUB and evolution in W. pigra. The analysis of codon usage patterns in W. pigra aids in understanding its evolution and genetic architecture.

The complete mitochondrial genome of a grasshopper endemic to the Qinghai-Tibet plateau, Uvaroviola multispinosa (Acrididae: Oedipodinae).

The grasshopper Uvaroviola multispinosa (Acrididae: Oedipodinae), is an agricultural pest to pasture and limitedly distributed in Qinghai-Tibet plateau of China. The complete mitochondrial genome of U. multispinosa is 15,620 bp long, which comprises of 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA (tRNA) genes and a putative non-coding control region (GenBank accession ID: MK829651). These genes are unequally distributed on different DNA chains, which 23 are located on the majority Chain. The nucleotide composition shows evidently bias (A, C, G, and T was 43.8, 14.8, 10.1, and 31.3%, respectively) with an overall AT content of 75.2%. All PCGs are initiated by ATN codons, among them ATG is the most preferred. Eleven PCGs use a common stop codon of TAA or TAG, whereas the remaining two were terminated with single T as incomplete stop codon. The phylogenetic relationships based on Bayes method showed that U. multispinosa is closely related to Compsorhipis davidiana, which is in accordance with its traditional morphological classification.

Draft Genome of a Blister Beetle Mylabris aulica.

Mylabris aulica is a widely distributed blister beetle of the Meloidae family. It has the ability to synthesize a potent defensive secretion that includes cantharidin, a toxic compound used to treat many major illnesses. However, owing to the lack of genetic studies on cantharidin biosynthesis in M. aulica, the commercial use of this species is less extensive than that of other blister beetle species in China. This study reports a draft assembly and possible genes and pathways related to cantharidin biosynthesis for the M. aulica blister beetle using nanopore sequencing data. The draft genome assembly size was 288.5 Mb with a 467.8 Kb N50, and a repeat content of 50.62%. An integrated gene finding pipeline performed for assembly obtained 16,500 protein coding genes. Benchmarking universal single-copy orthologs assessment showed that this gene set included 94.4% complete Insecta universal single-copy orthologs. Over 99% of these genes were assigned functional annotations in the gene ontology, Kyoto Encyclopedia of Genes and Genomes, or Genbank non-redundant databases. Comparative genomic analysis showed that the completeness and continuity of our assembly was better than those of Hycleus cichorii and Hycleus phaleratus blister beetle genomes. The analysis of homologous orthologous genes and inference from evolutionary history imply that the Mylabris and Hycleus genera are genetically close, have a similar genetic background, and have differentiated within one million years. This M. aulica genome assembly provides a valuable resource for future blister beetle studies and will contribute to cantharidin biosynthesis.

Draft Genome of the Asian Buffalo Leech Hirudinaria manillensis.

The Asian Buffalo leech, Hirudinaria manillensis, is an aquatic sanguivorous species distributed widely in Southeast Asia. H. manillensis has long been used clinically for bloodletting and other medical purposes. Recent studies have focused on artificial culturing, strain optimization, and the identification and development new drugs based on the anticoagulant effects of H. manillensis bites; however, data regarding its genome remain unclear. This study aimed to determine the genome sequence of an adult Asian Buffalo leech. We generated a draft assembly of 151.8 Mb and a N50 scaffold of 2.28 Mb. Predictions indicated that the assembled genome contained 21,005 protein-coding genes. Up to 17,865 genes were annotated in multiple databases including Gene Ontology. Sixteen anticoagulant proteins with a Hirudin or Antistasin domain were identified. This study is the first to report the whole-genome sequence of the Asian Buffalo leech, an important sanguivorous leech of clinical significance. The quality of the assembly is comparable to those of other annelids. These data will help further the current understanding of the biological mechanisms and genetic characteristics of leeches and serve as a valuable resource for future studies.

A Checklist of Short-horned Grasshopper Species (Orthoptera: Caelifera) from Littoral Region of Cameroon with description of a new species of the genus Hemierianthus Saussure, 1903 (Orthoptera: Chorotypidae).

The Littoral Region of Cameroon is one of the most highly disturbed regions in Cameroon that remain poorly studied and where no investigation regarding the grasshopper fauna has been conducted. Hence a survey was conducted from August 2015 to November 2018 in 11 localities of this region. In each locality, grasshoppers were sampled in three vegetation types (forests, fallows and cultivated farms) using sweep nets sampling method. We collected in total 51 short-horned grasshopper species belonging to the families Acrididae, Chorotypidae, Euschmidtiidae, Pyrgomorphidae and Thericleidae. Among these species, 43 belong to Acrididae, 5 to Pyrgomorphidae and only one belongs to the Chorotypidae, Euschmidtiidae and Thericleidae respectively. The Chorotypidae species Hemierianthus mbongueensis Yetchom Xu sp.nov. is also described in this study. Twelve grasshopper species were widely distributed as they occurred in all study sites and 11 were recorded exclusively from one study site. This study showed that there is high species richness of short-horned grasshoppers in the Littoral Region of Cameroon. However, more sampling effort in different ecosystems of this region and other regions of Cameroon will certainly improve our understanding on grasshopper taxonomy, ecology and biogeography.

The Complete Chloroplast Genomes of Two Lespedeza Species: Insights into Codon Usage Bias, RNA Editing Sites, and Phylogenetic Relationships in Desmodieae (Fabaceae: Papilionoideae).

The genus Lespedeza (tribe: Desmodieae) consists of about 40 species that have high medicinal and economic value. However, in this genus, using morphological characters, the species identification is quite complicated, which can be solved by the analysis of the complete chloroplast genomes. As primary organelle genomes, the complete genome sequences of chloroplasts (cp) provide unique molecular information to study the divergence of species, RNA editing, and phylogeny. Therefore, to the best of our knowledge, for the first time, we sequenced the complete cp genomes of two representative Lespedeza species: Lespedeza davurica and Lespedeza cuneata. The cp genomes of both the species were found to be 149,010 bp in length, exhibiting the typical angiosperm chloroplast structure containing four regions. The Lespedeza cp genomes showed similar conserved gene contents, order, and orientations with a total GC content of 35.0%. A total of 128 genes, including 83 protein-coding genes, 37 tRNAs, and eight rRNAs, were identified from each genome. Unique molecular features of the two Lespedeza cp genome sequences were obtained by performing the analysis of repeats, sequence divergence, codon usage, and predicting the RNA editing sites in addition to phylogenetic analysis with other key genera in tribe Desmodieae. Using the two datasets, the phylogenetic relationship of Lespedeza species among Deasmodieae was discovered, suggesting that whole cp genomes provided useful information for phylogenetic studies of these species.