A chloroplast DNA mutational hotspot and gene conversion in a noncoding region near rbcL in the grass family (Poaceae)

The tortoise and the hare: choosing between noncoding plastome and nuclear Adh sequences for phylogeny reconstruction in a recently diverged plant group

Phylogenetic resolution is often low within groups of recently diverged taxa due to a paucity of phylogenetically informative characters. We tested the relative utility of seven noncoding cpDNA regions and a pair of homoeologous nuclear genes for resolving recent divergences, using tetraploid cottons (Gossypium) as a model system. The five tetraploid species of Gossypium are a monophyletic assemblage derived from an allopolyploidization event that probably occurred within the last 0.5-2 million years. Previous analysis of cpDNA restriction site data provided only partial resolution within this clade despite a large number of enzymes employed. We sequenced three cpDNA introns (rpl16, rpoC1, ndhA) and four cpDNA spacers (accD-psaI, trnL-trnF, trnT-trnL, atpB-rbcL) for a total of over 7 kb of sequence per taxon, yet obtained only four informative nucleotide substitutions (0.05%) resulting in incomplete phylogenetic resolution. In addition, we sequenced a 1.65-kb region of a homoeologous pair of nuclear-encoded alcohol dehydrogenase (Adh) genes. In contrast with the cpDNA sequence data, the Adh homoeologues yielded 25 informative characters (0.76%) and provided a robust and completely resolved topology that is concordant with previous cladistic and phenetic analyses. The enhanced resolution obtained using the nuclear genes reflects an approximately three- to sixfold increase in nucleotide substitution rate relative to the plastome spacers and introns.

Molecular evolution and phylogenetic utility of the chloroplast rpl16 intron in Chusquea and the Bambusoideae (Poaceae)

Phylogenetic relationships within Chusquea, a diverse genus of neotropical woody bamboos, and among selected members of the Bambusoideae were explored using rpl16 intron sequence data from the chloroplast genome. Mechanisms of mutation, including slipped-strand mispairing, secondary structure, minute inversions, and base substitutions, were examined within the rpl16 intron, and their effects on sequence alignment and phylogenetic analysis were investigated. Thirty-five bamboo sequences were generated and two separate matrices were analyzed using maximum parsimony. In the first, 23 sequences from Chusquea, 1 of Neurolepis, and 3 outgroups were included. Neurolepis was supported as sister to Chusquea, Chusquea was strongly supported as a monophyletic lineage, and three species of Chusquea subg. Rettbergia were resolved as the most basal clade within the genus. In the second analysis, 15 sequences, 14 from across the subfamily and 1 outgroup, were included. A Bambusoideae clade was recovered with the Olyreae/Parianeae (herbaceous bamboos) and the Bambuseae (woody bamboos) each supported as monophyletic. Two clades corresponding to temperate and tropical woody bamboos were derived within the Bambuseae and the tropical taxa were further split into New World and Old World clades. The rpl16 intron in bamboos was found to be susceptible to frequent length mutations of multiple origins, nonindependent character evolution, and regions of high mutability, all of which created difficulties in alignment and phylogenetic analysis; nonetheless the rpl16 intron is phylogenetically informative at the inter- and intrageneric levels in bamboos.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

Analysis of a hotspot for deletion formation within the intron of the chloroplast trnI gene

The chloroplast genomes of higher plants encode several tRNA genes that contain highly conserved type II introns. Using primers specific to conserved 5' and 3' regions within the introns of the genes trnA (tRNA-ala) and trnI (tRNA-ile) we have PCR amplified parts of these introns from 36 plant species representing a wide range of plant families. Deletions were found in the introns of both tRNA genes. Fourteen species had detectable deletions within the intron of trnI and four species within the intron of trnA. The occurrence of these deletions among the various plant families suggests that the events leading to the formation of these deletions occurred independently many times during the evolution of higher plants. Analysis of the amplified PCR products from the trnI intron suggests that these independent deletions may not be random but appear to fall into two size classes. Several members of each class were cloned and sequenced and the end points of the deletions were mapped. The 3' ends of all deletions studied terminate within the same short region. The 5' ends of the deletions map to two different regions, giving rise to the two size classes. These two 5' deletion endpoint regions show some sequence similarity. Only two of the identified deletions contain directly repeated sequences at the deletion endpoints, a feature associated with homologous recombination. Our results suggest that within the trnI intron, there are preferred sites or "hotspots" for deletion formation involving a novel imprecise recombination mechanism. The significance of these sequences and possible mechanisms for deletion formation are discussed.

Ebb and flow of the chloroplast inverted repeat

The endpoints of the large inverted repeat (IR) of chloroplast DNA in flowering plants differ by small amounts between species. To quantify the extent of this movement and define a possible mechanism for IR expansion, DNA sequences across the IR-large single-copy (IR-LSC) junctions were compared among 13 Nicotiana species and other dicots. In most Nicotiana species the IR terminates just upstream of, or somewhere within, the 5' portion of the rps19 gene. The truncated copy of this gene, rps19', varies in length even between closely related species but is of constant size within a single species. In Nicotiana, six different rps19' structures were found. A phylogenetic tree of Nicotiana species based on restriction site data shows that the IR has both expanded and contracted during the evolution of this genus. Gene conversion is proposed to account for these small and apparently random IR expansions. A large IR expansion of over 12 kb has occurred in Nicotiana acuminata. The new IR-LSC junction in this species lies within intron 1 of the clpP gene. This rearrangement occurred via a double-strand DNA break and recombination between poly (A) tracts in clpP intron 1 and upstream of rps19. Nicotiana acuminata chloroplast DNA contains a "molecular fossil' of the IR-LSC junction that existed prior to this dramatic rearrangement.

Phylogenetic affinities of the grasses to other monocots as revealed by molecular analysis of chloroplast DNA

The distribution of structural alterations of the chloroplast genome found in grass chloroplast (cp) DNA in comparison with that of tobacco was systematically surveyed in the cpDNAs of monocots. Southern hybridization and/or PCR analyses for the detection of (1) three inversions in the large single-copy region, (2) loss of an intron in the rpoC1 gene, (3) an extra-sequence insertion in the rpoC2 gene, (4) the deletion of ORF2280, (5) rearrangements of the accD (ORF512) gene, and (6) non-reciprocal translocation of the rpl23 gene, were carried out on cpDNAs isolated from 58 species, 22 families, and 11 orders, which covered almost all families of monocots. These structural alterations of cpDNA mostly occurred at the family level. However, only part of the Restionaceae possessed the inversion that characterizes the lineage of grass differentiation. The order of mutational events made it possible to reconstruct grass phylogeny in monocots. Since no variations in structural alterations of the cpDNA were found among the Poaceae, grass plants were inferred to have originated from an ancestor harboring these structural alterations of the chloroplast genome. These phylogenetic relationships were supported by the sequence data of rbcL.

Phylogenetic relationships within Hevea brasiliensis as deduced from a polymorphic mitochondrial DNA region

We have cloned a 4.5-kb mtDNA fragment showing a high RFLP polymorphism between various Hevea genotypes. Subcloning and sequencing of a 1.4-kb segment of this clone allowed us to design PCR amplification primers to isolate homologous mtDNA segments of about 0.9 kb from 23 representative genotypes of Hevea. Complete sequences from 4 genotypes showed between 6.7% and 20.2% of nucleotide diversity, suggesting the presence of a hypervariable, or hotspot, region. A sequence of 345 nucleotides within this region was determined for the 23 genotypes. The phylogenetic relationships inferred from the sequence comparison are in general agreement with the results obtained from mtDNA RFLP analysis, indicating that this polymorphic mtDNA region is a useful molecular marker for phylogenetic analysis within Hevea.

Neighboring base composition and transversion/transition bias in a comparison of rice and maize chloroplast noncoding regions

The correspondence between the transversion/transition ratio and the neighboring base composition in chloroplast DNA is examined. For 18 noncoding regions of the chloroplast genome, alignments between rice (Oryza sativa) and maize (Zea mays) were generated by two different methods. Difficulties of aligning noncoding DNA are discussed, and the alignments are analyzed in a manner that reduces alignment artifacts. Sequence divergence is < 10%, so multiple substitutions at a site are assumed to be rare. Observed substitutions were analyzed with respect to the A+T content of the two immediately flanking bases. It is shown that as this content increases, the proportion of transversions also increases. When both the 5'- and 3'-flanking nucleotides are G or C (A+T content of 0), only 25% of the observed substitutions are transversions. However, when both the 5'- and 3'-flanking nucleotides are A or T (A+T content of 2), 57% of the observed substitutions are transversions. Therefore, the influence of flanking base composition on substitutions, previously reported for a single noncoding region, is a general feature of the chloroplast genome.

Extensive sequence divergence in the 3' inverted repeat of the chloroplast rbcL gene in non-flowering land plants and algae

A stem-loop region is present at the 3' terminus of the chloroplast rbcL mRNA in all taxa surveyed to date. In spinach, this structure has been shown by others to be involved in modulating transcript stability and correct 3' terminus processing, and is a conserved feature of other flowering plant rbcL mRNAs. In Chlamydomonas reinhardtii, an analogous structure has been shown by others to serve as a transcription terminator. Our sequencing data have shown that this region is highly divergent in several non-flowering land plants, as evidenced by representatives from the ferns, conifers, 'fern-allies' and liverworts. To extend our analysis, a computer-assisted survey of the stem-loop region of the 3' flanking region of published chloroplast rbcL genes was undertaken. The flowering plant rbcL inverted repeats (IR) were remarkably conserved in sequence, allowing for precise multiple alignments of both monocot and dicot sequences within a single matrix. Surprisingly, sequences obtained from non-flowering land plants, algae, photosynthetic protists and photosynthetic prokaryotes were extremely variant, in terms of both sequence composition and thermodynamic parameters.

Rates and patterns of chloroplast DNA evolution

The chloroplast genome (cpDNA) of plants has been a focus of research in plant molecular evolution and systematics. Several features of this genome have facilitated molecular evolutionary analyses. First, the genome is small and constitutes an abundant component of cellular DNA. Second, the chloroplast genome has been extensively characterized at the molecular level providing the basic information to support comparative evolutionary research. And third, rates of nucleotide substitution are relatively slow and therefore provide the appropriate window of resolution to study plant phylogeny at deep levels of evolution. Despite a conservative rate of evolution and a relatively stable gene content, comparative molecular analyses reveal complex patterns of mutational changes. Non-coding regions of cpDNA diverge through insertion/deletion changes that are sometimes site dependent. Coding genes exhibit different patterns of codon bias that appear to violate the equilibrium assumptions of some evolutionary models. Rates of molecular change often vary among plant families and orders in a manner that violates the assumption of a simple molecular clock. Finally, protein-coding genes exhibit patterns of amino acid change that appear to depend on protein structure, and these patterns may reveal subtle aspects of structure/function relationships. Only comparative studies of molecular sequences have the resolution to reveal this underlying complexity. A complete description of the complexity of molecular change is essential to a full understanding of the mechanisms of evolutionary change and in the formulation of realistic models of mutational processes.