Secondary structure models of 18S and 28S rRNAs of the true bugs based on complete rDNA sequences of Eurydema maracandica Oshanin, 1871 (Heteroptera, Pentatomidae)

Phylogenomics of the superfamily Lygaeoidea with proposals on taxonomic system (Hemiptera: Heteroptera)

Lygaeoidea, one of the largest superfamilies of the true bugs (Hemiptera: Heteroptera), currently comprises 16 families and exhibits abundant diversification in habitats, body plans, and feeding habits. The relationships within this superfamily are complex, and comprehensive molecular phylogenetic studies, particularly those focusing on subfamilies and tribes within Rhyparochromidae, have been lacking. In the present study, we sampled 125 lygaeoid species representing all 16 families of Lygaeoidea, with a focus on two subfamilies and 12 tribes within Rhyparochromidae. A phylogenetic dataset with 102 genes was assembled, including two nuclear rRNA genes (18S rDNA, 28S rDNA), two mitochondrial rRNA genes (12S rDNA, 16S rDNA), 13 mitochondrial protein-coding genes (PCGs) and 85 nuclear PCGs generated from the low-coverage genomes. Our inferences indicate that Rhyparochromidae is not monophyletic, with one subfamily, one tribe, and one genus deserving elevation to the rank of family. Additionally, this phylogenetic result is also supported by corresponding morphological evidence. Besides, the transfer of the Heissothignus from Heterogastridae to Meschiidae is supported by molecular evidence in this study.

Phylogenetic placement and comparative analysis of the mitochondrial genomes of Idiostoloidea (Hemiptera: Heteroptera)

The classification system and the higher level phylogenetic relationships of Pentatomomorpha, the second largest infraorder of Heteroptera (Insecta: Hemiptera), have been debated and remain controversial over decades. In particular, the placement and phylogenetic relationship of Idiostoloidea are not well resolved, which hampers a better understanding of the evolutionary history of Pentatomomorpha. In this study, for the first time, we reported the complete mitochondrial genome for two narrowly distributed families of Idiostoloidea (including Idiostolidae and Henicocoridae), respectively. The length of the mitochondrial genome of Monteithocoris hirsutus and Henicocoris sp. is 16,632 and 16,013 bp, respectively. The content of AT is ranging from 75.15% to 80.48%. The mitogenomic structure of Idiostoloidea is highly conservative and there are no gene arrangements. By using the Bayesian inference, maximum likelihood, and Bayesian site-heterogeneous mixture model, we inferred the phylogenetic relationships within Pentatomomorpha and estimated their divergence times based on concatenated mitogenomes and nuclear ribosomal genes. Our results support the classification system of six superfamilies within Pentatomomorpha and confirm the monophyletic groups of each superfamily, with the following phylogenetic relationships: (Aradoidea + (Pentatomoidea + (Idiostoloidea + (Coreoidea + (Pyrrhocoroidea + Lygaeoidea))))). Furthermore, estimated divergence times revealed that most pentatomomorphan superfamilies and families diverged during the Late Jurassic to Early Cretaceous, which coincides with the explosive radiation of angiosperms.

Inconsistencies in the Classification of the Family Cydnidae (Hemiptera: Heteroptera: Pentatomoidea) Revealed by Molecular Apomorphies in the Secondary and Tertiary Structures of 18S rRNA Length-Variable Region L (LVR L)

The SSU nuclear rDNA (encoding 18S ribosomal RNA) is one of the most frequently sequenced genes in the molecular analysis of insects. Molecular apomorphies in the secondary and tertiary structures of several 18S rRNA length-variable regions (LVRs) located within the V2, V4, and V7 hypervariable regions can be good indicators for recovering monophyletic groups within some heteropteran families. Among the LVRs that have been analysed, the LVR L in the V4 hypervariable region is the longest and most crucial for such assessments. We analysed the 18S rRNA V4 hypervariable region sequences of 45 species from the family Cydnidae, including all 6 subfamilies (Amaurocorinae, Amnestinae, Cephalocteinae, Cydninae, Garsauriinae, and Sehirinae) and three pentatomoid families (Parastrachiidae, Thaumastellidae, and Thyreocoridae), which have often been included in the broadly defined Cydnidae family. This is the first time that representatives of all Cydnidae subfamilies have been included in a molecular analysis. Only taxa from two subfamilies, Sehirinae and Cydninae, have been used in previous molecular studies. The secondary and tertiary structures of the LVR L were predicted for each species using the two-step procedure already accepted for such analyses to recover any molecular apomorphy essential for determining monophyly. The results of our comparative studies contradict the current understanding of the relationships among burrowing bugs and the current family classification.

Different roles of host and habitat in determining the microbial communities of plant-feeding true bugs

BACKGROUND

The true bugs (Heteroptera) occupy nearly all of the known ecological niches of insects. Among them, as a group containing more than 30,000 species, the phytophagous true bugs are making increasing impacts on agricultural and forestry ecosystems. Previous studies proved that symbiotic bacteria play important roles in these insects in fitting various habitats. However, it is still obscure about the evolutionary and ecological patterns of the microorganisms of phytophagous true bugs as a whole with comprehensive taxon sampling.

RESULTS

Here, in order to explore the symbiotic patterns between plant-feeding true bugs and their symbiotic microorganisms, 209 species belonging to 32 families of 9 superfamilies had been sampled, which covered all the major phytophagous families of true bugs. The symbiotic microbial communities were surveyed by full-length 16S rRNA gene and ITS amplicons respectively for bacteria and fungi using the PacBio platform. We revealed that hosts mainly affect the dominant bacteria of symbiotic microbial communities, while habitats generally influence the subordinate ones. Thereafter, we carried out the ancestral state reconstruction of the dominant bacteria and found that dramatic replacements of dominant bacteria occurred in the early Cretaceous and formed newly stable symbiotic relationships accompanying the radiation of insect families. In contrast, the symbiotic fungi were revealed to be horizontally transmitted, which makes fungal communities distinctive in different habitats but not significantly related to hosts.

CONCLUSIONS

Host and habitat determine microbial communities of plant-feeding true bugs in different roles. The symbiotic bacterial communities are both shaped by host and habitat but in different ways. Nevertheless, the symbiotic fungal communities are mainly influenced by habitat but not host. These findings shed light on a general framework for future microbiome research of phytophagous insects. Video Abstract.

Tribe Acalyptaini (Hemiptera: Tingidae: Tinginae) Revisited: Can Apomorphies in Secondary and Tertiary Structures of 18S rRNA Length-Variable Regions (LVRs) Support Tribe Validity?

The lace bug tribe Acalyptaini (Tingidae: Tinginae) includes five genera, Acalypta, Derephysia, Dictyonota, Kalama, and Recaredus, and it was recently resurrected based on morphological and karyological characters. We aimed to validate the distinctiveness of this tribe using 18S rDNA sequences, which have not been used in previous Tingidae phylogenomic studies. Our results confirmed the monophyly of the tribe. Moreover, the monophyly of the subfamily Cantacaderinae and its basal position within the family Tingidae were indicated, as well as the position of the tribe Litadeini as sister to all other Tinginae. In addition, we attempted to determine the apomorphic morpho-molecular characters in the secondary and tertiary structures of length-variable regions of the 18S rRNA sequences of the analysed species. The results showed that two LVRs (LVR X and LVR L) of the hypervariable region V4 exhibited significant variability in the number of nucleotides and could be considered for apomorphic recognition.

Phylogeny and diversification of planthoppers (Hemiptera Fulgoromorpha) based on a comprehensive molecular dataset and large taxon sampling

Our understanding of the evolution of Fulgoromorpha (Insects, Hemiptera) has relied on molecular studies that have only considered either a limited number of taxa where all the families were not represented simultaneously, or a reduced number of genes.The absence of a global analysis comparing all the available data has thus led to significant biases in the analyzes, as evidenced by the incongruence of the results reported for planthopper phylogeny. Here we provide a phylogenetic and dating analysis of the Fulgoromorpha with a large sampling of 531 ingroup taxa, representing about 80% of the currently described suprageneric taxonomic diversity in this group. This study is based on most of the molecular sequences available to date and duly verified, for a set of nuclear and mitochondrial genes from a taxonomic sampling as complete as possible. The most significant results of our study are: (1) the unexpected paraphyly of Delphacidae whose Protodelphacida seem more related to Cixiidae than to other Delphacidae;(2) the group Meenoplidae-Kinnaridae recovered sister to the remaining Fulgoroidea families; (3) the early branching node of Tettigometridae sister of all the other families;(4) the Achilidae-Derbidae clade with Achilidae Plectoderini including Achilixiidae recovered as monophyletic as well as theFulgoridae-Dictyopharidae clade; and (5) the Tropiduchidae placed sister to the other so called 'higher' families (sec. Shcherbakov, 2006).Our divergence times analysis, calibrated with a set of duly verified fossils, suggests that the first diversification of planthoppers occurred in the Early Triassic around 240 Mya and those of the superfamilies Delphacoidea and Fulgoroidea in the Middle-Late Triassic around 210 Mya and 230 Mya, respectively. By the end of the Jurassic, all major planthopper lineages were originated, and all families, around 125 Mya, might havebeen driven in their distribution and evolution (in their first subfamilial divisions) by the geographical constraints of the Gondwanan break-up.Rapid evolutionary radiations occurred particularly in Fulgoridae around 125-130 Mya. Our results stress the importance of the good quality of the sequences used in the molecular analyzes and the primordial importance of a large sampling when analyzing the phylogeny of the group.

Molecular Apomorphies in the Secondary and Tertiary Structures of Length-Variable Regions (LVRs) of 18S rRNA Shed Light on the Systematic Position of the Family Thaumastellidae (Hemiptera: Heteroptera: Pentatomoidea)

The SSU nrDNA, a small subunit of the nuclear ribosomal DNA (coding 18S rRNA), is one of the most frequently sequenced genes in molecular studies in Hexapoda. In insects, including true bugs (Hemiptera: Heteroptera), only its primary structures (i.e., aligned sequences) are predominantly used in phylogenetic reconstructions. It is known that including RNA secondary structures in the alignment procedure is essential for improving accuracy and robustness in phylogenetic tree reconstruction. Moreover, local plasticity in rRNAs might impact their tertiary structures and corresponding functions. To determine the systematic position of Thaumastellidae within the superfamily Pentatomoidea, the secondary and-for the first time among all Hexapoda-tertiary structures of 18S rRNAs in twelve pentatomoid families were compared and analysed. Results indicate that the shapes of the secondary and tertiary structures of the length-variable regions (LVRs) in the 18S rRNA are phylogenetically highly informative. Based on these results, it is suggested that the Thaumastellidae is maintained as an independent family within the superfamily Pentatomoidea, rather than as a part of the family Cydnidae. Moreover, the analyses indicate a close relationship between Sehirinae and Parastrachiidae, expressed in morpho-molecular synapomorphies in the predicted secondary and tertiary structures of the length-variable region L (LVR L).

Symbiotic Microorganisms and Their Different Association Types in Aquatic and Semiaquatic Bugs

True bugs (Hemiptera, suborder Heteroptera) constitute the largest suborder of nonholometabolous insects and occupy a wide range of habitats various from terrestrial to semiaquatic to aquatic niches. The transition and occupation of these diverse habitats impose various challenges to true bugs, including access to oxygen for the aquatic species and plant defense for the terrestrial phytophagans. Although numerous studies have demonstrated that microorganisms can provide multiple benefits to terrestrial host insects, a systematic study with comprehensive higher taxa sampling that represents aquatic and semiaquatic habitats is still lacking. To explore the role of symbiotic microorganisms in true bug adaptations, 204 samples belonging to all seven infraorders of Heteroptera were investigated, representing approximately 85% of its superfamilies and almost all known habitats. The symbiotic microbial communities of these insects were analyzed based on the full-length amplicons of the bacterial 16S rRNA gene and fungal ITS region. Bacterial communities varied among hosts inhabiting terrestrial, semiaquatic, and aquatic habitats, while fungal communities were more related to the geographical distribution of the hosts. Interestingly, co-occurrence networks showed that species inhabiting similar habitats shared symbiotic microorganism association types. Moreover, functional prediction analyses showed that the symbiotic bacterial community of aquatic species displayed richer amino acid and lipid metabolism pathways, while plant-feeding true bugs benefited more from the symbiont-provided xenobiotics biodegradation pathway. These results deepened the recognition that symbiotic microorganisms were likely to help heteropterans occupy diverse ecological habitats and provided a reference framework for further studies on how microorganisms affect host insects living in various habitats. IMPORTANCE Symbiotic bacteria and fungi generally colonize insects and provide various benefits for hosts. Although numerous studies have investigated symbionts in terrestrial plant-feeding insects, explorations of symbiotic bacterial and fungal communities in aquatic and semiaquatic insects are rare. In this study, the symbiotic microorganisms of 204 aquatic, semiaquatic, and terrestrial true bugs were explored. This comprehensive taxon sampling covers ~85% of the superfamilies of true bugs and most insect habitats. Analyses of the diversity of symbionts demonstrated that the symbiotic microbial diversities of true bugs were mainly affected by host habitats. Co-occurrence networks showed that true bugs inhabiting similar habitats shared symbiotic microbial association types. These correlations between symbionts and hosts together with the functions of bacterial communities indicated that symbiotic microbial communities may help true bugs adapt to (semi)aquatic habitats.

Diversification of the phytophagous lineages of true bugs (Insecta: Hemiptera: Heteroptera) shortly after that of the flowering plants

More than 95% of phytophagous true bug (Hemiptera: Heteroptera) species belong to four superfamilies: Miroidea (Cimicomorpha), Pentatomoidea, Coreoidea, and Lygaeoidea (all Pentatomomorpha). These iconic groups of highly diverse, overwhelmingly phytophagous insects include several economically prominent agricultural and silvicultural pest species, though their evolutionary history has not yet been well resolved. In particular, superfamily- and family-level phylogenetic relationships of these four lineages have remained controversial, and the divergence times of some crucial nodes for phytophagous true bugs have hitherto been little known, which hampers a better understanding of the evolutionary processes and patterns of phytophagous insects. In the present study, we used 150 species and concatenated nuclear and mitochondrial protein-coding genes and rRNA genes to infer the phylogenetic relationships within the Terheteroptera (Cimicomorpha + Pentatomomorpha) and estimated their divergence times. Our results support the monophyly of Cimicomorpha, Pentatomomorpha, Miroidea, Pentatomoidea, Pyrrhocoroidea, Coreoidea, and Lygaeoidea. The phylogenetic relationships across phytophagous lineages are largely congruent at deep nodes across the analyses based on different datasets and tree-reconstructing methods with just a few exceptions. Estimated divergence times and ancestral state reconstructions for feeding habit indicate that phytophagous true bugs explosively radiated in the Early Cretaceous-shortly after the angiosperm radiation-with the subsequent diversification of the most speciose clades (Mirinae, Pentatomidae, Coreinae, and Rhyparochromidae) in the Late Cretaceous.

Mitochondrial Genomes from Two Specialized Subfamilies of Reduviidae (Insecta: Hemiptera) Reveal Novel Gene Rearrangements of True Bugs

Reduviidae, a hyper-diverse family, comprise 25 subfamilies with nearly 7000 species and include many natural enemies of crop pests and vectors of human disease. To date, 75 mitochondrial genomes (mitogenomes) of assassin bugs from only 11 subfamilies have been reported. The limited sampling of mitogenome at higher categories hinders a deep understanding of mitogenome evolution and reduviid phylogeny. In this study, the first mitogenomes of Holoptilinae (Ptilocnemus lemur) and Emesinae (Ischnobaenella hainana) were sequenced. Two novel gene orders were detected in the newly sequenced mitogenomes. Combined 421 heteropteran mitogenomes, we identified 21 different gene orders and six gene rearrangement units located in three gene blocks. Comparative analyses of the diversity of gene order for each unit reveal that the tRNA gene cluster trnI-trnQ-trnM is the hotspot of heteropteran gene rearrangement. Furthermore, combined analyses of the gene rearrangement richness of each unit and the whole mitogenome among heteropteran lineages confirm Reduviidae as a 'hot-spot group' of gene rearrangement in Heteroptera. The phylogenetic analyses corroborate the current view of phylogenetic relationships between basal groups of Reduviidae with high support values. Our study provides deeper insights into the evolution of mitochondrial gene arrangement in Heteroptera and the early divergence of reduviids.

The SSU rRNA secondary structures of the Plagiorchiida species (Digenea), its applications in systematics and evolutionary inferences

The small subunit ribosomal RNA (SSU rRNA) is widely used phylogenetic marker in broad groups of organisms and its secondary structure increasingly attracts the attention of researchers as supplementary tool in sequence alignment and advanced phylogenetic studies. Its comparative analysis provides a great contribution to evolutionary biology, allowing find out how the SSU rRNA secondary structure originated, developed and evolved. Herein, we provide the first data on the putative SSU rRNA secondary structures of the Plagiorchiida species. The structures were found to be quite conserved across broad range of species studied, well compatible with those of others eukaryotic SSU rRNA and possessed some peculiarities: cross-shaped structure of the ES6b, additional shortened ES6c2 helix, and elongated ES6a helix and h39 + ES9 region. The secondary structures of variable regions ES3 and ES7 appeared to be tissue-specific while ES6 and ES9 were specific at a family level allowing considering them as promising markers for digenean systematics. Their uniqueness more depends on the length than on the nucleotide diversity of primary sequences which evolutionary rates well differ. The findings have important implications for understanding rRNA evolution, developing molecular taxonomy and systematics of Plagiorchiida as well as for constructing new anthelmintic drugs.

Characterization of the complete mitochondrial genome and phylogenetic implications for Eurydema maracandica (Hemiptera: Pentatomidae)

The complete mitochondrial genome of Eurydema maracandica was sequenced and was 15,391 bp long with 76.82% A + T. There were 37 typical genes including 13 protein-coding genes/PCGs, 22 transfer RNA/tRNAs, and 2 ribosomal RNA/rRNAs. The 690 bp D-loop region was located between 12S rRNA and trnI-Ile. All PCGs started with ATN codons, except COI, ATP8, and ND1, and ended with TAA, except COI, COII, and ND5. Phylogenetic analyses indicated highly supported monophyly for each family, and Pentatomidae species formed a solid monophyletic group. Eurydema maracandica and Eurydema gebleri were clustered sibling clades, and the genus Halyomorpha was close to Eurydema.

Sequencing of the mitochondrial genome of the avocado lace bug Pseudacysta perseae (Heteroptera, Tingidae) using a genome skimming approach

Lace bugs (Tingidae) are a family of phytophagous heteropterans, some of which are important agricultural and forestry pests. They currently comprise around 2500 species distributed worldwide, for which only one mitochondrial genome has been described so far. We sequenced the complete mitochondrial genome and the nuclear ribosomal gene segment of the avocado lace bug Pseudacysta perseae using a genome skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including the one of the sycamore lace bug Corythucha ciliata, were retrieved to allow for comparisons and phylogenetic analyses. P. perseae mitochondrial genome was determined to be 15,850 bp long, and presented the typical organisation of insect mitogenomes. The phylogenetic analysis placed P. perseae as a sister to C. ciliata but did not confirm the monophyly of Miroidae including Tingidae. Our results contradicted widely accepted phylogenetic hypothesis, which highlights the limits of analyses based on mitochondrial data only. Shotgun sequencing approaches should provide substantial improvements in harmonizing mitochondrial and nuclear databases.

Shotgun assembly of the assassin bug Brontostoma colossus mitochondrial genome (Heteroptera, Reduviidae)

The complete mitochondrial genome of the assassin bug Brontostoma colossus (Distant, 1902) (Heteroptera: Reduviidae) has been sequenced using a genome-skimming approach on an Illumina Hiseq 2000 platform. Fifty-four additional heteropteran mitogenomes, including five assassin bug species, were retrieved to allow for comparisons and phylogenetic analyses. The mitochondrial genome of B. colossus was determined to be 16,625 bp long, and consists of 13 protein-coding genes (PCGs), 23 transfer-RNA genes (tRNAs), two ribosomal-RNA genes (rRNAs), and one control region. The nucleotide composition is biased toward adenine and thymine (A+T=73.4%). Overall, architecture, nucleotide composition and genome asymmetry are similar among all available assassin bug mitogenomes. All PCGs have usual start-codons (Met and Ile). Three T and two TA incomplete termination codons were identified adjacent to tRNAs, which was consistent with the punctuation model for primary transcripts processing followed by 3' polyadenylation of mature mRNA. All tRNAs exhibit the classic clover-leaf secondary structure except for tRNASer(AGN) in which the DHU arm forms a simple loop. Two notable features are present in the B. colossus mitogenome: (i) a 131 bp duplicated unit including the complete tRNAArg gene, resulting in 23 potentially functional tRNAs in total, and (ii) a 857 bp duplicated region comprising 277 bp of the srRNA gene and 580 bp of the control region. A phylogenetic analysis based on 55 true bug mitogenomes confirmed that B. colossus belongs to Reduviidae, but contradicted a widely accepted hypothesis. This highlights the limits of phylogenetic analyses based on mitochondrial data only.