Ecological niche modelling and nDNA sequencing support a new, morphologically cryptic beetle species unveiled by DNA barcoding

Background

DNA sequencing techniques used to estimate biodiversity, such as DNA barcoding, may reveal cryptic species. However, disagreements between barcoding and morphological data have already led to controversy. Species delimitation should therefore not be based on mtDNA alone. Here, we explore the use of nDNA and bioclimatic modelling in a new species of aquatic beetle revealed by mtDNA sequence data.

Methodology/Principal Findings

The aquatic beetle fauna of Australia is characterised by high degrees of endemism, including local radiations such as the genus Antiporus. Antiporus femoralis was previously considered to exist in two disjunct, but morphologically indistinguishable populations in south-western and south-eastern Australia. We constructed a phylogeny of Antiporus and detected a deep split between these populations. Diagnostic characters from the highly variable nuclear protein encoding arginine kinase gene confirmed the presence of two isolated populations. We then used ecological niche modelling to examine the climatic niche characteristics of the two populations. All results support the status of the two populations as distinct species. We describe the south-western species as Antiporus occidentalis sp.n.

Conclusion/Significance

In addition to nDNA sequence data and extended use of mitochondrial sequences, ecological niche modelling has great potential for delineating morphologically cryptic species.     

Mitogenomic phylogeny of the common long-tailed macaque (Macaca fascicularis fascicularis)

Background

Long-tailed macaques (Macaca fascicularis) are an important model species in biomedical research and reliable knowledge about their evolutionary history is essential for biomedical inferences. Ten subspecies have been recognized, of which most are restricted to small islands of Southeast Asia. In contrast, the common long-tailed macaque (M. f. fascicularis) is distributed over large parts of the Southeast Asian mainland and the Sundaland region. To shed more light on the phylogeny of M. f. fascicularis, we sequenced complete mitochondrial (mtDNA) genomes of 40 individuals from all over the taxon’s range, either by classical PCR-amplification and Sanger sequencing or by DNA-capture and high-throughput sequencing.

Results

Both laboratory approaches yielded complete mtDNA genomes from M. f. fascicularis with high accuracy and/or coverage. According to our phylogenetic reconstructions, M. f. fascicularis initially diverged into two clades 1.70 million years ago (Ma), with one including haplotypes from mainland Southeast Asia, the Malay Peninsula and North Sumatra (Clade A) and the other, haplotypes from the islands of Bangka, Java, Borneo, Timor, and the Philippines (Clade B). The three geographical populations of Clade A appear as paraphyletic groups, while local populations of Clade B form monophyletic clades with the exception of a Philippine individual which is nested within the Borneo clade. Further, in Clade B the branching pattern among main clades/lineages remains largely unresolved, most likely due to their relatively rapid diversification 0.93-0.84 Ma.

Conclusions

Both laboratory methods have proven to be powerful to generate complete mtDNA genome data with similarly high accuracy, with the DNA-capture and high-throughput sequencing approach as the most promising and only practical option to obtain such data from highly degraded DNA, in time and with relatively low costs. The application of complete mtDNA genomes yields new insights into the evolutionary history of M. f. fascicularis by providing a more robust phylogeny and more reliable divergence age estimations than earlier studies.

Electronic supplementary material

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

Delimiting Species without Nuclear Monophyly in Madagascar's Mouse Lemurs

Background

Speciation begins when populations become genetically separated through a substantial reduction in gene flow, and it is at this point that a genetically cohesive set of populations attain the sole property of species: the independent evolution of a population-level lineage. The comprehensive delimitation of species within biodiversity hotspots, regardless of their level of divergence, is important for understanding the factors that drive the diversification of biota and for identifying them as targets for conservation. However, delimiting recently diverged species is challenging due to insufficient time for the differential evolution of characters—including morphological differences, reproductive isolation, and gene tree monophyly—that are typically used as evidence for separately evolving lineages.

Methodology

In this study, we assembled multiple lines of evidence from the analysis of mtDNA and nDNA sequence data for the delimitation of a high diversity of cryptically diverged population-level mouse lemur lineages across the island of Madagascar. Our study uses a multi-faceted approach that applies phylogenetic, population genetic, and genealogical analysis for recognizing lineage diversity and presents the most thoroughly sampled species delimitation of mouse lemur ever performed.

Conclusions

The resolution of a large number of geographically defined clades in the mtDNA gene tree provides strong initial evidence for recognizing a high diversity of population-level lineages in mouse lemurs. We find additional support for lineage recognition in the striking concordance between mtDNA clades and patterns of nuclear population structure. Lineages identified using these two sources of evidence also exhibit patterns of population divergence according to genealogical exclusivity estimates. Mouse lemur lineage diversity is reflected in both a geographically fine-scaled pattern of population divergence within established and geographically widespread taxa, as well as newly resolved patterns of micro-endemism revealed through expanded field sampling into previously poorly and well-sampled regions.     

Localisation of NORs and counterstain-enhanced fluorescence studies in Salmo gairdneri and Salmo trutta (Pisces,Salmonidae)

Summary The karyotypes of the rainbow trout (Salmo gairdneri R.) and the brown trout (Salmo trutta L.) were analyzed by means of silver staining and the chromomycin A₃/distamycin A/DAPI fluorescence banding technique. The nucleolus organizer regions (NORs) were localized at the secondary constrictions of chromosome no. 14 in S. gairdneri and of chromosome no. 10 in S. trutta. Additional silver positive dots were observed at or close to several centromeres in S. gairdneri. Brilliant chromomycin A₃ (CMA₃) fluorescence heterochromatin blocks were localized on both sides of the nucleolar constrictions in S. gairdneri. A polymorphic CMA₃ positive band was detected close to the NORs of S. trutta. No distamycin A/DAPI intense heterochromatin blocks were detected in the genomes of the two Salmo species investigated.     

Phylogenetic analysis of the genus Sorghum based on combined sequence data from cpDNA regions and ITS generate well-supported trees with two major lineages

Background and Aims

Wild Sorghum species provide novel traits for both biotic and abiotic stress resistance and yield for the improvement of cultivated sorghum. A better understanding of the phylogeny in the genus Sorghum will enhance use of the valuable agronomic traits found in wild sorghum.

Methods

Four regions of chloroplast DNA (cpDNA; psbZ-trnG, trnY-trnD, trnY-psbM and trnT-trnL) and the internal transcribed spacer (ITS) of nuclear ribosomal DNA were used to analyse the phylogeny of sorghum based on maximum-parsimony analyses.

Key Results

Parsimony analyses of the ITS and cpDNA regions as separate or combined sequence datasets formed trees with strong bootstrap support with two lineages: the Eu-sorghum species S. laxiflorum and S. macrospermum in one and Stiposorghum and Para-sorghum in the other. Within Eu-sorghum, S. bicolor-3, -11 and -14 originating from southern Africa form a distinct clade. S. bicolor-2, originally from Yemen, is distantly related to other S. bicolor accessions.

Conclusions

Eu-sorghum species are more closely related to S. macrospermum and S. laxiflorum than to any other Australian wild Sorghum species. S. macrospermum and S. laxiflorum are so closely related that it is inappropriate to classify them in separate sections. S. almum is closely associated with S. bicolor, suggesting that the latter is the maternal parent of the former given that cpDNA is maternally inherited in angiosperms. S. bicolor-3, -11 and -14, from southern Africa, are closely related to each other but distantly related to S. bicolor-2.     

The genome of the truffle-parasite Tolypocladium ophioglossoides and the evolution of antifungal peptaibiotics

Background

Two major mycoparasitic lineages, the family Hypocreaceae and the genus Tolypocladium, exist within the fungal order, Hypocreales. Peptaibiotics are a group of secondary metabolites almost exclusively described from Trichoderma species of Hypocreaceae. Peptaibiotics are produced by nonribosomal peptide synthetases (NRPSs) and have antibiotic and antifungal activities. Tolypocladium species are mainly truffle parasites, but a few species are insect pathogens.

Results

The draft genome sequence of the truffle parasite Tolypocladium ophioglossoides was generated and numerous secondary metabolite clusters were discovered, many of which have no known putative product. However, three large peptaibiotic gene clusters were identified using phylogenetic analyses. Peptaibiotic genes are absent from the predominantly plant and insect pathogenic lineages of Hypocreales, and are therefore exclusive to the largely mycoparasitic lineages. Using NRPS adenylation domain phylogenies and reconciliation of the domain tree with the organismal phylogeny, it is demonstrated that the distribution of these domains is likely not the product of horizontal gene transfer between mycoparasitic lineages, but represents independent losses in insect pathogenic lineages. Peptaibiotic genes are less conserved between species of Tolypocladium and are the product of complex patterns of lineage sorting and module duplication. In contrast, these genes are more conserved within the genus Trichoderma and consistent with diversification through speciation.

Conclusions

Peptaibiotic NRPS genes are restricted to mycoparasitic lineages of Hypocreales, based on current sampling. Phylogenomics and comparative genomics can provide insights into the evolution of secondary metabolite genes, their distribution across a broader range of taxa, and their possible function related to host specificity.

Electronic supplementary material

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

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.     

Evolutionary analyses of entire genomes do not support the association of mtDNA mutations with Ras/MAPK pathway syndromes

Background

There are several known autosomal genes responsible for Ras/MAPK pathway syndromes, including Noonan syndrome (NS) and related disorders (such as LEOPARD, neurofibromatosis type 1), although mutations of these genes do not explain all cases. Due to the important role played by the mitochondrion in the energetic metabolism of cardiac muscle, it was recently proposed that variation in the mitochondrial DNA (mtDNA) genome could be a risk factor in the Noonan phenotype and in hypertrophic cardiomyopathy (HCM), which is a common clinical feature in Ras/MAPK pathway syndromes. In order to test these hypotheses, we sequenced entire mtDNA genomes in the largest series of patients suffering from Ras/MAPK pathway syndromes analyzed to date (n = 45), most of them classified as NS patients (n = 42).

Methods/Principal Findings

The results indicate that the observed mtDNA lineages were mostly of European ancestry, reproducing in a nutshell the expected haplogroup (hg) patterns of a typical Iberian dataset (including hgs H, T, J, and U). Three new branches of the mtDNA phylogeny (H1j1, U5b1e, and L2a5) are described for the first time, but none of these are likely to be related to NS or Ras/MAPK pathway syndromes when observed under an evolutionary perspective. Patterns of variation in tRNA and protein genes, as well as redundant, private and heteroplasmic variants, in the mtDNA genomes of patients were as expected when compared with the patterns inferred from a worldwide mtDNA phylogeny based on more than 8700 entire genomes. Moreover, most of the mtDNA variants found in patients had already been reported in healthy individuals and constitute common polymorphisms in human population groups.

Conclusions/Significance

As a whole, the observed mtDNA genome variation in the NS patients was difficult to reconcile with previous findings that indicated a pathogenic role of mtDNA variants in NS.     

Corolla morphology influences diversification rates in bifid toadflaxes (Linaria sect. Versicolores)

Background and Aims

The role of flower specialization in plant speciation and evolution remains controversial. In this study the evolution of flower traits restricting access to pollinators was analysed in the bifid toadflaxes (Linaria sect. Versicolores), a monophyletic group of ∼30 species and subspecies with highly specialized corollas.

Methods

A time-calibrated phylogeny based on both nuclear and plastid DNA sequences was obtained using a coalescent-based method, and flower morphology was characterized by means of morphometric analyses. Directional trends in flower shape evolution and trait-dependent diversification rates were jointly analysed using recently developed methods, and morphological shifts were reconstructed along the phylogeny. Pollinator surveys were conducted for a representative sample of species.

Key Results

A restrictive character state (narrow corolla tube) was reconstructed in the most recent common ancestor of Linaria sect. Versicolores. After its early loss in the most species-rich clade, this character state has been convergently reacquired in multiple lineages of this clade in recent times, yet it seems to have exerted a negative influence on diversification rates. Comparative analyses and pollinator surveys suggest that the narrow- and broad-tubed flowers are evolutionary optima representing divergent strategies of pollen placement on nectar-feeding insects.

Conclusions

The results confirm that different forms of floral specialization can lead to dissimilar evolutionary success in terms of diversification. It is additionally suggested that opposing individual-level and species-level selection pressures may have driven the evolution of pollinator-restrictive traits in bifid toadflaxes.     

Characterizing genetic diversity of contemporary pacific chickens using mitochondrial DNA analyses

Background

Mitochondrial DNA (mtDNA) hypervariable region (HVR) sequences of prehistoric Polynesian chicken samples reflect dispersal of two haplogroups—D and E—by the settlers of the Pacific. The distribution of these chicken haplogroups has been used as an indicator of human movement. Recent analyses suggested similarities between prehistoric Pacific and South American chicken samples, perhaps reflecting prehistoric Polynesian introduction of the chicken into South America. These analyses have been heavily debated. The current distribution of the D and E lineages among contemporary chicken populations in the Western Pacific is unclear, but might ultimately help to inform debates about the movements of humans that carried them.

Objectives

We sought to characterize contemporary mtDNA diversity among chickens in two of the earliest settled archipelagoes of Remote Oceania, the Marianas and Vanuatu.

Methods

We generated HVR sequences for 43 chickens from four islands in Vanuatu, and for 5 chickens from Guam in the Marianas.

Results

Forty samples from Vanuatu and three from Guam were assigned to haplogroup D, supporting this as a Pacific chicken haplogroup that persists in the Western Pacific. Two haplogroup E lineages were observed in Guam and two in Vanuatu. Of the E lineages in Vanuatu, one was identical to prehistoric Vanuatu and Polynesian samples and the other differed by one polymorphism. Contrary to our expectations, we observed few globally distributed domesticate lineages not associated with Pacific chicken dispersal. This might suggest less European introgression of chickens into Vanuatu than expected. If so, the E lineages might represent lineages maintained from ancient Pacific chicken introductions. The Vanuatu sample might thus provide an opportunity to distinguish between maintained ancestral Pacific chicken lineages and replacement by global domesticates through genomic analyses, which could resolve questions of contemporary haplogroup E chicken relationships and inform interpretations of debated sequences from archaeological samples.     

Genomic Characterisation of Invasive Non-Typhoidal Salmonella enterica Subspecies enterica Serovar Bovismorbificans Isolates from Malawi

Background

Invasive Non-typhoidal Salmonella (iNTS) are an important cause of bacteraemia in children and HIV-infected adults in sub-Saharan Africa. Previous research has shown that iNTS strains exhibit a pattern of gene loss that resembles that of host adapted serovars such as Salmonella Typhi and Paratyphi A. Salmonella enterica serovar Bovismorbificans was a common serovar in Malawi between 1997 and 2004.

Methodology

We sequenced the genomes of 14 Malawian bacteraemia and four veterinary isolates from the UK, to identify genomic variations and signs of host adaptation in the Malawian strains.

Principal Findings

Whole genome phylogeny of invasive and veterinary S. Bovismorbificans isolates showed that the isolates are highly related, belonging to the most common international S. Bovismorbificans Sequence Type, ST142, in contrast to the findings for S. Typhimurium, where a distinct Sequence Type, ST313, is associated with invasive disease in sub-Saharan Africa. Although genome degradation through pseudogene formation was observed in ST142 isolates, there were no clear overlaps with the patterns of gene loss seen in iNTS ST313 isolates previously described from Malawi, and no clear distinction between S. Bovismorbificans isolates from Malawi and the UK.The only defining differences between S. Bovismorbificans bacteraemia and veterinary isolates were prophage-related regions and the carriage of a S. Bovismorbificans virulence plasmid (pVIRBov).

Conclusions

iNTS S. Bovismorbificans isolates, unlike iNTS S. Typhiumrium isolates, are only distinguished from those circulating elsewhere by differences in the mobile genome. It is likely that these strains have entered a susceptible population and are able to take advantage of this niche. There are tentative signs of convergent evolution to a more human adapted iNTS variant. Considering its importance in causing disease in this region, S. Bovismorbificans may be at the beginning of this process, providing a reference against which to compare changes that may become fixed in future lineages in sub-Saharan Africa.     

Combining FISH and model-based predictions to understand chromosome evolution in Typhonium (Araceae)

Background and Aims

Since the advent of molecular phylogenetics, numerous attempts have been made to infer the evolutionary trajectories of chromosome numbers on DNA phylogenies. Ideally, such inferences should be evaluated against cytogenetic data. Towards this goal, we carried out phylogenetic modelling of chromosome number change and fluorescence in situ hybridization (FISH) in a medium sized genus of Araceae to elucidate if data from chromosomal markers would support maximum likelihood-inferred changes in chromosome numbers among close relatives. Typhonium, the focal genus, includes species with 2n = 65 and 2n = 8, the lowest known count in the family.

Methods

A phylogeny from nuclear and plastid sequences (96 taxa, 4252 nucleotides) and counts for all included species (15 of them first reported here) were used to model chromosome number evolution, assuming discrete events, such as polyploidization and descending or ascending dysploidy, occurring at different rates. FISH with three probes (5S rDNA, 45S rDNA and Arabidopsis-like telomeres) was performed on ten species with 2n = 8 to 2n = 24.

Key Results

The best-fitting models assume numerous past chromosome number reductions. Of the species analysed with FISH, the two with the lowest chromosome numbers contained interstitial telomeric signals (Its), which together with the phylogeny and modelling indicates decreasing dysploidy as an explanation for the low numbers. A model-inferred polyploidization in another species is matched by an increase in rDNA sites.

Conclusions

The combination of a densely sampled phylogeny, ancestral state modelling and FISH revealed that the species with n = 4 is highly derived, with the FISH data pointing to a Robertsonian fusion-like chromosome rearrangement in the ancestor of this species.     

Paleogene radiation of a plant pathogenic mushroom

Background

The global movement and speciation of fungal plant pathogens is important, especially because of the economic losses they cause and the ease with which they are able to spread across large areas. Understanding the biogeography and origin of these plant pathogens can provide insights regarding their dispersal and current day distribution. We tested the hypothesis of a Gondwanan origin of the plant pathogenic mushroom genus Armillaria and the currently accepted premise that vicariance accounts for the extant distribution of the species.

Methods

The phylogeny of a selection of Armillaria species was reconstructed based on Maximum Parsimony (MP), Maximum Likelihood (ML) and Bayesian Inference (BI). A timeline was then placed on the divergence of lineages using a Bayesian relaxed molecular clock approach.

Results

Phylogenetic analyses of sequenced data for three combined nuclear regions provided strong support for three major geographically defined clades: Holarctic, South American-Australasian and African. Molecular dating placed the initial radiation of the genus at 54 million years ago within the Early Paleogene, postdating the tectonic break-up of Gondwana.

Conclusions

The distribution of extant Armillaria species is the result of ancient long-distance dispersal rather than vicariance due to continental drift. As these finding are contrary to most prior vicariance hypotheses for fungi, our results highlight the important role of long-distance dispersal in the radiation of fungal pathogens from the Southern Hemisphere.     

Does polyploidy facilitate long-distance dispersal?

Background and Aims

The ability of plant lineages to reach all continents contributes substantially to their evolutionary success. This is exemplified by the Poaceae, one of the most successful angiosperm families, in which most higher taxa (tribes, subfamilies) have global distributions. Due to the old age of the ocean basins relative to the major angiosperm radiations, this is only possible by means of long-distance dispersal (LDD), yet the attributes of lineages with successful LDD remain obscure. Polyploid species are over-represented in invasive floras and in the previously glaciated Arctic regions, and often have wider ecological tolerances than diploids; thus polyploidy is a candidate attribute of successful LDD.

Methods

The link between polyploidy and LDD was explored in the globally distributed grass subfamily Danthonioideae. An almost completely sampled and well-resolved species-level phylogeny of the danthonioids was used, and the available cytological information was assembled. The cytological evolution in the clade was inferred using maximum likelihood (ML) as implemented in ChromEvol. The biogeographical evolution in the clade was reconstructed using ML and Bayesian approaches.

Key Results

Numerous increases in ploidy level are demonstrated. A Late Miocene–Pliocene cycle of polyploidy is associated with LDD, and in two cases (the Australian Rytidosperma and the American Danthonia) led to secondary polyploidy. While it is demonstrated that successful LDD is more likely in polyploid than in diploid lineages, a link between polyploidization events and LDD is not demonstrated.

Conclusions

The results suggest that polyploids are more successful at LDD than diploids, and that the frequent polyploidy in the grasses might have facilitated the extensive dispersal among continents in the family, thus contributing to their evolutionary success.     

Baltic salmon, Salmo salar, from Swedish river Lule älv is more resistant to furunculosis compared to rainbow trout

Background

Furunculosis, caused by Aeromonas salmonicida, continues to be a major health problem for the growing salmonid aquaculture. Despite effective vaccination programs regular outbreaks occur at the fish farms calling for repeated antibiotic treatment. We hypothesized that a difference in natural susceptibility to this disease might exist between Baltic salmon and the widely used rainbow trout.

Study Design

A cohabitation challenge model was applied to investigate the relative susceptibility to infection with A. salmonicida in rainbow trout and Baltic salmon. The course of infection was monitored daily over a 30-day period post challenge and the results were summarized in mortality curves.

Results

A. salmonicida was recovered from mortalities during the entire test period. At day 30 the survival was 6.2% and 34.0% for rainbow trout and Baltic salmon, respectively. Significant differences in susceptibility to A. salmonicida were demonstrated between the two salmonids and hazard ratio estimation between rainbow trout and Baltic salmon showed a 3.36 higher risk of dying from the infection in the former.

Conclusion

The finding that Baltic salmon carries a high level of natural resistance to furunculosis might raise new possibilities for salmonid aquaculture in terms of minimizing disease outbreaks and the use of antibiotics.     

Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus

Background

The mariner family of transposable elements is one of the most widespread in the Metazoa. It is subdivided into several subfamilies that do not mirror the phylogeny of these species, suggesting an ancient diversification. Previous hybridization and PCR studies allowed a partial survey of mariner diversity in the Metazoa. In this work, we used a comparative genomics approach to access the genus-wide diversity and evolution of mariner transposable elements in twenty Drosophila sequenced genomes.

Results

We identified 36 different mariner lineages belonging to six distinct subfamilies, including a subfamily not described previously. Wide variation in lineage abundance and copy number were observed among species and among mariner lineages, suggesting continuous turn-over. Most mariner lineages are inactive and contain a high proportion of damaged copies. We showed that, in addition to substitutions that rapidly inactivate copies, internal deletion is a major mechanism contributing to element decay and the generation of non-autonomous sublineages. Hence, 23% of copies correspond to several Miniature Inverted-repeat Transposable Elements (MITE) sublineages, the first ever described in Drosophila for mariner. In the most successful MITEs, internal deletion is often associated with internal rearrangement, which sheds light on the process of MITE origin. The estimation of the transposition rates over time revealed that all lineages followed a similar progression consisting of a rapid amplification burst followed by a rapid decrease in transposition. We detected some instances of multiple or ongoing transposition bursts. Different amplification times were observed for mariner lineages shared by different species, a finding best explained by either horizontal transmission or a reactivation process. Different lineages within one species have also amplified at different times, corresponding to successive invasions. Finally, we detected a preference for insertion into short TA-rich regions, which appears to be specific to some subfamilies.

Conclusions

This analysis is the first comprehensive survey of this family of transposable elements at a genus scale. It provides precise measures of the different evolutionary processes that were hypothesized previously for this family based on PCR data analysis. mariner lineages were observed at almost all “life cycle” stages: recent amplification, subsequent decay and potential (re)-invasion or invasion of genomes.

Electronic supplementary material

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

Control region variability of haplogroup C1d and the tempo of the peopling of the Americas

Background

Among the founding mitochondrial haplogroups involved in the peopling of the Americas, haplogroup C1d has been viewed as problematic because of its phylogeny and because of the estimates of its antiquity, apparently being much younger than other founding haplogroups. Several recent analyses, based on data from the entire mitochondrial genome, have contributed to an advance in the resolution of these problems. The aim of our analysis is to compare the conclusions drawn from the available HVR-I and HVR-II data for haplogroup C1d with the ones based on whole mitochondrial genomes.

Methodology/Principal Findings

HVR-I and HVR-II sequences defined as belonging to haplogroup C1d by standard criteria were gathered from the literature as well as from population studies carried out in Uruguay. Sequence phylogeny was reconstructed using median-joining networks, geographic distribution of lineages was analyzed and the age of the most recent common ancestor estimated using the ρ-statistic and two different mutation rates. The putative ancestral forms of the haplogroup were found to be more widespread than the derived lineages, and the lineages defined by np 194 were found to be widely distributed and of equivalent age.

Conclusions/Significance

The analysis of control region sequences is found to still harbor great potential in tracing microevolutionary phenomena, especially those found to have occurred in more recent times. Based on the geographic distributions of the alleles of np 7697 and np 194, both discussed as possible basal mutations of the C1d phylogeny, we suggest that both alleles were part of the variability of the haplogroup at the time of its entrance. Moreover, based on the mutation rates of the different sites stated to be diagnostic, it is possible that the anomalies found when analyzing the haplogroup are due to paraphyly.