Protein and DNA sequence determinants of thermophilic adaptation

There have been considerable attempts in the past to relate phenotypic trait—habitat temperature of organisms—to their genotypes, most importantly compositions of their genomes and proteomes. However, despite accumulation of anecdotal evidence, an exact and conclusive relationship between the former and the latter has been elusive. We present an exhaustive study of the relationship between amino acid composition of proteomes, nucleotide composition of DNA, and optimal growth temperature (OGT) of prokaryotes. Based on 204 complete proteomes of archaea and bacteria spanning the temperature range from −10 °C to 110 °C, we performed an exhaustive enumeration of all possible sets of amino acids and found a set of amino acids whose total fraction in a proteome is correlated, to a remarkable extent, with the OGT. The universal set is Ile, Val, Tyr, Trp, Arg, Glu, Leu (IVYWREL), and the correlation coefficient is as high as 0.93. We also found that the G + C content in 204 complete genomes does not exhibit a significant correlation with OGT (R = −0.10). On the other hand, the fraction of A + G in coding DNA is correlated with temperature, to a considerable extent, due to codon patterns of IVYWREL amino acids. Further, we found strong and independent correlation between OGT and the frequency with which pairs of A and G nucleotides appear as nearest neighbors in genome sequences. This adaptation is achieved via codon bias. These findings present a direct link between principles of proteins structure and stability and evolutionary mechanisms of thermophylic adaptation. On the nucleotide level, the analysis provides an example of how nature utilizes codon bias for evolutionary adaptation to extreme conditions. Together these results provide a complete picture of how compositions of proteomes and genomes in prokaryotes adjust to the extreme conditions of the environment.     

Nucleotide Composition Bias Affects Amino Acid Content in Proteins Coded by Animal Mitochondria

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content. Received: 3 July 1996 / Accepted: 6 November 1996     

Chloroplast phylogenomics of liverworts: a reappraisal of the backbone phylogeny of liverworts with emphasis on Ptilidiales

As one of the four main lineages diverging from the early diversification of land plants, the phylogeny of liverworts holds the information about nearly 500 Myr of independent adaptation to changing environments. Thus, resolving the phylogenetic history of liverworts will provide unique insights into the successful diversification of early land plants in terrestrial ecosystems. However, the deep diverging events of this group remain incompletely resolved, such as the definite position of Ptilidiales. Here, we aimed to reconstruct the backbone relationships of liverworts using 84 protein-coding chloroplast genes, a dataset comprising 35 representatives from all major lineages of liverworts, and three phylogenetic analyses, namely maximum parsimony, maximum likelihood and Bayesian inference. To test the impact of composition biases, the phylogenetic analyses were carried out using three alignments representing the same dataset either as: (i) nucleotides, (ii) amino acids, or (iii) recoded nucleotides applying ambiguity base code. Chloroplast genome data consistently supported the monophyletic origin of three major lineages in liverworts, as well as the majority of backbone relationships. Ptilidiales were found to be sister to Jungermanniales. The rapid accumulation of G/C tracks as a consequence of increased GC content is an important cause for the long branches inferred in this group. Our study not only provides empirical evidence to support the significance of plastid genome sequencing to reconstruct the phylogeny of this important plant lineage, but also suggests that the GC content has played a critical role in the evolutionary dynamics of plastid genomes in land plants.     

Molecular phylogeny and divergence times of deuterostome animals

The phylogenetic relationships among deuterostome animals have been debated for many years, and a diversity of hypotheses have been proposed based on both morphological and molecular data. Here we have assembled sequences of 217 nuclear-encoded proteins to address specific questions concerning their relationships and times of origin. We recovered significant support for urochordates as the closest relative of vertebrates with an analysis of 59 proteins (17,400 amino acids) and suggest that the basal position of urochordates found in previous molecular studies may have been the result of long-branch attraction biases. Our results also support Ambulacraria, the pairing of hemichordates with echinoderms (nine proteins; 2,382 amino acids), and Cyclostomata, the pairing of lampreys with hagfish (25 proteins; 6,895 amino acids). In addition, 325 shared proteins (102,110 amino acids) were obtained from the complete genomes of six vertebrates and a urochordate for phylogenetic analysis and divergence time estimation. An evolutionary timescale was estimated using a local (Bayesian) molecular clock method. We found that most major lineages of deuterostomes arose prior to the Cambrian Explosion of fossils (approximately 520 MYA) and that several lineages had originated before periods of global glaciation in the Precambrian.     

Adaptation to environmental temperature is a major determinant of molecular evolutionary rates in archaea

Methods to infer the ancestral conditions of life are commonly based on geological and paleontological analyses. Recently, several studies used genome sequences to gain information about past ecological conditions taking advantage of the property that the G+C and amino acid contents of bacterial and archaeal ribosomal DNA genes and proteins, respectively, are strongly influenced by the environmental temperature. The adaptation to optimal growth temperature (OGT) since the Last Universal Common Ancestor (LUCA) over the universal tree of life was examined, and it was concluded that LUCA was likely to have been a mesophilic organism and that a parallel adaptation to high temperature occurred independently along the two lineages leading to the ancestors of Bacteria on one side and of Archaea and Eukarya on the other side. Here, we focus on Archaea to gain a precise view of the adaptation to OGT over time in this domain. It has been often proposed on the basis of indirect evidence that the last archaeal common ancestor was a hyperthermophilic organism. Moreover, many results showed the influence of environmental temperature on the evolutionary dynamics of archaeal genomes: Thermophilic organisms generally display lower evolutionary rates than mesophiles. However, to our knowledge, no study tried to explain the differences of evolutionary rates for the entire archaeal domain and to investigate the evolution of substitution rates over time. A comprehensive archaeal phylogeny and a non homogeneous model of the molecular evolutionary process allowed us to estimate ancestral base and amino acid compositions and OGTs at each internal node of the archaeal phylogenetic tree. The last archaeal common ancestor is predicted to have been hyperthermophilic and adaptations to cooler environments can be observed for extant mesophilic species. Furthermore, mesophilic species present both long branches and high variation of nucleotide and amino acid compositions since the last archaeal common ancestor. The increase of substitution rates observed in mesophilic lineages along all their branches can be interpreted as an ongoing adaptation to colder temperatures and to new metabolisms. We conclude that environmental temperature is a major factor that governs evolutionary rates in Archaea.     

Assessing Diversity of DNA Structure-Related Sequence Features in Prokaryotic Genomes

Prokaryotic genomes are diverse in terms of their nucleotide and oligonucleotide composition as well as presence of various sequence features that can affect physical properties of the DNA molecule. We present a survey of local sequence patterns which have a potential to promote non-canonical DNA conformations (i.e. different from standard B-DNA double helix) and interpret the results in terms of relationships with organisms'' habitats, phylogenetic classifications, and other characteristics. Our present work differs from earlier similar surveys not only by investigating a wider range of sequence patterns in a large number of genomes but also by using a more realistic null model to assess significant deviations. Our results show that simple sequence repeats and Z-DNA-promoting patterns are generally suppressed in prokaryotic genomes, whereas palindromes and inverted repeats are over-represented. Representation of patterns that promote Z-DNA and intrinsic DNA curvature increases with increasing optimal growth temperature (OGT), and decreases with increasing oxygen requirement. Additionally, representations of close direct repeats, palindromes and inverted repeats exhibit clear negative trends with increasing OGT. The observed relationships with environmental characteristics, particularly OGT, suggest possible evolutionary scenarios of structural adaptation of DNA to particular environmental niches.     

基于线粒体基因组序列的鞘翅目肉食亚目水生类群系统发育分析

【目的】明确肉食亚目(Adephaga)水生类群线粒体基因组的基本特征,并基于线粒体基因组序列分析肉食亚目水生类群的系统发育关系。【方法】基于Illumina HiSeq X Ten测序技术测定了圆鞘隐盾豉甲Dineutus mellyi和齿缘龙虱Eretes sticticus的线粒体全基因组序列,对其进行了基因注释,并对其tRNA基因二级结构进行了预测分析。加上已公布的鞘翅目(Coleoptera)肉食亚目水生类群17个种的线粒体基因组序列,对该类群共19个种线粒体的蛋白质编码基因(protein-coding genes, PCGs)开展了比较基因组学分析,包括AT含量、密码子偏好性、选择压力等。基于13个PCGs的氨基酸序列和核苷酸序列,利用最大似然法(ML)和贝叶斯法(BI)分别构建鞘翅目肉食亚目水生类群的系统发育关系,并通过FcLM分析进一步评估伪龙虱科(Noteridae)和瀑甲科(Meruidae)的系统发育位置。【结果】圆鞘隐盾豉甲和齿缘龙虱的线粒体基因组全长分别为16 123 bp(GenBank登录号: MN781126)和16 196 bp(GenBank登录号: MN781132),都包含13个PCGs、22个tRNA基因、2个rRNA基因和1个D-loop区(控制区)。19个肉食亚目水生类群线粒体基因组PCGs的碱基组成都呈现A+T偏好性,在密码子使用上也都偏向于使用富含A+T的密码子;在进化过程中13个PCGs的进化模式相同,都受到纯化选择。基于线粒体基因组13个PCGs的氨基酸序列的肉食亚目水生类群的系统发育关系为(豉甲科Gyrinidae+(沼梭甲科Haliplidae+((壁甲科Aspidytidae+(两栖甲科Amphizoidae+龙虱科Dytiscidae))+(水甲科Hygrobiidae+(瀑甲科Meruidae+伪龙虱科Noteridae)))))。【结论】研究结果表明,豉甲科是肉食亚目水生类群的基部类群,接下来是沼梭甲科和龙虱总科;伪龙虱科和瀑甲科互为姐妹群,并一起作为龙虱总科内部的一个分支;两栖甲科与龙虱科具有更近的亲缘关系。     

Compositional changes in RNA,DNA and proteins for bacterial adaptation to higher and lower temperatures

It is known that in thermophiles the G+C content of ribosomal RNA linearly correlates with growth temperature, while that of genomic DNA does not. Although the G+C contents (singlet) of the genomic DNAs of thermophiles and methophiles do not differ significantly, the dinucleotide (doublet) compositions of the two bacterial groups clearly do. The average amino acid compositions of proteins of the two groups are also distinct. Based on these facts, we here analyzed the DNA and protein compositions of various bacteria in terms of the optimal growth temperature (OGT). Regression analyses of the sequence data for thermophilic, mesophilic and psychrophilic bacteria revealed good linear relationships between OGT and the dinucleotide compositions of DNA, and between OGT and the amino acid compositions of proteins. Together with the above-mentioned linear relationship between ribosomal RNA and OGT, the DNA and protein compositions can be regarded as thermostability measures for RNA, DNA and proteins, covering a wide range of temperatures. Both the DNA and proteins of psychrophiles apparently exhibit characteristics diametrically opposite to those of thermophiles. The physicochemical parameters of dinucleotides suggested that supercoiling of DNA is relevant to its thermostability. Protein stability in thermophiles is realized primarily through global changes that increase charged residues (i.e., Glu, Arg, and Lys) on the molecular surface of all proteins. This kind of global change is attainable through a change in the amino acid composition coupled with alterations in the DNA base composition. The general strategies of thermophiles and psychrophiles for adaptation to higher and lower temperatures, respectively, that are suggested by the present study are discussed.     

Does the growth temperature of a prokaryote influence the purine content of its mRNAs?

The formation and breaking of hydrogen bonds between nucleic acid bases are dependent on temperature. The high G+C content of organisms was surmised to be an adaptation for high temperature survival because of the thermal stability of G:C pairs. However, a survey of genomic GC% and optimum growth temperature (OGT) of several prokaryotes revoked any direct relation between them. Significantly high purine (R=A or G) content in mRNAs is also seen as a selective response for survival among thermophiles. Nevertheless, the biological relevance of thermophiles loading their unstable mRNAs with excess purines (purine-loading or R-loading) is not persuasive. Here, we analysed the mRNA sequences from the genomes of 168 prokaryotes (as obtained from NCBI Genome database) with their OGTs ranging from -5 °C to 100 °C to verify the relation between R-loading and OGT. Our analysis fails to demonstrate any correlation between R-loading of the mRNA pool and OGT of a prokaryote. The percentage of purine-loaded mRNAs in prokaryotes is found to be in a rough negative correlation with the genomic GC% (r(2)=0.655, slope=-1.478, P<000.1). We conclude that genomic GC% and bias against certain combinations of nucleotides drive the mRNA-synonymous (sense) strands of DNA towards variations in R-loading.     

HGT-DB: a database of putative horizontally transferred genes in prokaryotic complete genomes

The Horizontal Gene Transfer DataBase (HGT-DB) is a genomic database that includes statistical parameters such as G+C content, codon and amino-acid usage, as well as information about which genes deviate in these parameters for prokaryotic complete genomes. Under the hypothesis that genes from distantly related species have different nucleotide compositions, these deviated genes may have been acquired by horizontal gene transfer. The current version of the database contains 88 bacterial and archaeal complete genomes, including multiple chromosomes and strains. For each genome, the database provides statistical parameters for all the genes, as well as averages and standard deviations of G+C content, codon usage, relative synonymous codon usage and amino-acid content. It also provides information about correspondence analyses of the codon usage, plus lists of extraneous group of genes in terms of G+C content and lists of putatively acquired genes. With this information, researchers can explore the G+C content and codon usage of a gene when they find incongruities in sequence-based phylogenetic trees. A search engine that allows searches for gene names or keywords for a specific organism is also available. HGT-DB is freely accessible at http://www.fut.es/~debb/HGT.     

Two Aspects of DNA Base Composition: G+C Content and Translation-Coupled Deviation from Intra-Strand Rule of A=T and G=C

The relative contribution of mutation and selection to the G+C content of DNA was analyzed in bacterial species having widely different G+C contents. The analysis used two methods that were developed previously. The first method was to plot the average G+C content of a set of nucleotides against the G+C content of the third codon position for each gene. This method was used to present the G+C distribution of the third codon position and to assess the relative neutrality of a set of nucleotides to that of the G+C content of the third codon position. The second method was to plot the intrastrand bias of the third codon position from Parity Rule 2 (PR2), where A=T and G=C. It was found that whereas intragenomic distributions of the DNA G+C content of these bacteria are narrow in the majority of species, in some species the G+C content of the minor class of genes distributes over wider ranges than the major class of genes. On the other hand, ubiquitous PR2 biases are amino acid specific and independent of the G+C content of DNA, so that when averaged over the amino acids, the biases are small and not correlated with the DNA G+C content. Therefore, translation coupled PR2-biases are unlikely to explain the wide range of G+C contents among different species. Considering all data available, it was concluded that the amino acid-specific PR2 bias has only a minor effect, if any, on the average G+C content. In addition, PR2 bias patterns of different species show phylogenetic relationships, and the pattern can be as a taxal fingerprint. Received: 5 November 1998 / Accepted: 1 March 1999     

Near Homogeneity of PR2-Bias Fingerprints in the Human Genome and Their Implications in Phylogenetic Analyses

Genes of a multicellular organism are heterogeneous in the G+C content, which is particularly true in the third codon position. The extent of deviation from intra-strand equality rule of A = T and G = C (Parity Rule 2, or PR2) is specific for individual amino acids and has been expressed as the PR2-bias fingerprint. Previous results suggested that the PR2-bias fingerprints tend to be similar among the genes of an organism, and the fingerprint of the organism is specific for different taxa, reflecting phylogenetic relationships of organisms. In this study, using coding sequences of a large number of human genes, we examined the intragenomic heterogeneity of their PR2-bias fingerprints in relation to the G+C content of the third codon position (P ₃ ). Result shows that the PR2-bias fingerprint is similar in the wide range of the G+C content at the third codon position (0.30–0.80). This range covers approximately 89% of the genes, and further analysis of the high G+C range (0.80–1.00), where genes with normal PR2-bias fingerprints and those with anomalous fingerprints are mixed, shows that the total of 95% of genes have the similar finger prints. The result indicates that the PR2-bias fingerprint is a unique property of an organism and represents the overall characteristics of the genome. Combined with the previous results that the evolutionary change of the PR2-bias fingerprint is a slow process, PR2-bias fingerprints may be used for the phylogenetic analyses to supplement and augment the conventional methods that use the differences of the sequences of orthologous proteins and nucleic acids. Potential advantages and disadvantages of the PR2-bias fingerprint analysis are discussed. Received: 21 December 2000 / Accepted: 16 February 2001     

DNA G+C content of the third codon position and codon usage biases of human genes

The human genome, as in other eukaryotes, has a wide heterogeneity in the DNA base composition. The evolutionary basis for this heterogeneity has been unknown. A previous study of the human genome (846 genes analyzed) has shown that, in the major range of the G+C content in the third codon position (0.25-0.75), biases from the Parity Rule 2 (PR2) among the synonymous codons of the four-codon amino acids are similar except in the highest G+C range (Sueoka, N., 1999. Translation-coupled violation of Parity Rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third codon position. Gene 238, 53-58.). PR2 is an intra-strand rule where A=T and G=C are expected when there are no biases between the two complementary strands of DNA in mutation and selection rates (substitution rates). In this study, 14,026 human genes were analyzed. In addition, the third codon positions of two-codon amino acids were analyzed. New results show the following: (a) The G+C contents of the third codon position of human genes are scattered in the G+C range of 0.22-0.96 in the third codon position. (b) The PR2 biases are similar in the range of 0.25-0.75, whereas, in the high G+C range (0.75-0.96; 13% of the genes), the PR2-bias fingerprints are different from those of the major range. (c) Unlike the PR2 biases, the G+C contents of the third codon position for both four-codon and two-codon amino acids are all correlated almost perfectly with the G+C content of the third codon position over the total G+C ranges. These results support the notion that the directional mutation pressure, rather than the directional selection pressure, is mainly responsible for the heterogeneity of the G+C content of the third codon position.     

烫伤大鼠早期口服谷氨酰胺对血浆氨基酸代谢的影响

研究烫伤早期口服谷氨酰胺 (GLN)后GLN及其相关氨基酸的代谢变化。以Wistar大鼠为烫伤模型 ,烫伤后早期口服GLN ,并以 83 5-50型氨基酸自动分析仪的生理体液法测定血浆游离氨基酸。结果烫伤后饲标准饲料(C)组血浆GLN在烫伤后 2d和 5d降低 ,与GLN代谢相关的丙氨酸和氨无显著性变化 ;GLN饲料 (G)组各时相点GLN均增加 ,1 0h和 8d增加显著 ,与GLN代谢相关的丙氨酸和血氨增加显著 ;GLN +精氨酸 (G +A)组GLN在 2d降低 ,血氨升高显著。与C组比 ,G和G +A组血浆总氨基酸、支琏氨基酸、GLN、丙氨酸、r-氨基丁酸和氨均较C组高。提示 ,烫伤后早期口服GLN能提高血GLN和与之代谢相关的氨基酸浓度。     

Effect of taxon sampling on recovering the phylogeny of squamate reptiles based on complete mitochondrial genome and nuclear gene sequence data

The complete nucleotide sequences of the mitochondrial (mt) genomes of three species of squamate lizards: Blanus cinereus (Amphisbaenidae), Anguis fragilis (Anguidae), and Tarentola mauritanica (Geckkonidae) were determined anew. The deduced amino acid sequences of all 13 mt protein-coding genes were combined into a single data set and phylogenetic relationships among main squamate lineages were analyzed under maximum likelihood (ML) and Bayesian Inference (BI). Within Squamata, the monophyly of Iguanidae, Anguimorpha, Amphisbaenia, Gekkota, Serpentes, and Acrodonta received high statistical support with both methods. It is particularly striking that this is the first molecular analysis that recovers the monophyly of Scincomorpha (including Scincidae, Xantusiidae, Cordylidae, and Lacertidae), although it is only supported in the Bayesian analysis, and it is sensitive to changes in the outgroup (see below). Phylogenetic relationships among the main squamate lineages could not be resolved with ML but received strong support with BI (above 95%). The newly reconstructed phylogeny of squamates does not support the Iguania-Scleroglossa split. Acrodonta and Serpentes form a clade, which is the sister group of the remaining squamate lineages. Within these, Gekkota were the first branching out, followed by Amphisbaenia, and a clade including Anguimorpha as sister group of Scincomorpha + Iguanidae. The recovered topology differed substantially from previously reported hypotheses on squamate relationships, and the relative effect of using different outgroups, genes, and taxon samplings were explored. The sister group relationship of Serpentes + Acrodonta, and their relative basal position within Squamata could be due to a long-branch attraction artifact. Phylogenetic relationships among Scincomorpha, Amphisbaenia, and Anguimorpha were found to be rather unresolved. Future improving of squamate phylogenetic relationships would rely on finding snake and acrodont species with slower mt evolutionary rates, ensuring thorough taxon coverage of squamate diversity, and incorporating more nuclear genes with appropriate evolutionary rates.     

Mutational bias affects protein evolution in flowering plants

Amino acid sequences from several thousand homologous gene pairs were compared for two plant genomes, Oryza sativa and Arabidopsis thaliana. The Arabidopsis genes all have similar G+C (guanine plus cytosine) contents, whereas their homologs in rice span a wide range of G+C levels. The results show that those rice genes that display increased divergence in their nucleotide composition (specifically, increased G+C content) showed a corresponding, predictable change in the amino acid compositions of the encoded proteins relative to their Arabidopsis homologs. This trend was not seen in a "control" set of rice genes that had nucleotide contents closer to their Arabidopsis homologs. In addition to showing an overall difference in the amino acid composition of the homologous proteins, we were also able to investigate the biased patterns of amino acid substitution since the divergence of these two species. We found that the amino acid exchange matrix was highly asymmetric when comparing the High G+C rice genes with their Arabidopsis homologs. Finally, we investigated the possible causes of this biased pattern of sequence evolution. Our results indicate that the biased pattern of protein evolution is the consequence, rather than the cause, of the corresponding changes in nucleotide content. In fact, there is an even more marked asymmetry in the patterns of substitution at synonymous nucleotide sites. Surprisingly, there is a very strong negative correlation between the level of nucleotide bias and the length of the coding sequences within the rice genome. This difference in gene length may provide important clues about the underlying mechanisms.     

Phylogenetic analysis of nitrite, nitric oxide, and nitrous oxide respiratory enzymes reveal a complex evolutionary history for denitrification

Denitrification is a facultative respiratory pathway in which nitrite (NO2(-)), nitric oxide (NO), and nitrous oxide (N2O) are successively reduced to nitrogen gas (N(2)), effectively closing the nitrogen cycle. The ability to denitrify is widely dispersed among prokaryotes, and this polyphyletic distribution has raised the possibility of horizontal gene transfer (HGT) having a substantial role in the evolution of denitrification. Comparisons of 16S rRNA and denitrification gene phylogenies in recent studies support this possibility; however, these results remain speculative as they are based on visual comparisons of phylogenies from partial sequences. We reanalyzed publicly available nirS, nirK, norB, and nosZ partial sequences using Bayesian and maximum likelihood phylogenetic inference. Concomitant analysis of denitrification genes with 16S rRNA sequences from the same organisms showed substantial differences between the trees, which were supported by examining the posterior probability of monophyletic constraints at different taxonomic levels. Although these differences suggest HGT of denitrification genes, the presence of structural variants for nirK, norB, and nosZ makes it difficult to determine HGT from other evolutionary events. Additional analysis using phylogenetic networks and likelihood ratio tests of phylogenies based on full-length sequences retrieved from genomes also revealed significant differences in tree topologies among denitrification and 16S rRNA gene phylogenies, with the exception of the nosZ gene phylogeny within the data set of the nirK-harboring genomes. However, inspection of codon usage and G + C content plots from complete genomes gave no evidence for recent HGT. Instead, the close proximity of denitrification gene copies in the genomes of several denitrifying bacteria suggests duplication. Although HGT cannot be ruled out as a factor in the evolution of denitrification genes, our analysis suggests that other phenomena, such gene duplication/divergence and lineage sorting, may have differently influenced the evolution of each denitrification gene.     

最大似然法和贝叶斯推论与似然比检验探讨泽泻目系统发育关系

应用最大似然法(ML)、贝叶斯推论(BI)、邻接法(NJ)和似然比检验(hLRTs)进行泽泻目分子系统学研究。所用的rbcL基因序列代表了泽泻目14科46属以及作为外类群的6相关属。研究结果表明,*等级制似然比检验表明泽泻目rbcL序列最适合的DNA进化模型为GTR+I+G,最大似然法、贝叶斯法和邻接法构建的系统发育树拓扑结构相似,没有显著的差异,但贝叶斯树支持率较高;泽泻目为一单系类群,由两个主要谱系分支构成,深层分布格局由5个主要分支构成。基于分子系统发育树,文中对泽泻目科间、水鳖科+茨藻科、泽泻科+花蔺科+黄花蔺科、和"Cymodoeaceae complex"的系统发育关系进行了讨论。研究结果还表明,泽泻目系统发育关系可能还需要更多的证据进一步的澄清。