ITS secondary structure derived from comparative analysis: implications for sequence alignment and phylogeny of the Asteraceae

An RNA secondary structure model is presented for the nuclear ribosomal internal transcribed spacers (ITS) based on comparative analysis of 340 sequences from the angiosperm family Asteraceae. The model based on covariation analysis agrees with structural features proposed in previous studies using mainly thermodynamic criteria and provides evidence for additional structural motifs within ITS1 and ITS2. The minimum structure model suggests that at least 20% of ITS1 and 38% of ITS2 nucleotide positions are involved in base pairing to form helices. The sequence alignment enabled by conserved structural features provides a framework for broadscale molecular evolutionary studies and the first family-level phylogeny of the Asteraceae based on nuclear DNA data. The phylogeny based on ITS sequence data is very well resolved and shows considerable congruence with relationships among major lineages of the family suggested by chloroplast DNA studies, including a monophyletic subfamily Asteroideae and a paraphyletic subfamily Cichorioideae. Combined analyses of ndhF and ITS sequences provide additional resolution and support for relationships in the family.     

Fifteen novel universal primer pairs for sequencing whole chloroplast genomes and a primer pair for nuclear ribosomal DNAs

Chloroplast genome information helps improve the phylogenetic resolution and can act as organelle-scale barcodes in recently radiated plant groups. Previously we reported that nine universal primer pairs could amplify angiosperm whole chloroplast genomes by long-range polymerase chain reaction and using next-generation sequencing. Although these primers show high universality and efficiency for sequencing whole chloroplast genomes in angiosperms, they did not fully resolve the following two issues surrounding sequencing angiosperm chloroplast genomes: (i) approximately 30% of angiosperms cannot be amplified successfully; and (ii) only fresh leaves can be applied. In this study, we designed another set of 15 universal primer pairs for amplifying angiosperm whole chloroplast genomes to complement the original nine primer pairs. Furthermore, we designed a primer pair for nuclear ribosomal DNAs (nrDNAs). To validate the functionality of the primers, we tested 44 species with silica gel-dried leaves and 15 species with fresh leaves that have been shown to not be amplified with the original nine primer pairs. The result showed that, in 65.9% and 88.6% of the 44 species with silica gel-dried leaves, the whole chloroplast genome and nrDNAs could be amplified, respectively. In addition, all 15 fresh leaf samples could have the whole chloroplast genome successfully amplified. The nrDNAs comprise partial sequences of 18S and 26S, along with the complete sequence of 5.8S and the internal transcribed spacers ITS1 and ITS2. The mean size of nrDNA was 5800 bp. This study shows that the 15 universal primer set is an indispensable tool for amplifying whole chloroplast genomes in angiosperms, and these are an important supplement to the nine reported primer pairs.     

The D4 set: primers that target highly variable intron loops in plant chloroplast genomes

Chloroplast group II introns offer high-quality, rapidly evolving single-copy loci for comparative sequence analysis. These introns feature diagnostic secondary structures with loops that are among the least evolutionarily constrained sequence in plastomes. We exploited these structures to develop universal primers that amplify and sequence the large Domain IV (D4) loop in several angiosperm introns. With a single sequence read, we recover 300-600 nucleotides of highly variable sequence across angiosperms, with rates of change that are equal to or higher than many of the best known intergenic spacers in plant chloroplast genomes.     

Lengths and nucleotide sequences of the internal spacers of nuclear ribosomal DNA in gymnosperms and pteridophytes

The nucleotide sequences of the internal transcribed spacers (ITS1 and ITS2) of the nuclear ribosomal DNA were analyzed in species belonging to gymnosperms and pteridophytes. The lengths of the ITSs of sixteen species of gymnosperms and seven species of pteridophytes were estimated. The gymnosperms have ITS1 regions larger than those observed in the pteridophytes and angiosperms (ca. 610–3100 bp versus 159–360 bp). On the other hand, the ITS2 regions appear to be of a conserved length (182–370 bp). We have determined the complete nucleotide sequences of ITS regions from four gymnosperm species and five pteridophyte species by cloning the PCR products. Sequence analysis showed the presence of three short tandem arranged subrepeats of about 70 bp in the 1112 bp ITS1 ofEphedra fragilis. Pyrimidine rich (about 90%) DNA segments of 40–50 bp were observed in the ITS1 ofGinkgo biloba. A highly conserved 16 bp long sequence known to be present in the ITS1 of the angiosperm species has been also found in the ITS1 ofCycas revoluta, Taxus baccata andEphedra fragilis. Dedicated to Prof.Emilio Battaglia.     

基于核rDNA的ITS序列在种子植物系统发育研究中的应用

种子植物核rDNA是高度重复的串联序列,由于同步进化的力量．大多数物种中这些重复单位间已发生纯合或接近纯合。5．8S rDNA把核rDNA的内转录间隔区分为ITS1和ITS2两部分．在被子植物中ITS1的长度为165～298bp,ITS2的长度为177～266bp,而在裸子植物中ITS片段较长。且其长度变化主要由ITS1的长度变异所致。可对这两个片段PCR产物进行直接测序或克隆测序。由于ITS序列变异较快．能够提供较丰富的变异位点和信息位点,已成为被子植物较低分类阶元的系统发育和分类研究中的重要分子标记,为探讨多倍体复合体网状进化关系,异源多倍体的起源提供了重要的系统学信息．但它一般不适合科以上水平的系统学研究。裸子植物中ITS片段较长,重复序列间的纯合程度不同,测序比较困难．因此对探讨裸子植物系统发育和分类受到了一定的限制,但近年来有所发展。     

Functional analysis of internal transcribed spacer 2 of Saccharomyces cerevisiae ribosomal DNA.

Using the previously described "tagged ribosome" (pORCS) system for in vivo mutational analysis of yeast rDNA, we show that small deletions in the 5'-terminal portion of ITS2 completely block maturation of 26 S rRNA at the level of the 29 SB precursor (5.8 S rRNA-ITS2-26 S rRNA). Various deletions in the 3'-terminal part, although severely reducing the efficiency of processing, still allow some mature 26 S rRNA to be formed. On the other hand, none of the ITS2 deletions affect the production of mature 17 S rRNA. Since all of the deletions severely disturb the recently proposed secondary structure of ITS2, these findings suggest an important role for higher order structure of ITS2 in processing. Analysis of the effect of complete or partial replacement of S. cerevisiae ITS2 with its counterpart sequences from Saccharomyces rosei or Hansenula wingei, points to helix V of the secondary structure model as an important element for correct and efficient processing. Direct mutational analysis shows that disruption of base-pairing in the middle of helix V does not detectably affect 26 S rRNA formation. In contrast, introduction of clustered point mutations at the apical end of helix V that both disrupt base-pairing and change the sequence of the loop, severely reduces processing. Since a mutant containing only point mutations in the sequence of the loop produces normal amounts of mature 26 S rRNA, we conclude that the precise (secondary and/or primary) structure at the lower end of helix V, but excluding the loop, is of crucial importance for efficient removal of ITS2.     

Deep-level diagnostic value of the rDNA-ITS region

The similarity of certain reported angiosperm rDNA internal transcribed spacer (ITS) region sequences to those of green algae prompted our analysis of the deep-level phylogenetic signal in the highly conserved but short 5.8S and hypervariable ITS2 sequences. We found that 5.8S sequences yield phylogenetic trees similar to but less well supported than those generated by a ca. 10-fold longer alignment from rDNA-18S sequences, as well as independent evidence. We attribute this result to our finding that, compared to 18S, the 5.8S has a higher proportion of sites subject to vary and greater among-site substitution rate homogeneity. We also determined that our phylogenetic results are not likely affected by intramolecular compensatory mutation to maintain RNA secondary structure nor by evident systematic biases in base composition. Despite historical homology, there appears to be no ITS2 primary sequence similarity shared sufficient similarity to cluster correctly on the basis of alignability. Our results indicate that groups, however, share sufficient similarity to cluster correctly on the basis of alignability. Our results indicate that ITS region sequences can diagnose organismal origins and phylogenetic relationships at many phylogenetic levels and provide a useful paradigm for molecular evolutionary study.      

Organization of the root apical meristem in angiosperms

Although flowers, leaves, and stems of the angiosperms have understandably received more attention than roots, the growing root tips, or root apical meristems (RAMs), are organs that could provide insight into angiosperm evolution. We studied RAM organization across a broad spectrum of angiosperms (45 orders and 132 families of basal angiosperms, monocots, and eudicots) to characterize angiosperm RAMs and cortex development related to RAMs. Types of RAM organization in root tips of flowering plants include open RAMs without boundaries between some tissues in the growing tip and closed RAMs with distinct boundaries between apical regions. Epidermis origin is associated with the cortex in some basal angiosperms and monocots and with the lateral rootcap in eudicots and other basal angiosperms. In most angiosperm RAMs, initials for the central region of the rootcap, or columella, are distinct from the lateral rootcap and its initials. Slightly more angiosperm families have exclusively closed RAMs than exclusively open RAMs, but many families have representatives with both open and closed RAMs. Root tips with open RAMs are generally found in angiosperm families considered sister to other families; certain open RAMs may be ancestral in angiosperms.     

Molecular Evolution and Phylogenetic Utility of the Internal Transcribed Spacer 2 (ITS2) in Calyptratae (Diptera: Brachycera)

The resolution potential of internal transcribed spacer 2 (ITS2) at deeper levels remains controversial. In this study, 105 ITS2 sequences of 55 species in Calyptratae were analyzed to examine the phylogenetic utility of the spacer above the subfamily level and to further understand its evolutionary characteristics. We predicted the secondary structure of each sequence using the minimum-energy algorithm and constructed two data matrixes for phylogenetic analysis. The ITS2 regions of Calyptratae display strong A-T bias and slight variation in length. The tandem and dispersed repeats embedded in the spacers possibly resulted from replication slippage or transposition. Most foldings conformed to the four-domain model. Sequence comparison in combination with the secondary structures revealed six conserved motifs. Covariation analysis from the conserved motifs indicated that the secondary structure restrains the sequence evolution of the spacer. The deep-level phylogeny derived from the ITS2 data largely agreed with the phylogenetic hypotheses from morphologic and other molecular evidence. Our analyses suggest that the accordant resolutions generated from different analyses can be used to infer deep-level phylogenetic relations.     

Evolution of seed storage protein genes: Legumin genes of Ginkgo biloba

Legumin-like seed storage proteins have been intensively studied in crop plants. However, little is known about the molecular evolution of these proteins and their genes and it was assumed that they originated from an ancestral gene that already existed at the beginning of angiosperm evolution. We have evidence for the ubiquitous occurrence of homologous proteins in gymnosperms as well. We have characterized the major seed storage globulin from Ginkgo biloba by amino acid sequencing, which reveals clear homology to legumin-like proteins from angiosperms. The Ginkgo legumin is encoded by a gene family; we describe two of its members. The promoter regions contain sequence motifs which are known to function as regulatory elements involved in seed-specific expression of angiosperm legumins, although the tissues concerned are different in gymnosperms and angiosperms. The Ginkgo legumin gene structure is divergent from that of angiosperms and suggests that the evolution of legumin genes implicated loss of introns. From our data and from functional approaches recently described it becomes obvious that the posttranslational processing site of legumin precursors is less conserved than hitherto assumed. Finally, we present a phylogenetic analysis of legumin encoding sequences and discuss their utility as molecular markers for the reconstruction of seed plant evolution.Correspondence to: K.-P. Häger     

Multigene analyses identify the three earliest lineages of extant flowering plants

Flowering plants (angiosperms) are by far the largest, most diverse, and most important group of land plants, with over 250,000 species and a dominating presence in most terrestrial ecosystems. Understanding the origin and early diversification of angiosperms has posed a long-standing botanical challenge [1]. Numerous morphological and molecular systematic studies have attempted to reconstruct the early history of this group, including identifying the root of the angiosperm tree. There is considerable disagreement among these studies, with various groups of putatively basal angiosperms from the subclass Magnoliidae having been placed at the root of the angiosperm tree (reviewed in [2-4]). We investigated the early evolution of angiosperms by conducting combined phylogenetic analyses of five genes that represent all three plant genomes from a broad sampling of angiosperms. Amborella, a monotypic, vessel-less dioecious shrub from New Caledonia, was clearly identified as the first branch of angiosperm evolution, followed by the Nymphaeales (water lillies), and then a clade of woody vines comprising Schisandraceae and Austrobaileyaceae. These findings are remarkably congruent with those from several concurrent molecular studies [5-7] and have important implications for whether or not the first angiosperms were woody and contained vessels, for interpreting the evolution of other key characteristics of basal angiosperms, and for understanding the timing and pattern of angiosperm origin and diversification.     

Structural evolution of nrDNA ITS in Pinaceae and its phylogenetic implications

Nuclear ribosomal DNA (nrDNA) has been considered as an important tool for inferring phylogenetic relationships at many taxonomic levels. In comparison with its fast concerted evolution in angiosperms, nrDNA is symbolized by slow concerted evolution and substantial ITS region length variation in gymnosperms, particularly in Pinaceae. Here we studied structure characteristics, including subrepeat composition, size, GC content and secondary structure, of nrDNA ITS regions of all Pinaceae genera. The results showed that the ITS regions of all taxa studied contained subrepeat units, ranging from 2 to 9 in number, and these units could be divided into two types, longer subrepeat (LSR) without the motif (5'-GGCCACCCTAGTC) and shorter subrepeat (SSR) with the motif. Phylogenetic analyses indicate that the homology of some SSRs still can be recognized, providing important informations for the evolutionary history of nrDNA ITS and phylogeny of Pinaceae. In particular, the adjacent tandem SSRs are not more closely related to one another than they are to remote SSRs in some genera, which may imply that multiple structure variations such as recombination have occurred in the ITS1 region of these groups. This study also found that GC content in the ITS1 region is relevant to its sequence length and subrepeat number, and could provide some phylogenetic information, especially supporting the close relationships among Picea, Pinus, and Cathaya. Moreover, several characteristics of the secondary structure of Pinaceae ITS1 were found as follows: (1) the structure is dominated by several extended hairpins; (2) the configuration complexity is positively correlated with subrepeat number; (3) paired subrepeats often partially overlap at the conserved motif (5'-GGCCACCCTAGTC), and form a long stem, while other subrepeats fold onto itself, leaving part of the conserved motif exposed in hairpin loops.     

The age and diversification of the angiosperms re-revisited

? Premise of the study: It has been 8 years since the last comprehensive analysis of divergence times across the angiosperms. Given recent methodological improvements in estimating divergence times, refined understanding of relationships among major angiosperm lineages, and the immense interest in using large angiosperm phylogenies to investigate questions in ecology and comparative biology, new estimates of the ages of the major clades are badly needed. Improved estimations of divergence times will concomitantly improve our understanding of both the evolutionary history of the angiosperms and the patterns and processes that have led to this highly diverse clade. ? Methods: We simultaneously estimated the age of the angiosperms and the divergence times of key angiosperm lineages, using 36 calibration points for 567 taxa and a "relaxed clock" methodology that does not assume any correlation between rates, thus allowing for lineage-specific rate heterogeneity. ? Key results: Based on the analysis for which we set fossils to fit lognormal priors, we obtained an estimated age of the angiosperms of 167-199 Ma and the following age estimates for major angiosperm clades: Mesangiospermae (139-156 Ma); Gunneridae (109-139 Ma); Rosidae (108-121 Ma); Asteridae (101-119 Ma). ? Conclusions: With the exception of the age of the angiosperms themselves, these age estimates are generally younger than other recent molecular estimates and very close to dates inferred from the fossil record. We also provide dates for all major angiosperm clades (including 45 orders and 335 families [208 stem group age only, 127 both stem and crown group ages], sensu APG III). Our analyses provide a new comprehensive source of reference dates for major angiosperm clades that we hope will be of broad utility.     

Secondary structure and phylogenetic analysis of the internal transcribed spacers 1 and 2 of bush crickets (Orthoptera: Tettigoniidae: Barbitistini)

We inferred secondary structure models of the internal transcribed spacers (ITS) 1 and 2 of bush crickets using a combined comparative and thermodynamic approach. The inferred secondary structure models were used to account for interdependency of interacting nucleotides in a phylogenetic analysis of the bush cricket genus Poecilimon. Our analysis indicates that the two previously reported conformational structures (i.e., hairpin and ring) of ITS2 are likely to fold in bush crickets as well and that both predicted structures are similar to those proposed for other eukaryotes. Comparing predicted ITS1 secondary structure models proved to be difficult because of substantial variation in their nucleotide sequence length. Our study revealed that the phylogenetic signal of ITS1 and ITS2 is largely congruent with that preserved in the mitochondrial genes 16S rRNA, tRNA‐Val and 12S rRNA. The phylogenetic signal in both the nuclear and the mitochondrial genome question the monophyly of the genus Poecilimon: species of the genera Poecilimonella, Parapoecilimon, Polysarcus and Phonochorion consistently cluster within Poecilimon.     

Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.     

The four-celled female gametophyte of Illicium (Illiciaceae; Austrobaileyales): implications for understanding the origin and early evolution of monocots, eumagnoliids,and eudicots

The recent consensus that Amborellaceae, Nymphaeales, and Austrobaileyales form the three earliest-diverging lineages of angiosperms has led comparative biologists to reconsider the origin and early developmental evolution of the angiosperm seven-celled/eight-nucleate (Polygonum-type) female gametophyte. Illicium mexicanum (Illiciaceae; Austrobaileyales) develops a four-celled/four-nucleate female gametophyte. The ontogenetic sequence of the Illicium female gametophyte is consistent with that of all other Austrobaileyales and also with all Nymphaeales and is likely a plesiomorphy of angiosperms. A character analysis based on more than 250 embryological studies indicates that a transition from an ancestrally four-celled/four-nucleate Illicium-like female gametophyte to a seven-celled/eight-nucleate female gametophyte occurred in the common ancestor of the sister group to Austrobaileyales (a clade that includes monocots, eumagnoliids, and eudicots). Comparative analysis of reconstructed ancestral female gametophyte ontogenies identifies specific early stages of ontogeny that were modified during this transition. These modifications generated two important angiosperm novelties-a set of three persistent antipodal cells and a binucleate central cell, which upon fertilization yields a triploid endosperm. Early angiosperms are anatomically quite diverse in these two features, although triploid endosperm, composed of one paternal genome and two maternal genomes, is a conserved feature of the overwhelming majority of angiosperms.     

Characterization of matteuccin, the 2.2S storage protein of the ostrich fern. Evolutionary relationship to angiosperm seed storage proteins

The 2.2S spore storage protein (matteuccin) of the ostrich fern, Matteuccia struthiopteris, has been isolated and characterized. It is a small basic protein consisting of two disulfide-linked polypeptides with approximate molecular masses of 3.0 kDa and 8.0 kDa. At least four different isoforms exist where two of the forms differ from the other by having a slightly smaller heavy chain. Amino acid analysis reveals that the 2.2S protein is rich in arginine. Almost complete amino acid sequence information was obtained for the light chain and a partial sequence for the heavy chain. Amino acid sequence comparison reveals that this protein shows a high similarity to seed storage proteins in different angiosperm species in spite of the fact that the common ancestor of ferns and angiosperms lived more than 300 million years ago.     

The mitochondrial genome of the moss Physcomitrella patens sheds new light on mitochondrial evolution in land plants

The phylogenetic positions of bryophytes and charophytes, together with their genome features, are important for understanding early land plant evolution. Here we report the complete nucleotide sequence (105,340 bp) of the circular-mapping mitochondrial DNA of the moss Physcomitrella patens. Available evidence suggests that the multipartite structure of the mitochondrial genome in flowering plants does not occur in Physcomitrella. It contains genes for 3 rRNAs (rnl, rns, and rrn5), 24 tRNAs, and 42 conserved mitochondrial proteins (14 ribosomal proteins, 4 ccm proteins, 9 nicotinamide adenine dinucleotide dehydrogenase subunits, 5 ATPase subunits, 2 succinate dehydrogenase subunits, apocytochrome b, 3 cytochrome oxidase subunits, and 4 other proteins). We estimate that 5 tRNA genes are missing that might be encoded by the nuclear genome. The overall mitochondrial genome structure is similar in Physcomitrella, Chara vulgaris, Chaetosphaeridium globosum, and Marchantia polymorpha, with easily identifiable inversions and translocations. Significant synteny with angiosperm and chlorophyte mitochondrial genomes was not detected. Phylogenetic analysis of 18 conserved proteins suggests that the moss-liverwort clade is sister to angiosperms, which is consistent with a previous analysis of chloroplast genes but is not consistent with some analyses using mitochondrial sequences. In Physcomitrella, 27 introns are present within 16 genes. Nine of its intron positions are shared with angiosperms and 4 with Marchantia, which in turn shares only one intron position with angiosperms. The phylogenetic analysis as well as the syntenic structure suggest that the mitochondrial genomes of Physcomitrella and Marchantia retain prototype features among land plant mitochondrial genomes.