Implications of the Plastid Genome Sequence of Typha (Typhaceae, Poales) for Understanding Genome Evolution in Poaceae

Plastid genomes of the grasses (Poaceae) are unusual in their organization and rates of sequence evolution. There has been a recent surge in the availability of grass plastid genome sequences, but a comprehensive comparative analysis of genome evolution has not been performed that includes any related families in the Poales. We report on the plastid genome of Typha latifolia, the first non-grass Poales sequenced to date, and we present comparisons of genome organization and sequence evolution within Poales. Our results confirm that grass plastid genomes exhibit acceleration in both genomic rearrangements and nucleotide substitutions. Poaceae have multiple structural rearrangements, including three inversions, three genes losses (accD, ycf1, ycf2), intron losses in two genes (clpP, rpoC1), and expansion of the inverted repeat (IR) into both large and small single-copy regions. These rearrangements are restricted to the Poaceae, and IR expansion into the small single-copy region correlates with the phylogeny of the family. Comparisons of 73 protein-coding genes for 47 angiosperms including nine Poaceae genera confirm that the branch leading to Poaceae has significantly accelerated rates of change relative to other monocots and angiosperms. Furthermore, rates of sequence evolution within grasses are lower, indicating a deceleration during diversification of the family. Overall there is a strong correlation between accelerated rates of genomic rearrangements and nucleotide substitutions in Poaceae, a phenomenon that has been noted recently throughout angiosperms. The cause of the correlation is unknown, but faulty DNA repair has been suggested in other systems including bacterial and animal mitochondrial genomes.     

Phylogenetic structure in the grass family (Poaceae): evidence from the nuclear gene phytochrome B

Phylogenetic analyses of partial phytochrome B (PHYB) nuclear DNA sequences provide unambiguous resolution of evolutionary relationships within Poaceae. Analysis of PHYB nucleotides from 51 taxa representing seven traditionally recognized subfamilies clearly distinguishes three early-diverging herbaceous "bambusoid" lineages. First and most basal are Anomochloa and Streptochaeta, second is Pharus, and third is Puelia. The remaining grasses occur in two principal, highly supported clades. The first comprises bambusoid, oryzoid, and pooid genera (the BOP clade); the second comprises panicoid, arundinoid, chloridoid, and centothecoid genera (the PACC clade). The PHYB phylogeny is the first nuclear gene tree to address comprehensively phylogenetic relationships among grasses. It corroborates several inferences made from chloroplast gene trees, including the PACC clade, and the basal position of the herbaceous bamboos Anomochloa, Streptochaeta, and Pharus. However, the clear resolution of the sister group relationship among bambusoids, oryzoids, and pooids in the PHYB tree is novel; the relationship is only weakly supported in ndhF trees and is nonexistent in rbcL and plastid restriction site trees. Nuclear PHYB data support Anomochlooideae, Pharoideae, Pooideae sensu lato, Oryzoideae, Panicoideae, and Chloridoideae, and concur in the polyphyly of both Arundinoideae and Bambusoideae.     

Distinctive Architecture of the Chloroplast Genome in the Chlorodendrophycean Green Algae Scherffelia dubia and Tetraselmis sp. CCMP 881

The Chlorodendrophyceae is a small class of green algae belonging to the core Chlorophyta, an assemblage that also comprises the Pedinophyceae, Trebouxiophyceae, Ulvophyceae and Chlorophyceae. Here we describe for the first time the chloroplast genomes of chlorodendrophycean algae (Scherffelia dubia, 137,161 bp; Tetraselmis sp. CCMP 881, 100,264 bp). Characterized by a very small single-copy (SSC) region devoid of any gene and an unusually large inverted repeat (IR), the quadripartite structures of the Scherffelia and Tetraselmis genomes are unique among all core chlorophytes examined thus far. The lack of genes in the SSC region is offset by the rich and atypical gene complement of the IR, which includes genes from the SSC and large single-copy regions of prasinophyte and streptophyte chloroplast genomes having retained an ancestral quadripartite structure. Remarkably, seven of the atypical IR-encoded genes have also been observed in the IRs of pedinophycean and trebouxiophycean chloroplast genomes, suggesting that they were already present in the IR of the common ancestor of all core chlorophytes. Considering that the relationships among the main lineages of the core Chlorophyta are still unresolved, we evaluated the impact of including the Chlorodendrophyceae in chloroplast phylogenomic analyses. The trees we inferred using data sets of 79 and 108 genes from 71 chlorophytes indicate that the Chlorodendrophyceae is a deep-diverging lineage of the core Chlorophyta, although the placement of this class relative to the Pedinophyceae remains ambiguous. Interestingly, some of our phylogenomic trees together with our comparative analysis of gene order data support the monophyly of the Trebouxiophyceae, thus offering further evidence that the previously observed affiliation between the Chlorellales and Pedinophyceae is the result of systematic errors in phylogenetic reconstruction.     

Small single-copy region of plastid DNA in the non-photosynthetic angiosperm Epifagus virginiana contains only two genes. Differences among dicots, monocots and bryophytes in gene organization at a non-bioenergetic locus.

We have determined the nucleotide sequence of a 7 kb (1 kb = 10(3) base-pairs) region that includes the entire small single-copy region (SSC) of the plastid genome of Epifagus virginiana, a non-photosynthetic, parasitic flowering plant. The SSC (4.8 kb) is considerably smaller than those of photosynthetic plants due to the complete deletion of all photosynthetic, chlororespiratory and ribosomal protein genes. This leaves only two genes: a protein gene of 1738 codons whose product is unlikely to be involved in bioenergetic processes and a leucine tRNA gene (trn(LUAG)). Both genes span junctions between the inverted repeat and the SSC, with the consequence that the terminal 20 base-pairs of the repeat is transcribed in both directions and functions both as the 3' end of the tRNA gene and as an internal segment of orf1738. We find that the region of tobacco plastid DNA homologous to Epifagus orf1738 contains a single open reading frame (ORF) of 1901 codons rather than the three ORFs of 1244, 273 and 228 codons originally reported. However, we confirm that the equivalent region of the bryophyte Marchantia contains two genes (1068 and 464 codons) corresponding to the N and C-terminal portions of the dicot protein. In contrast, rice plastid DNA contains a severely truncated pseudogene at this locus.     

The first complete plastid genome of Burmannia disticha L. from the mycoheterotrophic monocot family Burmanniaceae

Burmanniaceae is one major group within the monocot order Dioscoreales that has not had its plastome sequenced. Members of Burmanniaceae are mostly achlorophyllous, although the genus Burmannia also includes autotrophs. Here, we report sequencing and analysis of the first Burmanniaceae plastid genome from Burmannia disticha L.. This plastome is 157,480 bp and was assembled as a circular sequence with the typical quadripartite structure of plant plastid genomes. This plastome has a regular number of potentially functional genes with a total of 111, including 78 protein coding genes, 4 ribosomal RNA (rRNA) genes, and 29 tRNA genes. The ratio of the total length of genic:intergenic DNA is 1.58:1, and the mean length of intergenic regions is 398 bp, the longest being 1918 bp. The overall GC content of the B. disticha plastome is 34.90%, and the IR regions in B. disticha are more GC rich (39.50%) than the LSC (32.30%) and SSC (28.80%) regions. Phylogenetic analysis of protein-coding sequences from plastomes of related species in the order Dioscoreales support a clade comprising Burmanniaceae and Dioscoreaceae. This phylogenetic placement is congruent with previous findings based on nuclear and mitochondrial evidence.     

Complete Plastid Genome Sequence of the Brown Alga Undaria pinnatifida

In this study, we fully sequenced the circular plastid genome of a brown alga, Undaria pinnatifida. The genome is 130,383 base pairs (bp) in size; it contains a large single-copy (LSC, 76,598 bp) and a small single-copy region (SSC, 42,977 bp), separated by two inverted repeats (IRa and IRb: 5,404 bp). The genome contains 139 protein-coding, 28 tRNA, and 6 rRNA genes; none of these genes contains introns. Organization and gene contents of the U. pinnatifida plastid genome were similar to those of Saccharina japonica. There is a co-linear relationship between the plastid genome of U. pinnatifida and that of three previously sequenced large brown algal species. Phylogenetic analyses of 43 taxa based on 23 plastid protein-coding genes grouped all plastids into a red or green lineage. In the large brown algae branch, U. pinnatifida and S. japonica formed a sister clade with much closer relationship to Ectocarpus siliculosus than to Fucus vesiculosus. For the first time, the start codon ATT was identified in the plastid genome of large brown algae, in the atpA gene of U. pinnatifida. In addition, we found a gene-length change induced by a 3-bp repetitive DNA in ycf35 and ilvB genes of the U. pinnatifida plastid genome.     

虎杖叶绿体基因组结构与变异分析

虎杖(Reynoutria japonica Houtt.)为蓼科(Polygonaceae)蓼族(Polygoneae)虎杖属(Reynoutria Houtt.)植物,是一种传统的中草药,具有利湿退黄、清热解毒、散瘀止痛、止咳化痰的功效。本研究采取高通量测序技术获得5个虎杖品种的叶绿体全基因组序列,并与NCBI已公布的蓼族何首乌(Fallopia multiflora)和金线草(Antenoron filiforme)等植物的叶绿体全基因组序列进行了基因组学和系统发育分析。通过基因组学分析发现,5种虎杖的叶绿体基因组大小有163 376 bp和163 371 bp两种情况,并呈现出典型的环状四分体结构,85 784 bp的一条较长的单拷贝区(large single-copy region,LSC),18 616 bp的一条较短的单拷贝区(small single-copy region,SSC),还有两条长度一致的反向重复区,分别为IRa区和IRb区相间隔分布。通过注释得到161个基因,其中蛋白编码基因106个,rRNA编码基因10个,tRNA编码基因45个。总GC含量为36.7...     

Plastid genome structure and phylogenomics of Nymphaeales: conserved gene order and new insights into relationships

The plastid genomes of early-diverging angiosperms were among the first land plant plastomes investigated. Despite their importance to understanding angiosperm evolution, no investigation has so far compared gene content or gene synteny of these plastid genomes with a focus on the Nymphaeales. Here, we report an evaluation and comparison of gene content, gene synteny and inverted repeat length for a set of 15 plastid genomes of early-diverging angiosperms. Seven plastid genomes of the Nymphaeales were newly sequenced for this investigation. We compare gene order and inverted repeat (IR) length across all genomes, review the gene annotations of previously published genomes, generate a multi-gene alignment of 77 plastid-encoded genes and reconstruct the phylogenetic relationships of the taxa under study. Our results show that gene content and synteny are highly conserved across early-diverging angiosperms: All species analyzed display complete gene synteny when accounting for expansions and contractions of the IRs. This conservation was initially obscured by ambiguous and potentially incorrect gene annotations in previously published genomes. We also report the presence of intact open reading frames across all taxa analyzed. The multi-gene phylogeny displays maximum support for the families Cabombaceae and Hydatellaceae, but no support for a clade of all Nymphaeaceae. It further indicates that the genus Victoria is embedded within Nymphaea. Plastid genomes of Trithuria were found to deviate by numerous substitutions and length changes in the IRs. Phylogenetic analyses further indicate that a previously published plastome named Nymphaea mexicana falls into a clade of N. odorata and should be re-evaluated.     

Complete nucleotide sequence of the chloroplast genome from a leptosporangiate fern, Adiantum capillus-veneris L.

We determined the complete nucleotide sequence of the chloroplast genome of the leptosporangiate fern, Adiantum capillus-veneris L. (Pteridaceae). The circular genome is 150,568 bp, with a large single-copy region (LSC) of 82,282 bp, a small-single copy region (SSC) of 21,392 bp and inverted repeats (IR) of 23,447 bp each. We compared the sequence to other published chloroplast genomes to infer the location of putative genes. When the IR is considered only once, we assigned 118 genes, of which 85 encode proteins, 29 encode tRNAs and 4 encode rRNAs. Four protein-coding genes, all four rRNA genes and six tRNA genes occur in the IR. Most (57) putative protein-coding genes appear to start with an ATG codon, but we also detected five other possible start codons, some of which suggest tRNA editing. We also found 26 apparent stop codons in 18 putative genes, also suggestive of RNA editing. We found all but one of the tRNA genes necessary to encode the complete repertoire required for translation. The missing trnK gene appears to have been disrupted by a large inversion, relative to other published chloroplast genomes. We detected several structural rearrangements that may provide useful information for phylogenetic studies.     

Complete plastid genome sequence of Vaccinium macrocarpon: structure, gene content, and rearrangements revealed by next generation sequencing

The complete plastid genome sequence of the American cranberry (Vaccinium macrocarpon Ait.) was reconstructed using next-generation sequencing data by in silico procedures. We used Roche 454 shotgun sequence data to isolate cranberry plastid-specific sequences of “HyRed” via homology comparisons with complete sequences from several species available at the National Center for Biotechnology Information database. Eleven cranberry plastid contigs were selected for the construction of the plastid genome-based homologies and on raw reads flowing through contigs and connection information. We assembled and annotated a cranberry plastid genome (82,284 reads; 185x coverage) with a length of 176 kb and the typical structure found in plants, but with several structural rearrangements in the large single-copy region when compared to other plastid asterid genomes. To evaluate the reliability of the sequence data, phylogenetic analysis of 30 species outside the order Ericales (with 54 genes) showed Vaccinium inside the clade Asteridae, as reported in other studies using single genes. The cranberry plastid genome sequence will allow the accumulation of critical data useful for breeding and a suite of other genetic studies.     

Complete chloroplast genome sequences of <Emphasis Type="Italic">Hordeum vulgare</Emphasis>, <Emphasis Type="Italic">Sorghum bicolor</Emphasis> and <Emphasis Type="Italic">Agrostis stolonifera</Emphasis>, and comparative analyses with other grass genomes

Comparisons of complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera to six published grass chloroplast genomes reveal that gene content and order are similar but two microstructural changes have occurred. First, the expansion of the IR at the SSC/IRa boundary that duplicates a portion of the 5' end of ndhH is restricted to the three genera of the subfamily Pooideae (Agrostis, Hordeum and Triticum). Second, a 6 bp deletion in ndhK is shared by Agrostis, Hordeum, Oryza and Triticum, and this event supports the sister relationship between the subfamilies Erhartoideae and Pooideae. Repeat analysis identified 19-37 direct and inverted repeats 30 bp or longer with a sequence identity of at least 90%. Seventeen of the 26 shared repeats are found in all the grass chloroplast genomes examined and are located in the same genes or intergenic spacer (IGS) regions. Examination of simple sequence repeats (SSRs) identified 16-21 potential polymorphic SSRs. Five IGS regions have 100% sequence identity among Zea mays, Saccharum officinarum and Sorghum bicolor, whereas no spacer regions were identical among Oryza sativa, Triticum aestivum, H. vulgare and A. stolonifera despite their close phylogenetic relationship. Alignment of EST sequences and DNA coding sequences identified six C-U conversions in both Sorghum bicolor and H. vulgare but only one in A. stolonifera. Phylogenetic trees based on DNA sequences of 61 protein-coding genes of 38 taxa using both maximum parsimony and likelihood methods provide moderate support for a sister relationship between the subfamilies Erhartoideae and Pooideae.     

Characterization of the gene for prostate-specific antigen, a human glandular kallikrein

The gene for the human glandular kallikrein, prostate-specific antigen, has been cloned. The sequence of 7130 nucleotides encompassing the gene and 633 bp of 5' and 639 bp of 3' flanking DNA has been determined. The translation initiation site was slightly heterogeneous, yielding 5' non-translated leader sequences of 41 and 35 bp. The gene is divided into five exons, with introns located at positions identical with those found in other glandular kallikrein genes. The nucleotide sequence is very similar to that of the human kallikrein gene hGK-1, with 76 to 93% of the nucleotides being identical in the exons and 76 to 87% in the introns. The similarity also extends approximately 200 bp into the sequence flanking the 5' end of hGK-1 and several other, both human and rodent, glandular kallikrein genes.     

Nucleotide substitution rates in legume chloroplast DNA depend on the presence of the inverted repeat

The chloroplast genomes of some species of legumes lack the large inverted repeat (IR) that is a trademark of most land-plant chloroplasts. Our analysis of chloroplast genes in legume species that have an IR shows that the synonymous (silent) substitution rate in IR genes is 2.3-fold lower than in single-copy (SC) genes, which is largely in agreement with earlier findings. Given that all genes in species that lack the IR are single-copy, what level of synonymous substitution exists in these genes? We report a uniform substitution rate in IR-less genomes, and moreover, we find this rate to be at the level otherwise reserved for SC genes. In other words, the synonymous substitution rate has accelerated in the remaining copy of the duplicate region. We propose that this acceleration is a direct result of the decrease in the copy number of the sequence, rather than an intrinsic property of the genes normally located in the IR.     

Ubiquity of the St chloroplast genome in St-containing Triticeae polyploids.

Interspecific hybridization occurs between Tritceae species in the grass family (Poaceae) giving rise to allopolyploid species. To examine bias in cytoplasmic DNA inheritance in these hybridizations, the sequence of the 3' end of the chloroplast ndhF gene was compared among 29 allopolyploid Triticeae species containing the St nuclear genome in combination with the H, I, Ns, P, W, Y, and Xm nuclear genomes. These ndhF sequences were also compared with those from diploid or allotetraploid Triticeae species having the H, I, Ns, P, W, St, and Xm genomes. The cpDNA sequences were highly similar among diploid, allotetraploid, allohexaploid, and allooctoploid Triticeae accessions containing the St nuclear genome, with 0-6-nucleotide (nt) substitutions (0-0.8%) occurring between pairs of species. Neighbor-joining analysis of the sequences showed that the ndhF DNA sequences from species containing the St nuclear genome formed a strongly supported clade. The data indicated a strong preference for cpDNA inheritance from the St nuclear genome-containing parent in hybridizations between Triticeae species. This preference was independent of the presence of the H, I, Ns, P, W, and Xm nuclear genomes, the geographic distribution of the species, and the mode of reproduction. The data suggests that hybridizations having the St-containing parent as the female may be more successful.     

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

The Pharus genome fills an important genomic gap, providing numerous insights into how whole-genome duplication contributed to the origin and diversification of the grass family.