Identification of small juvenile scombrids from northwest tropical Australia using mitochondrial DNA cytochrome b sequences

Small juveniles of the nine species of scombrids in Australian waters are morphologically similar to one another and, consequently, difficult to identify to species level. We show that the sequence of the mitochondrial DNA cytochrome b gene region is a powerful tool for identification of these young fish. Using this method, we identified 50 juvenile scombrids collected from Exmouth Bay, Western Australia. Six species of scombrids were apparent in this sample of fish: narrow-barred Spanish mackerel (Scomberomorus commerson), Indian mackerel (Rastrelliger kanagurta), frigate tuna (Auxis thazard), bullet tuna (Auxis rochei), leaping bonito (Cybiosarda elegans), and kawakawa (Euthynnus affinis). The presence of Indian mackerel, frigate tuna, leaping bonito, and kawakawa is the first indication that coastal waters may be an important spawning habitat for these species, although offshore spawning may also occur. The occurrence of small juvenile S. commerson was predicted from the known spawning patterns of that species, but other mackerel species (Scomberomorus munroi, Scomberomorus queenslandicus, Scomberomorus semifasiciatus) likely to be spawning during the sampling period were not detected among the 50 small juveniles analyzed here.

Mitochondrial DNA supports the identification of two endangered river sharks (Glyphis glyphis and Glyphis garricki) across northern Australia

The river sharks (genus Glyphis) are a small group of poorly known sharks occurring in tropical rivers and estuarine waters across northern Australia, south-east Asia and the subcontinent. The taxonomy of the genus has long been unclear due to very few individuals having been caught and examined, resulting in a paucity of data regarding their distribution, biology and ecology. Only recently has attention focussed on the two Australian species, G. glyphis and G. garricki. This study is a result of a rare opportunity to collate the few samples that have been collected from these species and the bull shark Carcharhinus leucas, which shares an overlapping range. These samples were analysed using the DNA barcoding approach (cox1 mitochondrial gene), compared with six other species of carcharhinids and evaluated in light of the current taxonomic classification. Nine species-specific nucleotide differences were found between G. glyphis and G. garricki and no intra-specific variation provides strong support for the separation into distinct species. Significant differences were also observed at the inter-generic level, with Glyphis forming a distinct clade from Carcharhinus. This study provides the basis for future molecular studies required to better address conservation issues confronting G. glyphis and G. garricki in Australia.

Evidence for extensive population structure in the white-spotted eagle ray within the Indo-Pacific inferred from mitochondrial gene sequences

The white-spotted eagle ray Aetobatus narinari is a species complex that occurs circumglobally throughout warm-temperate waters. Aetobatus narinari is semi-pelagic and large (up to 300 cm disc width), suggesting high dispersal capabilities and gene flow on a wide spatial scale. Sequence data from two mitochondrial genes, cytochrome b (cytb) and NADH dehydrogenase subunit 4 (ND4), were used to determine the genetic variability within and among 18 sampling locations in the central Indo-Pacific biogeographical region. Populations in the Indo-Pacific were highly genetically structured with c. 70% of the total genetic variation found among three geographical regions (East China Sea, Southeast Asia and Australia). FST was 0.64 for cytb and 0.53 for ND4, with φST values being even larger, that is, 0.78 for cytb and 0.65 for ND4. This high-level genetic partitioning provides strong evidence against extensive gene flow in A. narinari. The degree of genetic population structuring in the Indo-Pacific was similar to that found on a global scale. Global FST was 0.63 for cytb and 0.57 for ND4, and global φST values were 0.94 for cytb and 0.82 for ND4. This suggests that the A. narinari complex may be more speciose than the two or three species proposed to date. Further sampling and genetic analyses are likely to uncover the ‘evolutionarily significant’ and ‘management’ units that are critical to determine the susceptibilities of individual populations to regional fishing pressures and to provide advice on management options. Network analyses showed a close genetic relationship between haplotypes from the central Indo-Pacific and South Africa, providing support for a proposed dispersal pathway from the possible centre of origin of the A. narinari species complex in the Indo-Pacific into the Atlantic Ocean.

Mitochondrial DNA diversity of Cleruchoides noackae (Hymenoptera: Mymaridae): a potential biological control agent for Thaumastocoris peregrinus (Hemiptera: Thaumastocoridae)

Carpintero and Dellap, (Hemiptera: Thaumastocoridae) is a native Australian sap-feeding insect that has become invasive and seriously damaging to commercially grown in the Southern Hemisphere. Lin and Huber (Hymenoptera: Mymaridae) was recently discovered as an egg parasitoid of the Thaumastocoridae in Australia. Mitochondrial DNA (mtDNA; cytochrome oxidase subunit I, COI) sequence diversity amongst 104 individuals from these native populations revealed 24 sequence haplotypes. The COI haplotypes of individuals collected from the Sydney and Southeast Queensland clustered in distinct groups, indicating limited spread of the insect between the regions. Individuals collected from Perth in Western Australia were represented by four COI haplotypes. Although this population is geographically more isolated from other populations, two COI haplotypes were identical to haplotypes found in the Sydney region. The results suggest that has recently been introduced into Perth, possibly from the Sydney area. The high mtDNA diversity and limited spread that is suggested for is in contrast to the lack of geographic associated mtDNA diversity and extensive spread of . If implemented as a biological control agent, this factor will need to be considered in collecting and releasing .

Hybridisation, paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies. © 2013.

The complete validated mitochondrial genome of the silver gemfish Rexea solandri (Cuvier, 1832) (Perciformes, Gempylidae)

The silver gemfish Rexea solandri is an important economic resource but vulnerable to overfishing in Australian waters. The complete mitochondrial genome sequence is described from 1.6 million reads obtained via next generation sequencing. The total length of the mitogenome is 16,350 bp comprising 2 rRNA, 13 protein-coding genes, 22 tRNA and 2 non-coding regions. The mitogenome sequence was validated against sequences of PCR fragments and BLAST queries of Genbank. Gene order was equivalent to that found in marine fishes.

A mitochondrial species identification assay for Australian blacktip sharks (Carcharhinus tilstoni, C. limbatus and C. amblyrhynchoides) using real-time PCR and high-resolution melt analysis

Tropical Australian shark fisheries target two morphologically indistinguishable blacktip sharks, the Australian blacktip (Carcharhinus tilstoni) and the common blacktip (C. limbatus). Their relative contributions to northern and eastern Australian coastal fisheries are unclear because of species identification difficulties. The two species differ in their number of precaudal vertebrae, which is difficult and time consuming to obtain in the field. But, the two species can be distinguished genetically with diagnostic mutations in their mitochondrial DNA ND4 gene. A third closely related sister species, the graceful shark C. amblyrhynchoides, can also be distinguished by species-specific mutations in this gene. DNA sequencing is an effective diagnostic tool, but is relatively expensive and time consuming. In contrast, real-time high-resolution melt (HRM) PCR assays are rapid and relatively inexpensive. These assays amplify regions of DNA with species-specific genetic mutations that result in PCR products with unique melt profiles. A real-time HRM PCR species-diagnostic assay (RT-HRM-PCR) has been developed based on the mtDNA ND4 gene for rapid typing of C. tilstoni, C. limbatus and C. amblyrhynchoides. The assay was developed using ND4 sequences from 66 C. tilstoni, 33. C. limbatus and five C. amblyrhynchoides collected from Indonesia and Australian states and territories; Western Australia, the Northern Territory, Queensland and New South Wales. The assay was shown to be 100% accurate on 160 unknown blacktip shark tissue samples by full mtDNA ND4 sequencing.

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.

DNA bar-coding reveals source and patterns of Thaumastocoris peregrinus invasions in South Africa and South America.

Thaumastocoris peregrinus is a recently introduced invertebrate pest of non-native Eucalyptus plantations in the Southern Hemisphere. It was first reported from South Africa in 2003 and in Argentina in 2005. Since then, populations have grown explosively and it has attained an almost ubiquitous distribution over several regions in South Africa on 26 Eucalyptus species. Here we address three key questions regarding this invasion, namely whether only one species has been introduced, whether there were single or multiple introductions into South Africa and South America and what the source of the introduction might have been. To answer these questions, bar-coding using mitochondrial DNA (COI) sequence diversity was used to characterise the populations of this insect from Australia, Argentina, Brazil, South Africa and Uruguay. Analyses revealed three cryptic species in Australia, of which only T. peregrinus is represented in South Africa and South America. Thaumastocoris peregrinus populations contained eight haplotypes, with a pairwise nucleotide distance of 0.2-0.9% from seventeen locations in Australia. Three of these haplotypes are shared with populations in South America and South Africa, but the latter regions do not share haplotypes. These data, together with the current distribution of the haplotypes and the known direction of original spread in these regions, suggest that at least three distinct introductions of the insect occurred in South Africa and South America before 2005. The two most common haplotypes in Sydney, one of which was also found in Brisbane, are shared with the non-native regions. Sydney populations of T. peregrinus, which have regularly reached outbreak levels in recent years, might thus have served as source of these three distinct introductions into other regions of the Southern Hemisphere.

The complete validated mitochondrial genome of the black marlin Istiompax indica (Cuvier, 1832)

Two complete mitochondrial genomes of the black marlin Istiompax indica were assembled from approximately 3.5 and 2.5 million reads produced by Ion Torrent next generation sequencing. The complete genomes were 16,531 bp and 16,532 bp in length consisting of 2 rRNA, 13 protein-coding genes, 22tRNA and 2 coding regions. They demonstrated a similar A + T base (52.6%) to other teleosts. Intraspecific sequence variation was 99.5% for three I. indica mitogenomes and 99.7% for X. gladius. A lower value (85%) was found for the I. platypterus mitogenomes from genbank and accredited to inadvertent inclusion of gene regions from a con-familial species in one record, highlighting the need for cautious downstream use of genbank data. © 2014 Informa UK Ltd.

Random Mating Between Two Widely Divergent Mitochondrial Lineages of Cryptolestes ferrugineus (Coleoptera: Laemophloeidae): A Test of Species Limits in a Phosphine-Resistant Stored Product Pest

Effective pest management relies on accurate delimitation of species and, beyond this, on accurate species identification. Mitochondrial COI sequences are useful for providing initial indications in delimiting species but, despite acknowledged limitations in the method, many studies involving COI sequences and species problems remain unresolved. Here we illustrate how such impasses can be resolved with microsatellite and nuclear sequence data, to assess more directly the amount of gene flow between divergent lineages. We use a population genetics approach to test for random mating between two 8 ± 2% divergent COI lineages of the rusty grain beetle, Cryptolestes ferrugineus (Stephens). This species has become strongly resistant to phosphine, a fumigant used worldwide for disinfesting grain. The possibility of cryptic species would have significant consequences for resistance management, especially if resistance was confined to one mitochondrial lineage. We find no evidence of restricted gene flow or nonrandom mating across the two COI lineages of these beetles, rather we hypothesize that historic population structure associated with early Pleistocene climate changes likely contributed to divergent lineages within this species.

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.

A universal test to determine the integrity of RNA, and its suitability for reverse transcription, in animal tissue laboratory specimens

Degradation of RNA in diagnostic specimens can cause false-negative test results and potential misdiagnosis when tests rely on the detection of specific RNA sequence. Current molecular methods of checking RNA integrity tend to be host species or group specific, necessitating libraries of primers and reaction conditions. The objective here was to develop a universal (multi-species) quality assurance tool for determining the integrity of RNA in animal tissues submitted to a laboratory for analyses. Ribosomal RNA (16S rRNA) transcribed from the mitochondrial 16S rDNA was used as template material for reverse transcription to cDNA and was amplified using polymerase chain reaction (PCR). As mitochondrial DNA has a high level of conservation, the primers used were shown to reverse transcribe and amplify RNA from every animal species tested. Deliberate degradation of rRNA template through temperature abuse of samples resulted in no reverse transcription and amplification. Samples spiked with viruses showed that single-stranded viral RNA and rRNA in the same sample degraded at similar rates, hence reverse transcription and PCR amplification of 16S rRNA could be used as a test of sample integrity and suitability for analysis that required the sample's RNA, including viral RNA. This test will be an invaluable quality assurance tool for determination of RNA integrity from tissue samples, thus avoiding erroneous test results that might occur if degraded target RNA is used unknowingly as template material for reverse transcription and subsequent PCR amplification.

Assessing the relationship between patch type and soil mites: A molecular approach

An urgent need exists for indicators of soil health and patch functionality in extensive rangelands that can be measured efficiently and at low cost. Soil mites are candidate indicators, but their identification and handling is so specialised and time-consuming that their inclusion in routine monitoring is unlikely. The aim of this study was to measure the relationship between patch type and mite assemblages using a conventional approach. An additional aim was to determine if a molecular approach traditionally used for soil microbes could be adapted for soil mites to overcome some of the bottlenecks associated with soil fauna diversity assessment. Soil mite species abundance and diversity were measured using conventional ecological methods in soil from patches with perennial grass and litter cover (PGL), and compared to soil from bare patches with annual grasses and/or litter cover (BAL). Soil mite assemblages were also assessed using a molecular method called terminal-restriction fragment length polymorphism (T-RFLP) analysis. The conventional data showed a relationship between patch type and mite assemblage. The Prostigmata and Oribatida were well represented in the PGL sites, particularly the Aphelacaridae (Oribatida). For T-RFLP analysis, the mite community was represented by a series of DNA fragment lengths that reflected mite sequence diversity. The T-RFLP data showed a distinct difference in the mite assemblage between the patch types. Where possible, T-RFLP peaks were matched to mite families using a reference 18S rDNA database, and the Aphelacaridae prevalent in the conventional samples at PGL sites were identified, as were prostigmatids and oribatids. We identified limits to the T-RFLP approach and this included an inability to distinguish some species whose DNA sequences were similar. Despite these limitations, the data still showed a clear difference between sites, and the molecular taxonomic inferences also compared well with the conventional ecological data. The results from this study indicated that the T-RFLP approach was effective in measuring mite assemblages in this system. The power of this technique lies in the fact that species diversity and abundance data can be obtained quickly because of the time taken to process hundreds of samples, from soil DNA extraction to data output on the gene analyser, can be as little as 4 days.

The genetic structure of Australasian green turtles (Chelonia mydas): Exploring the geographical scale of genetic exchange

Ecological and genetic studies of marine turtles generally support the hypothesis of natal homing, but leave open the question of the geographical scale of genetic exchange and the capacity of turtles to shift breeding sites. Here we combine analyses of mitochondrial DNA (mtDNA) variation and recapture data to assess the geographical scale of individual breeding populations and the distribution of such populations through Australasia. We conducted multiscale assessments of mtDNA variation among 714 samples from 27 green turtle rookeries and of adult female dispersal among nesting sites in eastern Australia. Many of these rookeries are on shelves that were flooded by rising sea levels less than 10 000 years (c. 450 generations) ago. Analyses of sequence variation among the mtDNA control region revealed 25 haplotypes, and their frequency distributions indicated 17 genetically distinct breeding stocks (Management Units) consisting either of individual rookeries or groups of rookeries in general that are separated by more than 500 km. The population structure inferred from mtDNA was consistent with the scale of movements observed in long-term mark-recapture studies of east Australian rookeries. Phylogenetic analysis of the haplotypes revealed five clades with significant partitioning of sequence diversity (Φ = 68.4) between Pacific Ocean and Southeast Asian/Indian Ocean rookeries. Isolation by distance was indicated for rookeries separated by up to 2000 km but explained only 12% of the genetic structure. The emerging general picture is one of dynamic population structure influenced by the capacity of females to relocate among proximal breeding sites, although this may be conditional on large population sizes as existed historically across this region.