Molecular phylogenetic systematics of twelve species of acipenseriformes based on mtDNA ND4L-ND4 gene sequence analysis

Acipenseriformes is an endangered primitive fish group, which occupies a special place in the history of ideas concerning fish evolution, even in vertebrate evolution. However, the classification and evolution of the fishes have been debated. The mitochondrial DNA (mtDNA) ND4L and partial ND4 genes were first sequenced in twelve species of the order Acipenseriformes, including endemic Chinese species. The following points were drawn from DNA sequences analysis: (i) the two species of Huso can be ascribed to Acipenser; (ii) A. dabryanus is the mostly closely related to A. sinensis, and most likely the landlocked form of A. sinensis; (iii) genus Acipenser in trans-Pacific region might have a common origin; (iv) mtDNA ND4L and ND4 genes are the ideal genetic markers for phylogenetic analysis of the order Acipenseriformes.

Generation and analysis of 5318 expressed sequence tags from the filamentous sporophyte of Porphyra haitanensis (Rhodophyta)

Porphyra haitanensis T. J. Chang et B. F. Zheng (Bangiales, Rhodophyta) is cultivated in China and widely consumed in Asia. To gain more insight into its physiological and biochemical properties, we generated 5318 expressed sequence tags (ESTs) from the sporophyte of P. haitanensis, and upon assembling into a nonredundant set, 2535 sequences were obtained, among which only 32.2% (816) shared certain similarity with published sequences (Nr and KOG). Functional classification of such ESTs revealed that most of the transcripts were related to its conservative biological metabolism, and P. haitanensis most likely possesses cyanide-resistant respiration and a C4-like carbon-fixation pathway, both of which have never been reported in a rhodophyte before. Twenty-eight percent of the nonredundant gene clusters exhibited significant similarity to those from P. yezoensis Ueda sporophytes, and 16 genes up-regulated in P. yezoensis sporophytes were also expressed abundantly in P. haitanensis. Codon usage analysis indicated that exposure to high GC pressure might occur during evolution of P. haitanensis. These findings represent the most extensive collection of ESTs from P. haitanensis to date, and all the ESTs in this study have been submitted to GenBank (accession nos. DN604790-DN608469, EG016226-EG018540).

Sequence analysis of the ITS region and 5.8S rDNA of Porphyra haitanensis

The sequences of the ITS (internal transcribed spacer) and 5.8S rDNA of three cultivated strains of Porphyra haitanensis thalli (NB, PT and ST) were amplified, sequenced and analyzed. In addition, the phylogenic relationships of the sequences identified in this study with those of other Porphyra retrieved from GenBank were evaluated. The results are as follows: the sequences of the ITS and 5.8S rDNA were essentially identical among the three strains. The sequences of ITS l were 331 by to 334 bp, while those of the 5.8S rDNA were 158 by and the sequences of ITS2 ranged from 673 by to 681 bp. The sequences of the ITS had a high level of homology (up to 99.5%) with that of P. haitanensis (DQ662228) retrieved from GenBank, but were only approximately 50% homologous with those of other species of Porphyra. The results obtained when a phylogenetic tree was constructed coincided with the results of the homology analysis. These results suggest that the three cultivated strains of P. haitanensis evolved conservatively and that the ITS showed evolutionary consistency. However, the sequences of the ITS and 5.8S rDNA of different Porphyra species showed great variations. Therefore, the relationship of Porphyra interspecies phyletic evolution could be judged, which provides the proof for Porphyra identification study. However, proper classifications of the subspecies and the populations of Porphyra should be determined through the use of other molecular techniques to determine the genetic variability and rational phylogenetic relationships.

The complete mitochondrial genome sequence of the cutlassfish Trichiurus japonicus (Perciformes: Trichiuridae) : Genome characterization and phylogenetic considerations

Mitochondrial genome sequence and structure analysis has become a powerful tool for studying molecular evolution and phylogenetic relationships. To understand the systematic status of Trichiurus japonicus in suborder Scombroidei, we determined the complete mitochondrial genome (mitogenome) sequence using the long-polymerase chain reaction (long-PCR) and shotgun sequencing method. The entire mitogenome is 16,796 by in length and has three unusual features, including (1) the absence of tRNA(Pro) gene, (2) the possibly nonfunctional light-strand replication origin (O-L) showing a shorter loop in secondary structure and no conserved motif (5'-GCCGG-3'), (3) two sets of the tandem repeats at the 5' and 3' ends of the control region. The three features seem common for Trichiurus mitogenomes, as we have confirmed them in other three T. japonicus individuals and in T nanhaiensis. Phylogenetic analysis does not support the monophyly of Trichiuridae, which is against the morphological result. T. japonicus is most closely related to those species of family Scombridae; they in turn have a sister relationship with Perciformes members including suborders Acanthuroidei, Caproidei, Notothenioidei, Zoarcoidei, Trachinoidei, and some species of Labroidei, based on the current dataset of complete mitogenome. T japonicus together with T. brevis, T lepturus and Aphanopus carbo form a clade distinct from Lepidopus caudatus in terms of the complete Cyt b sequences. T. japonicus mitogenome, as the first discovered complete mitogenome of Trichiuridae, should provide important information on both genomics and phylogenetics of Trichiuridae. (C) 2009 Elsevier B.V. All rights reserved.

Lineage-specific domain fusion in the evolution of purine nucleotide cyclases in cyanobacteria

Cyclic nucleotides (both cAMP and cGMP) play extremely important roles in cyanobacteria, such as regulating heterocyst formation, respiration, or gliding. Catalyzing the formation of cAMP and cGMP from ATP and GTP is a group of functionally important enzymes named adenylate cyclases and guanylate cyclases, respectively. To understand their evolutionary patterns, in this study, we presented a systematic analysis of all the cyclases in cyanobacterial genomes. We found that different cyanobacteria had various numbers of cyclases in view of their remarkable diversities in genome size and physiology. Most of these cyclases exhibited distinct domain architectures, which implies the versatile functions of cyanobacterial cyclases. Mapping the whole set of cyclase domain architectures from diverse prokaryotic organisms to their phylogenetic tree and detailed phylogenetic analysis of cyclase catalytic domains revealed that lineage-specific domain recruitment appeared to be the most prevailing pattern contributing to the great variability of cyanobacterial cyclase domain architectures. However, other scenarios, such as gene duplication, also occurred during the evolution of cyanobacterial cyclases. Sequence divergence seemed to contribute to the origin of putative guanylate cyclases which were found only in cyanobacteria. In conclusion, the comprehensive survey of cyclases in cyanobacteria provides novel insight into their potential evolutionary mechanisms and further functional implications.

Mesozoic-Cenozoic tectonic evolution of the Zhuanghai area, Bohai-Bay Basin, east China: the application of balanced cross-sections

The technique of balancing cross-sections, an important method for studying the tectonic history of sedimentary basins, has many applications. It enables one to compile charts for petroleum exploration and development, and growth sections of ancient structures can be restored so that the structural growth history can be studied. In order to study tectonic evolution in the Zhuanghai area of the Bohai-Bay basin, we selected two seismic profiles and compiled two structural growth sections. Based on the two balanced cross-sections, the evolution can be divided into four phases: the Triassic-Middle Jurassic phase, Late Jurassic - Cretaceous phase, Palaeogene extension phase, and Late Palaeogene-to-present phase. The whole area was uplifted during the Triassic-Middle Jurassic phase because of intense extrusion stress related to the Indo-China movement. During the Late Jurassic and Early Cretaceous, intense extension occurred in east China, and the whole area rifted, leading to the deposition of a thick sedimentary sequence. In the Late Cretaceous, the area suffered uplift and compression associated with the sinistral strike slip of the Tanlu fault. In the Palaeogene, a rifting basin developed in the area. Finally, it became stable and was placed in its present position by dextral strike-slip motion. In addition, some problems associated with compiling balanced cross-sections are discussed.

A close phylogenetic relationship between Sipuncula and Annelida evidenced from the complete mitochondrial genome sequence of Phascolosoma esculenta

Background: There are many advantages to the application of complete mitochondrial (mt) genomes in the accurate reconstruction of phylogenetic relationships in Metazoa. Although over one thousand metazoan genomes have been sequenced, the taxonomic sampling is highly biased, left with many phyla without a single representative of complete mitochondrial genome. Sipuncula (peanut worms or star worms) is a small taxon of worm-like marine organisms with an uncertain phylogenetic position. In this report, we present the mitochondrial genome sequence of Phascolosoma esculenta, the first complete mitochondrial genome of the phylum. Results: The mitochondrial genome of P. esculenta is 15,494 bp in length. The coding strand consists of 32.1% A, 21.5% C, 13.0% G, and 33.4% T bases (AT = 65.5%; AT skew = -0.019; GC skew = -0.248). It contains thirteen protein-coding genes (PCGs) with 3,709 codons in total, twenty-two transfer RNA genes, two ribosomal RNA genes and a non-coding AT-rich region (AT = 74.2%). All of the 37 identified genes are transcribed from the same DNA strand. Compared with the typical set of metazoan mt genomes, sipunculid lacks trnR but has an additional trnM. Maximum Likelihood and Bayesian analyses of the protein sequences show that Myzostomida, Sipuncula and Annelida (including echiurans and pogonophorans) form a monophyletic group, which supports a closer relationship between Sipuncula and Annelida than with Mollusca, Brachiopoda, and some other lophotrochozoan groups. Conclusion: This is the first report of a complete mitochondrial genome as a representative within the phylum Sipuncula. It shares many more similar features with the four known annelid and one echiuran mtDNAs. Firstly, sipunculans and annelids share quite similar gene order in the mitochondrial genome, with all 37 genes located on the same strand; secondly, phylogenetic analyses based on the concatenated protein sequences also strongly support the sipunculan + annelid clade (including echiurans and pogonophorans). Hence annelid "key-characters" including segmentation may be more labile than previously assumed.

Chromosomal rearrangement in Pectinidae revealed by rRNA loci and implications for bivalve evolution

Karyotype and chromosomal localization of major (18-5.8-28S) and minor (5S) ribosomal RNA genes were studied in two species of Pectinidae, zhikong (Chlamys farreri) and bay (Argopecten irradians irradians) scallops. using fluorescence in situ hybridization (FISH). C. farreri had a haploid number of 19 with a karyotype of 3m + 4sm + 7sm-st + 4st + 1st-t, and A. i. irradians had a haploid number of 16 with a karyotype of 5st + 11t. In C. farreri, the major and minor rRNA genes had one locus each and were mapped to the same chromosome-Chromosome 5. In A. i. irradians, the major rRNA genes had two loci, located on Chromosomes 4 and 8, and the 5S rRNA gene was found at a third chromosome-Chromosome 10. Results of this and other studies indicate that karyotype of A. i. irradians (n = 16, 21 arms) is secondary and derived from an ancestral karyotype similar to that of C. farreri (n = 19, 38 arms) through considerable chromosomal loss and rearrangements. The ability to tolerate significant chromosomal loss suggests that the modal karyotype of Pectinidae and possibly other bivalves with a haploid number of 19 is likely tetraploid; i.e., at least one genome duplication has occurred during the evolution of Bivalvia.

Molecular phylogeny, recent radiation and evolution of gross morphology of the rhubarb genus Rheum (Polygonaceae) inferred from chloroplast DNA trnL-F sequences

Background and Aims Rheum, a highly diversified genus with about 60 species, is mainly confined to the mountainous and desert regions of the Qinghai-Tibetan plateau and adjacent areas. This genus represents a good example of the extensive diversification of the temperate genera in the Qinghai-Tibetan plateau, in which the forces to drive diversification remain unknown. To date, the infrageneric classification of Rheum has been mainly based on morphological characters. However, it may have been subject to convergent evolution under habitat pressure, and the systematic position of some sections are unclear, especially Sect. Globulosa, which has globular inflorescences, and Sect. Nobilia, which has semi-translucent bracts. Recent palynological research has found substantial contradictions between exine patterns and the current classification of Rheum. Two specific objectives of this research were (1) to evaluate possible relationships of some ambiguous sections with a unique morphology, and (2) to examine possible occurrence of the radiative speciation with low genetic divergence across the total genus and the correlation between the extensive diversification time of Rheum and past geographical events, especially the recent large-scale uplifts of the Qinghai-Tibetan Plateau.Methods The chloroplast DNA trnL-F region of 29 individuals representing 26 species of Rheum, belonging to seven out of eight sections, was sequenced and compared. The phylogenetic relationships were further constructed based on the sequences obtained.Key Results Despite the highly diversified morphology, the genetic variation in this DNA fragment is relatively low. The molecular phylogeny is highly inconsistent with gross morphology, pollen exine patterns and traditional classifications, except for identifying all samples of Sect. Palmata, three species of Sect. Spiciformia and a few species of Sect. Rheum as corresponding monophyletic groups. The monotypic Sect. Globulosa showed a tentative position within the clade comprising five species of Sect. Rheum. All of the analyses revealed the paraphyly of R. nobile and R. alexandrae, the only two species of Sect. Nobilia circumscribed by the possession of large bracts. The crude calibration of lineages based on trnL-F sequence differentiation implied an extensive diversification of Rheum within approx. 7 million years.Conclusions Based on these results, it is suggested that the rich geological and ecological diversity caused by the recent large-scale uplifts of the Qinghai-Tibetan Plateau since the late Tertiary, coupled with the oscillating climate of the Quaternary stage, might have promoted rapid speciation in small and isolated populations, as well as allowing the fixation of unique or rare morphological characters in Rheum. Such a rapid radiation, combined with introgressive hybridization and reticulate evolution, may have caused the transfer of cpDNA haplotypes between morphologically dissimilar species, and might account for the inconsistency between morphological classification and molecular phylogeny reported here.

A new species of Saussurea (Asteraceae) from Tibet and its systematic position based on ITS sequence analysis

A new species of Saussurea, S. erecta S. W Liu, J. T Pan A J. Q. Liu sp. nov., is described from Tibet. It resembles S. kingii but may be distinguished by having distinct stems and glabrous achenes. Saussurea kingii was placed in sect. Pseudoeriocoryne of subgen. Eriocoryne; this section was circumscribed by acaulescence and an inflorescence with congested capitula surrounded by a rosette of leaves. The discovery of S. erecta with distinct stems, cauline leaves and corymbose capitula blurred the delimitation of sect. Pseudoeriocoryne and suggested that the section may be polyphyletic. Both the close relationship and the significant difference between S. erecta and S. kingii were confirmed by analyses of nrDNA ITS sequences. The resulting phylogenies based on ITS data further suggest that Saussurea sect. Pseudoeriocoryne, as traditionally defined, does not constitute a monophyletic group. The rapid radiation and speciation of Saussurea in the Qinghai-Tibetan Plateau, as inferred from ITS phylogeny, are discussed. (c) 2005 The Linnean Society of London.

A Phylogenetic Mixture Model for Detecting Pattern-Heterogeneity in Gene Sequence or Character-State Data

Mark Pagel, Andrew Meade (2004). A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data. Systematic Biology, 53(4), 571-581. RAE2008

Modeling the evolution of regulatory elements by simultaneous detection and alignment with phylogenetic pair HMMs.

The computational detection of regulatory elements in DNA is a difficult but important problem impacting our progress in understanding the complex nature of eukaryotic gene regulation. Attempts to utilize cross-species conservation for this task have been hampered both by evolutionary changes of functional sites and poor performance of general-purpose alignment programs when applied to non-coding sequence. We describe a new and flexible framework for modeling binding site evolution in multiple related genomes, based on phylogenetic pair hidden Markov models which explicitly model the gain and loss of binding sites along a phylogeny. We demonstrate the value of this framework for both the alignment of regulatory regions and the inference of precise binding-site locations within those regions. As the underlying formalism is a stochastic, generative model, it can also be used to simulate the evolution of regulatory elements. Our implementation is scalable in terms of numbers of species and sequence lengths and can produce alignments and binding-site predictions with accuracy rivaling or exceeding current systems that specialize in only alignment or only binding-site prediction. We demonstrate the validity and power of various model components on extensive simulations of realistic sequence data and apply a specific model to study Drosophila enhancers in as many as ten related genomes and in the presence of gain and loss of binding sites. Different models and modeling assumptions can be easily specified, thus providing an invaluable tool for the exploration of biological hypotheses that can drive improvements in our understanding of the mechanisms and evolution of gene regulation.

R5 clade C SHIV strains with tier 1 or 2 neutralization sensitivity: tools to dissect env evolution and to develop AIDS vaccines in primate models.

BACKGROUND: HIV-1 clade C (HIV-C) predominates worldwide, and anti-HIV-C vaccines are urgently needed. Neutralizing antibody (nAb) responses are considered important but have proved difficult to elicit. Although some current immunogens elicit antibodies that neutralize highly neutralization-sensitive (tier 1) HIV strains, most circulating HIVs exhibiting a less sensitive (tier 2) phenotype are not neutralized. Thus, both tier 1 and 2 viruses are needed for vaccine discovery in nonhuman primate models. METHODOLOGY/PRINCIPAL FINDINGS: We constructed a tier 1 simian-human immunodeficiency virus, SHIV-1157ipEL, by inserting an "early," recently transmitted HIV-C env into the SHIV-1157ipd3N4 backbone [1] encoding a "late" form of the same env, which had evolved in a SHIV-infected rhesus monkey (RM) with AIDS. SHIV-1157ipEL was rapidly passaged to yield SHIV-1157ipEL-p, which remained exclusively R5-tropic and had a tier 1 phenotype, in contrast to "late" SHIV-1157ipd3N4 (tier 2). After 5 weekly low-dose intrarectal exposures, SHIV-1157ipEL-p systemically infected 16 out of 17 RM with high peak viral RNA loads and depleted gut CD4+ T cells. SHIV-1157ipEL-p and SHIV-1157ipd3N4 env genes diverge mostly in V1/V2. Molecular modeling revealed a possible mechanism for the increased neutralization resistance of SHIV-1157ipd3N4 Env: V2 loops hindering access to the CD4 binding site, shown experimentally with nAb b12. Similar mutations have been linked to decreased neutralization sensitivity in HIV-C strains isolated from humans over time, indicating parallel HIV-C Env evolution in humans and RM. CONCLUSIONS/SIGNIFICANCE: SHIV-1157ipEL-p, the first tier 1 R5 clade C SHIV, and SHIV-1157ipd3N4, its tier 2 counterpart, represent biologically relevant tools for anti-HIV-C vaccine development in primates.

Functional evolution of mammalian odorant receptors.

The mammalian odorant receptor (OR) repertoire is an attractive model to study evolution, because ORs have been subjected to rapid evolution between species, presumably caused by changes of the olfactory system to adapt to the environment. However, functional assessment of ORs in related species remains largely untested. Here we investigated the functional properties of primate and rodent ORs to determine how well evolutionary distance predicts functional characteristics. Using human and mouse ORs with previously identified ligands, we cloned 18 OR orthologs from chimpanzee and rhesus macaque and 17 mouse-rat orthologous pairs that are broadly representative of the OR repertoire. We functionally characterized the in vitro responses of ORs to a wide panel of odors and found similar ligand selectivity but dramatic differences in response magnitude. 87% of human-primate orthologs and 94% of mouse-rat orthologs showed differences in receptor potency (EC50) and/or efficacy (dynamic range) to an individual ligand. Notably dN/dS ratio, an indication of selective pressure during evolution, does not predict functional similarities between orthologs. Additionally, we found that orthologs responded to a common ligand 82% of the time, while human OR paralogs of the same subfamily responded to the common ligand only 33% of the time. Our results suggest that, while OR orthologs tend to show conserved ligand selectivity, their potency and/or efficacy dynamically change during evolution, even in closely related species. These functional changes in orthologs provide a platform for examining how the evolution of ORs can meet species-specific demands.

Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits.

Alewife, Alosa pseudoharengus, populations occur in two discrete life-history variants, an anadromous form and a landlocked (freshwater resident) form. Landlocked populations display a consistent pattern of life-history divergence from anadromous populations, including earlier age at maturity, smaller adult body size, and reduced fecundity. In Connecticut (USA), dams constructed on coastal streams separate anadromous spawning runs from lake-resident landlocked populations. Here, we used sequence data from the mtDNA control region and allele frequency data from five microsatellite loci to ask whether coastal Connecticut landlocked alewife populations are independently evolved from anadromous populations or whether they share a common freshwater ancestor. We then used microsatellite data to estimate the timing of the divergence between anadromous and landlocked populations. Finally, we examined anadromous and landlocked populations for divergence in foraging morphology and used divergence time estimates to calculate the rate of evolution for foraging traits. Our results indicate that landlocked populations have evolved multiple times independently. Tests of population divergence and estimates of gene flow show that landlocked populations are genetically isolated, whereas anadromous populations exchange genes. These results support a 'phylogenetic raceme' model of landlocked alewife divergence, with anadromous populations forming an ancestral core from which landlocked populations independently diverged. Divergence time estimates suggest that landlocked populations diverged from a common anadromous ancestor no longer than 5000 years ago and perhaps as recently as 300 years ago, depending on the microsatellite mutation rate assumed. Examination of foraging traits reveals landlocked populations to have significantly narrower gapes and smaller gill raker spacings than anadromous populations, suggesting that they are adapted to foraging on smaller prey items. Estimates of evolutionary rates (in haldanes) indicate rapid evolution of foraging traits, possibly in response to changes in available resources.