Application of RNA-seq for mitogenome reconstruction, and reconsideration of long-branch artifacts in Hemiptera phylogeny

Phylogenomic analysis of Syrphoidea (Diptera: Syrphidae, Pipunculidae) based on the expanded mitogenomic data

The mitochondrial genome has become the most widely used genomic resource in resolving the insect phylogenetic relationships. In this study, we assess the interrelationships among the syrphid and pipunculid members of Syrphoidea using mitochondrial genome sequences of 152 taxa, 9 of which are newly reported and three are assembled from the existing transcriptome data. The Pipunculidae was found to be deeply nested members of Schizophora, which resulted in a nonmonophyletic Syrphoidea. In the monophyletic Syrphidae, unequivocal robust support was found for Microdontinae as the sister group of all other Syrphidae. The subfamily Eristalinae was nonmonophyletic. The Pipizinae was recovered as the sister group to the Syrphinae, albeit with strong support. As a whole, our results are concord with previously established hypotheses on Syrphoidea from the genome scale data. The mitochondrial genomes were successful in producing a robustly supported phylogenetic framework for the Syrphoidea.

Expansions and contractions in gene families of independently-evolved blood-feeding insects

BACKGROUND

The blood-feeding behavior evolved multiple times in Insecta lineages and it represents an excellent opportunity to study patterns of convergent molecular evolution regarding this habit. In insects the expansion of some gene families is linked with blood-feeding behavior, but a wide study comparing the evolution of these gene families among different lineages is still missing. Here we gathered genomic data from six independently-evolved hematophagous lineages, aiming to identify convergent expansions and/or contractions of gene families in hematophagous lineages of insects.

RESULTS

We found four rapidly evolving gene families shared by at least two hematophagous independently-evolved lineages, including a heat-shock and a chemosensory protein. On the expression of these four rapidly evolving gene families we found more genes expressed in mated individuals compared with virgin individuals in rapidly-expanded families and more genes expressed in non-blood-feeding individuals compared with blood-feeding individuals in rapidly-contracted families.

CONCLUSION

Our results reveal a new set of candidate genes to be explored in further analysis to help the development of new strategies to deal with blood-feeding vectors and also presents a new perspective to study the evolution of hematophagy identifying convergent molecular patterns.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

Description of soybean aphid (Aphis glycines Matsumura) mitochondrial genome and comparative mitogenomics of Aphididae (Hemiptera: Sternorrhyncha)

The complete mitochondrial genome of the soybean aphid (Aphis glycines Matsumura), a major agricultural pest in the world, is described for the first time, which consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, as well as a large repeat region between tRNA-Glu and tRNA-Phe, and an AT-rich control region. The 17,954 bp mtgenome is the largest reported from the family Aphididae, and its gene order follows the ancestral insect mtgenome except for the repeat region, which contains a 195 bp unit repeated 11.9 times, representing the highest reported repeats among the known aphid mtgenomes to date. A new molecular phylogeny of Aphidae is reconstructed based on all available aphid mtgenomes, and it is shown that the mtgenome data can robustly resolve relationships at the subfamily level, but do not have sufficient phylogenetic information to resolve deep relationships. A phylogeny-based comparative analysis of mtgenomes has been performed to investigate the evolution of the repeat region between tRNA-Glu and tRNA-Phe. So far, among aphids, 13 species are known to have this repeat region of variable lengths, and a phylogenetic analysis of the repeat region shows that a large proportion of the sequences are conserved across the phylogeny, suggesting that the repeat region evolved in the most recent common ancestor of Aphidinae and Eriosomatinae, and that it has gone through numerous episodes of lineage-specific losses and expansions. Combined together, this study provides novel insights into how the repeat regions have evolved within aphids.

The complete mt genomes of Lutzia halifaxia, Lt. fuscanus and Culex pallidothorax (Diptera: Culicidae) and comparative analysis of 16 Culex and Lutzia mt genome sequences

BACKGROUND

Despite the medical importance of the genus Culex, the mitochondrial genome (mt genome) characteristics of Culex spp. are not well understood. The phylogeny of the genus and particularly the generic status of the genus Lutzia and the subgenus Culiciomyia remain unclear.

METHODS

The present study sequenced and analyzed the complete mt genomes of Lutzia halifaxia, Lutzia fuscanus and Cx. (Culiciomyia) pallidothorax and assessed the general characteristics and phylogenetics of all known 16 mt genome sequences for species in the genera Culex and Lutzia.

RESULTS

The complete mt genomes of Lt. halifaxia, Lt. fuscanus and Cx. pallidothorax are 15,744, 15,803 and 15,578 bp long, respectively, including 13 PCGs, 22 tRNAs, two tRNAs and a control region (CR). Length variations in the Culex and Lutzia mt genomes involved mainly the CR, and gene arrangements are the same as in other mosquitoes. We identified four types of repeat units in the CR sequences, and the poly-T stretch exists in all of these mt genomes. The repeat units of CR are conserved to different extent and provide information on their evolution. Phylogenetic analyses demonstrated that the Coronator and Sitiens groups are each monophyletic, whereas the monophyletic status of the Pipiens Group was not supported; Cx. pallidothorax is more closely related to the Sitiens and Pipiens groups; and both phylogenetics analysis and repeat unit features in CR show that Lutzia is a characteristic monophyletic entity, which should be an independent genus.

CONCLUSIONS

To our knowledge, this is the first comprehensive review of the mt genome sequences and taxonomic discussion based on the mt genomes of Culex spp. and Lutzia spp. The research provides general information on the mt genome of these two genera, and the phylogenetic and taxonomic status of Lutzia and Culiciomyia.

Insights into the phylogeny of Hemiptera from increased mitogenomic taxon sampling

Although reconstruction of the phylogeny of Hemiptera has progressed tremendously over the past two decades, some higher-level relationships remain poorly resolved. Here, we investigated the Hemiptera higher-level relationships using full mitochondrial genome data from 357 ingroup species, representing the most comprehensive sampling yet undertaken for reconstructing the phylogeny of this group. In this study, 92 mitochondrial genomes were newly determined. Various data treatment methods and substitution models were applied to tree reconstructions. Effects of compositional heterogeneity, rate heterogeneity, model adequacy and taxon sampling on support values and topological stability were explored. Phylogenetic analyses (1) confirmed the monophyly of Hemiptera under site-heterogeneous model, (2) placed Sternorrhyncha as sister to all other Hemiptera, (3) recovered Coccoidea as the sister taxon of Aphidoidea, followed successively by Aleyrodoidea and Psylloidea, and (4) indicated that the grouping of Coleorrhyncha and Fulgoromorpha was the result of long-branch attraction effect.