Comparative mitochondrial genomes of the Rhus gall aphid Kaburagia rhusicola subspecies with variable gall shapes

Abstract

Rhus gall aphids (Hemiptera: Aphididae: Eriosomatinae) stimulate the formation of galls on their primary host plants (sumacs: Rhus spp., Anacardiaceae). The shapes of galls are often used as an extended phenotype to identify the aphid species and subspecies. We collected four Rhus galls with conspicuously different shapes formed by Kaburagia rhusicola aphids, whose sequences of the complete mitochondrial genomes (mitogenomes) were obtained by high-throughput sequencing. Each mitogenome was assembled into a circular molecule containing 13 protein-coding genes, two rRNAs, 22 tRNAs and one control region. All the protein-coding genes had a typical ATN initiation codon and TAA termination codon except for cox1 and nad4, which had a single T as stop codon. All the tRNAs could be folded as a typical clover-leaf secondary structure, except for trnS1 lacking a dihydrouridine (DHU) arm. The relative synonymous codon usage and ratio of nonsynonymous to synonymous substitution rates showed that the four K. rhusicola samples were highly similar to the subspecies K. r. ovogallis. The phylogenetic analyses grouped these samples with K. r. ovogallis in a clade sister to K. r. rhusicola. All these molecular analyses demonstrated that our current samples represented one subspecies of Kaburagia rhusicola, i.e., K. r. ovogallis, and the gall shape was variable even at the subspecies level in Kaburagia gall aphids.

Keywords: Gall shape; Kaburagia; Mitochondrial genome; Phylogeny; Rhus gall aphid.