Christophe J. Douady

Universal trees based on large combined protein sequence data sets

Universal trees of life based on small-subunit (SSU) ribosomal RNA (rRNA) support the separate mono/holophyly of the domains Archaea (archaebacteria), Bacteria (eubacteria) and Eucarya (eukaryotes) and the placement of extreme thermophiles at the base of the Bacteria. The concept of universal tree reconstruction recently has been upset by protein trees that show intermixing of species from different domains. Such tree topologies have been attributed to either extensive horizontal gene transfer or degradation of phylogenetic signals because of saturation for amino acid substitutions. Here we use large combined alignments of 23 orthologous proteins conserved across 45 species from all domains to construct highly robust universal trees. Although individual protein trees are variable in their support of domain integrity, trees based on combined protein data sets strongly support separate monophyletic domains. Within the Bacteria, we placed spirochaetes as the earliest derived bacterial group. However, elimination from the combined protein alignment of nine protein data sets, which were likely candidates for horizontal gene transfer, resulted in trees showing thermophiles as the earliest evolved bacterial lineage. Thus, combined protein universal trees are highly congruent with SSU rRNA trees in their strong support for the separate monophyly of domains as well as the early evolution of thermophilic Bacteria.

Phylogeography of a subterranean amphipod reveals cryptic diversity and dynamic evolution in extreme environments

Extreme conditions in subsurface are suspected to be responsible for morphological convergences, and so to bias biodiversity assessment. Subterranean organisms are also considered as having poor dispersal abilities that in turn generate a large number of endemic species when habitat is fragmented. Here we test these general hypotheses using the subterranean amphipod Niphargus virei. All our phylogenetic analyses (Bayesian, maximum likelihood and distance), based on two independent genes (28S and COI), revealed the same tripartite structure. N. virei populations from Benelux, Jura region and the rest of France appeared as independent evolutionary units. Molecular rates estimated via global or Bayesian relaxed clock suggest that this split is at least 13 million years old and accredit the cryptic diversity hypothesis. Moreover, the geographical distribution of these lineages showed some evidence of recent dispersal through apparent vicariant barrier. In consequence, we argue that future analyses of evolution and biogeography in subsurface, or more generally in extreme environments, should consider dispersal ability as an evolving trait and morphology as a potentially biased marker.

Molecular phylogenetic evidence refuting the hypothesis of Batoidea (rays and skates) as derived sharks

Early morphological studies regarding the evolutionary history of elasmobranchs suggested sharks and batoids (skates and rays) were respectively monophyletic. More modern morphological cladistic studies, however, have tended to suggest that batoids are derived sharks, closely related to sawsharks and angelsharks, a phylogenetic arrangement known as the Hypnosqualea hypothesis. Very few molecular studies addressing interordinal relationships of elasmobranchs have been published; the few that do exist, are very limited in terms of both taxon representation and/or aligned sequence positions, and are insufficient to answer the question of whether batoids are derived sharks. The purpose of this study was to address this issue with more complete taxon representation, concomitant with a reasonable number of aligned sequence positions. The data set included a 2.4-kb segment of the mitochondrial 12S rRNA-tRNA valine-16S rRNA locus, and in terms of taxa, representatives of two orders of Batoidea, at least one representative of all orders of sharks, and as an outgroup, the widely recognized sister group to elasmobranchs-Holocephali. The results provide the first convincing molecular evidence for shark monophyly and the rejection of the Hypnosqualea hypothesis. Our phylogenetic placement of batoids as a basal elasmobranch lineage means that much of the current thinking regarding the evolution of morphological and life history characteristics in elasmobranchs needs to be re-evaluated.

The Sahara as a vicariant agent, and the role of Miocene climatic events, in the diversification of the mammalian order Macroscelidea (elephant shrews).

Although the Sahara is a major geographical feature of the African continent, its role in the diversification of animal species is not well understood. We present here a molecular phylogeny for members of the endemic African mammalian order Macroscelidea (elephant shrews) with molecular-clock calculations; this molecular phylogeny provides convincing evidence that the genus Elephantulus is diphyletic. Elephantulus rozeti, the only elephant shrew species that resides north of the Sahara, is the sister group of a species from a different genus (Petrodromus tetradactylus), which resides just south of the Sahara. The split between these taxa coincided with major Miocene climatic events, which triggered the cooling and aridification of midlatitude continental regions, and a shift in the Sahara from a tropical to an arid environment. Thus, the North African distribution of E. rozeti is not the result of dispersion from an eastern species of the genus, but instead the result of a vicariant event involving the formation of the Sahara. The splitting events involved with most Elephantulus species in our analysis appear to coincide with these climatic events. This coincidence suggests that the environmental consequences associated with this period played an important role in the radiation of this order of mammals. The strongly supported phylogeny provides compelling evidence for a complex history of mosaic evolution, including pronounced bradytelic morphological evolution in some lineages, accelerated morphological evolution in others, and a remarkably slow rate of evolution of the male reproductive structure.

Plant resistance to mechanical stress: evidence of an avoidance–tolerance trade-off

External mechanical forces resulting from the pressure exerted by wind or water movement are a major stress factor for plants and may cause regular disturbances in many ecosystems. A plants ability to resist these forces relies either on minimizing the forces encountered by the plant (avoidance strategy), or on maximizing its resistance to breakage (tolerance strategy). We investigated plant resistance strategies using aquatic vegetation as a model, and examined whether avoidance and tolerance are negatively correlated. We tested the avoidance-tolerance correlation across 28 species using a phylogenetically corrected analysis, after construction of a molecular phylogeny for the species considered. Different species demonstrated contrasting avoidance and tolerance and we demonstrated a significant negative relationship between the two strategies, which suggests an avoidance-tolerance trade-off. Negative relationships may result from costs that each strategy incurs or from constraints imposed by physical laws on plant tissues. The existence of such a trade-off has important ecological and evolutionary consequences. It would lead to constraints on the evolution and variation of both strategies, possibly limiting their evolution and may constrain many morphological, anatomical and architectural traits that underlie avoidance and tolerance.

Intragenomic Heterogeneity and Intergenomic Recombination among Haloarchaeal rRNA Genes

More than one copy of rRNA operons, which code for both the small-subunit (SSU) and large-subunit (LSU) rRNA, are often found in prokaryotes. It is generally assumed that all rRNA operons within a single cell are almost identical. A notable exception is the extremely halophilic archaeal genus Haloarcula, most species of which are known to harbor highly divergent rRNA operons that differ at approximately 5% of the nucleotide positions in the SSU gene and at 1 to 2% of the nucleotide positions in the LSU gene. We report that such intragenomic heterogeneity is not unique to Haloarcula, as high levels of intragenomic sequence variation have been observed for the SSU genes of two other genera of extreme halophiles, Halosimplex and Natrinema. To investigate this in detail, the two rRNA operons of Halosimplex carlsbadense and the four operons of Natrinema sp. strain XA3-1 were cloned and completely sequenced. The SSU and LSU genes of H. carlsbadense show the highest levels of intragenomic heterogeneity observed so far in archaea (6.7 and 2.6%). The operons of Natrinema sp. strain XA3-1 have additional unusual characteristics, such as identical internal transcribed spacers, while one of four SSU genes is 5% divergent and all LSU genes differ from each other by 0.9 to 1.9%. The heterogeneity among the Natrinema sp. strain XA3-1 LSU genes is localized in hot spots, and one of these regions is shown to be the result of a recombination event with a distantly related halophile. This is the first example of interspecies recombination between rRNA genes in archaea, and the recombination occurred over one of the largest phylogenetic distances ever reported for such an event. We suggest that intragenomic heterogeneity of rRNA operons is an ancient and stable trait in several lineages of the Halobacteriales. The impact of this phenomenon on the taxonomy of extremely halophilic archaea is discussed.

A streamlined, bi-organelle, multiplex PCR approach to species identification: Application to global conservation and trade monitoring of the great white shark, Carcharodon carcharias

The great white shark, Carcharodoncarcharias, is the most widely protectedelasmobranch in the world, and is classified asVulnerable by the IUCN and listed on AppendixIII of CITES. Monitoring of trade in whiteshark products and enforcement of harvest andtrade prohibitions is problematic, however, inlarge part due to difficulties in identifyingmarketed shark parts (e.g., dried fins, meatand processed carcasses) to species level. Toaddress these conservation and managementproblems, we have developed a rapid, moleculardiagnostic assay based on species-specific PCRprimer design for accurate identification ofwhite shark body parts, including dried fins. The assay is novel in several respects: Itemploys a multiplex PCR assay utilizing bothnuclear (ribosomal internal transcribed spacer2) and mitochondrial (cytochrome b) locisimultaneously to achieve a highly robustmeasure of diagnostic accuracy; it is verysensitive, detecting the presence of whiteshark DNA in a mixture of genomic DNAs from upto ten different commercially fished sharkspecies pooled together in a single PCR tube;and it successfully identifies white shark DNAfrom globally distributed animals. Inaddition to its utility for white shark trademonitoring and conservation applications, thishighly streamlined, bi-organelle, multiplex PCRassay may prove useful as a general model forthe design of genetic assays aimed at detectingbody parts from other protected and threatenedspecies.

Molecular Evidence for the Major Clades of Placental Mammals

Higher-level relationships among placental mammals, as well as the historical biogeography of this group against the backdrop of continental fragmentation and reassembly, remain poorly understood. Here, we analyze two independent molecular data sets that represent all placental orders. The first data set includes six genes (A2AB, IRBP, vWF, 12S rRNA, tRNA valine, 16S rRNA; total = 5.71 kb) for 26 placental taxa and two marsupials; the second data set includes 2.95 kb of exon 11 of the BRCA1 gene for 51 placental taxa and four marsupials. We also analyzed a concatenation of these data sets (8.66 kb) for 26 placentals and one marsupial. Unrooted and rooted analyses were performed with parsimony, distance methods, maximum likelihood, and a Bayesian approach. Unrooted analyses provide convincing support for a fundamental separation of placental orders into groups with southern and northern hemispheric origins according to the current fossil record. On rooted trees, one or both of these groups are monophyletic depending on the position of the root. Maximum likelihood and Bayesian analyses with the BRCA1 and combined 8.66 kb data sets provide strong support for the monophyly of the northern hemisphere group (Boreoeutheria). Boreoeutheria is divided into Laurasiatheria (Carnivora + Cetartiodactyla + Chiroptera + Eulipotyphla + Perissodactyla + Pholidota) and Euarchonta (Dermoptera + Primates + Scandentia) + Glires (Lagomorpha + Rodentia). The southern hemisphere group is either monophyletic or paraphyletic, depending on the method of analysis used. Within this group, Afrotheria (Proboscidea + Sirenia + Hyracoidea + Tubulidentata + Macroscelidea + Afrosoricida) is monophyletic. A unique nine base-pair deletion in exon 11 of the BRCA1 gene also supports Afrotheria monophyly. Given molecular dates that suggest that the southern hemisphere group and Boreoeutheria diverged in the Early Cretaceous, a single trans-hemispheric dispersal event may have been of fundamental importance in the early history of crown-group Eutheria. Parallel adaptive radiations have subsequently occurred in the four major groups: Laurasiatheria, Euarchonta + Glires, Afrotheria, and Xenarthra.

Phylogeography of a Nematode (Heligmosomoides polygyrus) in the Western Palearctic region : Persistence of Northern cryptic populations during ice ages?

This study establishes the continental phylogeographical pattern of a European nematode, Heligmosomoides polygyrus (Dujardin, 1845; Heligmosomoidea). We sequenced 687 base pairs of the mitochondrial DNA (mtDNA) cyt b gene for 136 individuals collected in 22 localities. The results revealed that H. polygyrus populations are separated into five major units corresponding to the Italian, northern European (Denmark and Ireland), Iberian, western European, and Balkan populations. Different subclades were also observed within the first two groups. Based on the rate of molecular evolution of H. polygyrus cyt b gene‐estimated to 3.5%–3.7% divergence per million years (Myr) in a previous study – the isolation time of the five clades was estimated between 2.5 ± 0.24 and 1.5 ± 0.23 million years ago. Moreover, H. polygyrus presents a higher genetic variability in the Mediterranean peninsulas as compared to northwestern Europe, highlighting the role of these regions as refuge areas. Like its specific host, the wood mouse Apodemus sylvaticus, H. polygyrus’ pattern of postglacial recolonization of northwestern Europe was initiated from Iberian populations, while Italian and Balkan populations did not expand to the north. The results also suggest the existence of forested and temperate refuges in the southern British Isles during the Quaternary. Finally, the genetic diversity as well as the level of genetic divergence between the lineages of H. polygyrus are compared to those observed in other vertebrate and invertebrate phylogeographical studies: the existence of highly differentiated lineages in H. polygyrus (5%−10% of genetic divergence) highlights that the effects of Pleistocene climate changes on free‐living organisms are also reflected in their obligate parasites.

Preventing the pollution of mitochondrial datasets with nuclear mitochondrial paralogs (numts)

Molecular tools have become prominent in ecology and evolution. A target of choice for molecular ecologists and evolutionists is mitochondrial DNA (mtDNA), whose many advantages have also convinced broad-scale, pragmatic programmes such as barcode initiatives. Of course, mtDNA is also of interest to human geneticists investigating mitochondrial diseases. Studies using mtDNA are however put at great risk by the inadvertent co-amplification or preferred amplification of nuclear pseudogenes (numts). A posteriori analysis of putative mtDNA sequences can help in removing numts but faces severe limitations (e.g. recently translocated numts will most of the time go unnoticed). Counter-measures taken a priori, i.e. explicitly designed for avoiding numt co-amplification or preferred amplification, are appealing but have never been properly assessed. Here we investigate the efficiency of four such measures (mtDNA enrichment, cDNA amplification, long-range amplification and pre-PCR dilution) on a common set of numt cases, showing that mtDNA enrichment is the worst performer while the use of pre-PCR dilution is a simple, yet robust method to prevent the pollution of putative mtDNA datasets with numts. Therefore, straightforward recommendations can be made that, if followed, will considerably increase the confidence in the mitochondrial origin of any mtDNA-like sequence.