Merging ancient and modern DNA: extinct seabird taxon rediscovered in the North Tasman Sea

Ancient DNA has revolutionized the way in which evolutionary biologists research both extinct and extant taxa, from the inference of evolutionary history to the resolution of taxonomy. Here, we present, to our knowledge, the first study to report the rediscovery of an ‘extinct’ avian taxon, the Tasman booby (Sula tasmani), using classical palaeontological data combined with ancient and modern DNA data. Contrary to earlier work, we show an overlap in size between fossil and modern birds in the North Tasman Sea (classified currently as S. tasmani and Sula dactylatra fullagari, respectively). In addition, we show that Holocene fossil birds have mitochondrial control region sequences that are identical to those found in modern birds. These results indicate that the Tasman booby is not an extinct taxon: S. dactylatra fullagari O''Brien & Davies, 1990 is therefore a junior synonym of Sula tasmani van Tets, Meredith, Fullagar & Davidson, 1988 and all North Tasman Sea boobies should be known as S. d. tasmani. In addition to reporting the rediscovery of an extinct avian taxon, our study highlights the need for researchers to be cognizant of multidisciplinary approaches to understanding taxonomy and past biodiversity.     

Optimization of Mitochondrial DNA-based Hybridization Assays to Diagnostics in Soil

Nucleic acid hybridization among root-knot nematode mitochondrial DNAs can be used to identify several Meloidogyne species. Research was initiated to optimize mitochondrial DNA-based molecular diagnostics for the demanding environments likely to be encountered in field isolates. DNA hybridization using reconstituted DNA-soil mixtures revealed a loss of assay sensitivity ranging from 34% to 92% with four agronomic soils tested. This problem was alleviated by the addition of exogenously added DNA. Variation in nematode egg lysis procedures also affected hybridization efficiency, with NaOC1 treatment most effective at disrupting Meloidogyne eggs. These optimized conditions permit detection of mtDNA released from one to five Meloidogyne eggs using standard nucleic acid hybridization procedures.     

A Cytoplasmic Suppressor of a Nuclear Mutation Affecting Mitochondrial Functions in Drosophila

Phenotypes relevant to oxidative phosphorylation (OXPHOS) in eukaryotes are jointly determined by nuclear and mitochondrial DNA (mtDNA). Thus, in humans, the variable clinical presentations of mitochondrial disease patients bearing the same primary mutation, whether in nuclear or mitochondrial DNA, have been attributed to putative genetic determinants carried in the “other” genome, though their identity and the molecular mechanism(s) by which they might act remain elusive. Here we demonstrate cytoplasmic suppression of the mitochondrial disease-like phenotype of the Drosophila melanogaster nuclear mutant tko^25t, which includes developmental delay, seizure sensitivity, and defective male courtship. The tko^25t strain carries a mutation in a mitoribosomal protein gene, causing OXPHOS deficiency due to defective intramitochondrial protein synthesis. Phenotypic suppression was associated with increased mtDNA copy number and increased mitochondrial biogenesis, as measured by the expression levels of porin voltage dependent anion channel and Spargel (PGC1α). Ubiquitous overexpression of Spargel in tko^25t flies phenocopied the suppressor, identifying it as a key mechanistic target thereof. Suppressor-strain mtDNAs differed from related nonsuppressor strain mtDNAs by several coding-region polymorphisms and by length and sequence variation in the noncoding region (NCR), in which the origin of mtDNA replication is located. Cytoplasm from four of five originally Wolbachia-infected strains showed the same suppressor effect, whereas that from neither of two uninfected strains did so, suggesting that the stress of chronic Wolbachia infection may provide evolutionary selection for improved mitochondrial fitness under metabolic stress. Our findings provide a paradigm for understanding the role of mtDNA genotype in human disease.     

Genetic analyses reveal independent domestication origins of Eurasian reindeer

Although there is little doubt that the domestication of mammals was instrumental for the modernization of human societies, even basic features of the path towards domestication remain largely unresolved for many species. Reindeer are considered to be in the early phase of domestication with wild and domestic herds still coexisting widely across Eurasia. This provides a unique model system for understanding how the early domestication process may have taken place. We analysed mitochondrial sequences and nuclear microsatellites in domestic and wild herds throughout Eurasia to address the origin of reindeer herding and domestication history. Our data demonstrate independent origins of domestic reindeer in Russia and Fennoscandia. This implies that the Saami people of Fennoscandia domesticated their own reindeer independently of the indigenous cultures in western Russia. We also found that augmentation of local reindeer herds by crossing with wild animals has been common. However, some wild reindeer populations have not contributed to the domestic gene pool, suggesting variation in domestication potential among populations. These differences may explain why geographically isolated indigenous groups have been able to make the technological shift from mobile hunting to large-scale reindeer pastoralism independently.     

Organellar genomes of the four-toothed moss,Tetraphis pellucida

Background

Mosses are the largest of the three extant clades of gametophyte-dominant land plants and remain poorly studied using comparative genomic methods. Major monophyletic moss lineages are characterised by different types of a spore dehiscence apparatus called the peristome, and the most important unsolved problem in higher-level moss systematics is the branching order of these peristomate clades. Organellar genome sequencing offers the potential to resolve this issue through the provision of both genomic structural characters and a greatly increased quantity of nucleotide substitution characters, as well as to elucidate organellar evolution in mosses. We publish and describe the chloroplast and mitochondrial genomes of Tetraphis pellucida, representative of the most phylogenetically intractable and morphologically isolated peristomate lineage.

Results

Assembly of reads from Illumina SBS and Pacific Biosciences RS sequencing reveals that the Tetraphis chloroplast genome comprises 127,489 bp and the mitochondrial genome 107,730 bp. Although genomic structures are similar to those of the small number of other known moss organellar genomes, the chloroplast lacks the petN gene (in common with Tortula ruralis) and the mitochondrion has only a non-functional pseudogenised remnant of nad7 (uniquely amongst known moss chondromes).

Conclusions

Structural genomic features exist with the potential to be informative for phylogenetic relationships amongst the peristomate moss lineages, and thus organellar genome sequences are urgently required for exemplars from other clades. The unique genomic and morphological features of Tetraphis confirm its importance for resolving one of the major questions in land plant phylogeny and for understanding the evolution of the peristome, a likely key innovation underlying the diversity of mosses. The functional loss of nad7 from the chondrome is now shown to have occurred independently in all three bryophyte clades as well as in the early-diverging tracheophyte Huperzia squarrosa.     

Coprolites as a source of information on the genome and diet of the cave hyena

We performed high-throughput sequencing of DNA from fossilized faeces to evaluate this material as a source of information on the genome and diet of Pleistocene carnivores. We analysed coprolites derived from the extinct cave hyena (Crocuta crocuta spelaea), and sequenced 90 million DNA fragments from two specimens. The DNA reads enabled a reconstruction of the cave hyena mitochondrial genome with up to a 158-fold coverage. This genome, and those sequenced from extant spotted (Crocuta crocuta) and striped (Hyaena hyaena) hyena specimens, allows for the establishment of a robust phylogeny that supports a close relationship between the cave and the spotted hyena. We also demonstrate that high-throughput sequencing yields data for cave hyena multi-copy and single-copy nuclear genes, and that about 50 per cent of the coprolite DNA can be ascribed to this species. Analysing the data for additional species to indicate the cave hyena diet, we retrieved abundant sequences for the red deer (Cervus elaphus), and characterized its mitochondrial genome with up to a 3.8-fold coverage. In conclusion, we have demonstrated the presence of abundant ancient DNA in the coprolites surveyed. Shotgun sequencing of this material yielded a wealth of DNA sequences for a Pleistocene carnivore and allowed unbiased identification of diet.     

Complete Mitochondrial Genome of Echinostoma hortense (Digenea: Echinostomatidae)

Echinostoma hortense (Digenea: Echinostomatidae) is one of the intestinal flukes with medical importance in humans. However, the mitochondrial (mt) genome of this fluke has not been known yet. The present study has determined the complete mt genome sequences of E. hortense and assessed the phylogenetic relationships with other digenean species for which the complete mt genome sequences are available in GenBank using concatenated amino acid sequences inferred from 12 protein-coding genes. The mt genome of E. hortense contained 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 non-coding region. The length of the mt genome of E. hortense was 14,994 bp, which was somewhat smaller than those of other trematode species. Phylogenetic analyses based on concatenated nucleotide sequence datasets for all 12 protein-coding genes using maximum parsimony (MP) method showed that E. hortense and Hypoderaeum conoideum gathered together, and they were closer to each other than to Fasciolidae and other echinostomatid trematodes. The availability of the complete mt genome sequences of E. hortense provides important genetic markers for diagnostics, population genetics, and evolutionary studies of digeneans.     

Isolates of Meloidogyne hapla with Distinct Mitochondrial Genomes

Because two conflicting reports of the structure of the Meloidogyne hapla mitochondrial genome exist, we compared the mitochondrial DNA (mtDNA) purified from two isolates of M. hapla: one from San Bernardino County in southern California (BRDO) and the other from England. The authenticity of the BRDO isolate in particular was confirmed by examination of morphological characters, isoenzyme analysis, and differential host range tests. Restriction analysis revealed that mtDNA from the BRDO and English isolates corresponded to only the structure first reported, although significant differences between the two isolates were apparent. Southern blots probed with cloned, cytochrome oxidase I (cox-l) DNA from Romanomermis culicivorax mtDNA confirmed that the analyzed DNA was of mitochondrial origin. Thus, M. hapla has at least two distinct but presumably related mitchondrial genomes, plus at least one very different structure. These data are discussed with reference to recent molecular diagnostic and phylogenetic analyses of Meloidogyne.     

Complete Mitochondrial Genome Sequence Data Provides Genetic Evidence That the Brown Dog Tick Rhipicephalus sanguineus (Acari: Ixodidae) Represents a Species Complex

Ticks are blood-sucking ectoparasites of great medical and veterinary significance that can transmit bacteria, protozoa, fungi and viruses, and cause a variety of human and animal diseases worldwide. In the present study, we sequenced the complete mitochondrial (mt) genome of Rhipicephalus sanguineus from China (RSC) and compared with that of R. sanguineus from USA (RSU). Nucleotide sequence difference in the full mt genome was 11.23% between RSC and RSU. For the 13 protein-coding genes, comparison revealed sequence divergences at both the nucleotide (9.34-15.65%) and amino acid (2.54-19.23%) levels between RSC and RSU. In addition, sequence comparison of the conserved mt cox1 and cytb genes among multiple individual R. sanguineus revealed substantial nucleotide differences between RSC and RSU but limited sequence variation within RSC. Phylogenetic analysis of ticks based on the amino acid sequence data of 13 protein-coding genes revealed that R. sanguineus from China and R. sanguineus from USA represent sister taxa (likely separate species). Taken together, the findings support the recently proposal that R. sanguineus tick may represents a species complex of at least two closely related species.     

Assessing Nuclear and Mitochondrial DNA Sequence Variation Within Steinernema (Rhabditida: Steinernematidae)

DNA sequence analysis was used to characterize the nuclear ribosomal DNA ITS1 region and a portion of the COII and 16S rDNA genes of the mitochondrial genome from Steinernema entomopathogenic nematodes. Nuclear ITS1 nucleotide divergence among seven Steinernema spp. ranged from 6 to 22%, and mtDNA divergence among five species ranged from 12 to 20%. No intraspecific variation was observed among three S. feltiae strains. Phylogenetic analysis of both nuclear and mitochondrial DNA sequences confirms the existing morphological relationships of several Steinernema species. Both the rDNA ITS1 and mtDNA sequences were useful for resolving relationships among Steinernema taxa.     

Elbow joint adductor moment arm as an indicator of forelimb posture in extinct quadrupedal tetrapods

Forelimb posture has been a controversial aspect of reconstructing locomotor behaviour in extinct quadrupedal tetrapods. This is partly owing to the qualitative and subjective nature of typical methods, which focus on bony articulations that are often ambiguous and unvalidated postural indicators. Here we outline a new, quantitatively based forelimb posture index that is applicable to a majority of extant tetrapods. By determining the degree of elbow joint adduction/abduction mobility in several tetrapods, the carpal flexor muscles were determined to also play a role as elbow adductors. Such adduction may play a major role during the stance phase in sprawling postures. This role is different from those of upright/sagittal and sloth-like creeping postures, which, respectively, depend more on elbow extensors and flexors. Our measurements of elbow muscle moment arms in 318 extant tetrapod skeletons (Lissamphibia, Synapsida and Reptilia: 33 major clades and 263 genera) revealed that sprawling, sagittal and creeping tetrapods, respectively, emphasize elbow adductor, extensor and flexor muscles. Furthermore, scansorial and non-scansorial taxa, respectively, emphasize flexors and extensors. Thus, forelimb postures of extinct tetrapods can be qualitatively classified based on our quantitative index. Using this method, we find that Triceratops (Ceratopsidae), Anhanguera (Pterosauria) and desmostylian mammals are categorized as upright/sagittally locomoting taxa.     

Ancient Mitochondrial DNA Analyses of Ascaris Eggs Discovered in Coprolites from Joseon Tomb

Analysis of ancient DNA (aDNA) extracted from Ascaris is very important for understanding the phylogenetic lineage of the parasite species. When aDNAs obtained from a Joseon tomb (SN2-19-1) coprolite in which Ascaris eggs were identified were amplified with primers for cytochrome b (cyt b) and 18S small subunit ribosomal RNA (18S rRNA) gene, the outcome exhibited Ascaris specific amplicon bands. By cloning, sequencing, and analysis of the amplified DNA, we obtained information valuable for comprehending genetic lineage of Ascaris prevalent among pre-modern Joseon peoples.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-606) contains supplementary material, which is available to authorized users.     

A Case of Diphyllobothrium nihonkaiense Infection as Confirmed by Mitochondrial COX1 Gene Sequence Analysis

Diphyllobothrium nihonkaiense has been reported in Korea as Diphyllobothrium latum because of their close morphologic resemblance. We have identified a human case of D. nihonkaiense infection using the mitochondrial cytochrome c oxidase subunit I (cox1) gene sequence analysis. On 18 February 2012, a patient who had consumed raw fish a month earlier visited our outpatient clinic with a long tapeworm parasite excreted in the feces. The body of the segmented worm was 2 m long and divided into the scolex (head) and proglottids. It was morphologically close to D. nihonkaiense and D. latum. The cox1 gene analysis showed 99.4% (340/342 bp) homology with D. nihonkaiense but only 91.8% (314/342 bp) homology with D. latum. The present study suggested that the Diphyllobothrium spp. infection in Korea should be analyzed with specific DNA sequence for an accurate species identification.     

Complete Mitochondrial Genome of Anoplocephala magna Solidifying the Species

The 2 species of the genus Anoplocephala (Anoplocephalidae), A. perfoliata and A. magna, are among the most important equine cestode parasites. However, there is little information about their differences at the molecular level. The present study revealed that the mitochondrial (mt) genome of A. magna was 13,759 bp in size and 700 bp shorter than that of A. perfoliata. The 2 species includes 2 rRNA, 22 tRNA, and 12 protein-coding genes each. The size of each of the 36 genes was the same as that of A. perfoliata, except for cox1, rrnL, trnC, trnS2(UCN), trnG, trnH, trnQ, and trnP. In the full mitochondrial genome, the sequence similarity was 87.1%. The divergence in the nucleotide and amino acid sequences of individual protein-coding genes ranged from 11.1% to 16% and 6.8% to 16.4%, respectively. The 2 noncoding regions of the mt genome of A. magna were 199 bp and 271 bp in length, while the equivalent regions in A. perfoliata were 875 bp and 276 bp, respectively. The results of this study support the proposal that A. magna and A. perfoliata are separate species, consistent with previous morphological analyses.     

Mitochondrial DNA as effective molecular markers for the genetic variation and phylogeny of the family Osteoglossidae

The present study examined the genetic variation of the family Osteoglossidae from different geographical locations based on the mitochondrial NADH dehydrogenase subunit 2 (ND2) and ATPase subunit 6 (ATPase6) genes; we then re-constructed the phylogenetic relationships using the two sequences in combination. The results showed that the partial sequences of mitochondrial ND2 and ATPase6 of the family Osteoglossidae were 813 bp and 669 bp, respectively. A total of 42 species-specific nucleotide positions of the family Osteoglossidae were found to be useful for molecular identification. The sequence variation showed greater differences (8.3% ~ 28.1% for the combined sequences, 8.3% ~ 26.7% for the ND2 gene, and 9.3% ~ 28.7% for the ATPase6 gene) among the different species of Osteoglossidae, and there was a significant association between the genetic difference and geographical location. Phylogenetic analyses using neighbor-joining, Bayesian inference, and maximum parsimony (MP) methods based on the combined sequences of the two genes were able to distinguish the different species and were in agreement with the existing taxonomy based on morphological characters and in association with the geographical distribution among seven species of the family Osteoglossidae.     

Making teeth to order: conserved genes reveal an ancient molecular pattern in paddlefish (Actinopterygii)

Ray-finned fishes (Actinopterygii) are the dominant vertebrate group today (+30 000 species, predominantly teleosts), with great morphological diversity, including their dentitions. How dental morphological variation evolved is best addressed by considering a range of taxa across actinopterygian phylogeny; here we examine the dentition of Polyodon spathula (American paddlefish), assigned to the basal group Acipenseriformes. Although teeth are present and functional in young individuals of Polyodon, they are completely absent in adults. Our current understanding of developmental genes operating in the dentition is primarily restricted to teleosts; we show that shh and bmp4, as highly conserved epithelial and mesenchymal genes for gnathostome tooth development, are similarly expressed at Polyodon tooth loci, thus extending this conserved developmental pattern within the Actinopterygii. These genes map spatio-temporal tooth initiation in Polyodon larvae and provide new data in both oral and pharyngeal tooth sites. Variation in cellular intensity of shh maps timing of tooth morphogenesis, revealing a second odontogenic wave as alternate sites within tooth rows, a dental pattern also present in more derived actinopterygians. Developmental timing for each tooth field in Polyodon follows a gradient, from rostral to caudal and ventral to dorsal, repeated during subsequent loss of teeth. The transitory Polyodon dentition is modified by cessation of tooth addition and loss. As such, Polyodon represents a basal actinopterygian model for the evolution of developmental novelty: initial conservation, followed by tooth loss, accommodating the adult trophic modification to filter-feeding.     

Mitochondrial DNA variation in the malaria vector Anopheles minimus across China,Thailand and Vietnam: evolutionary hypothesis,population structure and population history

The effects of Pleistocene environmental fluctuations on the distribution and diversity of organisms in Southeast Asia are much less well known than in Europe and North America. In these regions, the combination of palaeoenvironmental reconstruction and inferences about population history from genetic data has been very powerful. In Southeast Asia, mosquitoes are good candidates for the genetic approach, with the added benefit that understanding the relative contributions of historical and current processes to population structure can inform management of vector species. Genetic variation among populations of Anopheles minimus was examined using 144 mtDNA COII sequences from 23 sites in China, Thailand and Vietnam. Haplotype diversity was high, with two distinct lineages that have a sequence divergence of over 2% and exhibit different geographical distributions. We compare alternative hypotheses concerning the origin of this pattern. The observed data deviate from the expectations based on a single-panmictic population with or without growth, or a stable but spatially structured population. However, they can be readily accommodated by a model of past fragmentation into eastern and western refugia, followed by growth and range expansion. This is consistent with the palaeoenvironmental reconstructions currently available for the region.