Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Numerous insertions of mitochondrial DNA in the genome of the northern mole vole, Ellobius talpinus

BACKGROUND

Ellobius talpinus is a subterranean rodent representing an attractive model in population ecology studies due to its highly special lifestyle and sociality. In such studies, mitochondrial DNA (mtDNA) is widely used. However, if nuclear copies of mtDNA, aka NUMTs, are present, they may co-amplify with the target mtDNA fragment, generating misleading results. The aim of this study was to determine whether NUMTs are present in E. talpinus.

METHODS AND RESULTS

PCR amplification of the putative mtDNA CytB-D-loop fragment using 'universal' primers from 56 E. talpinus samples produced multiple double peaks in 90% of the sequencing chromatograms. To reveal NUMTs, molecular cloning and sequencing of PCR products of three specimens was conducted, followed by phylogenetic analysis. The pseudogene nature of three out of the seven detected haplotypes was confirmed by their basal positions in relation to other Ellobius haplotypes in the phylogenetic tree. Additionally, 'haplotype B' was basal in relation to other E. talpinus haplotypes and found present in very distant sampling sites. BLASTN search revealed 195 NUMTs in the E. talpinus nuclear genome, including fragments of all four PCR amplified pseudogenes. Although the majority of the NUMTs studied were short, the entire mtDNA had copies in the nuclear genome. The most numerous NUMTs were found for rrnL, COXI, and D-loop.

CONCLUSIONS

Numerous NUMTs are present in E. talpinus and can be difficult to discriminate against mtDNA sequences. Thus, in future population or phylogenetic studies in E. talpinus, the possibility of cryptic NUMTs amplification should always be taken into account.

Morphological and Genetics Support for a Hitherto Undescribed Spotted Cat Species (Genus Leopardus; Felidae, Carnivora) from the Southern Colombian Andes

In 1989, a skin of a small spotted cat, from the Galeras Volcano in southern Colombia (Nariño Department), was donated to the Instituto Alexander von Humboldt (identification, ID 5857) at Villa de Leyva (Boyacá Department, Colombia). Although originally classified as Leopardus tigrinus, its distinctiveness merits a new taxonomic designation. The skin is distinct from all known L. tigrinus holotypes as well as from other Leopardus species. Analysis of the complete mitochondrial genomes from 44 felid specimens (including 18 L. tigrinus and all the current known species of the genus Leopardus), the mtND5 gene from 84 felid specimens (including 30 L. tigrinus and all the species of the genus Leopardus), and six nuclear DNA microsatellites (113 felid specimens of all the current known species of the genus Leopardus) indicate that this specimen does not belong to any previously recognized Leopardus taxon. The mtND5 gene suggests this new lineage (the Nariño cat as we name it) is a sister taxon of Leopardus colocola. The mitogenomic and nuclear DNA microsatellite analyses suggest that this new lineage is the sister taxon to a clade formed by Central American and trans-Andean L. tigrinus + (Leopardus geoffroyi + Leopardus guigna). The temporal split between the ancestor of this new possible species and the most recent ancestor within Leopardus was dated to 1.2-1.9 million years ago. We consider that this new unique lineage is a new species, and we propose the scientific name Leopardus narinensis.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Comparative mitogenome phylogeography of two anteater genera (Tamandua and Myrmecophaga; Myrmecophagidae, Xenarthra): Evidence of discrepant evolutionary traits

The species within Xenarthra (sloths, anteaters, and armadillos) are quintessential South American mammals. Of the three groups, Vermilingua (anteaters) contains the fewest extant and paleontological species. Here, we sampled and sequenced the entire mitochondrial genomes (mitogenomes) of two Tamandua species (Tamandua tetradactyla and Tamandua mexicana) (n=74) from Central and South America, as well as Myrmecophaga tridactyla (n=41) from South America. Within Tamandua, we detected three different haplogroups. The oldest (THI) contained many specimens with the T. tetradactyla morphotype (but also several with the T. mexicana morphotype) and originated in southeastern South America (currently Uruguay) before moving towards northern South America, where the THII haplogroup originated. THII primarily contained specimens with the T. mexicana morphotype (but also several with the T. tetradactyla morphotype) and was distributed in Central America, Colombia, and Ecuador. THI and THII yielded a genetic distance of 4%. THII originated in either northern South America or "in situ" in Central America with haplogroup THIII, which consisted of ~50% T. mexicana and 50% T. tetradactyla phenotypes. THIII was mostly located in the same areas as THII, i.e., Central America, Ecuador, and Colombia, though mainly in the latter. The three haplogroups overlapped in Colombia and Ecuador. Thus, T. tetradactyla and T. mexicana were not reciprocally monophyletic. For this reason, we considered that a unique species of Tamandua likely exists, i.e., T. tetradactyla. In contrast to Tamandua, M. tridactyla did not show different morphotypes throughout its geographical range in the Neotropics. However, two very divergent genetic haplogroups (MHI and MHII), with a genetic distance of ~10%, were detected. The basal haplogroup, MHI, originated in northwestern South America, whereas the more geographically derived haplogroup, MHII, overlapped with MHI, but also expanded into central and southern South America. Thus, Tamandua migrated from south to north whereas Myrmecophaga migrated from north to south. Our results also showed that temporal mitochondrial diversification for Tamandua began during the Late Pliocene and Upper Pleistocene, but for Myrmecophaga began during the Late Miocene. Furthermore, both taxa showed elevated levels of mitochondrial genetic diversity. Tamandua showed more evidence of female population expansion than Myrmecophaga. Tamandua experienced population expansion ~0.6-0.17 million years ago (Mya), whereas Myrmecophaga showed possible population expansion ~0.3-0.2 Mya. However, both taxa experienced a conspicuous female decline in the last 10 000-20 000 years. Our results also showed little spatial genetic structure for both taxa. However, several analyses revealed higher spatial structure in Tamandua than in Myrmecophaga. Therefore, Tamandua and Myrmecophaga were not subjected to the same biogeographical, geological, or climatological events in shaping their genetic structures.

Phylogeny of the supertribe Nebriitae (Coleoptera, Carabidae) based on analyses of DNA sequence data

The phylogeny of the carabid beetle supertribe Nebriitae is inferred from analyses of DNA sequence data from eight gene fragments including one nuclear ribosomal gene (28S), four nuclear-protein coding genes (CAD, topoisomerase 1, PEPCK, and wingless), and three mitochondrial gene fragments (16S + tRNA-Leu + ND1, COI ("barcode" region) and COI ("Pat/Jer" region)). Our taxon sample included 264 exemplars representing 241 species and subspecies (25% of the known nebriite fauna), 39 of 41 currently accepted genera and subgenera (all except Notiokasis and Archileistobrius), and eight outgroup taxa. Separate maximum likelihood (ML) analyses of individual genes, combined ML analyses of nuclear, nuclear protein-coding, and mitochondrial genes, and combined ML and Bayesian analyses of the eight-gene-fragment matrix resulted in a well-resolved phylogeny of the supertribe, with most nodes in the tree strongly supported. Within Nebriitae, 167 internal nodes of the tree (out of the maximum possible 255) are supported by maximum-likelihood bootstrap values of 90% or more. The tribes Notiophilini, Opisthiini, Pelophilini, and Nebriini are well supported as monophyletic but relationships among these are not well resolved. Nippononebria is a distinct genus more closely related to Leistus than Nebria. Archastes, Oreonebria, Spelaeonebria, and Eurynebria, previously treated as distinct genera by some authors, are all nested within a monophyletic genus Nebria. Within Nebria, four major clades are recognized: (1) the Oreonebria Series, including eight subgenera arrayed in two subgeneric complexes (the Eonebria and Oreonebria Complexes); (2) the Nebriola Series, including only subgenus Nebriola; (3) the Nebria Series, including ten subgenera arrayed in two subgeneric complexes, the Boreonebria and Nebria Complexes, with the latter further subdivided into three subgeneric subcomplexes (the Nebria, Epinebriola, and Eunebria Subcomplexes)); and (4) the Catonebria Series, including seven subgenera arrayed in two subgeneric complexes (the Reductonebria and Catonebria Complexes). A strong concordance of biogeography with the inferred phylogeny is noted and some evident vicariance patterns are highlighted. A revised classification, mainly within the Nebriini, is proposed to reflect the inferred phylogeny. Three genus-group taxa (Nippononebria, Vancouveria and Archastes) are given revised status and seven are recognized as new synonymies (Nebriorites Jeannel, 1941 and Marggia Huber, 2014 = Oreonebria Daniel, 1903; Pseudonebriola Ledoux & Roux, 1989 = Boreonebria Jeannel, 1937; Patrobonebria Bänninger, 1923, Paranebria Jeannel, 1937 and Barbonebriola Huber & Schmidt, 2017 = Epinebriola Daniel & Daniel, 1904; and Asionebria Shilenkov, 1982 = Psilonebria Andrewes, 1923). Six new subgenera are proposed and described for newly recognized clades: Parepinebriola Kavanaugh subgen. nov. (type species: Nebria delicata Huber & Schmidt, 2017), Insulanebria Kavanaugh subgen. nov. (type species: Nebria carbonaria Eschscholtz, 1829), Erwinebria Kavanaugh subgen. nov. (type species Nebria sahlbergii Fischer von Waldheim, 1828), Nivalonebria Kavanaugh subgen. nov. (type species: Nebria paradisi Darlington, 1931), Neaptenonebria Kavanaugh subgen. nov. (type species: Nebria ovipennis LeConte, 1878), and Palaptenonebria Kavanaugh subgen. nov. (type species: Nebria mellyi Gebler, 1847). Future efforts to better understand relationships within the supertribe should aim to expand the taxon sampling of DNA sequence data, particularly within subgenera Leistus and Evanoleistus of genus Leistus and the Nebria Complex of genus Nebria.

A preliminary phylogeny and review of the genus Tasmanitachoides, with descriptions of two new species (Coleoptera, Carabidae, Bembidarenini)

The genus Tasmanitachoides Erwin, a genus of very small carabid beetle endemic to Australia, is reviewed. Although uncommon in collections, they can be abundant and diverse on banks of fine gravel or coarse sand next to bodies of fresh water; samples from southeastern Australia suggest numerous undescribed species. An initial phylogenetic hypothesis for the genus is presented, including 19 of the 32 known species. The inferred phylogeny, based upon one mitochondrial and four nuclear genes, shows the kingi group to be sister to remaining Tasmanitachoides, with the wattsensis group and T. lutus (Darlington) also being phylogenetically isolated. Two new species are described: T. baehri sp. nov., from the Australian Capital Territory, is a member of the kingi group; T. erwini sp. nov., from Tasmania, is a member of the wattsensis group. Identification tools for described and some undescribed species are presented, including photographs of all known species.

Profile hidden Markov model sequence analysis can help remove putative pseudogenes from DNA barcoding and metabarcoding datasets

BACKGROUND

Pseudogenes are non-functional copies of protein coding genes that typically follow a different molecular evolutionary path as compared to functional genes. The inclusion of pseudogene sequences in DNA barcoding and metabarcoding analysis can lead to misleading results. None of the most widely used bioinformatic pipelines used to process marker gene (metabarcode) high throughput sequencing data specifically accounts for the presence of pseudogenes in protein-coding marker genes. The purpose of this study is to develop a method to screen for nuclear mitochondrial DNA segments (nuMTs) in large COI datasets. We do this by: (1) describing gene and nuMT characteristics from an artificial COI barcode dataset, (2) show the impact of two different pseudogene removal methods on perturbed community datasets with simulated nuMTs, and (3) incorporate a pseudogene filtering step in a bioinformatic pipeline that can be used to process Illumina paired-end COI metabarcode sequences. Open reading frame length and sequence bit scores from hidden Markov model (HMM) profile analysis were used to detect pseudogenes.

RESULTS

Our simulations showed that it was more difficult to identify nuMTs from shorter amplicon sequences such as those typically used in metabarcoding compared with full length DNA barcodes that are used in the construction of barcode libraries. It was also more difficult to identify nuMTs in datasets where there is a high percentage of nuMTs. Existing bioinformatic pipelines used to process metabarcode sequences already remove some nuMTs, especially in the rare sequence removal step, but the addition of a pseudogene filtering step can remove up to 5% of sequences even when other filtering steps are in place.

CONCLUSIONS

Open reading frame length filtering alone or combined with hidden Markov model profile analysis can be used to effectively screen out apparent pseudogenes from large datasets. There is more to learn from COI nuMTs such as their frequency in DNA barcoding and metabarcoding studies, their taxonomic distribution, and evolution. Thus, we encourage the submission of verified COI nuMTs to public databases to facilitate future studies.

Shards, sequences, and shorelines: two new species of Bembidion from North America (Coleoptera, Carabidae)

Two new species of Bembidion are described from river shores in North America. One, Bembidionmimbressp. nov., from the Gila River watershed in the lands of the Mimbres culture in New Mexico and Arizona, is closely related to the widespread Bembidionlevigatum. DNA sequences from several linkage groups and morphology provide evidence of the distinctiveness of B.mimbres. The second, Bembidioncorgenomasp. nov., has been the subject of recent genomic and transcriptomic studies. It belongs in the Bembidiontransversale subgroup, and occurs from California north to British Columbia, east to Montana and Nevada. The B.transversale subgroup as a whole is reviewed, and morphological characters that distinguish B.corgenoma from the similar and sympatric B.transversale and B.erosum are described and illustrated. DNA sequences of these three species show no consistent differences in 28S, COI, CAD, and Topoisomerase, and a coalescent species delimitation analysis reveals no notable structure within the complex. This is the first known trio of species within Bembidion for which those genes provide no clear signal of species boundaries. A neotype is designated for the one name in the group that lacks a primary type, Bembidiumhaplogonum Chaudoir. Chromosomes of the new species and their relatives are as is typical for Bembidion, with eleven pairs of autosomes and an XY/XX sex chromosome system.

Maximizing the acquisition of unique reads in noninvasive capture sequencing experiments

Noninvasive samples as a source of DNA are gaining interest in genomic studies of endangered species. However, their complex nature and low endogenous DNA content hamper the recovery of good quality data. Target capture has become a productive method to enrich the endogenous fraction of noninvasive samples, such as faeces, but its sensitivity has not yet been extensively studied. Coping with faecal samples with an endogenous DNA content below 1% is a common problem when prior selection of samples from a large collection is not possible. However, samples classified as unfavourable for target capture sequencing might be the only representatives of unique specific geographical locations, or to answer the question of interest. To explore how library complexity may be increased without repeating DNA extractions and generating new libraries, in this study we captured the exome of 60 chimpanzees (Pan troglodytes) using faecal samples with very low proportions of endogenous content (<1%). Our results indicate that by performing additional hybridizations of the same libraries, the molecular complexity can be maintained to achieve higher coverage. Also, whenever possible, the starting DNA material for capture should be increased. Finally, we specifically calculated the sequencing effort needed to avoid exhausting the library complexity of enriched faecal samples with low endogenous DNA content. This study provides guidelines, schemes and tools for laboratories facing the challenges of working with noninvasive samples containing extremely low amounts of endogenous DNA.

Mitochondrial genomes and polymorphic regions of Gonimbrasia belina and Gynanisa maja (Lepidoptera: Saturniidae), two important edible caterpillars of Southern Africa

Mopane worms are the vernacular designation for the edible caterpillars of the African emperor moths Gonimbrasia belina and Gynanisa maja. Both species, particularly G. belina, are widely harvested in Southern Africa, and their populations are declining. Despite their commercial, nutritional, and cultural importance, their genetic data are currently unavailable. We sequenced two complete mitogenomes from each species using Ion Torrent technology, and identified informative markers in the complete mitogenomes of the two species for use in future studies. Comparing the conspecific mitogenomes allowed the identification of regions with high nucleotide diversity in ATP6, ND1, ND4, ND5, ND6, and CYTB genes. The final panels of markers will allow for the survey of 3117 bp in G. belina, and 3990 bp in Gy. maja. Phylogenetic reconstruction within the family Saturniidae recovered the tribe Bunaeini as monophyletic and basal to Saturniidae, and the tribe Attacini as a monophyletic clade nested within the tribe Saturniini. The G. belina and Gy. maja mitogenomes are the first representatives of African Saturniidae, a taxonomic group with relevance as a food resource on the continent. This study represents the first step towards assessing the genetic diversity, population structure, and phylogeography of African edible caterpillars.

Categorization of species as native or nonnative using DNA sequence signatures without a complete reference library

New genetic diagnostic approaches have greatly aided efforts to document global biodiversity and improve biosecurity. This is especially true for organismal groups in which species diversity has been underestimated historically due to difficulties associated with sampling, the lack of clear morphological characteristics, and/or limited availability of taxonomic expertise. Among these methods, DNA sequence barcoding (also known as "DNA barcoding") and by extension, meta-barcoding for biological communities, has emerged as one of the most frequently utilized methods for DNA-based species identifications. Unfortunately, the use of DNA barcoding is limited by the availability of complete reference libraries (i.e., a collection of DNA sequences from morphologically identified species), and by the fact that the vast majority of species do not have sequences present in reference databases. Such conditions are critical especially in tropical locations that are simultaneously biodiversity rich and suffer from a lack of exploration and DNA characterization by trained taxonomic specialists. To facilitate efforts to document biodiversity in regions lacking complete reference libraries, we developed a novel statistical approach that categorizes unidentified species as being either likely native or likely nonnative based solely on measures of nucleotide diversity. We demonstrate the utility of this approach by categorizing a large sample of specimens of terrestrial insects and spiders (collected as part of the Moorea BioCode project) using a generalized linear mixed model (GLMM). Using a training data set of known endemic (n = 45) and known introduced species (n = 102), we then estimated the likely native/nonnative status for 4,663 specimens representing an estimated 1,288 species (412 identified species), including both those specimens that were either unidentified or whose endemic/introduced status was uncertain. Using this approach, we were able to increase the number of categorized specimens by a factor of 4.4 (from 794 to 3,497), and the number of categorized species by a factor of 4.8 from (147 to 707) at a rate much greater than chance (77.6% accuracy). The study identifies phylogenetic signatures of both native and nonnative species and suggests several practical applications for this approach including monitoring biodiversity and facilitating biosecurity.

DNA barcoding of the rodent genus Oligoryzomys (Cricetidae: Sigmodontinae): mitogenomic-anchored database and identification of nuclear mitochondrial translocations (Numts)

DNA barcoding has become a standard method for species identification in taxonomically complex groups. An important step of the barcoding process is the construction of a library of voucher-based material that was properly identified by independent methods, free of inaccurate identification, and paralogs. We provide here a cytochrome oxidase I (mt-Co1) DNA barcode database for species of the genus Oligoryzomys, based on type material and karyotyped specimens, and anchored on the mitochondrial genome of one species of Oligoryzomys, O. stramineus. To evaluate the taxonomic determination of new COI sequences, we assessed species intra/interspecific genetic distances (barcode gap), performed the General Mixed Yule Coalescent method (GMYC) for lineages' delimitation, and identified diagnostic nucleotides for each species of Oligoryzomys. Phylogenetic analyses of Oligoryzomys were performed on 2 datasets including 14 of the 23 recognized species of this genus: a mt-Co1 only matrix, and a concatenated matrix including mt-Co1, cytochrome b (mt-Cytb), and intron 7 of the nuclear fibrinogen beta chain gene (i7Fgb). We recovered nuclear-mitochondrial translocated (Numts) pseudogenes on our samples and identified several published sequences that are cases of Numts. We analyzed the rate of non-synonymous and synonymous substitution, which were higher in Numts in comparison to mtDNA sequences. GMYC delimitations and DNA barcode gap results highlight the need for further work that integrate molecular, karyotypic, and morphological analyses, as well as additional sampling, to tackle persistent problems in the taxonomy of Oligoryzomys.

The mystery of the origins of Cebus albifrons malitiosus and Cebus albifrons hypoleucus: mitogenomics and microsatellite analyses revealed an amazing evolutionary history of the Northern Colombian white-fronted capuchins

Knowledge of the genetic units of species is fundamental to the conservation of biodiversity. This is true for all regions, including the Neotropics where the Earth has its greatest diversity, including roughly 34% of primate species, a group that has almost 60% of its taxa threatened with extinction. The untufted (gracile) capuchins are medium-sized Neotropical primates, traditionally classified in four species: Cebus albifrons, C. capucinus, C. olivaceus, and C. kaapori. They have a very confusing intra-specific systematics with a large number of fragmented and isolated populations throughout their geographical distributions. We sequenced a large sample of gracile capuchins, including all of the recognized species, to offset the paucity of phylogenic and phylogeographic data regarding this group and to try to understand their phylogeny and evolution. A set of 189 gracile and robust capuchins were sequenced for their mitogenomes whereas another set of 394 gracile and robust capuchins were sequenced at two individual mitochondrial genes (mtCOI-COII). Additionally, 41 Colombian gracile capuchins were geno typified at eight nuclear DNA microsatellites. Our main findings are as follows: (1) Nineteen different groups of gracile capuchin were detected with the mitogenomics data set and more than twenty significant groups and sub-groups were identified with the mtCOI-COII genes; (2) The temporal splits of the older gracile capuchin haplogroups expanded between 2 and 4 million years ago (MYA), during the Pliocene; (3) The two most northern taxa of Colombian C. albifrons (malitiosus and hypoleucus) are the same taxon (C. a. hypoleucus) as was claimed by Cabrera. This taxon represents an old colonization event from the Amazon to current northern Colombia. It is intensely hybridized (evidence from both mitochondrial and nuclear genes) with a haplogroup of C. capucinus (H3) and also has an influx of robust capuchins; (4) Three different and independent migrations of C. albifrons from the Amazon arrived to northern Colombia giving rise to C. a. hypoleucus (including malitiosus), C. a versicolor (including leucocephalus, cesarae, and pleei), and C. a. adustus; (5) On the Caribbean island of Trinidad, two different gracile capuchin taxa exist, one autochthonous, which could correspond to a fourth migration into northwestern South America (C. a. trinitatis) and probably another one, introduced more recently (C. olivaceus brunneus); (6) The values of the genetic distance analyses, the inexistence of reciprocal mitochondrial monophylia for many clades of gracile capuchins and the strong hybridization detected with nuclear microsatellites, especially among hypoleucus (malitiosus), C. capucinus-H3, versicolor, and cesarae, support that all the gracile capuchins belong to one unique superspecies: C. capucinus (senior name for all the gracile capuchins).

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.

Ancient DNA reveals the timing and persistence of organellar genetic bottlenecks over 3,000 years of sunflower domestication and improvement

Here, we report a comprehensive paleogenomic study of archaeological and ethnographic sunflower remains that provides significant new insights into the process of domestication of this important crop. DNA from both ancient and historic contexts yielded high proportions of endogenous DNA, and although archaeological DNA was found to be highly degraded, it still provided sufficient coverage to analyze genetic changes over time. Shotgun sequencing data from specimens from the Eden's Bluff archaeological site in Arkansas yielded organellar DNA sequence from specimens up to 3,100 years old. Their sequences match those of modern cultivated sunflowers and are consistent with an early domestication bottleneck in this species. Our findings also suggest that recent breeding of sunflowers has led to a loss of genetic diversity that was present only a century ago in Native American landraces. These breeding episodes also left a profound signature on the mitochondrial and plastid haplotypes in cultivars, as two types were intentionally introduced from other Helianthus species for crop improvement. These findings gained from ancient and historic sunflower specimens underscore how future in-depth gene-based analyses can advance our understanding of the pace and targets of selection during the domestication of sunflower and other crop species.

An integrative approach to discovering cryptic species within the Bemisia tabaci whitefly species complex

Bemisia tabaci is a cryptic whitefly-species complex that includes some of the most damaging pests and plant-virus vectors of a diverse range of food and fibre crops worldwide. We combine experimental evidence of: (i) differences in reproductive compatibility, (ii) hybrid verification using a specific nuclear DNA marker and hybrid fertility confirmation and (iii) high-throughput sequencing-derived mitogenomes, to show that the "Mediterranean" (MED) B. tabaci comprises at least two distinct biological species; the globally invasive MED from the Mediterranean Basin and the "African silver-leafing" (ASL) from sub-Saharan Africa, which has no associated invasion records. We demonstrate that, contrary to its common name, the "ASL" does not induce squash silver-leafing symptoms and show that species delimitation based on the widely applied 3.5% partial mtCOI gene sequence divergence threshold produces discordant results, depending on the mtCOI region selected. Of the 292 published mtCOI sequences from MED/ASL groups, 158 (54%) are low quality and/or potential pseudogenes. We demonstrate fundamental deficiencies in delimiting cryptic B. tabaci species, based solely on partial sequences of a mitochondrial barcoding gene. We advocate an integrative approach to reveal the true species richness within cryptic species complexes, which is integral to the deployment of effective pest and disease management strategies.

Population genetic structure and evolutionary history of Bale monkeys (Chlorocebus djamdjamensis) in the southern Ethiopian Highlands

Background

Species with a restricted geographic distribution, and highly specialized habitat and dietary requirements, are particularly vulnerable to extinction. The Bale monkey (Chlorocebus djamdjamensis) is a little-known arboreal, bamboo-specialist primate endemic to the southern Ethiopian Highlands. While most Bale monkeys inhabit montane forests dominated by bamboo, some occupy forest fragments where bamboo is much less abundant. We used mitochondrial DNA (mtDNA) sequences to analyse the genetic structure and evolutionary history of Bale monkeys covering the majority of their remaining distribution range. We analysed 119 faecal samples from their two main habitats, continuous forest (CF) and fragmented forests (FF), and sequenced 735 bp of the hypervariable region I (HVI) of the control region. We added 12 orthologous sequences from congeneric vervets (C. pygerythrus) and grivets (C. aethiops) as well as animals identified as hybrids, previously collected in southern Ethiopia.

Results

We found strong genetic differentiation (with no shared mtDNA haplotypes) between Bale monkey populations from CF and FF. Phylogenetic analyses revealed two distinct and highly diverged clades: a Bale monkey clade containing only Bale monkeys from CF and a green monkey clade where Bale monkeys from FF cluster with grivets and vervets. Analyses of demographic history revealed that Bale monkey populations (CF and FF) have had stable population sizes over an extended period, but have all recently experienced population declines.

Conclusions

The pronounced genetic structure and deep mtDNA divergence between Bale monkey populations inhabiting CF and FF are likely to be the results of hybridization and introgression of the FF population with parapatric Chlorocebus species, in contrast to the CF population, which was most likely not impacted by hybridization. Hybridization in the FF population was probably enhanced by an alteration of the bamboo forest habitat towards a more open woodland habitat, which enabled the parapatric Chlorocebus species to invade the Bale monkey's range and introgress the FF population. We therefore propose that the CF and FF Bale monkey populations should be managed as separate units when developing conservation strategies for this threatened species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1217-y) contains supplementary material, which is available to authorized users.

Using a multi-gene approach to infer the complicated phylogeny and evolutionary history of lorises (Order Primates: Family Lorisidae)

Extensive phylogenetic studies have found robust phylogenies are modeled by using a multi-gene approach and sampling from the majority of the taxa of interest. Yet, molecular studies focused on the lorises, a cryptic primate family, have often relied on one gene, or just mitochondrial DNA, and many were unable to include all four genera in the analyses, resulting in inconclusive phylogenies. Past phylogenetic loris studies resulted in lorises being monophyletic, paraphyletic, or an unresolvable trichotomy with the closely related galagos. The purpose of our study is to improve our understanding of loris phylogeny and evolutionary history by using a multi-gene approach. We used the mitochondrial genes cytochrome b, and cytochrome c oxidase subunit 1, along with a nuclear intron (recombination activating gene 2) and nuclear exon (the melanocortin 1 receptor). Maximum Likelihood and Bayesian phylogenetic analyses were conducted based on data from each locus, as well as on the concatenated sequences. The robust, concatenated results found lorises to be a monophyletic family (Lorisidae) (PP ≥ 0.99) with two distinct subfamilies: the African Perodictinae (PP ≥ 0.99) and the Asian Lorisinae (PP ≥ 0.99). Additionally, from these analyses all four genera were all recovered as monophyletic (PP ≥ 0.99). Some of our single-gene analyses recovered monophyly, but many had discordances, with some showing paraphyly or a deep-trichotomy. Bayesian partitioned analyses inferred the most recent common ancestors of lorises emerged ∼42 ± 6 million years ago (mya), the Asian Lorisinae separated ∼30 ± 9 mya, and Perodictinae arose ∼26 ± 10 mya. These times fit well with known historical tectonic shifts of the area, as well as with the sparse loris fossil record. Additionally, our results agree with previous multi-gene studies on Lorisidae which found lorises to be monophyletic and arising ∼40 mya (Perelman et al., 2011; Pozzi et al., 2014). By taking a multi-gene approach, we were able to recover a well-supported, monophyletic loris phylogeny and inferred the evolutionary history of this cryptic family.

The impact of endogenous content, replicates and pooling on genome capture from faecal samples

Target-capture approach has improved over the past years, proving to be very efficient tool for selectively sequencing genetic regions of interest. These methods have also allowed the use of noninvasive samples such as faeces (characterized by their low quantity and quality of endogenous DNA) to be used in conservation genomic, evolution and population genetic studies. Here we aim to test different protocols and strategies for exome capture using the Roche SeqCap EZ Developer kit (57.5 Mb). First, we captured a complex pool of DNA libraries. Second, we assessed the influence of using more than one faecal sample, extract and/or library from the same individual, to evaluate its effect on the molecular complexity of the experiment. We validated our experiments with 18 chimpanzee faecal samples collected from two field sites as a part of the Pan African Programme: The Cultured Chimpanzee. Those two field sites are in Kibale National Park, Uganda (N = 9) and Loango National Park, Gabon (N = 9). We demonstrate that at least 16 libraries can be pooled, target enriched through hybridization, and sequenced allowing for the genotyping of 951,949 exome markers for population genetic analyses. Further, we observe that molecule richness, and thus, data acquisition, increase when using multiple libraries from the same extract or multiple extracts from the same sample. Finally, repeated captures significantly decrease the proportion of off-target reads from 34.15% after one capture round to 7.83% after two capture rounds, supporting our conclusion that two rounds of target enrichment are advisable when using complex faecal samples.

Small spotted bodies with multiple specific mitochondrial DNAs: existence of diverse and differentiated tigrina lineages or species (Leopardus spp: Felidae, Mammalia) throughout Latin America

We analysed two sets of mitochondrial (mt) DNA data from tigrinas (traditionally, Leopardus tigrinus) we sampled in Costa Rica, Venezuela, Colombia, Ecuador, Peru, Bolivia, northwestern and northeastern Argentina and southern Brazil. Additionally, the analysis included some GenBank sequences from southern, central and northeastern Brazil. The first mt set (mt ATP8+mt 16S rRNA with 41 tigrina) revealed the existence of seven different tigrina-like haplogroups. They could represent, at least, 4-6 different tigrina species following the Phylogenetic Species Concept (PSC). In the second mt set (mitogenomics with 18 tigrinas), we detected six different tigrina-like haplogroups. They could represent 4-5 different tigrina species - including a possible full new species, which has gone previously unnoticed to the world of science both morphologic and molecularly. Coat patterns of several of these different tigrinas support the molecular differences. We also detected intense hybridization in many Andean tigrina with margays (Leopardus wiedii) and ocelots (Leopardus pardalis) as well as hybridization of one Bolivian tigrina with Leopardus geoffroyi. Similar hybridization was found for many of the southern Brazilian tigrina (Leopardus guttulus). All of the temporal split estimates for these tigrina haplogroups, together with those of the Leopardus species recognized to date, began in the late Pliocene but mostly occurred during the Pleistocene. In agreement with the existence of multiple species within the traditional L. tigrinus species, we detected strong and significant spatial structure in the two mt data sets. There were clear circular clines. A major part of the analyses detected more genetic resemblance between the Central American + trans Andean Colombian and Ecuadorian tigrina (L. oncilla) with the most geographically distant tigrina from central and southern Brazil (L. guttulus; pure individuals not hybridized with L. geoffroyi). In comparison, the Andean tigrina taxa had intermediate geographical origins but were highly genetically differentiated both from the Central American + trans Andean Colombian-Ecuadorian tigrina and from the central and southern Brazilian tigrina.

Quantifying the Number of Independent Organelle DNA Insertions in Genome Evolution and Human Health

Fragments of organelle genomes are often found as insertions in nuclear DNA. These fragments of mitochondrial DNA (numts) and plastid DNA (nupts) are ubiquitous components of eukaryotic genomes. They are, however, often edited out during the genome assembly process, leading to systematic underestimation of their frequency. Numts and nupts, once inserted, can become further fragmented through subsequent insertion of mobile elements or other recombinational events that disrupt the continuity of the inserted sequence relative to the genuine organelle DNA copy. Because numts and nupts are typically identified through sequence comparison tools such as BLAST, disruption of insertions into smaller fragments can lead to systematic overestimation of numt and nupt frequencies. Accurate identification of numts and nupts is important, however, both for better understanding of their role during evolution, and for monitoring their increasingly evident role in human disease. Human populations are polymorphic for 141 numt loci, five numts are causal to genetic disease, and cancer genomic studies are revealing an abundance of numts associated with tumor progression. Here, we report investigation of salient parameters involved in obtaining accurate estimates of numt and nupt numbers in genome sequence data. Numts and nupts from 44 sequenced eukaryotic genomes reveal lineage-specific differences in the number, relative age and frequency of insertional events as well as lineage-specific dynamics of their postinsertional fragmentation. Our findings outline the main technical parameters influencing accurate identification and frequency estimation of numts in genomic studies pertinent to both evolution and human health.

DNA barcoding and traditional taxonomy: an integrated approach for biodiversity conservation

Biological diversity is depleting at an alarming rate. Additionally, a vast amount of biodiversity still remains undiscovered. Taxonomy has been serving the purpose of describing, naming, and classifying species for more than 250 years. DNA taxonomy and barcoding have accelerated the rate of this process, thereby providing a tool for conservation practice. DNA barcoding and traditional taxonomy have their own inherent merits and demerits. The synergistic use of both methods, in the form of integrative taxonomy, has the potential to contribute to biodiversity conservation in a pragmatic timeframe and overcome their individual drawbacks. In this review, we discuss the basics of both these methods of biological identification (traditional taxonomy and DNA barcoding), the technical advances in integrative taxonomy, and future trends. We also present a comprehensive compilation of published examples of integrative taxonomy that refer to nine topics within biodiversity conservation. Morphological and molecular species limits were observed to be congruent in ∼41% of the 58 source studies. The majority of the studies highlighted the description of cryptic diversity through the use of molecular data, whereas research areas like endemism, biological invasion, and threatened species were less discussed in the literature.

Comparing mitogenomic timetrees for two African savannah primate genera (Chlorocebus and Papio)

Complete mitochondrial (mtDNA) genomes have proved to be useful in reconstructing primate phylogenies with higher resolution and confidence compared to reconstructions based on partial mtDNA sequences. Here, we analyse complete mtDNA genomes of African green monkeys (genus Chlorocebus), a widely distributed primate genus in Africa representing an interesting phylogeographical model for the evolution of savannah species. Previous studies on partial mtDNA sequences revealed nine major clades, suggesting several cases of para- and polyphyly among Chlorocebus species. However, in these studies, phylogenetic relationships among several clades were not resolved, and divergence times were not estimated. We analysed complete mtDNA genomes for ten Chlorocebus samples representing major mtDNA clades to find stronger statistical support in the phylogenetic reconstruction than in the previous studies and to estimate divergence times. Our results confirmed para- and polyphyletic relationships of most Chlorocebus species, while the support for the phylogenetic relationships between the mtDNA clades increased compared to the previous studies. Our results indicate an initial west–east division in the northern part of the Chlorocebus range with subsequent divergence into north-eastern and southern clades. This phylogeographic scenario contrasts with that for another widespread African savannah primate genus, the baboons (Papio), for which a dispersal from southern Africa into East and West Africa was suggested.

Low genetic diversity and strong population structure shaped by anthropogenic habitat fragmentation in a critically endangered primate, Trachypithecus leucocephalus

Habitat fragmentation may strongly impact population genetic structure and reduce the genetic diversity and viability of small and isolated populations. The white-headed langur (Trachypithecus leucocephalus) is a critically endangered primate species living in a highly fragmented and human-modified habitat in southern China. We examined the population genetic structure and genetic diversity of the species and investigated the environmental and anthropogenic factors that may have shaped its population structure. We used 214 unique multi-locus genotypes from 41 social groups across the main distribution area of T. leucocephalus, and found strong genetic structure and significant genetic differentiation among local populations. Our landscape genetic analyses using a causal modelling framework suggest that a large habitat gap and geographical distance represent the primary landscape elements shaping genetic structure, yet high levels of genetic differentiation also exist between patches separated by a small habitat gap or road. This is the first comprehensive study that has evaluated the population genetic structure and diversity of T. leucocephalus using nuclear markers. Our results indicate strong negative impacts of anthropogenic land modifications and habitat fragmentation on primate genetic connectivity between forest patches. Our analyses suggest that two management units of the species could be defined, and indicate that habitat continuity should be enforced and restored to reduce genetic isolation and enhance population viability.

Mitogenomics phylogenetic relationships of the current sloth's genera and species (Bradypodidae and Megalonychidae)

We sequenced the complete mitogenome of 39 sloths (19 Bradypus variegatus, 4 B. tridactylus, 1 B. pygmaeus, 1 B. torquatus, 4 Choloepus didactylus, and 10 C. hoffmanni). A Bayesian tree (BI) indicated a temporal split between Bradypus and Choloepus around 31 million years ago (MYA, Oligocene) and the other major splits within each genera during the Miocene and Pliocene. A haplotype network (MJN) estimated a lower temporal split between the sloth genera (around 23.5 MYA). Both methods detected the ancestor of B. torquatus as the first to diverge within Bradypus (21 for BI and 19 MJN), followed by that of the ancestor of B. tridactylus. The split of B. pygmaeus from the common ancestor with B. variegatus was around 12 MYA (BI) or 4.3 MYA (MJN). The splits among the previous populations of B. variegatus began around 8 MYA (BI) or 3.6 MYA (MJN). The trans-Andean population was the first to diverge from the remaining cis-Andean populations of B. variegatus. The genetic differentiation of the trans-Andean B. variegatus population relative to the cis-Andean B. variegatus is similar to that found for different species of sloths. The mitogenomic analysis resolved the differentiation of C. hoffmanni from the C. didactylus individuals of the Guiana Shield. However, one C. didactylus from the Colombian Amazon specimen was inside the C. hoffmanni clade. This could be the first example of possible natural hybridization in the Amazon of both Choloepus taxa or the existence of un-differentiable phenotypes of these two species in some Amazonian areas.

Unbiased Mitoproteome Analyses Confirm Non-canonical RNA, Expanded Codon Translations

Proteomic MS/MS mass spectrometry detections are usually biased towards peptides cleaved by experimentally added digestion enzyme(s). Hence peptides resulting from spontaneous degradation and natural proteolysis usually remain undetected. Previous analyses of tryptic human proteome data (cleavage after K, R) detected non-canonical tryptic peptides translated according to tetra- and pentacodons (codons expanded by silent mono- and dinucleotides), and from transcripts systematically (a) deleting mono-, dinucleotides after trinucleotides (delRNAs), (b) exchanging nucleotides according to 23 bijective transformations. Nine symmetric and fourteen asymmetric nucleotide exchanges (X ↔ Y, e.g. A ↔ C; and X → Y → Z → X, e.g. A → C → G → A) produce swinger RNAs. Here unbiased reanalyses of these proteomic data detect preferentially non-canonical tryptic peptides despite assuming random cleavage. Unbiased analyses couldn't reconstruct experimental tryptic digestion if most detected non-canonical peptides were false positives. Detected non-tryptic non-canonical peptides map preferentially on corresponding, previously described non-canonical transcripts, as for tryptic non-canonical peptides. Hence unbiased analyses independently confirm previous trypsin-biased analyses that showed translations of del- and swinger RNA and expanded codons. Accounting for natural proteolysis completes trypsin-biased mitopeptidome analyses, independently confirms non-canonical transcriptions and translations.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

Reintroduction of confiscated and displaced mammals risks outbreeding and introgression in natural populations, as evidenced by orang-utans of divergent subspecies

Confiscated and displaced mammals are often taken to sanctuaries, where the explicit goal may be reintroduction to the wild. By inadvertently collecting animals from different source populations, however, such efforts risk reintroducing individuals that have not been in genetic contact for significant periods of time. Using genetic analyses and 44 years of data from Camp Leakey, an orang-utan rehabilitation site on Borneo, we determined the minimum extent to which orang-utans representing non-native, geographically and reproductively isolated taxa were reintroduced into the surrounding wild population. We found two reintroduced females were from a non-native subspecies, and have since produced at least 22 hybridized and introgressed descendants to date, of which at least 15 are living. Given that Bornean orang-utan subspecies are thought to have diverged from a common ancestor around 176,000 years ago, with marked differentiation over the last 80,000 years, we highlight the need for further evaluation of the effects of hybridizing orang-utans of different taxa--particularly in light of the ~1500 displaced orang-utans awaiting urgent reintroduction. As endangered mammals are increasing in number in sanctuaries worldwide, we stress the need for re-examination of historical reintroductions, to assess the extent and effects of de facto translocations in the past.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale–multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8–15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera)

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

Population structure of the Indonesian giant tiger shrimp Penaeus monodon: a window into evolutionary similarities between paralogous mitochondrial DNA sequences and their genomes

Here we used both microsatellites and mtCR (mitochondrial DNA control region) sequences as genetic markers to examine the genetic diversity and population structure of Penaeus monodon shrimp from six Indonesian regions. The microsatellite data showed that shrimp from the Indian and the Pacific Ocean were genetically distinct from each other. It has been reported previously that P. monodon mtCR sequences from the Indo-Pacific group into two major paralogous clades of unclear origin. Here we show that the population structure inferred from mtCR sequences matches the microsatellite-based population structure for one of these clades. This is consistent with the notion that this mtCR clade shares evolutionary history with nuclear DNA and may thus represent nuclear mitochondrial pseudogenes (Numts).

Odintifier--A computational method for identifying insertions of organellar origin from modern and ancient high-throughput sequencing data based on haplotype phasing

Background

Cellular organelles with genomes of their own (e.g. plastids and mitochondria) can pass genetic sequences to other organellar genomes within the cell in many species across the eukaryote phylogeny. The extent of the occurrence of these organellar-derived inserted sequences (odins) is still unknown, but if not accounted for in genomic and phylogenetic studies, they can be a source of error. However, if correctly identified, these inserted sequences can be used for evolutionary and comparative genomic studies. Although such insertions can be detected using various laboratory and bioinformatic strategies, there is currently no straightforward way to apply them as a standard organellar genome assembly on next-generation sequencing data. Furthermore, most current methods for identification of such insertions are unsuitable for use on non-model organisms or ancient DNA datasets.

Results

We present a bioinformatic method that uses phasing algorithms to reconstruct both source and inserted organelle sequences. The method was tested in different shotgun and organellar-enriched DNA high-throughput sequencing (HTS) datasets from ancient and modern samples. Specifically, we used datasets from lions (Panthera leo ssp. and Panthera leo leo) to characterize insertions from mitochondrial origin, and from common grapevine (Vitis vinifera) and bugle (Ajuga reptans) to characterize insertions derived from plastid genomes. Comparison of the results against other available organelle genome assembly methods demonstrated that our new method provides an improvement in the sequence assembly.

Conclusion

Using datasets from a wide range of species and different levels of complexity we showed that our novel bioinformatic method based on phasing algorithms can be used to achieve the next two goals: i) reference-guided assembly of chloroplast/mitochondrial genomes from HTS data and ii) identification and simultaneous assembly of odins. This method represents the first application of haplotype phasing for automatic detection of odins and reference-based organellar genome assembly.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0682-1) contains supplementary material, which is available to authorized users.

Analysis of the complete mitochondrial genome and characterization of diverse NUMTs of Macaca leonina

As a non-human primate, the pig-tailed macaque has received wide attention because it can be infected by HIV-1. In this study, we determined the complete mtDNA sequence of the northern pig-tailed macaque (Macaca leonina). Unexpectedly, during the amplification of the mtDNA control region (D-loop region) we observed several D-loop-like sequences, which were NUMTs (nuclear mitochondrial sequences) and a total of 14 D-loop-like NUMT haplotypes were later identified in five individuals. The neighbor-joining tree and estimated divergence time based on these D-loop-like NUMT sequences of M. leonina provide some insights into the understanding of the evolutionary history of NUMTs. D-loop-like haplotypes G and H, which also exist in the nuclear genome of mulatta, appear to have been translocated into the nuclear genome before the divergence of M. mulatta and M. leonina. The other D-loop-like NUMT haplotypes were translocated into the nuclear genome of M. leonina after the divergence of the two species. Later sequence conversion was predicted to occur among these 14 D-loop-like NUMT haplotypes. The overall structure of the mtDNA of M. leonina was found to be similar to that seen in other mammalian mitochondrial genomes. Phylogenetic analysis based on the maximum likelihood method shows M. leonina clustered with Macaca silenus among the analyzed mammalian species.

DNA barcoding to fishes: current status and future directions

DNA barcoding appears to be a promising approach for taxonomic identification, characterization, and discovery of newer species, facilitating biodiversity studies. It helps researchers to appreciate genetic and evolutionary associations by collection of molecular, morphological, and distributional data. Fish DNA barcoding, based on the sequencing of a uniform area of Cytochrome C Oxidase type I (COI) gene, has received significant interest as an accurate tool for species identification, authentication, and phylogenetic analysis. The aim of this review article was to investigate recent global status, approaches, and future direction of DNA barcoding in fisheries sectors. We have tried to highlight its possible impacts, complications, and validation issues at species levels for biodiversity analysis. Moreover, an effort has been put forward to understand issues related to various marker genes associated with barcode process as primer sequences and have concluded barcode promotion as an indispensable tool of molecular biology for the development of taxonomic support systems.

Low Genetic Diversity and Strong Geographical Structure of the Critically Endangered White-Headed Langur (Trachypithecus leucocephalus) Inferred from Mitochondrial DNA Control Region Sequences

Many Asian colobine monkey species are suffering from habitat destruction and population size decline. There is a great need to understand their genetic diversity, population structure and demographic history for effective species conservation. The white-headed langur (Trachypithecus leucocephalus) is a Critically Endangered colobine species endemic to the limestone karst forests in southwestern China. We analyzed the mitochondrial DNA (mtDNA) control region sequences of 390 fecal samples from 40 social groups across the main distribution areas, which represented one-third of the total extant population. Only nine haplotypes and 10 polymorphic sites were identified, indicating remarkably low genetic diversity in the species. Using a subset of 77 samples from different individuals, we evaluated genetic variation, population structure, and population demographic history. We found very low values of haplotype diversity (h = 0.570 ± 0.056) and nucleotide diversity (π = 0.00323 ± 0.00044) in the hypervariable region I (HVRI) of the mtDNA control region. Distribution of haplotypes displayed marked geographical pattern, with one population (Chongzuo, CZ) showing a complete lack of genetic diversity (having only one haplotype), whereas the other population (Fusui, FS) having all nine haplotypes. We detected strong population genetic structure among habit patches (Φ_ST = 0.375, P < 0.001). In addition, the Mantel test showed a significant correlation between the pairwise genetic distances and geographical distances among social groups in FS (correlation coefficient = 0.267, P = 0.003), indicting isolation-by-distance pattern of genetic divergence in the mtDNA sequences. Analyses of demographic history suggested an overall stable historical population size and modest population expansion in the last 2,000 years. Our results indicate different genetic diversity and possibly distinct population history for different local populations, and suggest that CZ and FS should be considered as one evolutionarily significant unit (ESU) and two management units (MUs) pending further investigation using nuclear markers.

Pleistocene aridification cycles shaped the contemporary genetic architecture of Southern African baboons

Plio-Pleistocene environmental change influenced the evolutionary history of many animal lineages in Africa, highlighting key roles for both climate and tectonics in the evolution of Africa’s faunal diversity. Here, we explore diversification in the southern African chacma baboon Papio ursinus sensu lato and reveal a dominant role for increasingly arid landscapes during past glacial cycles in shaping contemporary genetic structure. Recent work on baboons (Papio spp.) supports complex lineage structuring with a dominant pulse of diversification occurring 1-2Ma, and yet the link to palaeoenvironmental change remains largely untested. Phylogeographic reconstruction based on mitochondrial DNA sequence data supports a scenario where chacma baboon populations were likely restricted to refugia during periods of regional cooling and drying through the Late Pleistocene. The two lineages of chacma baboon, ursinus and griseipes, are strongly geographically structured, and demographic reconstruction together with spatial analysis of genetic variation point to possible climate-driven isolating events where baboons may have retreated to more optimum conditions during cooler, drier periods. Our analysis highlights a period of continuous population growth beginning in the Middle to Late Pleistocene in both the ursinus and the PG2 griseipes lineages. All three clades identified in the study then enter a state of declining population size (Nef) through to the Holocene; this is particularly marked in the last 20,000 years, most likely coincident with the Last Glacial Maximum. The pattern recovered here conforms to expectations based on the dynamic regional climate trends in southern Africa through the Pleistocene and provides further support for complex patterns of diversification in the region’s biodiversity.