Evolutionary rates and expression level in Chlamydomonas

In many biological systems, especially bacteria and unicellular eukaryotes, rates of synonymous and nonsynonymous nucleotide divergence are negatively correlated with the level of gene expression, a phenomenon that has been attributed to natural selection. Surprisingly, this relationship has not been examined in many important groups, including the unicellular model organism Chlamydomonas reinhardtii. Prior to this study, comparative data on protein-coding sequences from C. reinhardtii and its close noninterfertile relative C. incerta were very limited. We compiled and analyzed protein-coding sequences for 67 nuclear genes from these taxa; the sequences were mostly obtained from the C. reinhardtii EST database and our C. incerta EST data. Compositional and synonymous codon usage biases varied among genes within each species but were highly correlated between the orthologous genes of the two species. Relative rates of synonymous and nonsynonymous substitution across genes varied widely and showed a strong negative correlation with the level of gene expression estimated by the codon adaptation index. Our comparative analysis of substitution rates in introns of lowly and highly expressed genes suggests that natural selection has a larger contribution than mutation to the observed correlation between evolutionary rates and gene expression level in Chlamydomonas.     

Independent terrestrial origins of the Halosphaeriales (marine Ascomycota)

A phylogenetic study of marine ascomycetes was initiated to test and refine evolutionary hypotheses of marine-terrestrial transitions among ascomycetes. Taxon sampling focused on the Halosphaeriales, the largest order of marine ascomycetes. Approximately 1050 base pairs (bp) of the gene that codes for the nuclear small subunit (SSU) and 600 bp of the gene that codes for the nuclear large subunit (LSU) ribosomal RNAs (rDNA) were sequenced for 15 halosphaerialean taxa and integrated into a data set of homologous sequences from terrestrial ascomycetes. An initial set of phylogenetic analyses of the SSU rDNA from 38 taxa representing 15 major orders of the phylum Ascomycota confirmed a close phylogenetic relationship of the halosphaerialean species with several other orders of perithecial ascomycetes. A second set of analyses, which involved more intensive taxon sampling of perithecial ascomycetes, was performed using the SSU and LSU rDNA data in combined analyses. These second analyses included 15 halosphaerialean taxa, 26 terrestrial perithecial fungi from eight orders, and five outgroup taxa from the Pezizales. In these analyses the Halosphaeriales were polyphyletic and comprised two distinct lineages. One clade of Halosphaeriales comprised 12 taxa from 11 genera and was most closely related to terrestrial fungi of the Microascales. The second clade of halosphaerialean fungi comprised taxa from the genera Lulworthia and Lindra and was an isolated lineage among the perithecial fungi. Both the main clade of Halosphaeriales and the Lulworthia/Lindra clade are supported by the data as being independently derived from terrestrial ancestors.     

Rate acceleration and long-branch attraction in a conserved gene of cryptic daphniid (Crustacea) species.

The nuclear large subunit (LSU) rRNA gene is a rich source of phylogenetic characters because of its large size, mosaic of slowly and rapidly evolving regions, and complex secondary structure variation. Nevertheless, many studies have indicated that inconsistency, bias, and gene-specific error (e.g., within-individual gene family variation, cryptic sequence simplicity, and sequence coevolution) can complicate animal phylogenies based on LSU rDNA sequences. However, most of these studies sampled small gene fragments from expansion segments--among animals only five nonchordate complete LSU sequences are published. In this study, we sequenced near-complete nuclear LSU genes from 11 representative daphniids (Crustacea). The daphniid expansion segment V6 was larger and showed more length variation (90-351 bp) than is found in all other reported LSU V6 sequences. Daphniid LSU (without the V6 region) phylogenies generally agreed with the existing phylogenies based on morphology and mtDNA sequences. Nevertheless, a major disagreement between the LSU and the expected trees involved a positively misleading association between the two taxa with the longest branches, Daphnia laevis and D. occidentalis. Both maximum parsimony (MP) and maximum likelihood (ML) optimality criteria recovered this association, but parametric simulations indicated that MP was markedly more sensitive to this bias than ML. Examination of data partitions indicated that the inconsistency was caused by increased nucleotide substitution rates in the branches leading to D. laevis and D. occidentalis rather than among-taxon differences in base composition or distribution of sites that are free to vary. These results suggest that lineage-specific rate acceleration can lead to long-branch attraction even in the conserved genes of animal species that are almost morphologically indistinguishable.     

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (< 1-2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222-705 ka), Neandertals (108 ka; 70-156 ka), and modern humans (76 ka; 47-110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.     

Phylogenetic relationships of Gomphillaceae and Asterothyriaceae: evidence from a combined Bayesian analysis of nuclear and mitochondrial sequences

The phylogeny and systematic position of Gomphillaceae was reconstructed using a combined Bayesian analysis of nuclear LSU rDNA and mitochondrial SSU rDNA sequences. Twenty-four partial sequences of 12 taxa (11 Gomphillaceae and one Asterothyriaceae) plus two new sequences of Stictis radiata (Ostropales outgroup) were generated and aligned with the corresponding sequences retrieved from GenBank, resulting in an alignment of 82 taxa that was analyzed using a Bayesian approach with Markov chain Monte Carlo (B/MCMC) methods. Our results confirm Gomphillaceae sensu Vezda and Poelt plus Asterothyriaceae to be a monophyletic group, with an unresolved relationship between the two families. Placement of Gomphillaceae and Asterothyriaceae within Ostropales sensu Kauff and Lutzoni, as sister of Thelotremataceae, also is strongly supported. Alternative hypotheses placing Gomphillaceae in Lecanorales (Cladoniaceae), Agyriales (Baeomycetaceae) or within bitunicate Ascomycota (Arthoniomycetes, Chaetothyriomycetes, Dothideomycetes) were rejected with our dataset. After recent synonymization of Dimerella with Coenogonium (Ostropales: Coenogoniaceae), we propose the new combination Coenogonium pineti (one of our Ostropales outgroup taxa in this analysis).     

Evolutionary rates of mitochondrial genomes correspond to diversification rates and to contemporary species richness in birds and reptiles

Rates of biological diversification should ultimately correspond to rates of genome evolution. Recent studies have compared diversification rates with phylogenetic branch lengths, but incomplete phylogenies hamper such analyses for many taxa. Herein, we use pairwise comparisons of confamilial sauropsid (bird and reptile) mitochondrial DNA (mtDNA) genome sequences to estimate substitution rates. These molecular evolutionary rates are considered in light of the age and species richness of each taxonomic family, using a random-walk speciation–extinction process to estimate rates of diversification. We find the molecular clock ticks at disparate rates in different families and at different genes. For example, evolutionary rates are relatively fast in snakes and lizards, intermediate in crocodilians and slow in turtles and birds. There was also rate variation across genes, where non-synonymous substitution rates were fastest at ATP8 and slowest at CO3. Family-by-gene interactions were significant, indicating that local clocks vary substantially among sauropsids. Most importantly, we find evidence that mitochondrial genome evolutionary rates are positively correlated with speciation rates and with contemporary species richness. Nuclear sequences are poorly represented among reptiles, but the correlation between rates of molecular evolution and species diversification also extends to 18 avian nuclear genes we tested. Thus, the nuclear data buttress our mtDNA findings.     

Algal and Fungal Diversity in Antarctic Lichens

The composition of lichen ecosystems except mycobiont and photobiont has not been evaluated intensively. In addition, recent studies to identify algal genotypes have raised questions about the specific relationship between mycobiont and photobiont. In the current study, we analyzed algal and fungal community structures in lichen species from King George Island, Antarctica, by pyrosequencing of eukaryotic large subunit (LSU) and algal internal transcribed spacer (ITS) domains of the nuclear rRNA gene. The sequencing results of LSU and ITS regions indicated that each lichen thallus contained diverse algal species. The major algal operational taxonomic unit (OTU) defined at a 99% similarity cutoff of LSU sequences accounted for 78.7–100% of the total algal community in each sample. In several cases, the major OTUs defined by LSU sequences were represented by two closely related OTUs defined by 98% sequence similarity of ITS domain. The results of LSU sequences indicated that lichen‐associated fungi belonged to the Arthoniomycetes, Eurotiomycetes, Lecanoromycetes, Leotiomycetes, and Sordariomycetes of the Ascomycota, and Tremellomycetes and Cystobasidiomycetes of the Basidiomycota. The composition of major photobiont species and lichen‐associated fungal community were mostly related to the mycobiont species. The contribution of growth forms or substrates on composition of photobiont and lichen‐associated fungi was not evident.     

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.     

Assessment of soil fungal communities using pyrosequencing

Pyrosequencing, a non-electrophoretic method of DNA sequencing, was used to investigate the extensive fungal community in soils of three islands in the Yellow Sea of Korea, between Korea and China. Pyrosequencing was carried out on amplicons derived from the 5′ region of 18S rDNA. A total of 10,166 reads were obtained, with an average length of 103 bp. The maximum number of fungal phylotypes in soil predicted at 99% similarity was 3,334. The maximum numbers of phylotypes predicted at 97% and 95% similarities were 736 and 286, respectively. Through phylogenetic assignment using BLASTN, a total of 372 tentative taxa were identified. The majority of true fungal sequences recovered in this study belonged to the Ascomycota (182 tentative taxa in 2,708 reads) and Basidiomycota (172 tentative taxa in 6,837 reads). The predominant species of Ascomycota detected have been described as lichen-forming fungi, litter/wood decomposers, plant parasites, endophytes, and saprotrophs: Peltigera neopolydactyla (Lecanoromycetes), Paecilomyces sp. (Sordariomycetes), Phacopsis huuskonenii (Lecanoromycetes), and Raffaelea hennebertii (mitosporicAscomycota). The majority of sequences in the Basidiomycota matched ectomycorrhizal and wood rotting fungi, including species of the Agaricales and Aphyllophorales, respectively. A high number of sequences in the Thelephorales, Boletales, Stereales, Hymenochaetales, and Ceratobasidiomycetes were also detected. By applying high-throughput pyrosequencing, we observed a high diversity of soil fungi and found evidence that pyrosequencing is a reliable technique for investigating fungal communities in soils.     

Molecular phylogeny and paclitaxel screening of fungal endophytes from Taxus globosa

We studied the endophytic mycoflora associated with Taxus globosa, the Mexican yew. The study localities; Las Avispas (LA), San Gaspar (SG), and La Mina (LM) were three segments of cloud forest within the range of Sierra Gorda Biosphere Reserve, México. Overall, 245 endophytes were isolated and 105 representative Ascomycota (morphotaxons) were chosen for phylogenetic and genotypic characterization. Maximum likelihood analyses of large subunit of ribosomal RNA (LSU) rDNA showed well-supported clades of Dothideomycetes, Eurotiomycetes, Leotiomycetes, Pezizomycetes, and Sordariomycetes. Analyses of ITS rDNA groups showed 57 genotypes (95% sequence similarity), in general consistent with the phylogenetically delimitated taxa based on LSU rDNA sequences. The endophyte diversity measured by Fisher's α, Shanonn, and Simpson indices was ca. three-fold and ca. two-fold greater in LM than in LA and SG respectively. A screening for paclitaxel using a competitive inhibition enzyme immunoassay showed 16 positive isolates producing between 65 and 250 ng l(-1). The isolates included Acremonium, Botryosphaeria, Fusarium, Gyromitra, Nigrospora, Penicillium, three novel Pleosporales, and Xylaria.     

Fungal endophyte diversity in Sarracenia

Fungal endophytes were isolated from 4 species of the carnivorous pitcher plant genus Sarracenia: S. minor, S. oreophila, S. purpurea, and S. psittacina. Twelve taxa of fungi, 8 within the Ascomycota and 4 within the Basidiomycota, were identified based on PCR amplification and sequencing of the internal transcribed spacer sequences of nuclear ribosomal DNA (ITS rDNA) with taxonomic identity assigned using the NCBI nucleotide megablast search tool. Endophytes are known to produce a large number of metabolites, some of which may contribute to the protection and survival of the host. We speculate that endophyte-infected Sarracenia may benefit from their fungal associates by their influence on nutrient availability from within pitchers and, possibly, by directly influencing the biota within pitchers.     

Phylogenetic analysis of ribosomal RNA operons from uncultivated coastal marine bacterioplankton

Analyses of small subunit ribosomal RNA genes (SSU rDNAs) have significantly influenced our understanding of the composition of aquatic microbial assemblages. Unfortunately, SSU rDNA sequences often do not have sufficient resolving power to differentiate closely related species. To address this general problem for uncultivated bacterioplankton taxa, we analysed and compared sequences of polymerase chain reaction (PCR)-generated and bacterial artificial chromosome (BAC)-derived clones that contained most of the SSU rDNAs, the internal transcribed spacer (ITS) and the large subunit ribosomal RNA gene (LSU rDNA). The phylogenetic representation in the rRNA operon PCR library was similar to that reported previously in coastal bacterioplankton SSU rDNA libraries. We observed good concordance between the phylogenetic relationships among coastal bacterioplankton inferred from SSU or LSU rDNA sequences. ITS sequences confirmed the close intragroup relationships among members of the SAR11, SAR116 and SAR86 clades that were predicted by SSU and LSU rDNA sequence analyses. We also found strong support for homologous recombination between the ITS regions of operons from the SAR11 clade.     

megB1, a novel macroevolutionary genomic marker of the fungal phylum Basidiomycota

Molecular studies on the evolution and systematics of fungi have been established primarily based on the neutral theory by analyzing neutral mutations in some defined segments of housekeeping genes as genetic markers. Such an approach is, however, hardly applicable for analyzing ancient evolutionary radiations. In the present study, we looked for DNA sequences characterizing higher taxa, and discovered a unique macroevolutionary genomic marker, megB1, that specifies the phylum Basidiomycota. megB1 is an approximately 500-bp DNA element, which is defined by terminal sequences and five internal segments conserved throughout the phylum. megB1 resides on the rDNA intergenic spacer 1 (IGS1) from 27 species of 10 Basidiomycota genera examined. While megB1 was not found in IGS1 from the other 92 species of the 27 Basidiomycota genera, several genera representing them carry megB1 in some other genomic regions. No known taxonomic criteria fit into the classification on the basis of whether megB1 resides on rDNA. Neighbor-joining analysis of the megB1 sequence, however, properly assigned species to their respective genera. Thus far, megB1 has not been found in any genomic or genetic databases currently available for other phyla. These results suggest that megB1 may have emerged upon the occurrence of Basidiomycota, and that this phylum evolved thereafter leaving this element conserved throughout their further differentiation. megB1 may be a novel genomic marker useful in the analysis of ancient through the latest evolutionary radiation in Basidiomycota.     

Multigene molecular phylogenetics reveals true morels (Morchella) are especially species-rich in China

The phylogenetic diversity of true morels (Morchella) in China was estimated by initially analyzing nuclear ribosomal internal transcribed spacer (ITS) rDNA sequences from 361 specimens collected in 21 provinces during the 2003-2011 growing seasons, together with six collections obtained on loan from three Chinese herbaria. Based on the results of this preliminary screen, 40 Esculenta Clade (yellow morels) and 30 Elata Clade (black morels) were chosen to represent the full range of phylogenetic diversity sampled. To investigate their species limits, we generated DNA sequences from portions of three protein-coding genes (RPB1, RPB2 and EF-1α) and domains D1 and D2 of the nuclear large subunit (LSU) rDNA for all 70 collections. To fully assess evolutionary relationships, previously published multilocus DNA sequence data representing all known Morchella species was included in this study. Phylogenetic analyses employing maximum parsimony and maximum likelihood frameworks resolved 30 species in China compared with 22 in Europe and 19 within North America. Eleven novel phylogenetically distinct species were discovered in China, including two species within the Elata Clade and nine within the Esculenta Clade. Of the 30 species in China, 20 appear to be endemic, nine were also represented in Europe, and four putatively fire-adapted species have disjunct distributions in China, Europe and western North America. Although the diversification time estimates place the Esculenta Clade in China as early as the late Cretaceous and the Elata Clade by the early Oligocene, 27 of the 30 species evolved between the middle Miocene 12Mya and present.     

Heterogeneous Rates of Molecular Evolution Among Cryptic Species of the Ciliate Morphospecies Chilodonella uncinata

While molecular analyses have provided insight into the phylogeny of ciliates, the few studies assessing intraspecific variation have largely relied on just a single locus [e.g., nuclear small subunit rDNA (nSSU-rDNA) or mitochondrial cytochrome oxidase I]. In this study, we characterize the diversity of several nuclear protein-coding genes plus both nSSU-rDNA and mitochondrial small subunit rDNA (mtSSU-rDNA) of five isolates of the ciliate morphospecies Chilodonella uncinata. Although these isolates have nearly identical nSSU-rDNA sequences, they differ by up to 8.0% in mtSSU-rDNA. Comparative analyses of all loci, including β-tubulin paralogs, indicate a lack of recombination between strains, demonstrating that the morphospecies C. uncinata consists of multiple cryptic species. Further, there is considerable variation in substitution rates among loci as some protein-coding domains are nearly identical between isolates, while others differ by up to 13.2% at the amino acid level. Combining insights on macronuclear variation among isolates, the focus of this study, with published data from the micronucleus of two of these isolates, indicates that C. uncinata lineages are able to maintain both highly divergent and highly conserved genes within a rapidly evolving germline genome.     

Lineage-specific variations of congruent evolution among DNA sequences from three genomes,and relaxed selective constraints on <Emphasis Type="Italic">rbc</Emphasis>L in <Emphasis Type="Italic">Cryptomonas</Emphasis> (Cryptophyceae)

Background  

Plastid-bearing cryptophytes like Cryptomonas contain four genomes in a cell, the nucleus, the nucleomorph, the plastid genome and the mitochondrial genome. Comparative phylogenetic analyses encompassing DNA sequences from three different genomes were performed on nineteen photosynthetic and four colorless Cryptomonas strains. Twenty-three rbc L genes and fourteen nuclear SSU rDNA sequences were newly sequenced to examine the impact of photosynthesis loss on codon usage in the rbc L genes, and to compare the rbc L gene phylogeny in terms of tree topology and evolutionary rates with phylogenies inferred from nuclear ribosomal DNA (concatenated SSU rDNA, ITS2 and partial LSU rDNA), and nucleomorph SSU rDNA.     

A multigene molecular phylogenetic assessment of true morels (Morchella) in Turkey

A collection of 247 true morels (Morchella spp.) primarily from the Mediterranean and Aegean Regions of Southern Turkey, were analyzed for species diversity using partial RNA polymerase I (RPB1) and nuclear ribosomal large subunit (LSU) rDNA gene sequences. Based on the result of this initial screen, 62 collections representing the full range of genetic diversity sampled were subjected to multigene phylogenetic species recognition based on genealogical concordance (GCPSR). The 62-taxon dataset consisted of partial sequences from three nuclear protein-coding genes, RNA polymerase I (RPB1), RNA polymerase II (RPB2), translation elongation factor (EF1-α), and partial LSU rDNA gene sequences. Phylogenetic analyses of the individual and combined datasets, using maximum parsimony (MP) and maximum likelihood (ML), yielded nearly fully resolved phylogenies that were highly concordant topologically. GCPSR analysis of the 62-taxon dataset resolved 15 putative phylogenetically distinct species. The early diverging Elata (black morels) and Esculenta Clades (yellow morels) were represented, respectively, by 13 and two species. Because a Latin binomial can be applied with confidence to only one of the 15 species (Morchella semilibera), species were identified by clade (Mel for Elata and Mes for Esculenta) followed by a unique Arabic number for each species within these two clades. Eight of the species within the Elata Clade appear to be novel, including all seven species within the Mel-20-to-31 subclade and its sister designated Mel-25. Results of the present study provide essential data for ensuring the sustainability of morel harvests through the formulation of sound conservation policies.     

Phylogenetic analysis of the Metastrongyloidea (Nematoda: Strongylida) inferred from ribosomal RNA gene sequences

Phylogenetic relationships among nematodes of the strongylid superfamily Metastrongyloidea were analyzed using partial sequences from the large-subunit ribosomal RNA (LSU rRNA) and small-subunit ribosomal RNA (SSU rRNA) genes. Regions of nuclear ribosomal DNA (rDNA) were amplified by polymerase chain reaction, directly sequenced, aligned, and phylogenies inferred using maximum parsimony. Phylogenetic hypotheses inferred from the SSU rRNA gene supported the monophyly of representative taxa from each of the 7 currently accepted metastrongyloid families. Metastrongyloid taxa formed the sister group to representative trichostrongyloid sequences based on SSU data. Sequences from either the SSU or LSU RNA regions alone provided poor resolution for relationships within the Metastrongyloidea. However, a combined analysis using sequences from all rDNA regions yielded 3 equally parsimonious trees that represented the abursate Filaroididae as polyphyletic, Parafilaroides decorus as the sister species to the monophyletic Pseudaliidae, and a sister group relationship between Oslerus osleri and Metastrongylus salmi. Relationships among 3 members of the Crenosomatidae, and 1 representative of the Skrjabingylidae (Skrjabingylus chitwoodorum) were not resolved by these combined data. However, members of both these groups were consistently resolved as the sister group to the other metastrongyloid families. These relationships are inconsistent with traditional classifications of the Metastrongyloidea and existing hypotheses for their evolution.     

Advantages of a mechanistic codon substitution model for evolutionary analysis of protein-coding sequences

BACKGROUND: A mechanistic codon substitution model, in which each codon substitution rate is proportional to the product of a codon mutation rate and the average fixation probability depending on the type of amino acid replacement, has advantages over nucleotide, amino acid, and empirical codon substitution models in evolutionary analysis of protein-coding sequences. It can approximate a wide range of codon substitution processes. If no selection pressure on amino acids is taken into account, it will become equivalent to a nucleotide substitution model. If mutation rates are assumed not to depend on the codon type, then it will become essentially equivalent to an amino acid substitution model. Mutation at the nucleotide level and selection at the amino acid level can be separately evaluated. RESULTS: The present scheme for single nucleotide mutations is equivalent to the general time-reversible model, but multiple nucleotide changes in infinitesimal time are allowed. Selective constraints on the respective types of amino acid replacements are tailored to each gene in a linear function of a given estimate of selective constraints. Their good estimates are those calculated by maximizing the respective likelihoods of empirical amino acid or codon substitution frequency matrices. Akaike and Bayesian information criteria indicate that the present model performs far better than the other substitution models for all five phylogenetic trees of highly-divergent to highly-homologous sequences of chloroplast, mitochondrial, and nuclear genes. It is also shown that multiple nucleotide changes in infinitesimal time are significant in long branches, although they may be caused by compensatory substitutions or other mechanisms. The variation of selective constraint over sites fits the datasets significantly better than variable mutation rates, except for 10 slow-evolving nuclear genes of 10 mammals. An critical finding for phylogenetic analysis is that assuming variable mutation rates over sites lead to the overestimation of branch lengths.     

Suitability of chloroplast LSU rDNA and its diverse group I introns for species recognition and phylogenetic analyses of lichen-forming Trebouxia algae

To date, species identification of lichen photobionts has been performed principally on the basis of microscopic examinations and molecular data from nuclear-encoded genes. In plants, the chloroplast genome has been more readily exploited than the nuclear genome for systematic investigations. At the present time, very little information is available about the chloroplast genome of lichen-forming algae. For this reason, we have sequenced a portion of the gene encoding for the chloroplast large sub-unit rRNA (LSU rDNA) as a new molecular marker. Sequencing of the chloroplast LSU rDNAs revealed the existence of an unusual diversity of group I introns (a total of 31) within 15 analyzed Trebouxia species. The number, sequence and insertion site of these introns were very different among species, contributing to their recognition. A relatively large intron-free portion of the chloroplast LSU rDNA and part of the nuclear ribosomal cistron (18S–5.8S–26S) between the nuclear internal transcribed spacers (nrITS) were subjected to phylogenetic analyses. The obtained results indicate that data combination from both nuclear and chloroplast sequences can improve phylogenetic accuracy. Herein, we propose the suitability of both intronic and exonic sequences of the chloroplast LSU rDNA for species recognition, and an exonic sequence spanning from position 879 to 1837 in the Escherichia coli 23S rDNA for phylogenetic analyses of Trebouxia phycobionts.