Phylogenetic study of complete cytochrome b genes in musk deer (genus Moschus) using museum samples

As an endangered animal group, musk deer (genus Moschus) are not only a great concern of wildlife conservation, but also of special interest to evolutionary studies due to long-standing arguments on the taxonomic and phylogenetic associations in this group. Using museum samples, we sequenced complete mitochondrial cytochrome b genes (1140 bp) of all suggested species of musk deer in order to reconstruct their phylogenetic history through molecular information. Our results showed that the cytochrome b gene tree is rather robust and concurred for all the algorithms employed (parsimony, maximum likelihood, and distance methods). Further, the relative rate test indicated a constant sequence substitution rate among all the species, permitting the dating of divergence events by molecular clock. According to the molecular topology, M. moschiferus branched off the earliest from a common ancestor of musk deer (about 700,000 years ago); then followed the bifurcation forming the M. berezouskii lineage and the lineage clustering M. fuscus, M. chrysogaster, and M. leucogaster (around 370,000 years before present), interestingly the most recent speciation event in musk deer happened rather recently (140,000 years ago), which might have resulted from the diversified habitats and geographic barriers in southwest China caused by gigantic movements of the Qinghai-Tibetan Plateau in history. Combining the data of current distributions, fossil records, and molecular data of this study, we suggest that the historical dispersion of musk deer might be from north to south in China. Additionally, in our further analyses involving other pecora species, musk deer was strongly supported as a monophyletic group and a valid family in Artiodactyla, closely related to Cervidae. (C) 1999 Academic Press.

Phylogenetic relationships of five pika species (genus Ochotona) based on mitochondrial DNA restriction maps

Restriction site mapping of mitochondrial DNA (mtDNA) with 16 restriction endonucleases was used to examine the phylogenetic relationships of Ochotona cansus, O. huangensis, O. thibetana, O. curzoniae and O. erythrotis. A 1-kb length variation between O. erythrotis of subgenus Pika and other four species of subgenus Ochotona was observed, which may be a useful genetic marker for identifying the two subgenera. The phylogenetic tree constructed using PAUP based on 61 phylogenetically informative sites suggests that O. erythrotis diverged first, followed by O. cansus, while O. curzoniae and O. huangensis are sister taxa related to O. thibetana, The results indicate that both O. cansus and O. huangensis should be treated as independent species. If the base substitution rate of pikas mtDNA was 2% per million years, then the divergence time of the two subgenera, Pika and Ochotana, is about 8.8 Ma ago of late Miocence, middle Bao-dian of Chinese mammalian age, and the divergence of the four species in subgenus Ochotona would have occurred about 2.5 - 4.2 Ma ago, Yushean of Chinese mammalian age. This calculation appears to be substantiated by the fossil record.

Phylogeny and evolution of the Drosophila nasuta subgroup based on mitochondrial ND4 and ND4L gene sequences

The sequences of the mitochondrial ND4 gene (1339 bp) and the ND4L gene (290 bp) were determined for all the 14 extant taxa of the Drosophila nasuta subgroup The average A + T content of ND4 genes is 76.5% and that of ND4L genes is 83.5%. A total of 114 variable sites were scored. The ND4 gene sequence divergence ranged from 0 to 5.4% within the subgroup. The substitution rate of the ND4 gene is about 1.25% per million years. The base substitution of the genesis strongly transition biased. Neighbor-joining and parsimony were used to construct a phylogeny based on the resultant sequence data set. According to these trees, five, distinct mtDNA clades can be identified. D. niveifrons represents the most diverged lineage. D, sulfurigaster bilimbata and D. kepulauana form two independent lineages. The other two clades are the kohkoa complex and the albomicans complex. The Kohkoa complex consists of D. sulfurigaster sulfurigaster, D. pulaua, D. kohkoa, and Taxon-F. The albomicans complex can be divided into two groups: D. nasuta, D. sulfurigaster neonasuta, D. sulfurigaster albostrigata, and D.. albomicans from Chiangmai form one group; and D. pallidifrons, Taxon-I, Taxon-J, and D. albomicans from China form the other group. High genetic differentiation was found among D. albomicans populations. Based on our phylogenetic results, we hypothesize that D. niveifrons diverged first from the D, nasuta subgroup in Papua New Guinea about 3.5 Mya. The ancestral population spread to the north and when it reached Borneo, it diversified sequentially into the kohkoa complex, D. s. bilimbata, and D. kepulauana. About 1 Mya, another radiation occurred when the ancestral populations reached the Indo-China Peninsula, forming the albomicans complex. Discrepancy between morphological groupings and phylogenetic results suggests that the male morphological traits may not be orthologous. (C) 1999 Academic Press.

Episodic evolution of prolactin receptor gene in mammals: coevolution with its ligand

Divergence of proteins in signaling pathways requires ligand and receptor coevolution to maintain or improve binding affinity and/or specificity. In this paper we show a clear case of coevolution between the prolactin (PRL) gene and its receptor (prolactin receptor, PRLR) in mammals. First we observed episodic evolution of the extracellular and intracellular domains of the PRLR, which is closely consistent with that seen in PRL. Correlated evolution was demonstrated both between PRL and its receptor and between the two domains of the PRLR using Pearson's correlation coefficient. On comparing the ratio of the nonsynonymous substitution rate to synonymous substitution rate (omega=d(N)/d(S)) for each branch of the star phylogeny of mammalian PRLRs, separately for the extracellular domain (ECD) and the transmembrane domain/intracellular domain (TMD/ICD), we observed a lower omega ratio for ECD than TMD/ICD along those branches leading to pig, dog and rabbit but a higher ratio for ECD than TMD/ICD on the branches leading to primates, rodents and ruminants, on which bursts of rapid evolution were observed. These observations can be best explained by coevolution between PRL and its receptor and between the two domains of the PRLR.

Evolutionary implications of multiple SINE insertions in an intronic region from diverse mammals

An analysis of the nuclear beta-fibrinogen intron 7 locus from 30 taxa representing 12 placental orders of mammals reveals the enriched occurrences of short interspersed clement (SINE) insertion events. Mammalian-wide interspersed repeats (MIRs) are present at orthologous sites of all examined species except those in the order Rodentia. The higher substitution rate in mouse and a rare MIR deletion from rat account for the absence of MIR in the rodents. A minimum of five lineage-specific SINE sequences are also found to have independently inserted into this intron in Carnivora, Artiodactyla and Lagomorpha. In the case of Carnivora, the unique amplification pattern of order-specific CAN SINE provides important evidence for the "pan-carnivore" hypothesis of this repeat element and reveals that the CAN SINE family may still be active today. Particularly interesting is the finding that all identified lineage-specific SINE elements show a strong tendency to insert within or in very close proximity to the preexisting MIRs for their efficient integrations, suggesting that the MIR clement is a hot spot for successive insertions of other SINEs. The unexpected MIR excision as a result of a random deletion in the rat intron locus and the non-random site targeting detected by this study indicate that SINEs actually have a greater insertional flexibility and regional specificity than had previously been recognized. Implications for SINE sequence evolution upon and following integration, as well as the fascinating interactions between retroposons and the host genomes are discussed.

Episodic evolution of prolactin gene in primates

In the present study, we obtained exon 2-5 of prolactin (PRL) gene from four primate species by PCR and sequencing. Adding other genes available in GenBank, we calculate amino acid substitution rates for prolactin gene in primate. Comparison of nonsynonymous substitution rate to synonymous substitution rate ratios shows no evidence of positive selection for any lineage of primate prolactin gene. According to this and the facts that (i) no sites under positive selection are inferred by using maximum-likelihood method; (ii) among 32 amino acid replacement that occurred along the rapid evolutionary phase, only two are included in the 40 functionally important residues, indicating that amino acid replacement tends to occur in those functionally unimportant residues; (iii) partial of prolactin function is replaced by placental lactogen in primate at the rapid evolutionary phase of prolactin gene, we thus deem that it is relaxation of purifying selection to some extent rather than positive selection that enforces the rapid evolution of primate prolactin gene.

Origin and evolution of new exons in rodents

Gene number difference among organisms demonstrates that new gene origination is a fundamental biological process in evolution. Exon shuffling has been universally observed in the formation of new genes. Yet to be learned are the ways new exons originate and evolve, and how often new exons appear. To address these questions, we identified 2695 newly evolved exons in the mouse and rat by comparing the expressed sequences of 12,419 orthologous genes between human and mouse, using 743,856 pig ESTs as the outgroup. The new exon origination rate is about 2.71 x 10(-3) per gene per million years. These new exons have markedly accelerated rates both of nonsynonymous substitutions and of insertions/ deletions (indels). A much higher proportion of new exons have Kappa(a)/Kappa(s) ratios > 1 (where K-a is the nonsynonymous substitution rate and K-s is the synonymous substitution rate) than K do the old exons shared by human and mouse, implying a role of positive selection in the rapid evolution. The majority of these new exons have sequences unique in the genome, suggesting that most new exons might originate through "exonization" of intronic sequences. Most of the new exons appear to be alternative exons that are expressed at low levels.

Sox genes evolution in closely related young tetraploid cyprinid fishes and their diploid relative

Previous study and analysis of cytochrome b suggested that polyploidization event in the genus Tor occurred about 10 Mya ago. In order to understand evolutionary fates of Sox gene in the early stage of genome duplication at the nucleotide level, PCR surveys for Sox genes in three closely related cyprinid fishes T douronensis (2n = 100), T qiaojiensis (2n = ?), T sinensis (2n = 100) and their relative T brevifilis (2n = 50) were performed. Totally, 52 distinct Sox genes were obtained in these four species, representing SoxB, SoxC, and SoxE group. As expected, isoforms of some Sox genes correspond with the ploidy of species, such as two copies of Sox9a exist in tetraploid species. Analysis indicated that duplicated Sox gene pairs caused by polyploidization evolved independently of each other within polyploid species. Results of substitution rate showed nearly equal rate of nonsynonymous substitution of duplicated Sox orthologs among different polyploid species and their diploid relative orthologs, suggesting at the early stage of genome duplicated Sox orthologs are under similar selective constraints in different polyploidy species and their diploid relative at the amino acid level. All PCR fragments of Sox genes obtained in this study are not accompanied by obvious increase in mutations and pseudogene formation which means that they are under strong purifying selection, suggesting that they are functional at the DNA level. Cenealogical analysis revealed that T qiaojiensis was tetraploid, and T douronensis, T qiaojiensis as well as T sinensis had an allotetraploid ancestor. (C) 2009 Elsevier B.V. All rights reserved.

Phylogeny and biogeography of Chinese sisorid catfishes re-examined using mitochondrial cytochrome b and 16S rRNA gene sequences

The family Sisoridae is one of the largest and most diverse Asiatic catfish families, most species occurring in the water systems of the Qinhai-Tibetan Plateau and East Himalayas. To date published morphological and molecular phylogenetics hypotheses of sisorid catfishes are part congruent, and there are some areas of significant disagreement with respect to intergeneric relationships. We used mitochondrial cytochrome b and 16S rRNA gene sequences to clarify existing gaps in phylogenetics and to test conflicting vicariant and dispersal biogeographical hypotheses of Chinese sisorids using dispersal-vicariance analysis and weighted ancestral area analysis in combination with palaeogeographical data as well as molecular clock calibration. Our results suggest that: (1) Chinese sisorid catfishes form a monophyletic group with two distinct clades, one represented by (Gagata (Bagarius, Glyptothorax)) and the other by (glyptosternoids, Pseudecheneis); (2) the glyptosternoid is a monophyletic group and Glyptosternum, Glaridoglanis, and Exostoma are three basal species having a primitive position among it; (3) a hypothesis referring to Pseudecheneis as the sister group of the glyptosternoids, based on morphological evidence, is supported; (4) the genus Pareuchiloglanis, as presently defined, is not monophyletic; (5) congruent with previous hypotheses, the uplift of Qinghai-Tibetan Plateau played a primary role in the speciation and radiation of the Chinese sisorids; and (6) an evolutionary scenario combining aspects of both vicariance and dispersal theory is necessary to explain the distribution pattern of the glyptosternoids. In addition, using a cytochrome b substitution rate of 0.91% per million years and 0.23% for 16S rRNA, we tentatively date that the glyptosternoids most possibly originated in Oligocene-Miocene boundary (19-24Myr), and radiated from Miocene to Pleistocene, along with a center of origin in the Irrawaddy-Tsangpo drainages and several rapid speciation in a relatively short time. (c) 2005 Elsevier Inc. All rights reserved.

Mitochondrial 16S rRNA sequence variations and phylogeny of the Chinese sisorid catfishes

Partial sequences of mitochondrial 16S rRNA gene were obtained by PCR amplification for comparisons among nine species of glyptosternoid fishes and six species of non-glyptosternoids representing 10 sisorid genera. There are compositional biases in the A-rich impaired regions and G-rich paired regions. A-G transitions are primarily responsible for the Ts/Tv bias in impaired regions. The overall substitution rate in impaired regions is almost two times higher than that in the paired regions. Saturation plots at comparable levels of sequence divergence demonstrate no saturation effects. Phylogenetic analyses using both maximum likelihood and Bayesian methods support the monophyly of Sisoridae. Chinese sisorid catfishes are composed of two major lineages, one represented by (Gagata (Bagarius, Glyptothorax)) and the other by "glyptosternoids + Pseudecheneis". The glyptosternoids may not be a monophyletic group. A previous hypothesis referring to Pseudecheneis as the sister group of monophyletic glyptosternoids, based on morphological evidence, is not supported by the molecular data. Pseudecheneis is shown to be a sister taxon of Glaridoglanis. Pareuchiloglanis might be paraphyletic with Pseudexostoma and Euchiloglanis. Our results also support the hypothesis that Pareuchiloglanis anteanalis might be considered as the synonyms of Pareuchiloglanis sinensis, and genus Euchiloglanis might have only one valid species, Euchiloglanis davidi.

Phylogenetic relationships of five pika species (genus Ochotona) based on mitochondrial DNA restriction maps

Restriction site mapping of mitochondrial DNA (mtDNA) with 16 restriction endonucleases was used to examine the phylogenetic relationships of Ochotona cansus, O. huangensis, O. thibetana, O. curzoniae and O. erythrotis. A 1-kb length variation between O. erythrotis of subgenus Pika and other four species of subgenus Ochotona was observed, which may be a useful genetic marker for identifying the two subgenera. The phylogenetic tree constructed using PAUP based on 61 phylogenetically informative sites suggests that O. erythrotis diverged first, followed by O. cansus, while O. curzoniae and O. huangensis are sister taxa related to O. thibetana, The results indicate that both O. cansus and O. huangensis should be treated as independent species. If the base substitution rate of pikas mtDNA was 2% per million years, then the divergence time of the two subgenera, Pika and Ochotana, is about 8.8 Ma ago of late Miocence, middle Bao-dian of Chinese mammalian age, and the divergence of the four species in subgenus Ochotona would have occurred about 2.5 - 4.2 Ma ago, Yushean of Chinese mammalian age. This calculation appears to be substantiated by the fossil record.

The evolutionary origin of Indian Ocean tortoises (Dipsochelys).

Today, the only surviving wild population of giant tortoises in the Indian Ocean occurs on the island of Aldabra. However, giant tortoises once inhabited islands throughout the western Indian Ocean. Madagascar, Africa, and India have all been suggested as possible sources of colonization for these islands. To address the origin of Indian Ocean tortoises (Dipsochelys, formerly Geochelone gigantea), we sequenced the 12S, 16S, and cyt b genes of the mitochondrial DNA. Our phylogenetic analysis shows Dipsochelys to be embedded within the Malagasy lineage, providing evidence that Indian Ocean giant tortoises are derived from a common Malagasy ancestor. This result points to Madagascar as the source of colonization for western Indian Ocean islands by giant tortoises. Tortoises are known to survive long oceanic voyages by floating with ocean currents, and thus, currents flowing northward towards the Aldabra archipelago from the east coast of Madagascar would have provided means for the colonization of western Indian Ocean islands. Additionally, we found an accelerated rate of sequence evolution in the two Malagasy Pyxis species examined. This finding supports previous theories that shorter generation time and smaller body size are related to an increase in mitochondrial DNA substitution rate in vertebrates.

MODELING HETEROTACHY IN PHYLOGENETICS

Il a été démontré que l’hétérotachie, variation du taux de substitutions au cours du temps et entre les sites, est un phénomène fréquent au sein de données réelles. Échouer à modéliser l’hétérotachie peut potentiellement causer des artéfacts phylogénétiques. Actuellement, plusieurs modèles traitent l’hétérotachie : le modèle à mélange des longueurs de branche (MLB) ainsi que diverses formes du modèle covarion. Dans ce projet, notre but est de trouver un modèle qui prenne efficacement en compte les signaux hétérotaches présents dans les données, et ainsi améliorer l’inférence phylogénétique. Pour parvenir à nos fins, deux études ont été réalisées. Dans la première, nous comparons le modèle MLB avec le modèle covarion et le modèle homogène grâce aux test AIC et BIC, ainsi que par validation croisée. A partir de nos résultats, nous pouvons conclure que le modèle MLB n’est pas nécessaire pour les sites dont les longueurs de branche diffèrent sur l’ensemble de l’arbre, car, dans les données réelles, le signaux hétérotaches qui interfèrent avec l’inférence phylogénétique sont généralement concentrés dans une zone limitée de l’arbre. Dans la seconde étude, nous relaxons l’hypothèse que le modèle covarion est homogène entre les sites, et développons un modèle à mélanges basé sur un processus de Dirichlet. Afin d’évaluer différents modèles hétérogènes, nous définissons plusieurs tests de non-conformité par échantillonnage postérieur prédictif pour étudier divers aspects de l’évolution moléculaire à partir de cartographies stochastiques. Ces tests montrent que le modèle à mélanges covarion utilisé avec une loi gamma est capable de refléter adéquatement les variations de substitutions tant à l’intérieur d’un site qu’entre les sites. Notre recherche permet de décrire de façon détaillée l’hétérotachie dans des données réelles et donne des pistes à suivre pour de futurs modèles hétérotaches. Les tests de non conformité par échantillonnage postérieur prédictif fournissent des outils de diagnostic pour évaluer les modèles en détails. De plus, nos deux études révèlent la non spécificité des modèles hétérogènes et, en conséquence, la présence d’interactions entre différents modèles hétérogènes. Nos études suggèrent fortement que les données contiennent différents caractères hétérogènes qui devraient être pris en compte simultanément dans les analyses phylogénétiques.

Estimation des corrélations phylogénétiques entre paramètres d'évolution moléculaire et Traits d'histoire de vie

Depuis quelques années, l'évolution moléculaire cherche à caractériser les variations et l'intensité de la sélection grâce au rapport entre taux de substitution synonyme et taux de substitution non-synonyme (dN/dS). Cette mesure, dN/dS, a permis d'étudier l'histoire de la variation de l'intensité de la sélection au cours du temps ou de détecter des épisodes de la sélection positive. Les liens entre sélection et variation de taille efficace interfèrent cependant dans ces mesures. Les méthodes comparatives, quant a elle, permettent de mesurer les corrélations entre caractères quantitatifs le long d'une phylogénie. Elles sont également utilisées pour tester des hypothèses sur l'évolution corrélée des traits d'histoire de vie, mais pour être employées pour étudier les corrélations entre traits d'histoire de vie, masse, taux de substitution ou dN/dS. Nous proposons ici une approche combinant une méthode comparative basée sur le principe des contrastes indépendants et un modèle d'évolution moléculaire, dans un cadre probabiliste Bayésien. Intégrant, le long d'une phylogénie, sur les reconstructions ancestrales des traits et et de dN/dS nous estimons les covariances entre traits ainsi qu'entre traits et paramètres du modèle d'évolution moléculaire. Un modèle hiérarchique, a été implémenté dans le cadre du logiciel coevol, publié au cours de cette maitrise. Ce modèle permet l'analyse simultané de plusieurs gènes sans perdre la puissance donnée par l'ensemble de séquences. Un travail deparallélisation des calculs donne la liberté d'augmenter la taille du modèle jusqu'à l'échelle du génome. Nous étudions ici les placentaires, pour lesquels beaucoup de génomes complets et de mesures phénotypiques sont disponibles. À la lumière des théories sur les traits d'histoire de vie, notre méthode devrait permettre de caractériser l'implication de groupes de gènes dans les processus biologique liés aux phénotypes étudiés.

Cytochrome-b variation in Apis mellifera samples and its association with COI-COII patterns

Five Mbo I (Mbo-A, Mbo-M, Mbo-C(1), Mbo-C(2) and Mbo-C(3)) and Hinf I (Hinf-1 to Hinf-5) patterns were observed in Apis mellifera samples after restriction of a 485 bp fragment of the mitochondrial cytochrome-b (cyt-b) gene. Associating the cyt-b Restriction fragment length polymorphism (RFLP) pattern of each sample to its respective previously established COI-COII (Dra I sites) pattern, five restriction patterns (Mbo-C(1), Mbo-C(2), Mbo-C(3), Hinf-1 and Hinf-4) were observed in samples of maternal origin associated to the evolutionary branch C. No deletions or insertions were observed and the nucleotide substitution rate was estimated at 5.4%. Higher nucleotide diversity was observed among the branch C-haplotypes when compared with A and M lineages. Further studies are needed to confirm if the cyt-b + COI-COII haplotypes help to assign certain phylogeographic patterns to the branch C and to clarify phylogenetic relationships among A. mellifera subspecies.