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.

Adaptive evolution of primate TRIM5 alpha, a gene restricting HIV-1 infection

Recent studies showed that nonhuman primate TRIM5 alpha can efficiently block HIV-1 infection in human cell lines. It can also restrict other retroviruses, therefore, suggested as a general defender against retrovirus infection. Here, we present an evolutionary analysis of TRIM5 alpha in primates. Our results demonstrated that TRIM5a has been evolving rapidly in primates, which is likely caused by Darwinian positive selection. The SPRY domain of TRM5 alpha, which may be responsible for recognition of incoming viral capsids showed higher nonsynonymous/synonymous substitution ratios than the non-SPRY domain, indicating that the adaptive evolution of TRIM5a ill primates might be an innate strategy developed in defending retrovirus infection during primate evolution. In addition, the comparative protein sequence analysis suggested that the amino acid substitution pattern at a single site (344R/Q/P) located in the SPRY domain may explain the differences in Susceptibilities of HIV-1 infection in diverse primate species. (c) 2005 Elsevier B.V. All rights reserved.

Maximin 9, a novel free thiol containing antimicrobial peptide with antimycoplasma activity from frog Bombina maxima

Amphibian skin is a rich resource of antimicrobial peptides, like maximins and maximin Hs from frog Bombina maxima. Novel cDNA clones encoding a precursor protein, which comprises a novel maximin peptide (maximin 9) and reported maximin H3, were isolated from two constructed skin cDNA libraries of B. maxima. The predicted primary structure of maximin 9 is GIGRKFLGGVKTTFRCGVKDFASKHLY-NH2. A surprising substitution is at position 16, with a free cysteine in maximin 9 rather than usual conserved glycine in other reported maximins. Maximin 9, the homodimer form and its Cys(16) to Gly(16) mutant were synthesized and their antimicrobial activities were evaluated. Unlike previously reported maximin 3, the tested bacterial and fungal strains were resistant to maximin 9, its homodimer and the Cys(16) to Gly(16) mutant (with MICs > 100 mu M). On the other hand, interestingly, while eight clinical Mollicutes strains were generally resistant to maximin 9 homodimer and its Cys(16) to Gly(16) mutant, most of them are sensitive to maximin 9 at a peptide concentration of 30 mu M, especially in the presence of dithiothreitol. These results indicate that the presence of a reactive Cys residue in maximin 9 is important for its antimycoplasma activity. The diversity of antimicrobial peptide cDNA structures encountered in B. maxima skin cDNA libraries and the antimicrobial specificity differences of the peptides may reflect well the species' adaptation to the unique microbial environments. (c) 2005 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserved.

Variety of antimicrobial peptides in the Bombina maxima toad and evidence of their rapid diversification

Antimicrobial peptides secreted by the skin of many amphibians play an important role in innate immunity. From two skin cDNA libraries of two individuals of the Chinese red belly toad (Bombina maxima), we identified 56 different antimicrobial peptide cDNA sequences, each of which encodes a precursor peptide that can give rise to two kinds of antimicrobial peptides, maximin and maximin H. Among these cDNA, we found that the mean number of nucleotide substitution per non-synonymous site in both the maximin and maximin H domains significantly exceed the mean number of nucleotide substitution per synonymous site, whereas the same pattern was not observed in other structural regions, such as the signal and propiece peptide regions, suggesting that these antimicrobial peptide genes have been experiencing rapid diversification driven by Darwinian selection. We cloned and sequenced seven genes amplified from skin or liver genomic DNA. These genes have three exons and share the same gene structure, in which both maximin and maximin H are encoded by the third exon. This suggests that alternative splicing and somatic recombination are less likely to play a role in creating the diversity of maximins and maximin Hs. The gene trees based on different domain regions revealed that domain shuffling or gene conversion among these genes might have happened frequently.

Comparative study of three short-chain neurotoxins from the venom of Naja kaouthia (Yunnan, China)

Three short-chain neurotoxins named NT-I, NT-II, and NT-III were purified from the venom of Naja kaouthia, a snake distributed throughout the south of Yunnan province, China, by a series of chromatographic steps, including an FPLC Resource S column. Their molecular weights, determined by matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) MS, were 6952.19 Da, 6854.92 Da, and 6828.80 Da, respectively. NT-I consisted of 62 amino acid residues, and the other two consisted of 61 amino acid residues, including 8 cysteines. After hydrolysis by endoproteinase Glu-C, their primary sequences were determined. A test of their activities demonstrated that they effectively inhibited muscle contractions induced by electric stimulation. Furthermore, the extent of inhibition caused by NT-II and NT-III was less than that of NT-I. The IC(50)s were 0.04 mug/ml, 0.20 mug/ml, and 0.23 mug/ml for NT-I, NT-II, and NT-III, respectively. Compared with NT-II and NT-III, the higher activity of NT-I may be a result of the amino acid residue substitution Ile36 to Arg36.

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.

Sox and Zfx genes in giant panda

By using PCR cloning techniques, the DNA sequences of the HMG box regions of six Sox genes (pSox) and the zinc finger domains of two Zfx genes (pZfx) in the giant panda were identified. The giant panda Sox genes fell into two subfamilies, SOX-S1 and SOX-S2. The pSox and pZfx genes of the giant panda were highly homologous to the corresponding genes in mammals and revealed close substitution rates to those in the primates.

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.

Adaptive diversification of bitter taste receptor genes in mammalian evolution

The diversity and evolution of bitter taste perception in mammals is not well understood. Recent discoveries of bitter taste receptor (T2R) genes provide an opportunity for a genetic approach to this question. We here report the identification of 10 and 30 putative T2R genes from the draft human and mouse genome sequences, respectively, in addition to the 23 and 6 previously known T2R genes from the two species. A phylogenetic analysis of the T2R genes suggests that they can be classified into three main groups, which are designated A, B, and C. Interestingly, while the one-to-one gene orthology between the human and mouse is common to group B and C genes, group A genes show a pattern of species- or lineage-specific duplication. It is possible that group B and C genes are necessary for detecting bitter tastants common to both humans and mice, whereas group A genes are used for species-specific bitter tastants. The analysis also reveals that phylogenetically closely related T2R genes are close in their chromosomal locations, demonstrating tandem gene duplication as the primary source of new T2Rs. For closely related paralogous genes, a rate of nonsynonymous nucleotide substitution significantly higher than the rate of synonymous substitution was observed in the extracellular regions of T2Rs, which are presumably involved in tastant-binding. This suggests the role of positive selection in the diversification of newly duplicated T2R genes. Because many natural poisonous substances are bitter, we conjecture that the mammalian T2R genes are under diversifying selection for the ability to recognize a diverse array of poisons that the organisms may encounter in exploring new habitats and diets.

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.

Evolutionary implications of avian infectious Bronchitis Virus (AIBV) analysis

For developing efficient vaccines, it is essential to identify which amino acid changes are most important to the survival of the virus. We investigate the amino acid substitution features in the Avian Infectious Bronchitis Virus (AIBV) antigenic domain o

Phylogeny and biogeography of the Petaurista philippensis complex (Rodentia : Sciuridae), inter- and intraspecific relationships inferred from molecular and morphometric analysis

With modified DNA extraction and Purification protocols, the complete cytochrome b gene sequences (1140 bp) were determined from degraded museum specimens. Molecular analysis and morphological examination of cranial characteristics of the giant flying squirrels of Petaurista philippensis complex (P. grandis, P. hainana, and P. yunanensis) and other Petaurista species yielded new insights into long-standing controversies in the Petaurista systematics. Patterns of genetic variations and morphological differences observed in this study indicate that P. hainana, P. albiventer, and P. yunanensis can be recognized as distinct species, and P. grandis and P. petaurista are conspecific populations. Phylogenetic relationships reconstructed by using parsimony, likelihood, and Bayesian methods reveal that, with P. leucogenys as the basal branch, all Petaurista groups formed two distinct clades. Petaurista philippensis, P. hainana, P. yunanensis, and P. albiventer are clustered in the same clade, while P. grandis shows a close relationship to P. petaurista. Deduced divergence times based on Bayesian analysis and the transversional substitution at the third codon suggest that the retreating of glaciers and upheavals or movements of tectonic plates in the Pliocene-Pleistocene were the major factors responsible for the present geographical distributions of Petaurista groups. (c) 2005 Elsevier Inc. All rights reserved.

Evolution of cd59 gene in mammals

The CD59-coding sequences were obtained from 5 mammals by PCR and BLAST, and combined with the available sequences in GenBank, the nucleotide substitution rates of mammalian cd59 were calculated. Results of synonymous and nonsynonymous substitution rates