Combined mitochondrial and nuclear DNA sequences resolve the interrelations of the major Australasian marsupial radiations

Australasian marsupials include three major radiations, the insectivorous/carnivorous Dasyuromorphia, the omnivorous bandicoots (Peramelemorphia), and the largely herbivorous diprotodontians. Morphologists have generally considered the bandicoots and diprotodontians to be closely related, most prominently because they are both syndactylous (with the 2nd and 3rd pedal digits being fused). Molecular studies have been unable to confirm or reject this Syndactyla hypothesis. Here we present new mitochondrial (mt) genomes from a spiny bandicoot (Echymipera rufescens) and two dasyurids, a fat-tailed dunnart (Sminthopsis crassicaudata) and a northern quoll (Dasyurus hallucatus). By comparing trees derived from pairwise base-frequency differences between taxa with standard (absolute, uncorrected) distance trees, we infer that composition bias among mt protein-coding and RNA sequences is sufficient to mislead tree reconstruction. This can explain incongruence between trees obtained from mt and nuclear data sets. However, after excluding major sources of compositional heterogeneity, both the “reduced-bias” mt and nuclear data sets clearly favor a bandicoot plus dasyuromorphian association, as well as a grouping of kangaroos and possums (Phalangeriformes) among diprotodontians. Notably, alternatives to these groupings could only be confidently rejected by combining the mt and nuclear data. Elsewhere on the tree, Dromiciops appears to be sister to the monophyletic Australasian marsupials, whereas the placement of the marsupial mole (Notoryctes) remains problematic. More generally, we contend that it is desirable to combine mt genome and nuclear sequences for inferring vertebrate phylogeny, but as separately modeled process partitions. This strategy depends on detecting and excluding (or accounting for) major sources of nonhistorical signal, such as from compositional nonstationarity.

Phylogenetics of small horseshoe bats from east Asia based on mitochondrial DNA sequence variation

We undertook analyses of mitochondrial DNA gene sequences and echolocation calls to resolve phylogenetic relationships among the related bat taxa Rhinolophus pusillus (sampled across China), R. monoceros (Taiwan), R. cornutus (main islands of Japan), and R. c. pumilus (Okinawa, Japan), Phylogenetic trees and genetic divergence analyses were constructed by combining new complete mitochondrial cytochrome-b gene sequences and partial mitochondrial control region sequences with published sequences. Our work showed that these 4 taxa formed monophyletic groups in the phylogenetic tree. However, low levels of sequence divergence among the taxa, together with similarities in body size and overlapping echolocation call frequencies, point to a lack of taxonomic distinctiveness. We therefore suggest that these taxa are better considered as geographical subspecies rather than distinct species, although this should not diminish the conservation importance of these island populations, which are important evolutionarily significant units. Based on our findings, we suggest that the similarities in body size and echolocation call frequency in these rhinolophids result from their recent common ancestry, whereas similarities in body size and call frequency with R. hipposideros of Europe are the result of convergent evolution.

Isolation and characterisation of novel microsatellite and mitochondrial DNA markers for the Eastern Water Dragon (Physignathus lesueurii)

Habitat fragmentation as a result of urbanisation is a growing problem for native lizard species. The Eastern Water Dragon (Physignathus lesueurii) is a social arboreal agamid lizard, native to Australia. This species represents an ideal model species to investigate the effect of urbanisation because of their prominent abundance in the urban landscape. Here we describe the isolation and characterisation of a novel set of 74 di-, tri-, and tetramicrosatellites from which 18 were selected and optimised into two multiplexes. The 18 microsatellites generated a total 148 alleles across the two populations. The number of alleles per locus varied from 2 to 18 alleles and measures of Ho and He varied from 0.395 to 0.877 and from 0.441 to 0.880, respectively. We also present primers for four novel mitochondrial DNA (mtDNA) markers. The combined length of the four mtDNA marker pairs was 2,528 bp which included 15 nucleotides changes. In comparison to threatened species, which are generally characterised by small population sizes, the Eastern Water Dragon represents an ideal model species to investigate the effect of urbanisation on their behavioural ecology and connectivity patterns among populations.

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.

Australian lungfish (Neoceratodus forsteri: Dipnoi) have low genetic variation at allozyme and mitochondrial DNA loci : a conservation alert?

Genetic variation at allozyme and mitochondrial DNA loci was investigated in the Australian lungfish, Neoceratodus forsteri Krefft 1870. Tissue samples for genetic analysis were taken non-lethally from 278 individuals representing two spatially distinct endemic populations (Mary and Burnett rivers), as well as one population thought to be derived from an anthropogenic translocation in the 1890's (Brisbane river). Two of 24 allozyme loci resolved from muscle tissue were polymorphic. Mitochondrial DNA nucleotide sequence diversity estimated across 2,235 base pairs in each of 40 individuals ranged between 0.000423 and 0.001470 per river. Low genetic variation at allozyme and mitochondrial loci could be attributed to population bottlenecks, possibly induced by Pleistocene aridity. Limited genetic differentiation was detected among rivers using nuclear and mitochondrial markers suggesting that admixture may have occurred between the endemic Mary and Burnett populations during periods of low sea level when the drainages may have converged before reaching the ocean. Genetic data was consistent with the explanation that lungfish were introduced to the Brisbane river from the Mary river. Further research using more variable genetic loci is needed before the conservation status of populations can be determined, particularly as anthropogenic demands on lungfish habitat are increasing. In the interim we recommend a management strategy aimed at conserving existing genetic variation within and between rivers.

Restoration of native folding of single-stranded DNA sequences through reverse mutations: an indication of a new epigenetic mechanism

We used in vivo (biological), in silico (computational structure prediction), and in vitro (model sequence folding) analyses of single-stranded DNA sequences to show that nucleic acid folding conservation is the selective principle behind a high-frequency single-nucleotide reversion observed in a three-nucleotide mutated motif of the Maize streak virus replication associated protein (Rep) gene. In silico and in vitro studies showed that the three-nucleotide mutation adversely affected Rep nucleic acid folding, and that the single-nucleotide reversion [C(601)A] restored wild-type-like folding. In vivo support came from infecting maize with mutant viruses: those with Rep genes containing nucleotide changes predicted to restore a wild-type-like fold [A(601)/G(601)] preferentially accumulated over those predicted to fold differently [C(601)/T(601)], which frequently reverted to A(601) and displaced the original population. We propose that the selection of native nucleic acid folding is an epigenetic effect, which might have broad implications in the evolution of plants and their viruses.

Mitochondrial genome acquisition restores respiratory function and tumorigenic potential of cancer cells without mitochondrial DNA

We report that tumor cells devoid of their mitochondrial genome (mtDNA) show delayed tumor growth and that tumor formation is associated with acquisition of mtDNA from host cells. This leads to partial recovery of mitochondrial function in cells derived from primary tumors grown from cells without mtDNA and a shorter lag in tumor growth. Cell lines from circulating tumor cells showed further recovery of mitochondrial respiration and an intermediate lag to tumor growth, while cells from lung metastases exhibited full restoration of respiratory function and no lag in tumor growth. Stepwise assembly of mitochondrial respiratory supercomplexes was correlated with acquisition of respiratory function. Our findings indicate horizontal transfer of mtDNA from host cells in the tumor microenvironment to tumor cells with compromised respiratory function to re-establish respiration and tumor-initiating efficacy. These results suggest a novel pathophysiological process for overcoming mtDNA damage and support the notion of high plasticity of malignant cells.

Computational analysis of epigenetic information in human DNA sequences

Over the last few years, investigations of human epigenetic profiles have identified key elements of change to be Histone Modifications, stable and heritable DNA methylation and Chromatin remodeling. These factors determine gene expression levels and characterise conditions leading to disease. In order to extract information embedded in long DNA sequences, data mining and pattern recognition tools are widely used, but efforts have been limited to date with respect to analyzing epigenetic changes, and their role as catalysts in disease onset. Useful insight, however, can be gained by investigation of associated dinucleotide distributions. The focus of this paper is to explore specific dinucleotides frequencies across defined regions within the human genome, and to identify new patterns between epigenetic mechanisms and DNA content. Signal processing methods, including Fourier and Wavelet Transformations, are employed and principal results are reported.

Mitochondrial DNA sequence analysis from multiple gene fragments reveals genetic heterogeneity of Crassostrea ariakensis in East Asia

The native Asian oyster, Crassostrea ariakensis is one of the most common and important Crassostrea species that occur naturally along the coast of East Asia. Molecular species diagnosis is a prerequisite for population genetic analysis of wild oyster populations because oyster species cannot be discriminated reliably using external morphological characters alone due to character ambiguity. To date there have been few phylogeographic studies of natural edible oyster populations in East Asia, in particular this is true of the common species in Korea C. ariakensis. We therefore assessed the levels and patterns of molecular genetic variation in East Asian wild populations of C. ariakensis from Korea, Japan, and China using DNA sequence analysis of five concatenated mtDNA regions namely; 16S rRNA, cytochrome oxidase I, cytochrome oxidase II, cytochrome oxidase III, and cytochrome b. Two divergent C. ariakensis clades were identified between southern China and remaining sites from the northern region. In addition, hierarchical AMOVA and pairwise UST analyses showed that genetic diversity was discontinuous among wild populations of C. ariakensis in East Asia. Biogeographical and historical sea level changes are discussed as potential factors that may have influenced the genetic heterogeneity of wild C. ariakensis stocks across this region.

Arsenic trioxide promotes mitochondrial DNA mutation and cell apoptosis in primary APL cells and NB4 cell line

This study aimed to investigate the effects of arsenic trioxide (As2O3) on the mitochondrial DNA (mtDNA) of acute promyelocytic leukemia (APL) cells. The NB4 cell line was treated with 2.0 μmol/L As2O3in vitro, and the primary APL cells were treated with 2.0 μmol/L As2O3in vitro and 0.16 mg kg-1 d-1 As2O3in vivo. The mitochondrial DNA of all the cells above was amplified by PCR, directly sequenced and analyzed by Sequence Navigatore and Factura software. The apoptosis rates were assayed by flow cytometry. Mitochondrial DNA mutation in the D-loop region was found in NB4 and APL cells before As2O3 use, but the mutation spots were remarkably increased after As2O3 treatment, which was positively correlated to the rates of cellular apoptosis, the correlation coefficient: rNB4-As2O3=0.973818, and rAPL-As2O3=0.934703. The mutation types include transition, transversion, codon insertion or deletion, and the mutation spots in all samples were not constant and regular. It is revealed that As2O3 aggravates mtDNA mutation in the D-loop region of acute promyelocytic leukemia cells both in vitro and in vivo. Mitochondrial DNA might be one of the targets of As2O3 in APL treatment.

DNA sequences and austral taxa indicate generic synonymy of Paratrichocladius Santos-Abreu with Cricotopus Wulp (Diptera: Chironomidae)

In the century since the description of the orthoclad genus Paratrichocladius Santos-Abreu (Diptera: Chironomidae), separation in any life stage from the cosmopolitan, diverse Cricotopus Wulp has been problematic. Molecular analysis reveals the presence of two species in Australia that conform in morphology to Paratrichocladius and which form a well-supported clade including Paratrichocladius micans (Kieffer) from Africa and a distinct southern African larva. This clade clusters with taxa allied with Cricotopus albitibia (Walker), in turn nested within all other sampled Australian Cricotopus. Relevant nodes strongly support Cricotopus as nonmonophyletic without inclusion of Paratrichocladius. We synonymize Paratrichocladius with Cricotopus syn.n, treating Paratrichocladius as a subgenus. Cricotopus (Paratrichocladius) australiensis Cranston sp.n. is described for Trichocladius pluriserialis Freeman from Australia, which is not the same species under that name in New Zealand. Cricotopus (Paratrichocladius) bifenestrus Cranston sp.n. from Australia is described, also in all life stages. The many new combinations, listed in an Appendix, include three replacement names for new secondary homonyms, namely: Cricotopus (Paratrichocladius) sinobicinctus Cranston & Krosch nom.n. for Paratrichocladius bicinctus Fu, Sæther & Wang, Cricotopus draysoni Cranston & Krosch nom.n. for Cricotopus brevicornis Drayson, Krosch & Cranston, and Cricotopus (Paratrichocladius) sikhotealinus Makarchenko & Makarchenko nom.n. for Cricotopus orientalis Kieffer. We conclude with comments on wider issues in the taxonomy of Paratrichocladius, especially concerning New Zealand species.

Seasonal occurrence and distribution of Bryde's whales in the Hauraki Gulf, New Zealand

The Hauraki Gulf is a large, shallow embayment located north of Auckland City (36°51′S, 174°46′E), New Zealand. Bryde's whales (Balaenoptera edeni) are the most frequently observed balaenopterid in these waters. To assess the use of the Hauraki Gulf for this species, we examined the occurrence and distribution in relation to environmental parameters. Data were collected from a platform of opportunity during 674 daily surveys between March 2003 and February 2006. A total of 760 observations of Bryde's whales were recorded throughout the study period during 371 surveys. The number of Bryde's whales sighted/day was highest in winter, coinciding with the coolest median sea-surface temperature (14.6°C). Bryde's whales were recorded throughout the Hauraki Gulf in water depths ranging from 12.1–59.8 m (mean = 42.3, SD = 5.1). Cow–calf pairs were most frequently observed during the austral autumn in water depths of 29.9–53.9 m (mean = 40.8, SD = 5.2). Data from this study suggest Bryde's whales in the Hauraki Gulf exhibit a mix of both “inshore” and “offshore” characteristics from the Bryde's whales examined off the coast of South Africa. Based on complete mitochondrial DNA sequences, Sasaki et al. (2006) recognized two sister species of Bryde's whales: Balaenoptera brydei and B. edeni, with the latter including small-type, more coastal Bryde's whales from Japan, Hong Kong, and Australia. Their samples and samples in previous analyses of small-type whales, all originated from eastern and southeastern Asia. These authors did not include the forms of Bryde's whales that occur in other regions, e.g., in the Pacific off Peru (Valdivia et al. 1981), in the Atlantic off Brazil (Best 1977) and in the western Indian Ocean off South Africa (Best 1977). Recent genetic analysis using mtDNA from the “inshore” and “offshore” forms from South Africa confirms the offshore form is B. brydei, and establishes that the inshore form is more closely related to B. brydei than to B. edeni (Penry 2010). These different forms do vary considerably in their habitat use and ecology (refer to Table 1 for a detailed comparison between the South African inshore and offshore forms, as described by Best (1967, 1977) and the Bryde's whales from New Zealand (Wiseman 2008). Recent genetic analysis on the Bryde's whales in the Hauraki Gulf suggests they are B. brydei (Wiseman 2008). However, pending resolution of the uncertainty within and between species of this genus, we follow the Society of Marine Mammal's committee on taxonomy, who state that B. edeni applies to all Bryde's whales.

Bird evolution: testing the Metaves clade with six new mitochondrial genomes

Background Evolutionary biologists are often misled by convergence of morphology and this has been common in the study of bird evolution. However, the use of molecular data sets have their own problems and phylogenies based on short DNA sequences have the potential to mislead us too. The relationships among clades and timing of the evolution of modern birds (Neoaves) has not yet been well resolved. Evidence of convergence of morphology remain controversial. With six new bird mitochondrial genomes (hummingbird, swift, kagu, rail, flamingo and grebe) we test the proposed Metaves/Coronaves division within Neoaves and the parallel radiations in this primary avian clade. Results Our mitochondrial trees did not return the Metaves clade that had been proposed based on one nuclear intron sequence. We suggest that the high number of indels within the seventh intron of the β-fibrinogen gene at this phylogenetic level, which left a dataset with not a single site across the alignment shared by all taxa, resulted in artifacts during analysis. With respect to the overall avian tree, we find the flamingo and grebe are sister taxa and basal to the shorebirds (Charadriiformes). Using a novel site-stripping technique for noise-reduction we found this relationship to be stable. The hummingbird/swift clade is outside the large and very diverse group of raptors, shore and sea birds. Unexpectedly the kagu is not closely related to the rail in our analysis, but because neither the kagu nor the rail have close affinity to any taxa within this dataset of 41 birds, their placement is not yet resolved. Conclusion Our phylogenetic hypothesis based on 41 avian mitochondrial genomes (13,229 bp) rejects monophyly of seven Metaves species and we therefore conclude that the members of Metaves do not share a common evolutionary history within the Neoaves.

The complete mitochondrial genome of the eastern grey kangaroo (Macropus giganteus)

We present the complete mitochondrial genome (accession number: LK995454) of an iconic Australian species, the eastern grey kangaroo (Macropus giganteus). The mitogenomic organization is consistent with other marsupials, encoding 13 protein-coding genes, 22 tRNA genes, 2 ribosomal RNA genes, an origin of light strand replication and a control region or Dloop. No repetitive sequences were detected in the control region. The M. giganteus mitogenome exemplifies a combination of tRNA gene order and structural peculiarities that appear to be unique to marsupials. We present a maximum likelihood phylogeny based on complete mitochondrial protein and RNA coding sequences that confirms the phylogenetic position of the grey kangaroo among macropodids.