Utility of COI, CAD and morphological data for resolving relationships within the genus Sarcophaga (sensu lato) (Diptera: Sarcophagidae) : a preliminary study

Currently there are ~3000 known species of Sarcophagidae (Diptera), which are classified into 173 genera in three subfamilies. Almost 25% of sarcophagids belong to the genus Sarcophaga (sensu lato) however little is known about the validity of, and relationships between the ~150 (or more) subgenera of Sarcophaga s.l. In this preliminary study, we evaluated the usefulness of three sources of data for resolving relationships between 35 species from 14 Sarcophaga s.l. subgenera: the mitochondrial COI barcode region, ~800. bp of the nuclear gene CAD, and 110 morphological characters. Bayesian, maximum likelihood (ML) and maximum parsimony (MP) analyses were performed on the combined dataset. Much of the tree was only supported by the Bayesian and ML analyses, with the MP tree poorly resolved. The genus Sarcophaga s.l. was resolved as monophyletic in both the Bayesian and ML analyses and strong support was obtained at the species-level. Notably, the only subgenus consistently resolved as monophyletic was Liopygia. The monophyly of and relationships between the remaining Sarcophaga s.l. subgenera sampled remain questionable. We suggest that future phylogenetic studies on the genus Sarcophaga s.l. use combined datasets for analyses. We also advocate the use of additional data and a range of inference strategies to assist with resolving relationships within Sarcophaga s.l.

Deeply divergent mitochondrial lineages reveal patterns of local endemism in chironomids of the Australian Wet Tropics

The Wet Tropics bioregion of north-eastern Australia has been subject to extensive fluctuations in climate throughout the late Pliocene and Pleistocene. Cycles of rainforest contraction and expansion of dry sclerophyll forest associated with such climatic fluctuations are postulated to have played a major role in driving geographical endemism in terrestrial rainforest taxa. Consequences for the distributions of aquatic organisms, however, are poorly understood.The Australian non-biting midge species Echinocladius martini Cranston (Diptera: Chironomidae), although restricted to cool, well-forested freshwater streams, has been considered to be able to disperse among populations located in isolated rainforest pockets during periods of sclerophyllous forest expansion, potentially limiting the effect of climatic fluctuations on patterns of endemism. In this study, mitochondrial COI and 16S data were analysed for E. martini collected from eight sites spanning theWet Tropics bioregion to assess the scale and extent of phylogeographic structure. Analyses of genetic structure showed several highly divergent cryptic lineages with restricted geographical distributions. Within one of the identified lineages, strong genetic structure implied that dispersal among proximate (<1 km apart) streams was extremely restricted. The results suggest that vicariant processes, most likely due to the systemic drying of the Australian continent during the Plio-Pleistocene, might have fragmented historical E. martini populations and, hence, promoted divergence in allopatry.

Genetic resources in wild and cultured stocks of the Asian mud crab, Scylla paramamosain

The mud crab (Scylla spp.) aquaculture industry has expanded rapidly in recent years in many countries in the Indo - West Pacific (IWP) region as an alternative to marine shrimp culture because of significant disease outbreaks and associated failures of many shrimp culture industries in the region. Currently, practices used to produce and manage breeding crabs in hatcheries may compromise levels of genetic diversity, ultimately compromising growth rates, disease resistance and stock productivity. Therefore, to avoid “genetic pollution” and its harmful effects and to promote further development of mud crab aquaculture and fisheries in a sustainable way, a greater understanding of the genetic attributes of wild and cultured mud crab stocks is required. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations for multiple purposes including for commercial production, recreation and conservation and to increase profitability and sustainability of newly emerging crab culture industries. Phylogeographic patterns and the genetic structure of Asian mud crab populations across the IWP were assessed to determine if they were concordant with those of other widespread taxa possessing pelagic larvae of relatively long duration. A 597 bp fragment of the mitochondrial DNA COI gene was amplified and screened for variation in a total of 297 individuals of S. paramamosain from six sampling sites across the species’ natural geographical distribution in the IWP and 36 unique haplotypes were identified. Haplotype diversities per site ranged from 0.516 to 0.879. Nucleotide diversity estimates among haplotypes were 0.11% – 0.48%. Maximum divergence observed among S. paramamosain samples was 1.533% and samples formed essentially a single monophyletic group as no obvious clades were related to geographical location of sites. A weak positive relationship was observed however, between genetic distance and geographical distance among sites. Microsatellite markers were then used to assess contemporary gene flow and population structure in Asian mud crab populations sampled across their natural distribution in the IWP. Eight microsatellite loci were screened in sampled S. paramamosain populations and all showed high allelic diversity at all loci in sampled populations. In total, 344 individuals were analysed, and 304 microsatellite alleles were found across the 8 loci. The mean number of alleles per locus at each site ranged from 20.75 to 28.25. Mean allelic richness per site varied from 17.2 to 18.9. All sites showed high levels of heterozygosity as average expected heterozygosities for all loci ranged from 0.917 – 0.953 while mean observed heterozygosity ranged from 0.916 – 0.959. Allele diversities were similar at all sites and across all loci. The results did not show any evidence for major differences in allele frequencies among sites and patterns of allele frequencies were very similar in all populations across all loci. Estimates of population differentiation (FST) were relatively low and most probably largely reflect intra – individual variation for very highly variable loci. Results from nDNA analysis showed evidence for only very limited population genetic structure among sampled S. paramamosain, and a positive and significant association for genetic and geographical distance among sample sites. Microsatellite markers were then employed to determine if adequate levels of genetic diversity has been captured in crab hatcheries for the breeding cycle. The results showed that all microsatellite loci were polymorphic in hatchery samples. Culture populations were in general, highly genetically depauperate, compared with comparable wild populations, with only 3 to 8 alleles recorded for the same loci set per population. In contrast, very high numbers of alleles per locus were found in reference wild S. paramamosain populations, which ranged from 18 to 46 alleles per locus per population. In general, this translates into a 3 to 10 fold decline in mean allelic richness per locus in all culture stocks compared with wild reference counterparts. Furthermore, most loci in all cultured S. paramamosain samples showed departures from HWE equilibrium. Allele frequencies were very different in culture samples from that present in comparable wild reference samples and this in particular, was reflected in a large decline in allele diversity per locus. The pattern observed was best explained by significant impacts of breeding practices employed in hatcheries rather than natural differentiation among wild populations used as the source of brood stock. Recognition of current problems and management strategies for the species both for the medium and long-term development of the new culture industry are discussed. The priority research to be undertaken over the medium term for S. paramamosain should be to close the life cycle fully to allow individuals to be bred on demand and their offspring equalised to control broodstock reproductive contributions. Establishing a broodstock register and pedigree mating system will be required before any selection program is implemented. This will ensure that sufficient genetic variation will be available to allow genetic gains to be sustainably achieved in a future stock improvement program. A fundamental starting point to improve hatchery practices will be to encourage farmers and hatchery managers to spawn more females in their hatcheries as it will increase background genetic diversity in culture stocks. Combining crablet cohorts from multiple hatcheries into a single cohort for supply to farmers or rotation of breeding females regularly in hatcheries will help to address immediate genetic diversity problems in culture stocks. Application of these results can provide benefits for managing wild and cultured Asian mud crab populations more efficiently. Over the long-term, application of data on genetic diversity in wild and cultured stocks of Asian mud crab will contribute to development of sustainable and productive culture industries in Vietnam and other countries in the IWP and can contribute towards conservation of wild genetic resources.

Spatial population genetic structure reveals strong natal site fidelity in Echinocladius martini (Diptera: Chironomidae) in northeast Queensland, Australia

1. A diverse array of patterns has been reported regarding the spatial extent of population genetic structure and effective dispersal in freshwater macroinvertebrates. In river systems, the movements of many taxa can be restricted to varying degrees by the natural stream channel hierarchy. 2. In this study, we sampled populations of the non-biting freshwater midge Echinocladius martini in the Paluma bioregion of tropical northeast Queensland to investigate fine scale patterns of within- and among-stream dispersal and gene flow within a purported historical refuge. We amplified a 639 bp fragment of mitochondrial COI and analysed genetic structure using pairwise ΦST, hierarchical AMOVA, Mantel tests and a parsimony network. Genetic variation was partitioned among stream sections using Streamtree to investigate the effect of potential instream dispersal barriers. 3. The data revealed strong natal site fidelity and significant differentiation among neighbouring, geographically proximate streams. We found evidence for only episodic adult flight among sites on separate stream reaches. Overall, however, our data suggested that both larval and adult dispersal was largely limited to within a stream channel. 4. This may arise from a combination of the high density of riparian vegetation physically restricting dispersal and from the joint effects of habitat stability and large population sizes. Together these may mitigate the requirement for movement among streams to avoid inbreeding and local extinction due to habitat change and may thus enable persistence of upstream populations in the absence of regular compensatory upstream flight. Taken together, these data suggest that dispersal of E. martini is highly restricted, to the scale of only a few kilometres, and hence occurs predominantly within the natal stream.

Molecular phylogeny supports the paraphyletic nature of the genus Trogoderma (Coleoptera: Dermestidae) collected in the Australasian ecozone

To date, a molecular phylogenetic approach has not been used to investigate the evolutionary structure of Trogoderma and closely related genera. Using two mitochondrial genes, Cytochrome Oxidase I and Cytochrome B, and the nuclear gene, 18S, the reported polyphyletic positioning of Trogoderma was examined. Paraphyly in Trogoderma was observed, with one Australian Trogoderma species reconciled as sister to all Dermestidae and the Anthrenocerus genus deeply nested within the Australian Trogoderma clade. In addition, time to most recent common ancestor for a number of Dermestidae was calculated. Based on these estimations, the Dermestidae origin exceeded 175 million years, placing the origins of this family in Pangaea.

Phylogeography of Echinocladius martini Cranston (Diptera: Chironomidae) in closed forest streams of eastern Australia

The eastern Australian rainforests have experienced several cycles of range contraction and expansion since the late Miocene that are closely correlated with global glaciation events. Together with ongoing aridification of the continent, this has resulted in current distributions of native closed forest that are highly fragmented along the east coast. Several closed forest endemic taxa exhibit patterns of population genetic structure that are congruent with historical isolation of populations in discrete refugia and reflect evolutionary histories dramatically affected by vicariance. Currently, limited data are available regarding the impact of these past climatic fluctuations on freshwater invertebrate taxa. The non-biting midge species Echinocladius martini Cranston is distributed along the east coast and inhabits predominantly montane streams in closed forest habitat. Phylogeographic structure in E. martini was resolved here at a continental scale by incorporating data from a previous pilot study and expanding the sampling design to encompass populations in the Wet Tropics of north-eastern Queensland, south-east Queensland, New South Wales and Victoria. Patterns of phylogeographic structure revealed several deeply divergent mitochondrial lineages from central and south-eastern Australia that were previously unrecognised and were geographically endemic to closed forest refugia. Estimated divergence times were congruent with late Miocene onset of rainforest contractions across the east coast of Australia. This suggested that dispersal and gene flow among E. martini populations isolated in refugia has been highly restricted historically. Moreover, these data imply, in contrast to existing preconceptions about freshwater invertebrates, that this taxon may be acutely susceptible to habitat fragmentation.

A molecular phylogeny for the Tribe Dacini (Diptera: Tephritidae) : systematic and biogeographic implications

With well over 700 species, the Tribe Dacini is one of the most species-rich clades within the dipteran family Tephritidae, the true fruit flies. Nearly all Dacini belong to one of two very large genera, Dacus Fabricius and Bactrocera Macquart. The distribution of the genera overlap in or around the Indian subcontinent, but the greatest diversity of Dacus is in Africa and the greatest diversity of Bactrocera is in south-east Asia and the Pacific. The monophyly of these two genera has not been rigorously established, with previous phylogenies only including a small number of species and always heavily biased to one genus over the other. Moreover, the subgeneric taxonomy within both genera is complex and the monophyly of many subgenera has not been explicitly tested. Previous hypotheses about the biogeography of the Dacini based on morphological reviews and current distributions of taxa have invoked an out-of-India hypothesis; however this has not been tested in a phylogenetic framework. We attempted to resolve these issues with a dated, molecular phylogeny of 125 Dacini species generated using 16S, COI, COII and white eye genes. The phylogeny shows that Bactrocera is not monophyletic, but rather consists of two major clades: Bactrocera s.s. and the ‘Zeugodacus group of subgenera’ (a recognised, but informal taxonomic grouping of 15 Bactrocera subgenera). This ‘Zeugodacus’ clade is the sister group to Dacus, not Bactrocera and, based on current distributions, split from Dacus before that genus moved into Africa. We recommend that taxonomic consideration be given to raising Zeugodacus to genus level. Supportive of predictions following from the out-of-India hypothesis, the first common ancestor of the Dacini arose in the mid-Cretaceous approximately 80 mya. Major divergence events occurred during the Indian rafting period and diversification of Bactrocera apparently did not begin until after India docked with Eurasia (50–35 mya). In contrast, diversification in Dacus, at approximately 65 mya, apparently began much earlier than predicted by the out-of-India hypothesis, suggesting that, if the Dacini arose on the Indian plate, then ancestral Dacus may have left the plate in the mid to late Cretaceous via the well documented India–Madagascar–Africa migration route. We conclude that the phylogeny does not disprove the predictions of an out-of-India hypothesis for the Dacini, although modification of the original hypothesis is required.

Population structure of Bactrocera dorsalis s.s., B. papayae and B. philippinensis (Diptera: Tephritidae) in southeast Asia: evidence for a single species hypothesis using mitochondrial DNA and wingshape data

Background Bactrocera dorsalis s.s. is a pestiferous tephritid fruit fly distributed from Pakistan to the Pacific, with the Thai/Malay peninsula its southern limit. Sister pest taxa, B. papayae and B. philippinensis, occur in the southeast Asian archipelago and the Philippines, respectively. The relationship among these species is unclear due to their high molecular and morphological similarity. This study analysed population structure of these three species within a southeast Asian biogeographical context to assess potential dispersal patterns and the validity of their current taxonomic status. Results Geometric morphometric results generated from 15 landmarks for wings of 169 flies revealed significant differences in wing shape between almost all sites following canonical variate analysis. For the combined data set there was a greater isolation-by-distance (IBD) effect under a ‘non-Euclidean’ scenario which used geographical distances within a biogeographical ‘Sundaland context’ (r2 = 0.772, P < 0.0001) as compared to a ‘Euclidean’ scenario for which direct geographic distances between sample sites was used (r2 = 0.217, P < 0.01). COI sequence data were obtained for 156 individuals and yielded 83 unique haplotypes with no correlation to current taxonomic designations via a minimum spanning network. BEAST analysis provided a root age and location of 540kya in northern Thailand, with migration of B. dorsalis s.l. into Malaysia 470kya and Sumatra 270kya. Two migration events into the Philippines are inferred. Sequence data revealed a weak but significant IBD effect under the ‘non-Euclidean’ scenario (r2 = 0.110, P < 0.05), with no historical migration evident between Taiwan and the Philippines. Results are consistent with those expected at the intra-specific level. Conclusions Bactrocera dorsalis s.s., B. papayae and B. philippinensis likely represent one species structured around the South China Sea, having migrated from northern Thailand into the southeast Asian archipelago and across into the Philippines. No migration is apparent between the Philippines and Taiwan. This information has implications for quarantine, trade and pest management.

A fourier-transform infrared study of CO2 and CO2/H2 interactions with caesium-doped copper catalysts

FTIR spectra are reported of CO2 and COi/Hi on a silica-supported caesium-doped copper catalyst. Adsorption of COj on a "caesium"/silica surface resulted in the formation of COj and complexed CO species. Exposure of CO2 to' a caesium-doped reduced copper catalyst produced not only these species but also two forms of adsorbed carboxylate giving bands at 1550, 1510, 1365 and 1345 cm"1. Reaction of carboxylate species with hydrogen at 388 K gave formate species on copper and caesium oxide in addition to methoxy groups associated with caesium oxide. Methoxy species were not detected on undoped copper catalyst suggesting that caesium may be a promoter for the methanol synthesis reaction. Methanol decomposition on a caesium-doped copper catalyst produced a small number of formate species on copper and caesium oxide. Methoxy groups on caesium oxide decomposed to CO and U.2, and subsequent reaction between CO and adsorbed oxygen resulted in carboxylate formation. Methoxy species located at interfacial sites appeared to exhibit unusual adsorption properties.

The evolution and biogeography of the austral horse fly tribe Scionini (Diptera: Tabanidae: Pangoniinae) inferred from multiple mitochondrial and nuclear genes

Phylogenetic relationships within the Tabanidae are largely unknown, despite their considerable medical and ecological importance. The first robust phylogenetic hypothesis for the horse fly tribe Scionini is provided, completing the systematic placement of all tribes in the subfamily Pangoniinae. The Scionini consists of seven mostly southern hemisphere genera distributed in Australia, New Guinea, New Zealand and South America. A 5757. bp alignment of 6 genes, including mitochondrial (COI and COII), ribosomal (28S) and nuclear (AATS and CAD regions 1, 3 and 4) genes, was analysed for 176 taxa using both Bayesian and maximum likelihood approaches. Results indicate the Scionini are strongly monophyletic, with the exclusion of the only northern hemisphere genus Goniops. The South American genera Fidena, Pityocera and Scione were strongly monophyletic, corresponding to current morphology-based classification schemes. The most widespread genus Scaptia was paraphyletic and formed nine strongly supported monophyletic clades, each corresponding to either the current subgenera or several previously synonymised genera that should be formally resurrected. Molecular results also reveal a newly recognised genus endemic to New Zealand, formerly placed within Scaptia. Divergence time estimation was employed to assess the global biogeographical patterns in the Pangoniinae. These analyses demonstrated that the Scionini are a typical Gondwanan group whose diversification was influenced by the fragmentation of that ancient land mass. Furthermore, results indicate that the Scionini most likely originated in Australia and subsequently radiated to New Zealand and South American by both long distance dispersal and vicariance. The phylogenetic framework of the Scionini provided herein will be valuable for taxonomic revisions of the Tabanidae.

A molecular phylogeny of the checkered beetles and a description of Epiclininae a new subfamily (Coleoptera: Cleroidea: Cleridae)

We provide the first molecular phylogeny of the clerid lineage (Coleoptera: Cleridae, Thanerocleridae) within the superfamily Cleroidea to examine the two most recently-proposed hypotheses of higher-level classification. Phylogenetic relationships of checkered beetles were inferred from approximately ~5,000nt of both nuclear and mitochondrial rDNA (28S, 16S, and 12S) and the mitochondrial protein-coding gene COI. A worldwide sample of ~70 genera representing almost a quarter of generic diversity of the clerid lineage was included and phylogenies were reconstructed using Bayesian and Maximum Likelihood approaches. Results support the monophyly of many proposed subfamilies but were not entirely congruent with either current classification system. The subfamilial relationships within the Cleridae are resolved with support for three main lineages. Tillinae are supported as the sister group to all other subfamilies within the Cleridae, whereas Thaneroclerinae, Korynetinae and a new subfamily formally described here, Epiclininae subf. n, form a sister group to Clerinae + Hydnocerinae.

The complete mitochondrial genome of the yellow coaster, Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini): sequence, gene organization and a unique tRNA translocation event

In this paper, the complete mitochondrial genome of Acraea issoria (Lepidoptera: Nymphalidae: Heliconiinae: Acraeini) is reported; a circular molecule of 15,245 bp in size. For A. issoria, genes are arranged in the same order and orientation as the complete sequenced mitochondrial genomes of the other lepidopteran species, except for the presence of an extra copy of tRNAIle(AUR)b in the control region. All protein-coding genes of A. issoria mitogenome start with a typical ATN codon and terminate in the common stop codon TAA, except that COI gene uses TTG as its initial codon and terminates in a single T residue. All tRNA genes possess the typical clover leaf secondary structure except for tRNASer(AGN), which has a simple loop with the absence of the DHU stem. The sequence, organization and other features including nucleotide composition and codon usage of this mitochondrial genome were also reported and compared with those of other sequenced lepidopterans mitochondrial genomes. There are some short microsatellite-like repeat regions (e.g., (TA)9, polyA and polyT) scattered in the control region, however, the conspicuous macro-repeats units commonly found in other insect species are absent.

The complete mitochondrial genome of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae)

We determined the nucleotide sequence of the mitochondrial genome (mtgenome) of Spilonota lechriaspis Meyrick (Lepidoptera: Tortricidae). The entire closed circular molecule is 15,368 bp and contains 37 genes with the typical gene complement and order for lepidopteran mtgenomes. All tRNAs except tRNASer(AGN) can be folded into the typical cloverleaf secondary structures. The protein-coding genes (PCGs) have typical mitochondrial start codons, with the exception of COI, which uses the unusual CGA one as is found in all other Lepidoptera sequenced to date. In addition, six of 13 PCGs harbor the incomplete termination codons, a single T. The A+T-rich region contains some conserved structures that are similar to those found in other lepidopteran mtgenomes, including a structure combining the motif 'ATAGA', a 19-bp poly(T) stretch and three microsatellite (AT)n elements which are part of larger 122+ bp macrorepeats. This is the first report of macrorepeats in a lepidopteran mtgenome.

Towards a phylogeny of the Tenebrionoidea (Coleoptera)

The phylogenetic relationships of the beetle superfamily Tenebrionoidea are investigated using the most comprehensive genetic data set compiled to date. With ∼34,000 described species in approximately 1250 genera and 28 families, Tenebrionoidea represent one of the most diverse and species-rich superfamilies of beetles. The interfamilial relationships of the Tenebrionoidea are poorly known; previous morphological and molecular phylogenies recovered few well-supported and often conflicting relationships between families. Here we present a molecular phylogeny of Tenebrionoidea based on genes commonly used to resolve family and superfamily-level phylogenies of beetles (18S, 28S, 16S, 12S, tRNA Val and COI). The alignment spanned over 6.5 KB of DNA sequence and over 300 tenebrionoid genera from 24 of the 28 families were sampled. Maximum Likelihood and Bayesian analysis could not resolve deeper level divergences within the superfamily and very few relationships between families were supported. Increasing gene coverage in the alignment by removing taxa with missing data did not improve clade support but when rogue taxa were removed increased resolution was recovered. Investigation of signal strength suggested conflicting phylogenetic signal was present in the standard genes used for beetle phylogenetics, even when rogue taxa were removed. Our study of Tenebrionoidea highlights that even with relatively comprehensive taxon sampling within a lineage, this standard set of genes is unable to resolve relationships within this superfamily.

Comparison of intraspecific genetic structure among related chironomids (Diptera) from New Zealand and Patagonia : disparity between potential and realized dispersal

Population genetic studies of freshwater invertebrate taxa in New Zealand and South America are currently few despite the geologically and climatically dynamic histories of these regions. The focus of our study was a comparison of the influence on realized dispersal of 2 closely related nonbiting midges (Chironomidae) of population fragmentation on these separated austral land masses. We used a 734-base pair (bp) fragment of cytochrome c oxidase subunit I (COI) to investigate intraspecific genetic structure in Naonella forsythi Boothroyd in New Zealand and Ferringtonia patagonica Edwards in Patagonia. We proposed hypotheses about their potential dispersal and, hence, expected patterns of genetic structure in these 2 species based on published patterns for the closely related Australian taxon Echinocladius martini Cranston. Genetic structure revealed for both N. forsythi and F. patagonica was characterized by several highly divergent (2.0–10.5%) lineages of late Miocene–Pliocene age within each taxon that were not geographically localized. Many were distributed widely. This pattern differed greatly from population structure in E. martini, which was typified by much greater endemicity of divergent genetic lineages. Nevertheless, diversification of lineages in all 3 taxa appeared to be temporally congruent with the onset of late Miocene glaciations in the southern hemisphere that may have driven fragmentation of suitable habitat, promoting isolation of populations and divergence in allopatry. We argue that differences in realized dispersal post-isolation may be the result of differing availability of suitable habitat in interglacial periods.