How to sequence and annotate insect mitochondrial genomes for systematic and comparative genomics research

Over the past decade the mitochondrial (mt) genome has become the most widely used genomic resource available for systematic entomology. While the availability of other types of ‘–omics’ data – in particular transcriptomes – is increasing rapidly, mt genomes are still vastly cheaper to sequence and are far less demanding of high quality templates. Furthermore, almost all other ‘–omics’ approaches also sequence the mt genome, and so it can form a bridge between legacy and contemporary datasets. Mitochondrial genomes have now been sequenced for all insect orders, and in many instances representatives of each major lineage within orders (suborders, series or superfamilies depending on the group). They have also been applied to systematic questions at all taxonomic scales from resolving interordinal relationships (e.g. Cameron et al., 2009; Wan et al., 2012; Wang et al., 2012), through many intraordinal (e.g. Dowton et al., 2009; Timmermans et al., 2010; Zhao et al. 2013a) and family-level studies (e.g. Nelson et al., 2012; Zhao et al., 2013b) to population/biogeographic studies (e.g. Ma et al., 2012). Methodological issues around the use of mt genomes in insect phylogenetic analyses and the empirical results found to date have recently been reviewed by Cameron (2014), yet the technical aspects of sequencing and annotating mt genomes were not covered. Most papers which generate new mt genome report their methods in a simplified form which can be difficult to replicate without specific knowledge of the field. Published studies utilize a sufficiently wide range of approaches, usually without justification for the one chosen, that confusion about commonly used jargon such as ‘long PCR’ and ‘primer walking’ could be a serious barrier to entry. Furthermore, sequenced mt genomes have been annotated (gene locations defined) to wildly varying standards and improving data quality through consistent annotation procedures will benefit all downstream users of these datasets. The aims of this review are therefore to: 1. Describe in detail the various sequencing methods used on insect mt genomes; 2. Explore the strengths/weakness of different approaches; 3. Outline the procedures and software used for insect mt genome annotation, and; 4. Highlight quality control steps used for new annotations, and to improve the re-annotation of previously sequenced mt genomes used in systematic or comparative research.

One and the same: Integrative taxonomic evidence that Bactrocera invadens (Diptera: Tephritidae) is the same species as the Oriental fruit fly Bactrocera dorsalis

The invasive fruit fly Bactrocera invadens Drew, Tsuruta & White, and the Oriental fruit fly Bactrocera dorsalis (Hendel) are highly destructive horticultural pests of global significance. Bactrocera invadens originates from the Indian subcontinent and has recently invaded all of sub-Saharan Africa, while B. dorsalis principally occurs from the Indian subcontinent towards southern China and South-east Asia. High morphological and genetic similarity has cast doubt over whether B. invadens is a distinct species from B. dorsalis. Addressing this issue within an integrative taxonomic framework, we sampled from across the geographic distribution of both taxa and: (i) analysed morphological variation, including those characters considered diagnostic (scutum colour, length of aedeagus, width of postsutural lateral vittae, wing size, and wing shape); (ii) sequenced four loci (ITS1, ITS2, cox1 and nad4) for phylogenetic inference, and; (iii) generated a cox1 haplotype network to examine population structure. Molecular analyses included the closely related species, Bactrocera kandiensis Drew & Hancock. Scutum colour varies from red-brown to fully black for individuals from Africa and the Indian subcontinent. All individuals east of the Indian subcontinent are black except for a few red-brown individuals from China. The postsutural lateral vittae width of B. invadens is narrower than B. dorsalis from eastern Asia, but the variation is clinal, with subcontinent B. dorsalis populations intermediate in size. Aedeagus length, wing shape and wing size cannot discriminate between the two taxa. Phylogenetic analyses failed to resolve B. invadens from B. dorsalis, but did resolve B. kandiensis. Bactrocera dorsalis and B. invadens shared cox1 haplotypes, yet the haplotype network pattern does not reflect current taxonomy or patterns in thoracic colour. Some individuals of B. dorsalis/B. invadens possessed haplotypes more closely related to B. kandiensis than to conspecifics, suggestive of mitochondrial introgression between these species. The combined evidence fails to support the delimitation of B. dorsalis and B. invadens as separate biological species. Consequently, existing biological data for B. dorsalis may be applied to the invasive population in Africa. Our recommendation, in line with other recent publications, is that B. invadens be synonymized with B. dorsalis.

Leaf and inflorescence structure and phylogenetics of Tripogon and affiliated genera (Poaceae: Chloridoideae)

This thesis provides new knowledge on an understudied group of grasses, some of which are resurrection grasses (i.e. able to withstand extreme drought). The sole Australian species (Tripogon loliiformis) is morphologically diverse and could be more than one species. This study sought to determine how many species of Tripogon occur in Australia, their relationships to other species in the genus and to two other genera of resurrection grasses (Eragrostiella and Oropetium). Results of the research indicate there is not enough evidence, from DNA sequence data, to warrant splitting up T. loliiformis into multiple species. The extensive morphological diversity seems to be influenced by environmental conditions. The three genera are so closely related that they could be grouped into a single genus. This new knowledge opens up pathways for future investigations, including studying genes responsible for desiccation tolerance and the conservation of native grasses that occur in rocky habitats.

Morphology of scale lepidonts in the genusAphanius(Teleostei, Cyprinodontidae) using SEM

The microstructure of the anterior region of the scales in several species of the genus Aphanius was studied by SEM with the aim of determining whether scale morphology could be used to discriminate between the species of this genus. The characters examined concern the morphology of lepidonts, or “scale‐teeth”, their distribution and mode of implantation on the circuli. These characters were also subjected to UPGMA cluster analysis. Results from phenetic analysis of scale‐teeth characters agree overall with those of previously published morphological and biogeographical studies and in part with molecular analysis of the phylogenetic relationships between species of Aphanius. An affinity between A. danfordii and A. mento (found previously in studies based on osteological observations) was seen. The separation of A. apodus from the other species of the fasciatus group, which had also been noticed from molecular observations, was also observed, as well as the affinity of A. ginaonis with the group of A. dispar+A. sirhani. This study demonstrates that scale morphology can provide useful information on the relationships among species of the genus Aphanius encouraging the use of scale characters, combined with other traits, in phylogenetic analyses.

Structure and vascular tissue expression of duplicated TERMINAL EAR1-like paralogues in poplar

TERMINAL EAR1-like (TEL) genes encode putative RNA-binding proteins only found in land plants. Previous studies suggested that they may regulate tissue and organ initiation in Poaceae. Two TEL genes were identified in both Populus trichocarpa and the hybrid aspen Populus tremula × P. alba, named, respectively, PoptrTEL1-2 and PtaTEL1-2. The analysis of the organisation around the PoptrTEL genes in the P. trichocarpa genome and the estimation of the synonymous substitution rate for PtaTEL1-2 genes indicate that the paralogous link between these two Populus TEL genes probably results from the Salicoid large-scale gene-duplication event. Phylogenetic analyses confirmed their orthology link with the other TEL genes. The expression pattern of both PtaTEL genes appeared to be restricted to the mother cells of the plant body: leaf founder cells, leaf primordia, axillary buds and root differentiating tissues, as well as to mother cells of vascular tissues. Most interestingly, PtaTEL1-2 transcripts were found in differentiating cells of secondary xylem and phloem, but probably not in the cambium itself. Taken together, these results indicate specific expression of the TEL genes in differentiating cells controlling tissue and organ development in Populus (and other Angiosperm species).

Comparative genomic analysis of human Chlamydia pneumoniae isolates from respiratory, brain and cardiac tissues

Chlamydia pneumoniae is an obligate intracellular bacterium implicated in a wide range of human diseases including atherosclerosis and Alzheimer's disease. Efforts to understand the relationships between C. pneumoniae detected in these diseases have been hindered by the availability of sequence data for non-respiratory strains. In this study, we sequenced the whole genomes for C. pneumoniae isolates from atherosclerosis and Alzheimer's disease, and compared these to previously published C. pneumoniae genomes. Phylogenetic analyses of these new C. pneumoniae strains indicate two sub-groups within human C. pneumoniae, and suggest that both recombination and mutation events have driven the evolution of human C. pneumoniae. Further fine-detailed analyses of these new C. pneumoniae sequences show several genetically variable loci. This suggests that similar strains of C. pneumoniae are found in the brain, lungs and cardiovascular system and that only minor genetic differences may contribute to the adaptation of particular strains in human disease.

Nucleocapsid gene variability reveals two subgroups of Lettuce necrotic yellows virus

The complete nucleocapsid (N) genes of eight Australian isolates of Lettuce necrotic yellows virus (LNYV) were amplified by reverse transcription PCR, cloned and sequenced. Phylogenetic analyses of these sequences revealed two distinct subgroups of LNYV isolates. Nucleotide sequences within each subgroup were more than 96% identical but heterogeneity between groups was about 20% at the nucleotide sequence level. However, less than 4% heterogeneity was noted at the amino acid level, indicating mostly third nucleotide position changes and a strong conservation for N protein function. There was no obvious geographical or temporal separation of the subgroups in Australia.

Phylogeography of the pharaoh cuttle Sepia pharaonis based on partial mitochondrial 16S sequence data

The pharaoh cuttle Sepia pharaonis Ehrenberg, 1831 (Mollusca: Cephalopoda: Sepiida) is a broadly distributed species of substantial fisheries importance found from east Africa to southern Japan. Little is known about S. pharaonis phylogeography, but evidence from morphology and reproductive biology suggests that Sepia pharaonis is actually a complex of at least three species. To evaluate this possibility, we collected tissue samples from Sepia pharaonis from throughout its range. Phylogenetic analyses of partial mitochondrial 16S sequences from these samples reveal five distinct clades: a Gulf of Aden/Red Sea clade, a northern Australia clade, a Persian Gulf/Arabian Sea clade, a western Pacific clade (Gulf of Thailand and Taiwan) and an India/Andaman Sea clade. Phylogenetic analyses including several Sepia species show that S. pharaonis sensu lato may not be monophyletic. We suggest that "S. pharaonis" may consist of up to five species, but additional data will be required to fully clarify relationships within the S. pharaonis complex.

Identification of viral and non-viral reverse transcribing elements in pineapple ( Ananas comosus ), including members of two new badnavirus species

A previously published partial sequence of pineapple bacilliform virus was shown to be from a retrotransposon (family Metaviridae) and not from a badnavirus as previously thought. Two newly discovered sequence groups isolated from pineapple were associated with bacilliform virions and were transmitted by mealybugs. Phylogenetic analyses indicated that they were members of new badnavirus species. A third caulimovirid sequence was also amplified from pineapple, but available evidence suggests that this DNA is not encapsidated, but more likely derived from an endogenous virus.

The genetic diversity of ampeloviruses in Australian pineapples and their association with mealybug wilt disease

Pineapple mealybug wilt-associated virus 1 (PMWaV-1), 2 (PMWaV-2) and -3 (PMWaV-3) have been detected in Australian commercial pineapple crops, along with a previously undescribed ampelovirus, for which the name Pineapple mealybug wilt-associated virus 5 (PMWaV-5) is proposed. Partial sequences extending from open reading frame 1b through to the heat shock protein homologue were obtained for PMWaV-1, -3 and -5. Phylogenetic analyses of selected regions of these sequences indicated that PMWaV-5 is a distinct species and most closely related to PMWaV-1. The amino acid sequence variation observed in the RNA-dependent RNA polymerase region of PMWaV-1 isolates was 95.8–98.4% and of PMWaV-3 isolates was 92.2–99.5%. In surveys of mealybug wilt disease (MWD) affected crops, none of the four viruses was clearly associated with the disease at all survey sites. A statistically significant association (P < 0.001) between the presence of PMWaV-2 and symptoms was observed at one survey site (site 3), but the virus was at a low incidence at the remaining three survey sites. By contrast, although PMWaV-1 and -3 were equally distributed between symptomless and MWD-affected plants at site 3, there was a statistically significant (P < 0.001) association between each of these two viruses and MWD at sites 1 and 4. At site 2, there was a statistically significant (P < 0.001) association only between PMWaV-3 and MWD. PMWaV-1 was the most commonly found of the four viruses and conversely PMWaV-5 was only occasionally found. Australian isolates of PMWaV-1, -2 and -3 were transmitted by the mealybug species Dysmicoccus brevipes.

Evolutionary relationships between bat coronaviruses and their hosts

Recent studies have suggested that bats are the natural reservoir of a range of coronaviruses (CoVs), and that rhinolophid bats harbor viruses closely related to the severe acute respiratory syndrome (SARS) CoV, which caused an outbreak of respiratory illness in humans during 2002-2003. We examined the evolutionary relationships between bat CoVs and their hosts by using sequence data of the virus RNA-dependent RNA polymerase gene and the bat cytochrome b gene. Phylogenetic analyses showed multiple incongruent associations between the phylogenies of rhinolophid bats and their CoVs, which suggested that host shifts have occurred in the recent evolutionary history of this group. These shifts may be due to either virus biologic traits or host behavioral traits. This finding has implications for the emergence of SARS and for the potential future emergence of SARS-CoVs or related viruses.

Abaca bunchy top virus, a new member of the genus Babuvirus (family Nanoviridae)

Two isolates of a novel babuvirus causing "bunchy top" symptoms were characterised, one from abaca (Musa textilis) from the Philippines and one from banana (Musa sp.) from Sarawak (Malaysia). The name abacá bunchy top virus (ABTV) is proposed. Both isolates have a genome of six circular DNA components, each ca. 1.0-1.1 kb, analogous to those of isolates of Banana bunchy top virus (BBTV). However, unlike BBTV, both ABTV isolates lack an internal ORF in DNA-R, and the ORF in DNA-U3 found in some BBTV isolates is also absent. In all phylogenetic analyses of nanovirid isolates, ABTV and BBTV fall in the same clade, but on separate branches. However, ABTV and BBTV isolates shared only 79-81% amino acid sequence identity for the putative coat protein and 54-76% overall nucleotide sequence identity across all components. Stem-loop and major common regions were present in ABTV, but there was less than 60% identity with the major common region of BBTV. ABTV and BBTV were also shown to be serologically distinct, with only two out of ten BBTV-specific monoclonal antibodies reacting with ABTV. The two ABTV isolates may represent distinct strains of the species as they are less closely related to each other than are isolates of the two geographic subgroups (Asian and South Pacific) of BBTV.

Characterization of disease resistance gene candidates of the nucleotide binding site (NBS) type from banana and correlation of a transcriptional polymorphism with resistance to Fusarium oxysporum f.sp. cubense race 4

Most plant disease resistance (R) genes encode proteins with a nucleotide binding site and leucine-rich repeat structure (NBS-LRR). In this study, degenerate primers were used to amplify genomic NBS-type sequences from wild banana (Musa acuminata ssp. malaccensis) plants resistant to the fungal pathogen Fusarium oxysporum formae specialis (f. sp.) cubense (FOC) race 4. Five different classes of NBS-type sequences were identified and designated as resistance gene candidates (RGCs). The deduced amino acid sequences of the RGCs revealed the presence of motifs characteristic of the majority of known plant NBS-LRR resistance genes. Structural and phylogenetic analyses grouped the banana RGCs within the non-TIR (homology to Toll/interleukin-1 receptors) subclass of NBS sequences. Southern hybridization showed that each banana RGC is present in low copy number. The expression of the RGCs was assessed by RT-PCR in leaf and root tissues of plants resistant or susceptible to FOC race 4. RGC1, 3 and 5 showed a constitutive expression profile in both resistant and susceptible plants whereas no expression was detected for RGC4. Interestingly, RGC2 expression was found to be associated only to FOC race 4 resistant lines. This finding could assist in the identification of a FOC race 4 resistance gene.

Two novel mastreviruses from chickpea (Cicer arietinum) in Australia

Two novel mastreviruses (genus Mastrevirus; family Geminiviridae), with proposed names chickpea chlorosis virus (CpCV) and chickpea redleaf virus, are described from chickpea (Cicer arietinum) from eastern Australia. The viruses have genomes of 2,582 and 2,605 nucleotides, respectively, and share similar features and organisation with typical dicot-infecting mastreviruses. Two distinct strains of CpCV were suggested by phylogenetic analysis. Additionally, a partial mastrevirus Rep sequence from turnip weed (Rapistrum rugosum) indicated the presence of a distinct strain of Tobacco yellow dwarf virus (TYDV). In phylogenetic analyses, isolates of Bean yellow dwarf virus, Chickpea chlorotic dwarf Pakistan virus and Chickpea chlorotic dwarf Sudan virus from southern and northern Africa and south-central and western Asia clustered separately from these three viruses from Australia. An Australian, eastern Asian, or south-eastern Asian origin for the novel mastreviruses and TYDV is discussed.

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.