Skip to main content

Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication

Abstract

Background

The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. The bamboo-feeding leafhopper genus Abrus belong to the tribe Athysanini of subfamily Deltocephalinae, which currently comprises 19 valid described species, and are limited to the Oriental and Palaearctic regions in China. Although the taxonomy of Abrus are well updated, the references on comparative mitogenomic analyses of Abrus species are only known for a single species. In this study, we sequenced and analyzed the complete mitochondrial genomes (mitogenomes) of Abrus daozhenensis Chen, Yang & Li, 2012 (16,391bp) and A. yunshanensis Chen, Yang & Li, 2012 (15,768bp) (Athysanini), and compared with published mitogenome sequence of A. expansivus Xing & Li, 2014 (15,904bp).

Results

These Abrus species shared highly conserved mitogenomes with similar gene order to that of the putative ancestral insect with 37 typical genes and a non-coding A + T-rich region. The nucleotide composition of these genomes is highly biased toward A + T nucleotides (76.2%, 76.3%, and 74.7%), AT-skews (0.091 to 0.095, and 0.095), negative GC-skews (− 0.138, − 0.161, and − 0.138), and codon usage. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. Phylogenetic analyses based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus) based on limited available mitogenome sequence data of 37 species.

Conclusion

At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered monophylitic Athysanini, except Watanabella graminea (Athysanini) in Opsiini with high support values.

Peer Review reports

Background

Insect mitochondrial genome is a small double-stranded circular molecule with remarkable conservation in size ranging from 14–20 kb that encodes 37 genes: 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs genes, and a non-coding A + T-rich region (or control region) [1,2,3]. The representative mitochondrial (mt) genomes of almost all insect orders from higher-level to lower taxonomic ranks have extensively been utilized for studying phylogeny, population genetics, comparative and evolutionary genomics, molecular evolution, and identification at various taxonomic levels [3,4,5]. Due to their high genome copy numbers, multiple genome-level characteristics, relatively high evolutionary rate, and greater phylogenetic informativeness than single mitochondrial genes, mitogenome sequences have been widely used in various phylogenomic studies [3, 5, 6].

The leafhopper genus Abrus belonging to the tribe Athysanini of subfamily Deltocephalinae, was originally described by Dai & Zhang [7] with six new species (type species: A. hengshanensis Dai & Zhang, 2002) from China. After that, Li & Wang [8], Dai & Zhang [9], Li et al. [10], Chen et al. [11], Yang & Chen [12], and Xing & Li [13] further added 13 new species from the Oriental and Palaearctic parts of China. So far, Abrus is only restricted to China with 19 valid described species, which are quite similar in body coloration and difficult to distinguish based on external morphology, but the male genitalia with a unique structure of aedeagus are considerably differences among these species [13]. Among these, 19 species are widely distributed in the Oriental Region (Guizhou, Sichuan, Hunan, Hubei, Guangxi, Guangdong, Fujian, and Zhejiang), and only A. coneus is also distributed in the Palaearctic Region (China: Gansu) [13]. It belongs to the tribe Athysanini, which is one of the largest and most diverse tribes of the subfamily Deltocephalinae, which includes 228 genera and 1181 described species worldwide, with a majority of species serve as a vector of economically important plant diseases [14, 15]. The generic and species composition in Athysanini is not constant at present and is continuously changed with time due to ongoing taxonomic and systematic revisions, which include the discovery of new genera or transfer of described genera to other tribes [14,15,16,17]. At present, there are approximately 20 genera and 70 species of Athysanini in China [10, 13, 16, 17].

The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess [14, 15]. Numerous phylogenetic studies consistently supported for a paraphyletic Athysanini with a large group of polyphyletic genera based on morphology alone [18] or various combinations of morphological and molecular datasets: two nuclear gene fragments (28S ribosomal DNA and Histone H3) [14, 19, 20], or different combinations of 13-37 mitochondrial genes [21, 22]. Despite being the largest tribe with a cosmopolitan distribution, so far, only the representative three genera (Abrus, Norvellina, and Watanabella) and five species of Athysanini have been sequenced and analyzed (including two novel sequences in this study) from China [23, 24].

In this study, we sequenced and analyzed the complete mitogenomes of Abrus daozhenensis and Abrus yunshanensis, and compared them with the published sequence of Abrus expansivus to examine the phylogenetic relationships among the newly sequenced and published athysanine species based on mtDNA genome data, using the concatenated amino acid sequences and nucleotide sequences from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs).

Results

General features and gene orders in Abrus mitogenomes

The complete mitochondrial genomes of A. daozhenensis (GenBank: MZ274046), A. yunshanensis (GenBank: MZ274047), and A. expansivus (GenBank: MK033020) are 16,391 bp, 15,768 bp, and 15,904 bp in size, respectively (Table 1). The circular genome maps of these species are presented in Fig. 1A–C. All mitogenomes contained a typical set of 37 mitochondrial genes (13 PCGs, 22 tRNA genes, and two rRNA genes) and one control region (Supplementary Tables S1–S3). Gene order was invariant and identical to Drosophila yakuba and to other published deltocephaline mitogenomes [23, 25, 26].

Table 1 List of the mitochondrial genomes analyzed in the present study. List of taxa included in this study
Fig. 1
figure 1

The complete mitochondrial genome maps of three Abrus species. A Abrus daozhenensis; B Abrus yunshanensis; C Abrus expansivus

Base composition

The complete mitochondrial genome of three Abrus species exhibited heavy AT nucleotide bias, with 76.2%, 76.3%, and 74.7% in A. daozhenensis, A. yunshanensis, and A. expansivus, respectively. The A+T content of the CR (mean value = 82.4%) was always significantly higher than in other regions, while PCGs showed the lowest A+T content values (mean value = 74.0%) (Table 2). All three species had higher A+T content in rrnL than rrnS, with significantly different mean values (80.8% and 77.9%) (Table 2). All three mitogenomes showed positive AT-skews (0.091 to 0.097) and negative GC-skews (−0.138 to −0.161). Similarly, the control region showed positive AT-skewed (0.011 to 0.041) and negative GC-skewed (−0.103 to −0.119) (Table 2).

Table 2 Comparative nucleotide compositions, AT- skews and GC-skews in three Abrus species

Gene overlaps, intergenic spacers and non-coding regions

Gene overlaps are present in these three mitogenomes which ranged from 1 bp to 8 bp in length: A. daozhenensis (11 gene junctions, 36 bp overlaps), A. yunshanensis (12 gene junctions, 36 bp overlaps), and A. expansivus (12 gene junctions, 29 bp overlaps). A. daozhenensis shares 11 gene overlaps, whereas, A. yunshanensis and A. expansivus, each with 12 gene overlaps: Ile-Gln (3/3/1 bp), nad2-Trp (8/8/8 bp), Asp-atp8 (7/7/7 bp), atp8-atp6 (-/-/1 bp), Gly-nad3 (2/2/1 bp), Ala-Arg (1/1/1 bp), Arg-Asn (1/1/1 bp), Asn-Ser (4/4/3 bp), Ser-Glu (-/-/1 bp), nad5-His (1/1/- bp), His-nad4 (7/7/4 bp), Pro-nad6 (1/1/1 bp), and cob-Ser (1/1/1 bp) (Supplementary Tables S1–S3).

Intergenic spacers are present in these three mitogenomes which ranged from 1 bp to 39 bp in length: A. daozhenensis (7 intergenic spacers), A. yunshanensis (8 intergenic spacers), and A. expansivus (8 intergenic spacers). The longest intergenic spacers were present between Cys-Tyr with 38, 39, and 18 in A. daozhenensis, A. yunshanensis, and A. expansivus, respectively (Supplementary Tables S1–S3).

The putative control region, or A+T rich region, located between rrnS and trnI, was the most variable region in the whole mitogenome. The full lengths of CR in three Abrus mitogenomes were 2,035, 1,947 and 1,545 bp, respectively.

Transfer RNA and ribosomal RNA genes

For the 22 typical animal tRNA genes in each Abrus mitogenome, 14 tRNAs were encoded by the J-strand and the remaining eight were located on the N-strand, ranging from 61 to 76 bp in length. All tRNAs could be folded into the canonical cloverleaf secondary structure except for trnS1 (AGN), which lacks the dihydrouridine (DHU) arm and instead forms a loop in A. daozhenensis, A. yunshanensis and A. expansivus (Supplementary Figs. S1–S3). In addition to the lack of the dihydrouridine (DHU) arm in trnS1 (AGN), A. expansivus has a trnG (GGN) that has failed to form a typical clover secondary structure due to the lack of the TѰC arm (Supplementary Fig. S3). Besides the classic A-U and C-G pairs in the secondary structure, there were 35, 36, and 26 G-U base pairings in A. daozhenensis, A. yunshanensis and A. expansivus, respectively. Some other mismatched base pairs (U-U, A-A, C-U, A-G and A-C) were also found in the acceptor arm and anticodon arm (Supplementary Figs. S1–S3). The large rRNA subunit was located at a conserved position between trnL1 (CUN) and trnV, while the small rRNA subunits was between trnV and the control region (Fig. 1A–C). The 16S-rRNA (rrnL) gene length are 1,204, 1,203, and 1,204 in A. daozhenensis, A. yunshanensis and A. expansivus, respectively. Whereas 12S-rRNA (rrnS) genes with average lengths of 745 bp, 743 bp, and 756 bp, respectively. The mean A+T contents of the two rRNA genes (rrnL and rrnS) in all three Abrus mitogenomes were 80.1% and 77.8%, respectively (Table 2), and rrn genes were encoded on the N-strand.

Protein-coding genes (PCGs), codon usage, and relative synonymous codon usage (RSCU)

The total length of 13 PCG of A. daozhenensis, A. yunshanensis, and A. expansivus are 10,945bp, 10,945bp, and 10,936bp, respectively. Among the 13 protein-coding genes, nine are located on J-strand while the remaining four are on N-strand (Fig. 1A–C). The first codon position had a significantly higher A+T content than the second and third positions (78.1 versus 71.1% and 72.9%) (Supplementary Table 4).

All 13 PCGs started with the standard ATN codons. The starting codons of 13 PCGs are the same in A. daozhenensis and A. yunshanensis, among which the starting codons of cox1, atp6, cox3, nad4 and cob genes are all ATG, the starting codons of nad2, nad3, nad5, nad4l, and nad6 genes are all ATT, and the starting codons of cox2, atp8, and nad1 genes are all ATA. Except for cox1 and nad6, which started by ATA codon, the remaining starting codons in A. expansivus are the same as in A. daozhenensis and A. yunshanensis. Twelve of the thirteen PCGs in Abrus species are terminated with a TAA or TAG codon, except cox2 gene terminates with an incomplete T residue. Except for the cox3 gene that terminated with TAA in A. daozhenensis and TAG in A. yunshanensis, the stop codons of other genes were the same. Among them, the cox1, atp8, atp6, cox3, nad4, nad4l, nad6, cob, nad1 genes are terminated by standard TAA codon, whereas TAG terminates nad2, nad3, and nad5 genes, and the cox2 gene is terminated with incomplete T as the termination codon. Except for the stop codons of nad3 and nad4l genes in A. expansivus, the remaining stop codons are consistent with A. daozhenensis and A. yunshanensis (Table 3).

Table 3 Comparison of length, start and stop codons of 13 protein-coding genes (PCGs) among Abrus daozhenensis, Abrus yunshanrnsis and Abrus expansivus

After excluding the termination codons, the relative synonymous codon usage (RSCU) was calculated and summarized for A. daozhenensis, A. yunshanensis and A. expansivus in Fig. 2. The total numbers of non-stop codons were 3,636, 3,636 and 3,625 in A. daozhenensis, A. yunshanensis and A. expansivus respectively. The most frequently used amino acids: Isoleucine (Ile) and Methionine (Met), Asparagine (Asn), and Lysine (Lys) were the most frequently used amino acids.

Fig. 2
figure 2

Relative synonymous codon usage (RSCU) of mitochondrial genomes of Abrus daozhenensis, Abrus yunshanensis and Abrus expansivus. Codons missing in mitogenomes are shown at the top of columns

Phylogenetic relationship

Previous phylogenetic studies based on morphological and molecular datasets have largely elucidated the relationships among the deltocephaline tribes [14, 18, 19, 26, 29]. The phylogenetic relationships based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered Fieberiellini as sister to the remaining tribes followed Penthimiini and a clad (Selenocephalinithe + Cicadulini), (Scaphoideini + Drabescini + Macrostelini + Drabescini), (Drabescini + Athysanini + Opsiini), (Mukariini + Macrostelini + Eupelicini + Chiasmini), and (Athysanini + Deltocephalini + Paralimnini), in both BI and ML analyses with high support values (Fig. 3, Supplementary Figs. S4-6, S10-12). However, similar combinations of 37 genes with amino acid sequences formed a different clade (Supplementary Figs. S7-9, S13-15). The paraphyly of the following five tribes, Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini has not yet been resolved in recent phylogenetic studies based on morphological and molecular datasets. In this study, we recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus, see Zhang & Xing [16] and Yang et al. [24] based on limited available mitogenome sequence data of 37 species. So far, only 14 tribes, 32 genera, and 37 species (including two novel sequences in this study) have been sequenced and analyzed from China and are available on NCBI (https://www.ncbi.nlm.nih.gov/). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered Balclutha sp. (Macrostelini) in tribe Drabescini and Watanabella graminea (Athysanini) in Opsiini with high support values (Figs. S4-15).

Fig. 3
figure 3

Phylogenetic tree of subfamily Deltocephalinae using the Bayesian inference (BI) analysis method based on the concatenated nucleotide sequences of PCG123 + 2 rRNA + 22 tRNA datasets. Numbers on each node correspond to the posterior probability (PP) values

Discussion

The goal of this study was to investigate the comparative mitochondrial genome analysis of three Abrus species and its phylogenetic position in the tribe Athysanini. The phylogenetic position of Athysanini is still poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. Previous studies based on morphology and molecular datasets consistently recovered a paraphyletic Athysanini [14, 15, 21,22,23]. Similarly, the paraphyly of the following tribes, Athysanini, Opsiini, Penthimiini, Selenocephalini, and Scaphoideini has not yet been resolved in recent phylogenetic studies [14, 15, 21,22,23]. In the present study, we recovered the above mentioned tribes as monophylitic, with the exception of Athysanini, in which Watanabella graminea was consistently recovered in Opsiini with high strong support values. Further studies must include additional representative species from the above representative tribes to reconstruct the phylogenetic relationships among the tribes in the subfamily Deltocephalinae.

The length of complete mitochondrial genome in insects are remarkable conserved, typically ranging from 14–20 kb in size that encodes 37 genes: 13 PCGs, 2 rRNAs, 22 tRNA, and a non-coding A + T-rich region [3]. With a few exceptions, the length of nearly all known mitochondrial genome sequenced of the subfamily Deltocephalinae ranged from 14-17 kb in size. However, there are few exceptions, in which one or two protein-coding genes are missing (Dryadomorpha sp.: KX437736; Hishimonus phycitis: KX437727) or with relatively short length of of non-coding A + T-rich region (Phlogotettix sp. 2: KX437721; Cicadulina mbila: MK251127; Psammotettix sp.: KX437742). Herein, we found the consistent genome size in Abrus species with previous observations of Deltocephalinae: A. daozhenensis (GenBank: MZ274046), A. yunshanensis (GenBank: MZ274047), and A. expansivus (GenBank: MK033020) are 16,391 bp, 15,768 bp, and 15,904 bp in size, respectively (Table 1). The variation in mitogenome size among the different deltocephaline species is mainly due to the variable number of repeats in the control region. The mitogenome sequences of Abrus species were highly conserved in gene content, gene size, gene order, base composition, codon usage of PCGs and tRNA secondary structures.

Comparative studies of leafhopper mitochondrial genomes suggest that genome rearrangements are very rare, and was proposed as a genus-level trait [32]. All Abrus species from this study, and the other deltocephaline species sequenced so far are arranged in the putative ancestral insect gene arrangement, tRNA cluster trnW-trnC-trnY. So far, only three leafhopper species, Macrosteles quadrilineatus [31] and Macrosteles quadrimaculatus [32] shares the same tRNA rearrangement, that tRNA cluster trnW-trnC-trnY is rearranged to trnC-trnW-trnY. Species Japananus hyalinus are rearranged tRNA cluster trnY-trnW-trnC [29]. Complete mitogenome sequences for additional Deltocephaline specimens are needed to facilitate broader comparison and to identify features of the potential novel rearrangements in the gene order of the insects. The mitogenome of leafhoppers is a typical conserved circular, about 14.5–17 kb in length which shows the conserved gene structures, it contains double-stranded DNA molecule, 37 typical mitochondrial genes, and does not have introns. Although conservative, the novel rearrangements of the mitogenomes allowed us to identify with enough variation for population genetic or phylogeographic studies.

Conclusions

The present phylogenetic analyses based on 12 phylogenetic trees based on the various combination of concatenated datasets of 13 PCGs, 2 rRNA genes, and 22 tRNA genes along with different partition datasets and models strongly supported for paraphyletic Athysanini in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). In all analyses, Fieberiellini is consistently recovered as sister to the remaining tribes, however, the internal tribal relationships by using the PCGs and amino acid sequences formed a different clade. The unresolve paraphyly of the following five tribes, Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini have recovered the monophyletic, except Athysanini with Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus, consistently recovered in Opsiini with high support values.

In this study, we newly sequenced and analyzed the complete mitochondrial genomes of two Abrus species and compared them with the published mitogenome of A. expansivus Xing & Li, 2014. Comparative mitogenome analyses showed that the gene content, gene order, gene size, nucleotide base composition, codon usage of PCGs, and tRNA secondary structures are highly conserved in Abrus species. The complete mitochondrial genome in Abrus species is around 10,936-10,945 bp. The full lengths of CR in Abrus species range from 2,035bp in A. daozhenensis to 15,45bp in A. expansivus. All PCGs in Abrus species were initiated with ATN codons, however, 12 of the 13 PCGs are terminated with a TAA or TAG codon, except cox2 gene terminates with an incomplete T residue. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN), which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs).

Materials and methods

Taxon sampling

The adult specimens of Abrus daozhenensis and A. yunshanensis were collected with a hand-net from the bamboo plant at the Kuankuoshui National Nature Reserve, Guiyang, Guizhou, China, from June to August 2020. The permission was taken from Kuankuoshui National Nature Reserve to collect these leaf hopper species. The field-collected specimens were initially preserved in 100% ethanol which was later identified based on morphological characters following Chen et al. [11] and Xing & Li [13] and stored at − 20 °C until DNA extraction. Vouchers were deposited at the Institute of Entomology, Guizhou University, Guiyang, China.

DNA extraction and sequencing

The entire body of A. daozhenensis and A. yunshanensis were sent to Guangzhou Ruike Gene Technology Co. (China) for genomic DNA extraction and sequencing. Total genomic DNA was extracted from the thoracic muscle tissues and legs of each individual using the DNeasy© Tissue Kit (Qiagen) according to the manufacturer's protocol. For Illumina sequencing, genomic DNA was isolated using TIANamp Genomic DNA Kit (Tiangen, China). The Illumina sequencing library was generated using Truseq Nano DNA HT Sample Prep Kit (Illumina, USA). The complete mitogenome was sequenced using high-throughput sequencing on the Illumina Novaseq 6000 platform with an average insert size of 350 bp and a paired-end 150 bp (PE 150) sequencing strategy to generate sequencing data not less than 2 GB. Raw reads were trimmed of adapters using Trimmomatic [41].

Sequence assembly, annotation and analyses

In order to investigate the phylogenetic relationships among the deltocephaline tribes, we retrieved the GenBank file of almost all available mitogenomic data of deltocephaline tribes (14 out of 39 currently recognized tribes), including 37 species in 32 genera as ingroup. Outgroup sampling comprised eight species in three subfamilies: Eurymelinae, Coelidiinae, and Iassinae (Table 1). Consensus sequences of A. daozhenensis and A. yunshanensis were assembled with the complete mitochondrial genome of A. expansivus (NC_045238) as a reference in GENEIOUS v. 10.2.3 (https://www.geneious.com/) [42]. Both the assembled sequences were aligned and compared with the published sequence of Deltocephalinae to extract the 13 PCG and 2 rRNA genes by using MUSCLE [43] in MEGA 7.0 [44]. The 22 transfer RNA (tRNA) genes were annotated using the MITOS web server (http://mitos.bioinf.unileipzig.de/index.py) [45]. The secondary structure of tRNA was obtained from the MITOS web server and manually edited in Adobe Photoshop CS 6.0 (Figs. S1–S3). The graphical map of the circular genome and annotation was made using the CGView Server (http://stothard.afns.ualberta.ca/cgview_server/) [46]. MEGA Version 7.0 was used to analyze the base composition of the complete mitogenome. The strand asymmetry was calculated by using the formulas: GC-skew = [(G – C)/(G + C)] and AT-skew = [(A – T)/(A + T)] [47].

Phylogenetic analysis

For phylogenetic analyses, the concatenated nucleotide and amino acid sequences of all 13 PCGs, 22 tRNA genes, and 2 rRNA genes were extracted using PhyloSuite (v1.2.3) [48] and aligned with the codon and Normal mode using the invertebrate mitochondrial code and automatic strategy in MAFFT [49]. The aligned PCGs were then refined using MACSE and trimmed by GBlock. The aligned RNA genes and amino acid sequences were trimmed by trimAl. Aligned genes were concatenated to generate six datasets using PhyloSuite: (1) PCG123 (nucleotide data: 1,067bp); (2) PCG123 + 2 rRNA (nucleotide data: 12,521bp); (3) PCG123 + 22 tRNA + 2 rRNA (nucleotide data: 13,904bp); (4) PCG_AA (amino acid sequences: 3,572bp); (5) PCG_AA + 2 rRNA (amino acid sequences: 5,422bp); and PCG_AA + 2 rRNA + 22 tRNA (amino acid sequences: 6,805bp). We used PartitionFinder v.2.1.1 [50] to determine the best partitioning schemes for the datasets under the Bayesian Information Criterion (BIC). Detailed information on the best partitioning schemes and substitution model selection for ML and BI analyses using PartionFinder2 and ModelFinder are summarized in Supplementary Tables S5-S16. The phylogenetic trees were reconstructed using two methods: Bayesian inference (BI) and Maximum Likelihood (ML) based on various combinations of nucleotide and amino acid sequences. BI analyses were conducted using MrBayes v.3.2.7a [51] implemented in PhyloSuite software with various data partition schemes and best-fitting models determined by PartitionFinder, and Maximum Likelihood (ML) was performed on the IQ-tree using ultrafast bootstrap with 5000 replicates as implemented on the website server (http://iqtree.cibiv.univie.ac.at). The BI analyses contains four simultaneous Markov chain Monte Carlo (MCMC) runs of 2 million generations, and sampled every 1000 generations. The initial 25% of the sampled data were discarded as burn-in. Other parameters were kept at default settings. The finalized trees were visualized and edited with FIGTREE v1.3.1 [52] and the Interactive Tree of Life (iTOL: https://itol.embl.de) version 5 [53].

Availability of data and materials

The mitogenome sequences of Abrus daozhenensis Chen, Yang & Li, 2012 (GenBank accession no: MZ274046), and Abrus expansivus Xing & Li, 2014 (GenBank: MK033020) and A. yunshanensis Chen, Yang & Li, 2012 (GenBank accession no: MZ274047) have been deposited on GenBank.

References

  1. Wolstenholme DR. Animal mitochondrial DNA: structure and evolution. Int Rev Cytol. 1992;141:173–216. https://doi.org/10.1016/s0074-7696(08)62066-5.

    Article  CAS  PubMed  Google Scholar 

  2. Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80. https://doi.org/10.1093/nar/27.8.1767.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Cameron SL. Insect mitochondrial genomics: implications for evolution and phylogeny. Annu Rev Entomol. 2014;59:95–117. https://doi.org/10.1146/annurev-ento-011613-162007.

    Article  CAS  PubMed  Google Scholar 

  4. Wilson AC, Cann RL, Carr SM, George M, Gyllensten UB, Helm-Bychowski KM, Higuchi RG, Palumbi SR, Prager EM, Sage RD, Stoneking M. Mitochondrial DNA and two perspectives on evolutionary genetics. Biol J Linn Soc. 1985;26(4):375–400.

    Article  Google Scholar 

  5. Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach AT. Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu Rev Ecol Evol Syst. 2006;37:545–79. https://doi.org/10.1146/annurev.ecolsys.37.091305.110018.

    Article  Google Scholar 

  6. Morgan B, Wang TY, Chen TZ, Moctezuma V, Burgos O, Le MH, Huang JP. Long-read Sequencing Data Reveals Dynamic Evolution of Mitochondrial Genome Size and the Phylogenetic Utility of Mitochondrial DNA in Hercules Beetles (Dynastes; Scarabaeidae). Genome Biol Evol. 2022;14(10):evac147. https://doi.org/10.1093/gbe/evac147.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Dai W, Zhang YL. A new genus and six new species of Deltocephalinae from China (Homoptera: Cicadellidae). Acta Zootaxonomica Sinica. 2002;27:313–22.

    Google Scholar 

  8. Li ZZ, Wang LM. Descriptions of two new species of the genus Abrus from China (Hemiptera, Cicadellidae). Acta Zootaxonomica Sinica. 2006;31:840–2.

    Google Scholar 

  9. Dai W, Zhang YL. A review of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with description of one new species. Zootaxa. 2008;1688:37–53. https://doi.org/10.5281/zenodo.180583.

    Article  Google Scholar 

  10. Li ZZ, Dai RH, Xing JC: Deltocephalinae from China (Hemiptera: Cicadellidae). Popular Science Press, Beijing. 336 pp; 2011.

    Google Scholar 

  11. Chen XS, Yang L, Li ZZ. Bamboo-feeding leafhoppers in China. Beijing: China Forestry Publishing House; 2012. p. 218.

    Google Scholar 

  12. Yang L, Chen XS. Two new species of the bamboo-feeding leafhopper genus Abrus Dai & Zhang (Hemiptera, Cicadellidae, Deltocephalinae) from China. ZooKeys. 2013;318:81–9. https://doi.org/10.3897/zookeys.318.5799[14].

    Article  Google Scholar 

  13. Xing JC, Li ZZ. Two new species of the genus Abrus Dai & Zhang, 2002 (Hemiptera, Cicadellidae, Deltocephalinae) from China. ZooKeys. 2014;419:103–9. https://doi.org/10.3897/zookeys.419.7481.

    Article  Google Scholar 

  14. Zahniser JN, Dietrich CH. A review of the tribes of Deltocephalinae (Hemiptera: Auchenorrhyncha: Cicadellidae). Eur J Taxon. 2013;45:1–211. https://doi.org/10.5852/ejt.2013.45.

    Article  Google Scholar 

  15. Cao Y, Dietrich CH, Zahniser JN, Dmitriev DA. Dense sampling of taxa and characters improves phylogenetic resolution among deltocephaline leafhoppers (Hemiptera: Cicadellidae: Deltocephalinae). Syst Entomol. 2022;47:430–44. https://doi.org/10.13012/B2IDB-8842653_V2.

    Article  Google Scholar 

  16. Zhang L, Xing JC. Description of a new leafhopper species of the genus Watanabella (Hemiptera, Cicadellidae, Deltocephalinae) from China. Zootaxa. 2020;4808:539–46. https://doi.org/10.11646/zootaxa.4808.3.8.

    Article  Google Scholar 

  17. Yao BB, Dietrich CH, Xing JC. Description of a new species of the leafhopper genus Tambocerus (Hemiptera: Cicadellidae: Deltocephalinae) from Guizhou, China. Zootaxa. 2022;5169(1):87–92. https://doi.org/10.11646/zootaxa.5169.1.8.

    Article  PubMed  Google Scholar 

  18. Zahniser JN, Dietrich CH. Phylogeny of the leafhopper subfamily Deltocephalinae (Auchenorrhyncha: Cicadellidae) and related subfamilies based on morphology. System Biodivers. 2008;6:1–24. https://doi.org/10.1017/S1477200007002617.

    Article  Google Scholar 

  19. Zahniser JN, Dietrich CH. Phylogeny of the leafhopper subfamily Deltocephalinae based on molecular and morphological data with a revised family-group classification. Syst Entomol. 2010;35:489–511. https://doi.org/10.1111/j.1365-3113.2010.00522.x.

    Article  Google Scholar 

  20. Sunil Meshram. NM, Hashmi TR, Shashank PR: A new genus Vittaliana belonging to the tribe Opsiini (Hemiptera: Cicadellidae) from India and its molecular phylogeny. PeerJ. 2020;8:e9515. https://doi.org/10.7717/peerj.9515.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Luo H, Wang Y, Di XC, Shan LCY, Wang SS. The complete mitochondrial genome of Fieberiella septentrionalis (Hemiptera: Cicadellidae: Deltocephalinae). Mitochondrial DNA Part B. 2021;6(4):1426–8. https://doi.org/10.1080/23802359.2021.1912671.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Wu K, Yan M, Zhang Y, Dietrich CH, Duan Y. The complete mitochondrial genome of Aconurella prolixa (Lethierry 1885) (Deltocephalinae: Chiasmini). Mitochondrial DNA Part B. 2022;7:81–3. https://doi.org/10.1080/23802359.2021.2008834.

    Article  PubMed  Google Scholar 

  23. Wang JJ, Xing JC. Complete mitochondrial genome of Abrus expansivus (Hemiptera: Cicadellidae: Deltocephalinae) from China. Mitochondrial DNA Part B. 2019;4(1):197–8. https://doi.org/10.1080/23802359.2018.1545548.

    Article  Google Scholar 

  24. Yang W, Gao Y, Li C, Song Y. The complete mitochondrial genome of Chlorotettix nigromaculatus (Deltocephalinae) with phylogenetic consideration. Mitochondrial DNA Part B. 2019;4:624–5. https://doi.org/10.1080/23802359.2018.1564389.

    Article  Google Scholar 

  25. Lu JK, Wang JJ, Li DF, Wang XY, Dai RH. Description of the whole mitochondrial genome of Bhatia longiradiata (Hemiptera, Cicadellidae, Deltocephalinae: Drabescini) and phylogenetic relationship. Genes Genom. 2023;45:59–70. https://doi.org/10.1007/s13258-022-01338-6.

    Article  CAS  Google Scholar 

  26. Du Y, Dai W, Dietrich CH. Mitochondrial genomic variation and phylogenetic relationships of three groups in the genus Scaphoideus (Hemiptera: Cicadellidae: Deltocephalinae). Sci Rep. 2017;7:16908.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Xu TL, Dai RH. Complete mitogenome of Reticuluma hamata (Hemiptera: Cicadellidae: Deltocephalinae) from China. Mitochondrial DNA Part B. 2020;5(2):1437–8. https://doi.org/10.1080/23802359.2020.1721351.

    Article  Google Scholar 

  28. Song N, Cai W, Li H. Deep-level phylogeny of Cicadomorpha inferred from mitochondrial genomes sequenced by NGS. Sci Rep. 2017;7:10429.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Du Y, Zhang C, Dietrich CH, Zhang Y, Dai W: Characterization of the complete mitochondrial genomes of Maiestas dorsalis and Japananus hyalinus (Hemiptera: Cicadellidae) and comparison with other Membracoidea. Sci Rep 2017a,14197. https://doi.org/10.1038/s41598-017-14703-3.

  30. Song N, Cai W, Li H. Insufficient power of mitogenomic data in resolving the auchenorrhynchan monophyly. Zool J Linn Soc. 2018;183:776–90.

    Article  Google Scholar 

  31. Mao M, Yang X, Bennett G. The complete mitochondrial genome of Macrosteles quadrilineatus (Hemiptera: Cicadellidae). Mitochondrial DNA Part B. 2017;2:173–5.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Du Y, Dietrich CH, Dai W. Complete mitochondrial genome of Macrosteles quadrimaculatus (Matsumura) (Hemiptera: Cicadellidae: Deltocephalinae) with a shared tRNA rearrangement and its phylogenetic implications. Int J Biol Macromol. 2019;122:1027–34.

    Article  CAS  PubMed  Google Scholar 

  33. Yang HY, Dai RH. Complete mitochondrial genome of Mukaria splendida Distant (Hemiptera: Cicadellidae: Deltocephalinae: Mukariini) and phylogenetic analysis. Mitochondrial DNA Part B. 2021;6(2):622–3.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Xing JC, Wang JJ. Complete mitochondrial genome of Paralaevicephalus gracilipenis (Hemiptera: Cicadellidae: Deltocephalinae) from China. Mitochondrial DNA Part B. 2019;4:1372–3.

    Article  Google Scholar 

  35. Yu P, Wang M, Cui L, Chen X, Han B. The complete mitochondrial genome of Tambocerus sp. (Hemiptera: Cicadellidae). Mitochondrial DNA Part A. 2017;28:133–4. https://doi.org/10.3109/19401736.2015.1111357.

    Article  CAS  Google Scholar 

  36. Wu YF, Dai RH, Zhan HP, Qu L. Complete mitochondrial genome of Drabescoides nuchalis (Hemiptera: Cicadellidae). Mitochondrial DNA Part A. 2016;27(5):3626–7. https://doi.org/10.3109/19401736.2015.1079827.

    Article  CAS  Google Scholar 

  37. Xu DL, Yu TH, Zhang Y. Characterization of the complete mitochondrial genome of Drabescus ineffectus and Roxasellana stellata (Hemiptera: Cicadellidae: Deltocephalinae: Drabescini) and their phylogenetic implications. Insects. 2020;11(8):534.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Wang JJ, Wu YF, Yang MF, Dai RH. The Phylogenetic Implications of the Mitochondrial Genomes of Macropsis notata and Oncopsis nigrofasciata. Front Genet. 2020;11:443.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wang XY, Wang JJ, Dai RH. Mitogenomics of five Olidiana leafhoppers (Hemiptera: Cicadellidae: Coelidiinae) and their phylogenetic implications. PeerJ. 2021;9:e11086. https://doi.org/10.7717/peerj.11086.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Wang JJ, Wu YF, Dai RH, Yang MF. Comparative mitogenomes of six species in the subfamily Iassinae (Hemiptera: Cicadellidae) and phylogenetic analysis. Int J Biol Macromol. 2020;149:1294–303.

    Article  CAS  PubMed  Google Scholar 

  41. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinform. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.

    Article  CAS  Google Scholar 

  42. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9. https://doi.org/10.1093/bioinformatics/bts199.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7. https://doi.org/10.1093/nar/gkh340.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Phylogenet Evol. 2016;33:1870–4. https://doi.org/10.1093/molbev/msw054.

    Article  CAS  Google Scholar 

  45. Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Pütz J, Middendorf M, Stadler PF. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69:313–9. https://doi.org/10.1016/j.ympev.2012.08.023.

    Article  PubMed  Google Scholar 

  46. Perna NT, Kocher TD. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J Mol Evol. 1995;41:353–8.

    Article  CAS  PubMed  Google Scholar 

  47. Grant JR, Stothard P. The CGView Server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 2008;36(Suppl 2):W181. https://doi.org/10.1093/nar/gkn179.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Xiang CY, Gao FL, Jakovlić I, Lei HP, Hu Y, Zhang H, Zou H, Wang GT, Zhang D. Using PhyloSuite for molecular phylogeny and tree-based analyses. iMeta. 2023;e87:1–42. https://doi.org/10.1002/imt2.87.

    Article  Google Scholar 

  49. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80. https://doi.org/10.1093/molbev/mst010.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34:772–3.

    Google Scholar 

  51. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MRBAYES 3.2: efficient Bayesian phylogenetic inference and model selection across a large model space. Syst Biol. 2012;61:539–42. https://doi.org/10.1093/sysbio/sys029.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Rambaut A: FigTree Version 1.3. 1. Available in: http://tree.bio.ed.ac.uk/software/figtree/. Accessed 4 Jan 2011.

  53. Letunic I, Bork P: Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 2021,49(W1).

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (32060123, 31301909).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, X.J.C..; writing—original draft preparation, M.A.H. and Z.X.T.; writing—review and editing, X.J.C., M.A.H., Z.X.T. and R.R.S; supervision, X.J.C.; figures and tables, M.A.H. and Z.X.T.; funding acquisition, X.J.C. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Jichun Xing.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent to publication

Not applicable

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary Figure S1. Inferred secondary structures of 22 tRNA genes in the mitochondrial genome of Abrus daozhenensis. Watson-Crick base pairings are illustrated by lines (-), whereas GU base pairings are illustrated by dots (·).Supplementary Figure S2. Inferred secondary structures of 22 tRNA genes in the mitochondrial genome of Abrus yunshanensis. Watson-Crick base pairings are illustrated by lines (-), whereas GU base pairings are illustrated by dots (·).Supplementary Figure S3. Inferred secondary structures of 22 tRNA genes in the mitochondrial genome of Abrus expansivus. Watson-Crick base pairings are illustrated by lines (-), whereas GU base pairings are illustrated by dots (·).Supplementary Figure S4. Phylogenetic tree produced by Bayesian inference analysis of the PCG123 datasets. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S5. Phylogenetic tree produced by Bayesian inference analysis of the PCG123 + 2 rRNA datasets. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S6. Phylogenetic tree produced by Bayesian inference analysis of the PCG123 + 2 rRNA + 22 tRNA datasets. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S7. Phylogenetic tree produced by Bayesian inference analysis of the PCG123_AA dataset. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S8. Phylogenetic tree produced by Bayesian inference analysis of the PCG123_AA + 2 rRNA datasets. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S9. Phylogenetic tree produced by Bayesian inference analysis of the PCG123_AA + 2 rRNA + 22 tRNA datasets. Numbers at nodes are Bayesian posterior probability (BPP) support values. Supplementary Figure S10. Phylogenetic tree produced by maximum likelihood analyses based on PCG123 datasets. Numbers at nodes are bootstrap support values (BS). Supplementary Figure S11. Phylogenetic tree produced by maximum likelihood analyses based on PCG123+ 2 rRNA datasets. Numbers at nodes are bootstrap support values (BS). Supplementary Figure S12. Phylogenetic tree produced by maximum likelihood analyses based on PCG123+ 2 rRNA + 22 tRNA datasets. Numbers at nodes are bootstrap support values (BS). Supplementary Figure S13. The phylogenetic tree produced by maximum likelihood analyses based on PCG123_AA datasets. Numbers at nodes are bootstrap support values (BS). Supplementary Figure S14. The phylogenetic tree produced by maximum likelihood analyses based on PCG123_AA+ 2 rRNA datasets. Numbers at nodes are bootstrap support values (BS). Supplementary Figure S15. The phylogenetic tree produced by maximum likelihood analyses based on PCG123_AA+ 2 rRNA + 22 tRNA datasets. Numbers at nodes are bootstrap support values (BS).

Additional file 2:

 Supplementary Table S1. Organization of the Abrus daozhenensis mitochondrial genome. Supplementary Table S2. Organization of the Abrus yunshanensis mitochondrial genome. Supplementary Table S3. Organization of the Abrus expansivus mitochondrial genome. Supplementary Table S4. Overall A+T content in the first (P1), second (P2) and third (P3) codon positions. Supplementary Table S5. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG dataset selected by PartitionFinder. Supplementary Table S6. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG + 2 rRNA dataset selected by PartitionFinder. Supplementary Table S7. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG + 2 rRNA + 22 tRNA dataset selected by PartitionFinder. Supplementary Table S8. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG_AA dataset selected by PartitionFinder. Supplementary Table S9. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG_AA + 2 rRNA dataset selected by PartitionFinder. Supplementary Table S10. The best partitioning schemes and models for the Bayesian inference (BI) method based on 123PCG_AA + 2 rRNA + 22 tRNA dataset selected by PartitionFinder. Supplementary Table S11. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG dataset selected by PartitionFinder. Supplementary Table S12. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG + 2rRNA dataset selected by PartitionFinder. Supplementary Table S13. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG + 2rRNA + 22tRNA dataset selected by PartitionFinder. Supplementary Table S14. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG_AA dataset selected by PartitionFinder. Supplementary Table S15. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG_AA + 2rRNA dataset selected by PartitionFinder. Supplementary Table S16. The best partitioning schemes and models for maximum-likelihood (ML) analyses on 123PCG_AA + 2rRNA + 22tRNA dataset selected by PartitionFinder.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hassan, M.A., Tan, Z., Shen, R. et al. Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication. BMC Genomics 24, 714 (2023). https://doi.org/10.1186/s12864-023-09809-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09809-0

Keywords