Construction of a bacterial artificial chromosome library from the spikemoss Selaginella moellendorffii: a new resource for plant comparative genomics

Huperziacrassifolia (Lycopodiaceae), a new species from China based on morphological characters and molecular evidence

A new species of the firmoss from China, Huperziacrassifolia sp. nov., is described and illustrated based on morphological characters and molecular evidence. The new species resembles species associated with the H.javanica complex, in particular H.javanica based on leaf shape and serrations, but can be easily distinguished by elliptic lanceolate and thick coriaceous leaves, well differentiated seasonal constriction zones, and reflexed leaf margins when get dried. Phylogenomic reconstruction using whole chloroplast genome sequences recovered H.crassifolia as sister to H.sutchueniana and only distantly related to morphological similar species H.javanica, H.nanlingensis, and H.serrata. The genome size 2C = 17.2 pg indicated the new species to be a tetraploid, whereas diploid H.javanica had a genome size of 8.7 pg. Morphological characters, distribution, and conservation status of the new species are also presented.

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana-chosen as the model plant system for this reason-as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

Ammonium metabolism in Selaginella bryopteris in response to dehydration-rehydration and characterisation of desiccation tolerant, thermostable, cytosolic glutamine synthetase from plant

Water deficit (WD) has adverse effects on plant growth, and acclimation requires responses allowing primary metabolism to continue. Resurrection plants can serve as model system to gain insight into metabolic regulation during WD. We herein report the response of a resurrection lycophyte, Selaginella bryopteris, to dehydration-rehydration cycle with emphasis on ammonium metabolism. Dehydration of S. bryopteris fronds resulted in decrease of total protein and increase of free ammonium levels and the effect was reversed on rehydration. The proline content increased twice after 24 h of dehydration, which again recovered to background levels comparable to that at full turgor state. The specific activity of glutamine synthetase (GS) didn't change significantly till 6 h and then declined by 21% after 24 h of dehydration, whereas specific activities of glutamate synthase (GOGAT) and aminating glutamate dehydrogenase (GDH) were enhanced significantly during dehydration. The deaminating activity of GDH also increased during dehydration albeit at a slower rate. Immunoblot analysis indicated overexpression of GS and GDH polypeptides during dehydration and their levels declined on rehydration. The results suggested significant role of GDH along with GS/GOGAT in production of nitrogen-rich amino acids for desiccation tolerance. Unlike higher plants S. bryopteris expressed GS only in cytosol. The enzyme had pH and temperature optima of 5.5 and 60°C, respectively, and it retained 96% activity on preincubation at 60°C for 30 min indicating thermostability. Hence, like higher plants the cytosolic GS from S. bryopteris has a conserved role in stress tolerance.

Engineering the Erythromycin-Producing Strain Saccharopolyspora erythraea HOE107 for the Heterologous Production of Polyketide Antibiotics

Bacteria of the genus Saccharopolyspora produce important polyketide antibiotics, including erythromycin A (Sac. erythraea) and spinosad (Sac. spinosa). We herein report the development of an industrial erythromycin-producing strain, Sac. erythraea HOE107, into a host for the heterologous expression of polyketide biosynthetic gene clusters (BGCs) from other Saccharopolyspora species and related actinomycetes. To facilitate the integration of natural product BGCs and auxiliary genes beneficial for the production of natural products, the erythromycin polyketide synthase (ery) genes were replaced with two bacterial attB genomic integration sites associated with bacteriophages ϕC31 and ϕBT1. We also established a highly efficient conjugation protocol for the introduction of large bacterial artificial chromosome (BAC) clones into Sac. erythraea strains. Based on this optimized protocol, an arrayed BAC library was effectively transferred into Sac. erythraea. The large spinosad gene cluster from Sac. spinosa and the actinorhodin gene cluster from Streptomyces coelicolor were successfully expressed in the ery deletion mutant. Deletion of the endogenous giant polyketide synthase genes pkeA1-pkeA4, the product of which is not known, and the flaviolin gene cluster (rpp) from the bacterium increased the heterologous production of spinosad and actinorhodin. Furthermore, integration of pJTU6728 carrying additional beneficial genes dramatically improved the yield of actinorhodin in the engineered Sac. erythraea strains. Our study demonstrated that the engineered Sac. erythraea strains SLQ185, LJ161, and LJ162 are good hosts for the expression of heterologous antibiotics and should aid in expression-based genome-mining approaches for the discovery of new and cryptic antibiotics from Streptomyces and rare actinomycetes.

Cloning and Functional Characterization of Two 4-Coumarate: CoA Ligase Genes from Selaginella moellendorffii

Selaginella is an extant lycopodiophyte genus, which is representative of an ancient lineage of tracheophytes. The important evolutionary status makes it a valuable resource for the study of metabolic evolution in vascular plants. 4-coumarate: CoA ligase (4CL) is the pivotal enzyme that controls the flow of carbon through the phenylpropanoid metabolic pathway into the specific lignin, flavonoid, and wall-bound phenolics biosynthesis pathways. Although 4CLs have been extensively characterized in other vascular plants, little is known of their functions in Selaginella. Here, we isolated two 4CL genes (Sm4CL1 and Sm4CL2) from Selaginella moellendorffii. Based on the enzymatic activities of the recombinant proteins, both of these genes encoded bona fide 4CLs. The 4CL isoforms in S. moellendorffii have different activities: Sm4CL2 was more active than Sm4CL1. The enzymatic properties and gene expression patterns indicated that the 4CL genes have been conserved in the evolution of vascular plants.

Phosphoproteomic analysis of the non-seed vascular plant model Selaginella moellendorffii

BACKGROUND

Selaginella (Selaginella moellendorffii) is a lycophyte which diverged from other vascular plants approximately 410 million years ago. As the first reported non-seed vascular plant genome, Selaginella genome data allow comparative analysis of genetic changes that may be associated with land plant evolution. Proteomics investigations on this lycophyte model have not been extensively reported. Phosphorylation represents the most common post-translational modifications and it is a ubiquitous regulatory mechanism controlling the functional expression of proteins inside living organisms.

RESULTS

In this study, polyethylene glycol fractionation and immobilized metal ion affinity chromatography were employed to isolate phosphopeptides from wild-growing Selaginella. Using liquid chromatography-tandem mass spectrometry analysis, 1593 unique phosphopeptides spanning 1104 non-redundant phosphosites with confirmed localization on 716 phosphoproteins were identified. Analysis of the Selaginella dataset revealed features that are consistent with other plant phosphoproteomes, such as the relative proportions of phosphorylated Ser, Thr, and Tyr residues, the highest occurrence of phosphosites in the C-terminal regions of proteins, and the localization of phosphorylation events outside protein domains. In addition, a total of 97 highly conserved phosphosites in evolutionary conserved proteins were identified, indicating the conservation of phosphorylation-dependent regulatory mechanisms in phylogenetically distinct plant species. On the other hand, close examination of proteins involved in photosynthesis revealed phosphorylation events which may be unique to Selaginella evolution. Furthermore, phosphorylation motif analyses identified Pro-directed, acidic, and basic signatures which are recognized by typical protein kinases in plants. A group of Selaginella-specific phosphoproteins were found to be enriched in the Pro-directed motif class.

CONCLUSIONS

Our work provides the first large-scale atlas of phosphoproteins in Selaginella which occupies a unique position in the evolution of terrestrial plants. Future research into the functional roles of Selaginella-specific phosphorylation events in photosynthesis and other processes may offer insight into the molecular mechanisms leading to the distinct evolution of lycophytes.

DNA content variation and its significance in the evolution of the genus Micrasterias (Desmidiales, Streptophyta)

It is now clear that whole genome duplications have occurred in all eukaryotic evolutionary lineages, and that the vast majority of flowering plants have experienced polyploidisation in their evolutionary history. However, study of genome size variation in microalgae lags behind that of higher plants and seaweeds. In this study, we have addressed the question whether microalgal phylogeny is associated with DNA content variation in order to evaluate the evolutionary significance of polyploidy in the model genus Micrasterias. We applied flow-cytometric techniques of DNA quantification to microalgae and mapped the estimated DNA content along the phylogenetic tree. Correlations between DNA content and cell morphometric parameters were also tested using geometric morphometrics. In total, DNA content was successfully determined for 34 strains of the genus Micrasterias. The estimated absolute 2C nuclear DNA amount ranged from 2.1 to 64.7 pg; intraspecific variation being 17.4–30.7 pg in M. truncata and 32.0–64.7 pg in M. rotata. There were significant differences between DNA contents of related species. We found strong correlation between the absolute nuclear DNA content and chromosome numbers and significant positive correlation between the DNA content and both cell size and number of terminal lobes. Moreover, the results showed the importance of cell/life cycle studies for interpretation of DNA content measurements in microalgae.

Construction and characterization of a BAC library from the Coffea arabica genotype Timor Hybrid CIFC 832/2

Most of the world's coffee production originates from Coffea arabica, an allotetraploid species with low genetic diversity and for which few genomic resources are available. Genomic libraries with large DNA fragment inserts are useful tools for the study of plant genomes, including the production of physical maps, integration studies of physical and genetic maps, genome structure analysis and gene isolation by positional cloning. Here, we report the construction and characterization of a Bacterial Artificial Chromosome (BAC) library from C. arabica Timor Hybrid CIFC 832/2, a parental genotype for several modern coffee cultivars. The BAC library consists of 56,832 clones with an average insert size of 118 kb, which represents a dihaploid genome coverage of five to sixfold. The content of organellar DNA was estimated at 1.04 and 0.5 % for chloroplast and mitochondrial DNA, respectively. The BAC library was screened for the NADPH-dependent mannose-6-phosphate reductase gene (CaM6PR) with markers positioned on four linkage groups of a partial C. arabica genetic map. A mixed approach using PCR and membrane hybridization of BAC pools allowed for the discovery of nine BAC clones with the CaM6PR gene and 53 BAC clones that were anchored to the genetic map with simple sequence repeat markers. This library will be a useful tool for future studies on comparative genomics and the identification of genes and regulatory elements controlling major traits in this economically important crop species.

A role for more axillary growth1 (MAX1) in evolutionary diversity in strigolactone signaling upstream of MAX2

Strigolactones (SLs) are carotenoid-derived phytohormones with diverse roles. They are secreted from roots as attractants for arbuscular mycorrhizal fungi and have a wide range of endogenous functions, such as regulation of root and shoot system architecture. To date, six genes associated with SL synthesis and signaling have been molecularly identified using the shoot-branching mutants more axillary growth (max) of Arabidopsis (Arabidopsis thaliana) and dwarf (d) of rice (Oryza sativa). Here, we present a phylogenetic analysis of the MAX/D genes to clarify the relationships of each gene with its wider family and to allow the correlation of events in the evolution of the genes with the evolution of SL function. Our analysis suggests that the notion of a distinct SL pathway is inappropriate. Instead, there may be a diversity of SL-like compounds, the response to which requires a D14/D14-like protein. This ancestral system could have been refined toward distinct ligand-specific pathways channeled through MAX2, the most downstream known component of SL signaling. MAX2 is tightly conserved among land plants and is more diverged from its nearest sister clade than any other SL-related gene, suggesting a pivotal role in the evolution of SL signaling. By contrast, the evidence suggests much greater flexibility upstream of MAX2. The MAX1 gene is a particularly strong candidate for contributing to diversification of inputs upstream of MAX2. Our functional analysis of the MAX1 family demonstrates the early origin of its catalytic function and both redundancy and functional diversification associated with its duplication in angiosperm lineages.

Characterization of the basal angiosperm Aristolochia fimbriata: a potential experimental system for genetic studies

BACKGROUND

Previous studies in basal angiosperms have provided insight into the diversity within the angiosperm lineage and helped to polarize analyses of flowering plant evolution. However, there is still not an experimental system for genetic studies among basal angiosperms to facilitate comparative studies and functional investigation. It would be desirable to identify a basal angiosperm experimental system that possesses many of the features found in existing plant model systems (e.g., Arabidopsis and Oryza).

RESULTS

We have considered all basal angiosperm families for general characteristics important for experimental systems, including availability to the scientific community, growth habit, and membership in a large basal angiosperm group that displays a wide spectrum of phenotypic diversity. Most basal angiosperms are woody or aquatic, thus are not well-suited for large scale cultivation, and were excluded. We further investigated members of Aristolochiaceae for ease of culture, life cycle, genome size, and chromosome number. We demonstrated self-compatibility for Aristolochia elegans and A. fimbriata, and transformation with a GFP reporter construct for Saruma henryi and A. fimbriata. Furthermore, A. fimbriata was easily cultivated with a life cycle of just three months, could be regenerated in a tissue culture system, and had one of the smallest genomes among basal angiosperms. An extensive multi-tissue EST dataset was produced for A. fimbriata that includes over 3.8 million 454 sequence reads.

CONCLUSIONS

Aristolochia fimbriata has numerous features that facilitate genetic studies and is suggested as a potential model system for use with a wide variety of technologies. Emerging genetic and genomic tools for A. fimbriata and closely related species can aid the investigation of floral biology, developmental genetics, biochemical pathways important in plant-insect interactions as well as human health, and various other features present in early angiosperms.

Comparative metabolic profiling between desiccation-sensitive and desiccation-tolerant species of Selaginella reveals insights into the resurrection trait

Spike mosses (Selaginellaceae) represent an ancient lineage of vascular plants in which some species have evolved desiccation tolerance (DT). A sister-group contrast to reveal the metabolic basis of DT was conducted between a desiccation-tolerant species, Selaginella lepidophylla, and a desiccation-sensitive species, Selaginella moellendorffii, at 100% relative water content (RWC) and 50% RWC using non-biased, global metabolomics profiling technology, based on GC/MS and UHLC/MS/MS(2) platforms. A total of 301 metabolites, including 170 named (56.5%) and 131 (43.5%) unnamed compounds, were characterized across both species. S.  lepidophylla retained significantly higher abundances of sucrose, mono- and polysaccharides, and sugar alcohols than did S. moellendorffii. Aromatic amino acids, the well-known osmoprotectant betaine and flavonoids were also more abundant in S. lepidophylla. Notably, levels of γ-glutamyl amino acid, linked with glutathione metabolism in the detoxification of reactive oxygen species, and with possible nitrogen remobilization following rehydration, were markedly higher in S. lepidophylla. Markers for lipoxygenase activity were also greater in S. lepidophylla, especially at 50% RWC. S. moellendorffii contained more than twice the number of unnamed compounds, with only a slightly greater abundance than in S. lepidophylla. In contrast, S. lepidophylla contained 14 unnamed compounds of fivefold or greater abundance than in S. moellendorffii, suggesting that these compounds might play critical roles in DT. Overall, S. lepidophylla appears poised to tolerate desiccation in a constitutive manner using a wide range of metabolites with some inducible components, whereas S. moellendorffii mounts only limited metabolic responses to dehydration stress.

Comparative metabolic profiling between desiccation-sensitive and desiccation-tolerant species of Selaginella reveals insights into the resurrection trait

Spike mosses (Selaginellaceae) represent an ancient lineage of vascular plants in which some species have evolved desiccation tolerance (DT). A sister-group contrast to reveal the metabolic basis of DT was conducted between a desiccation-tolerant species, Selaginella lepidophylla, and a desiccation-sensitive species, Selaginella moellendorffii, at 100% relative water content (RWC) and 50% RWC using non-biased, global metabolomics profiling technology, based on GC/MS and UHLC/MS/MS(2) platforms. A total of 301 metabolites, including 170 named (56.5%) and 131 (43.5%) unnamed compounds, were characterized across both species. S.  lepidophylla retained significantly higher abundances of sucrose, mono- and polysaccharides, and sugar alcohols than did S. moellendorffii. Aromatic amino acids, the well-known osmoprotectant betaine and flavonoids were also more abundant in S. lepidophylla. Notably, levels of γ-glutamyl amino acid, linked with glutathione metabolism in the detoxification of reactive oxygen species, and with possible nitrogen remobilization following rehydration, were markedly higher in S. lepidophylla. Markers for lipoxygenase activity were also greater in S. lepidophylla, especially at 50% RWC. S. moellendorffii contained more than twice the number of unnamed compounds, with only a slightly greater abundance than in S. lepidophylla. In contrast, S. lepidophylla contained 14 unnamed compounds of fivefold or greater abundance than in S. moellendorffii, suggesting that these compounds might play critical roles in DT. Overall, S. lepidophylla appears poised to tolerate desiccation in a constitutive manner using a wide range of metabolites with some inducible components, whereas S. moellendorffii mounts only limited metabolic responses to dehydration stress.

The glycosyltransferase repertoire of the spikemoss Selaginella moellendorffii and a comparative study of its cell wall

Spike mosses are among the most basal vascular plants, and one species, Selaginella moellendorffii, was recently selected for full genome sequencing by the Joint Genome Institute (JGI). Glycosyltransferases (GTs) are involved in many aspects of a plant life, including cell wall biosynthesis, protein glycosylation, primary and secondary metabolism. Here, we present a comparative study of the S. moellendorffii genome across 92 GT families and an additional family (DUF266) likely to include GTs. The study encompasses the moss Physcomitrella patens, a non-vascular land plant, while rice and Arabidopsis represent commelinid and non-commelinid seed plants. Analysis of the subset of GT-families particularly relevant to cell wall polysaccharide biosynthesis was complemented by a detailed analysis of S. moellendorffii cell walls. The S. moellendorffii cell wall contains many of the same components as seed plant cell walls, but appears to differ somewhat in its detailed architecture. The S. moellendorffii genome encodes fewer GTs (287 GTs including DUF266s) than the reference genomes. In a few families, notably GT51 and GT78, S. moellendorffii GTs have no higher plant orthologs, but in most families S. moellendorffii GTs have clear orthologies with Arabidopsis and rice. A gene naming convention of GTs is proposed which takes orthologies and GT-family membership into account. The evolutionary significance of apparently modern and ancient traits in S. moellendorffii is discussed, as is its use as a reference organism for functional annotation of GTs.

Phylogenetics and evolution of Trx SET genes in fully sequenced land plants

Plant Trx SET proteins are involved in H3K4 methylation and play a key role in plant floral development. Genes encoding Trx SET proteins constitute a multigene family in which the copy number varies among plant species and functional divergence appears to have occurred repeatedly. To investigate the evolutionary history of the Trx SET gene family, we made a comprehensive evolutionary analysis on this gene family from 13 major representatives of green plants. A novel clustering (here named as cpTrx clade), which included the III-1, III-2, and III-4 orthologous groups, previously resolved was identified. Our analysis showed that plant Trx proteins possessed a variety of domain organizations and gene structures among paralogs. Additional domains such as PHD, PWWP, and FYR were early integrated into primordial SET-PostSET domain organization of cpTrx clade. We suggested that the PostSET domain was lost in some members of III-4 orthologous group during the evolution of land plants. At least four classes of gene structures had been formed at the early evolutionary stage of land plants. Three intronless orphan Trx SET genes from the Physcomitrella patens (moss) were identified, and supposedly, their parental genes have been eliminated from the genome. The structural differences among evolutionary groups of plant Trx SET genes with different functions were described, contributing to the design of further experimental studies.

Flow cytometry confirms reticulate evolution and reveals triploidy in Central European Diphasiastrum taxa (Lycopodiaceae, Lycophyta)

BACKGROUND AND AIMS

Interspecific Diphasiastrum hybrids have been assumed to be homoploid and to produce well-formed spores serving sexual reproduction. If this were the case, forms intermediate between hybrids and parents or hybrid swarms should be expected. The purpose of this study was: (1) to check whether homoploidy consistently applies to the three hybrids throughout their Central European range; (2) to examine whether their genome sizes confirm their parentage as assumed by morphology; and (3) to perform a screening for detection of ploidy levels other than diploid and variation in DNA content due to backcrossing.

METHODS

Flow cytometry was used first to measure the relative DNA values [with 4',6-diamidino-2-phenylindole (DAPI) staining] and ploidy level as a general screening, and secondly to determine the absolute DNA 2C values [with propidium iodide (PI) staining] in a number of selected samples with the main focus on the hybrids.

KEY RESULTS

A considerable variation of DNA 2C values (5·26-7·52 pg) was detected between the three European Diphasiastrum species. The values of the diploid hybrids are highly constant without significant variation between regions. They are also intermediate between their assumed parents and agree closely with those calculated from their putative parents. This confirms their hybrid origin, assumed parentage and homoploid status. Considerably higher DNA amounts (9·48-10·30 pg) were obtained for three populations, suggesting that these represent triploid hybrids, an interpretation that is strongly supported by their morphology.

CONCLUSIONS

Diploid hybrids have retained their genetic and morphological identites throughout their Central European range, and thus no indications for diploid backcrossing were found. The triploid hybrids have probably originated from backcrossing between a diploid gametophyte of a hybrid (derived from a diplospore) and a haploid gametophyte of a diploid parental species. By repeated crossing events, reticulate evolution patterns arise that are similar to those known for a number of ferns.

The Selaginella genome identifies genetic changes associated with the evolution of vascular plants

Vascular plants appeared ~410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).

Estimates of nuclear DNA content in 98 species of brown algae (Phaeophyta)

BACKGROUND AND AIMS

Brown algae are critical components of marine ecosystems around the world. However, the genome of only one species of the class has so far been sequenced. This contrasts with numerous sequences available for model organisms such as higher plants, flies or worms. The present communication expands our coverage of DNA content information to 98 species of brown algae with a view to facilitating further genomic investigations of the class.

METHODOLOGY

The DNA-localizing fluorochrome DAPI (4',6-diamidino-2-phenylindole) and the red blood cell (chicken erythrocyte) standard were used to estimate 2C values by static microspectrophotometry.

PRINCIPAL RESULTS

2C DNA contents are reported for 98 species of brown algae, almost doubling the number of estimates available for the class. The present results also expand the reported DNA content range to 0.2-3.6 pg, with several species of Fucales and Laminariales containing apparent polyploid genomes with 2C = 1.8-3.6 pg.

CONCLUSIONS

The data provide DNA content values for 12 of the 19 recognized orders of brown algae spanning the breadth of the class. Despite earlier contentions concerning DNA content and the presence of oogamy, the present results do not support a correlation between phylogenetic placement and genome size. The closest sister groups to the brown algae have genome sizes on the order of 0.3 pg (e.g. Schizocladiophyceae), suggesting that this may be the ancestral genome size. However, DNA content ranges widely across the class.

Evolution of the Cinnamyl/Sinapyl Alcohol Dehydrogenase (CAD/SAD) gene family: the emergence of real lignin is associated with the origin of Bona Fide CAD

Lignin plays a vital role in plant adaptation to terrestrial environments. The cinnamyl alcohol dehydrogenase (CAD) catalyzes the last step in monolignol biosynthesis and might have contributed to the lignin diversity in plants. To investigate the evolutionary history and functional differentiation of the CAD gene family, we made a comprehensive evolutionary analysis of this gene family from 52 species, including bacteria, early eukaryotes and green plants. The phylogenetic analysis, together with gene structure and function, indicates that all members of land plants, except two of moss, could be divided into three classes. Members of Class I (bona fide CAD), generally accepted as the primary genes involved in the monolignol biosynthesis, are all from vascular plants, and form a robustly supported monophyletic group with the lycophyte CADs at the basal position. This class is also conserved in the predicted three-dimensional structure and the residues constituting the substrate-binding pocket of the proteins. Given that Selaginella has real lignin, the above evidence strongly suggests that the earliest occurrence of the bona fide CAD in the lycophyte could be directly correlated with the origin of lignin. Class II comprises members more similar to the aspen sinapyl alcohol dehydrogenase gene, and includes three groups corresponding to lycophyte, gymnosperm, and angiosperm. Class III is conserved in land plants. The three classes differ in patterns of evolution and expression, implying that functional divergence has occurred among them. Our study also supports the hypothesis of convergent evolution of lignin biosynthesis between red algae and vascular plants.

Comparison of the chloroplast peroxidase system in the chlorophyte Chlamydomonas reinhardtii, the bryophyte Physcomitrella patens, the lycophyte Selaginella moellendorffii and the seed plant Arabidopsis thaliana

BACKGROUND

Oxygenic photosynthesis is accompanied by the formation of reactive oxygen species (ROS), which damage proteins, lipids, DNA and finally limit plant yield. The enzymes of the chloroplast antioxidant system are exclusively nuclear encoded. During evolution, plastid and mitochondrial genes were post-endosymbiotically transferred to the nucleus, adapted for eukaryotic gene expression and post-translational protein targeting and supplemented with genes of eukaryotic origin.

RESULTS

Here, the genomes of the green alga Chlamydomonas reinhardtii, the moss Physcomitrella patens, the lycophyte Selaginella moellendorffii and the seed plant Arabidopsis thaliana were screened for ORFs encoding chloroplast peroxidases. The identified genes were compared for their amino acid sequence similarities and gene structures. Stromal and thylakoid-bound ascorbate peroxidases (APx) share common splice sites demonstrating that they evolved from a common ancestral gene. In contrast to most cormophytes, our results predict that chloroplast APx activity is restricted to the stroma in Chlamydomonas and to thylakoids in Physcomitrella. The moss gene is of retrotransposonal origin.The exon-intron-structures of 2CP genes differ between chlorophytes and streptophytes indicating an independent evolution. According to amino acid sequence characteristics only the A-isoform of Chlamydomonas 2CP may be functionally equivalent to streptophyte 2CP, while the weakly expressed B- and C-isoforms show chlorophyte specific surfaces and amino acid sequence characteristics. The amino acid sequences of chloroplast PrxII are widely conserved between the investigated species. In the analyzed streptophytes, the genes are unspliced, but accumulated four introns in Chlamydomonas. A conserved splice site indicates also a common origin of chlorobiont PrxQ.The similarity of splice sites also demonstrates that streptophyte glutathione peroxidases (GPx) are of common origin. Besides a less related cysteine-type GPx, Chlamydomonas encodes two selenocysteine-type GPx. The latter were lost prior or during streptophyte evolution.

CONCLUSION

Throughout plant evolution, there was a strong selective pressure on maintaining the activity of all three investigated types of peroxidases in chloroplasts. APx evolved from a gene, which dates back to times before differentiation of chlorobionts into chlorophytes and streptophytes, while Prx and presumably also GPx gene patterns may have evolved independently in the streptophyte and chlorophyte branches.

Evolutionary genomics revealed interkingdom distribution of Tcn1-like chromodomain-containing Gypsy LTR retrotransposons among fungi and plants

Background

Chromodomain-containing Gypsy LTR retrotransposons or chromoviruses are widely distributed among eukaryotes and have been found in plants, fungi and vertebrates. The previous comprehensive survey of chromoviruses from mosses (Bryophyta) suggested that genomes of non-seed plants contain the clade which is closely related to the retrotransposons from fungi. The origin, distribution and evolutionary history of this clade remained unclear mainly due to the absence of information concerning the diversity and distribution of LTR retrotransposons in other groups of non-seed plants as well as in fungal genomes.

Results

In present study we preformed in silico analysis of chromodomain-containing LTR retrotransposons in 25 diverse fungi and a number of plant species including spikemoss Selaginella moellendorffii (Lycopodiophyta) coupled with an experimental survey of chromodomain-containing Gypsy LTR retrotransposons from diverse non-seed vascular plants (lycophytes, ferns, and horsetails). Our mining of Gypsy LTR retrotransposons in genomic sequences allowed identification of numerous families which have not been described previously in fungi. Two new well-supported clades, Galahad and Mordred, as well as several other previously unknown lineages of chromodomain-containing Gypsy LTR retrotransposons were described based on the results of PCR-mediated survey of LTR retrotransposon fragments from ferns, horsetails and lycophytes. It appeared that one of the clades, namely Tcn1 clade, was present in basidiomycetes and non-seed plants including mosses (Bryophyta) and lycophytes (genus Selaginella).

Conclusions

The interkingdom distribution is not typical for chromodomain-containing LTR retrotransposons clades which are usually very specific for a particular taxonomic group. Tcn1-like LTR retrotransposons from fungi and non-seed plants demonstrated high similarity to each other which can be explained by strong selective constraints and the 'retained' genes theory or by horizontal transmission.

Construction and characterization of twoCitrusBAC libraries and identification of clones containing the phytoene synthase gene

Two deep-coverage Bacterial Artificial Chromosome (BAC) libraries of Citrus sinensis (L.) Osbeck ‘Cara Cara’ navel orange and Citrus reticulata (L.) Blanco ‘Egan No. 1’ Ponkan mandarin, which belong to the two most important species of the Citrus genus, have been constructed and characterized to facilitate gene cloning and to analyze variety-specific genome composition. The C. sinensis BAC library consists of 36 000 clones with negligible false-positive clones and an estimated average insert size of 126 kb covering ~4.5 × 109 bp and thus providing an 11.8-fold coverage of haploid genome equivalents, whereas the C. reticulata library consists of 21 000 clones also with negligible false-positive clones and an estimated average of 120 kb covering ~2.5 × 109 bp representing a 6.6-fold coverage of haploid genome equivalents. Both libraries were evaluated for contamination with high-copy vector, empty pIndigoBAC536 vector, and organellar DNA sequences. Screening has been performed by Southern hybridization of BAC filters, which results in <0.5% chloroplast DNA contamination and no mitochondrial DNA contamination in both libraries. Eight and five positive clones harboring the gene encoding Phytoene synthase (Psy (EC 2.5.1.32)) were identified from the C. sinensis and C. reticulata libraries, respectively, using the filter hybridization procedure. These results suggest that the two BAC libraries are useful tools for the isolation of functional genes and advanced genomics research in the two important species C. sinensis and C. reticulata. Resources, high-density filters, individual clones, and whole libraries are available for public distribution and are accessible at the National Key Laboratory of Crop Genetic Improvement, Huazhong Agricultural University.

Evolution of phage-type RNA polymerases in higher plants: characterization of the single phage-type RNA polymerase gene from Selaginella moellendorffii

Selaginella moellendorfii (spikemoss) sequence trace data encoding a polypeptide highly similar to angiosperm and moss phage-type organelle RNA polymerases (RpoTs) were used to isolate a BAC clone containing the full-length gene SmRpoT as well as the corresponding cDNA. The SmRpoT mRNA comprises 3452 nt with an open reading frame of 3006 nt, encoding a putative protein of 1002 amino acids with a molecular mass of 113 kDa. The SmRpoT gene comprises 19 exons and 18 introns, conserved in their position with those of the angiosperm and Physcomitrella RpoT genes. In phylogenetic analyses, the Selaginella RpoT polymerase is in a sister position to all other phage-type polymerases of angiosperms. However, according to its conserved exon-intron structure, the Selaginella RpoT gene is representative of the molecular evolutionary lineage giving rise to the RpoT gene family of flowering plants. The N-terminal transit peptide of SmRpoT is shown to confer targeting of green fluorescent protein exclusively to mitochondria after transient expression in Arabidopsis and Selaginella protoplasts. Angiosperms and the moss P. patens possess small gene families encoding RpoTs, which include mitochondrial- and chloroplast-targeted RNA polymerases. In striking contrast, the Selaginella RpoT gene is shown to be single-copy, although Selaginella, as a lycophyte, has a phylogenetic position between Physcomitrella and angiosperms. Thus, there is no evidence that Selaginella may contain a nuclear-encoded phage-type chloroplast RNA polymerase.

Selaginella and 400 million years of separation

Selaginella (spikemoss) is an enigma in the plant kingdom. Although a fascination to botanists at the turn of the twentieth century, members of this genus are unremarkable in appearance, never flower, and are of no agronomic value. However, members of this genus are relicts from ancient times, and one has to marvel at how this genus has survived virtually unchanged in appearance for hundreds of millions of years. In light of the recent completion of the Selaginella moellendorffii genome sequence, this review is intended to survey what is known about Selaginella, with a special emphasis on recent inquiries into its unique biology and importance in understanding the early evolution of vascular plants.

Chapter 4. New model systems for the study of developmental evolution in plants

The number of genetically tractable plant model systems is rapidly increasing, thanks to the decreasing cost of sequencing and the wide amenability of plants to stable transformation and other functional approaches. In this chapter, I discuss emerging model systems from throughout the land plant phylogeny and consider how their unique attributes are contributing to our understanding of development, evolution, and ecology. These new models are being developed using two distinct strategies: in some cases, they are selected because of their close relationship to the established models, while in others, they are chosen with the explicit intention of exploring distantly related plant lineages. Such complementary approaches are yielding exciting new results that shed light on both micro- and macroevolutionary processes in the context of developmental evolution.

The highest-copy repeats are methylated in the small genome of the early divergent vascular plant Selaginella moellendorffii

Background

The lycophyte Selaginella moellendorffii is a vascular plant that diverged from the fern/seed plant lineage at least 400 million years ago. Although genomic information for S. moellendorffii is starting to be produced, little is known about basic aspects of its molecular biology. In order to provide the first glimpse to the epigenetic landscape of this early divergent vascular plant, we used the methylation filtration technique. Methylation filtration genomic libraries select unmethylated DNA clones due to the presence of the methylation-dependent restriction endonuclease McrBC in the bacterial host.

Results

We conducted a characterization of the DNA methylation patterns of the S. moellendorffii genome by sequencing a set of S. moellendorffii shotgun genomic clones, along with a set of methylation filtered clones. Chloroplast DNA, which is typically unmethylated, was enriched in the filtered library relative to the shotgun library, showing that there is DNA methylation in the extremely small S. moellendorffii genome. The filtered library also showed enrichment in expressed and gene-like sequences, while the highest-copy repeats were largely under-represented in this library. These results show that genes and repeats are differentially methylated in the S. moellendorffii genome, as occurs in other plants studied.

Conclusion

Our results shed light on the genome methylation pattern in a member of a relatively unexplored plant lineage. The DNA methylation data reported here will help understanding the involvement of this epigenetic mark in fundamental biological processes, as well as the evolutionary aspects of epigenetics in land plants.

Duplication and functional diversification of HAP3 genes leading to the origin of the seed-developmental regulatory gene, LEAFY COTYLEDON1 (LEC1), in nonseed plant genomes

The HAP3 gene encodes a subunit of the CCAAT-box-binding factor (CBF), a highly conserved trimeric activator that recognizes and binds the ubiquitous CCAAT promoter element with high affinity. Two types of HAP3 gene have been identified in plant genomes. The LEAFY COTYLEDON1 (LEC1)-type HAP3 genes encode a functionally specialized subunit of CBF, which is expressed specifically in developing seeds. In contrast, most non-LEC1-type HAP3 genes are expressed in various tissues. It has been proposed that the LEC1-type HAP3 genes originated from the duplication and functional divergence of non-LEC1-type HAP3 genes. However, it is not yet known when this duplication event took place or whether the LEC1-type HAP3 genes appeared at the same time as the origin of seed plants. Here we describe a comprehensive comparison of the duplication patterns of HAP3 genes in different plant genomes. We recognize a major expansion of the HAP3 gene family accompanying the origin and early diversification of land plants and postulate that retrotransposition and other mechanisms of gene duplication have been involved in the expansion of the plant HAP3 gene family. We provide evidence that the LEC1-type HAP3 genes originated in nonseed vascular plant genomes and demonstrate that they are inductively expressed under drought stress in nonseed plants. These genes, however, were recruited to a novel regulatory network in the early stages of seed plant evolution and steadily expressed during seed development and maturation.

Evolutionary conservation of plant gibberellin signalling pathway components

BACKGROUND

Gibberellins (GA) are plant hormones that can regulate germination, elongation growth, and sex determination. They ubiquitously occur in seed plants. The discovery of gibberellin receptors, together with advances in understanding the function of key components of GA signalling in Arabidopsis and rice, reveal a fairly short GA signal transduction route. The pathway essentially consists of GID1 gibberellin receptors that interact with F-box proteins, which in turn regulate degradation of downstream DELLA proteins, suppressors of GA-controlled responses.

RESULTS

Arabidopsis sequences of the gibberellin signalling compounds were used to screen databases from a variety of plants, including protists, for homologues, providing indications for the degree of conservation of the pathway. The pathway as such appears completely absent in protists, the moss Physcomitrella patens shares only a limited homology with the Arabidopsis proteins, thus lacking essential characteristics of the classical GA signalling pathway, while the lycophyte Selaginella moellendorffii contains a possible ortholog for each component. The occurrence of classical GA responses can as yet not be linked with the presence of homologues of the signalling pathway. Alignments and display in neighbour joining trees of the GA signalling components confirm the close relationship of gymnosperms, monocotyledonous and dicotyledonous plants, as suggested from previous studies.

CONCLUSION

Homologues of the GA-signalling pathway were mainly found in vascular plants. The GA signalling system may have its evolutionary molecular onset in Physcomitrella patens, where GAs at higher concentrations affect gravitropism and elongation growth.

Nuclear genome size in Selaginella

Estimates of nuclear genome size for 9 Selaginella species were obtained using flow cytometry, and measurements for 7 of these species are reported for the first time. Estimates range from 0.086 to 0.112 pg per holoploid genome (84-110 Mb). The data presented here agree with the previously published flow cytometric results for S. moellendorffii. Within the 9 species sampled here, chromosome number varies from 2n = 16 to 2n = 27. Nuclear genome size appears to be strongly correlated with chromosome number (Spearman's rank correlation; p = 0.00003725). Cultivated S. moellendorffii lacks sexual reproduction--manifest by the production of abortive megasporangia. Flow cytometric data generated from a herbarium specimen of a fertile wild-collected S. moellendorffii are virtually indistinguishable from the data generated from fresh material (0.088 vs. 0.089 pg/1C). Therefore, the limited fertility observed in cultivated plants is probably not the result of abnormal chromosome number (e.g., induced by interspecific hybridization).

Common functions for diverse small RNAs of land plants

Endogenous small RNAs, including microRNAs (miRNAs) and short interfering RNAs (siRNAs), are critical components of plant gene regulation. Some abundant miRNAs involved in developmental control are conserved between anciently diverged plants, while many other less-abundant miRNAs appear to have recently emerged in the Arabidopsis thaliana lineage. Using large-scale sequencing of small RNAs, we extended the known diversity of miRNAs in basal plants to include 88 confidently annotated miRNA families in the moss Physcomitrella patens and 44 in the lycopod Selaginella moellendorffii. Cleavage of 29 targets directed by 14 distinct P. patens miRNA families and a trans-acting siRNA (ta-siRNA) was experimentally confirmed. Despite a core set of 12 miRNA families also expressed in angiosperms, weakly expressed and apparently lineage-specific miRNAs accounted for the majority of miRNA diversity in both species. Nevertheless, the molecular functions of several of these lineage-specific small RNAs matched those of angiosperms, despite dissimilarities in the small RNA sequences themselves, including small RNAs that mediated negative feedback regulation of the miRNA pathway and miR390-dependent ta-siRNAs that guided the cleavage of AUXIN RESPONSE FACTOR mRNAs. Diverse, lineage-specific, small RNAs can therefore perform common biological functions in plants.

Nuclear DNA content estimates in green algal lineages: chlorophyta and streptophyta

BACKGROUND AND AIMS

Consensus higher-level molecular phylogenies present a compelling case that an ancient divergence separates eukaryotic green algae into two major monophyletic lineages, Chlorophyta and Streptophyta, and a residuum of green algae, which have been referred to prasinophytes or micromonadophytes. Nuclear DNA content estimates have been published for less than 1% of the described green algal members of Chlorophyta, which includes multicellular green marine algae and freshwater flagellates (e.g. Chlamydomonas and Volvox). The present investigation summarizes the state of our knowledge and adds substantially to our database of C-values, especially for the streptophyte charophycean lineage which is the sister group of the land plants. A recent list of 2C nuclear DNA contents for isolates and species of green algae is expanded by 72 to 157.

METHODS

The DNA-localizing fluorochrome DAPI (4',6-diamidino-2-phenylindole) and red blood cell (chicken erythrocytes) standard were used to estimate 2C values with static microspectrophotometry.

KEY RESULTS

In Chlorophyta, including Chlorophyceae, Prasinophyceae, Trebouxiophyceae and Ulvophyceae, 2C DNA estimates range from 0.01 to 5.8 pg. Nuclear DNA content variation trends are noted and discussed for specific problematic taxon pairs, including Ulotrichales-Ulvales, and Cladophorales-Siphonocladales. For Streptophyta, 2C nuclear DNA contents range from 0.2 to 6.4 pg, excluding the highly polyploid Charales and Desmidiales, which have genome sizes of up to 14.8 and 46.8 pg, respectively. Nuclear DNA content data for Streptophyta superimposed on a contemporary molecular phylogeny indicate that early diverging lineages, including some members of Chlorokybales, Coleochaetales and Klebsormidiales, have genomes as small as 0.1-0.5 pg. It is proposed that the streptophyte ancestral nuclear genome common to both the charophyte and the embryophyte lineages can be characterized as 1C = 0.2 pg and 1n = 6.

CONCLUSIONS

These data will help pre-screen candidate species for the on-going construction of bacterial artificial chromosome nuclear genome libraries for land plant ancestors. Data for the prasinophyte Mesostigma are of particular interest as this alga reportedly most closely resembles the 'ancestral green flagellate'. Both mechanistic and ecological processes are discussed that could have produced the observed C-value increase of >100-fold in the charophyte green algae whereas the ancestral genome was conserved in the embryophytes.

Purification of mitochondrial and plastid DNA

The size, structure and conformation of mitochondrial and plastid genomes differ dramatically among eukaryotes. Similarly, the yield and purity of extracted organelle DNA also vary, and are crucial factors for the success of restriction mapping and sequencing experiments. We describe here procedures for the purification of organelle DNA from a broad range of eukaryotes. By emphasizing the underlying principles, these procedures will facilitate the development of new species-specific protocols. The presented purification schemes involve either isolation of organelles and subsequent extraction of DNA from this subcellular fraction, or processing of whole-cell lysates followed by CsCl gradient centrifugation to separate nuclear and organelle DNAs according to their A + T content. We have successfully used the described procedures for organelle genome sequencing from diverse eukaryotes, including non-axenic protists. Procedures can be completed in 3-5 days, typically yielding a few micrograms of DNA-ample for sequencing complete genomes.

The fern as a model system to study photomorphogenesis

The fern gametophyte is a good model system for studying cell biological, physiological, and photobiological aspects of the fundamental processes of plant development and physiological phenomena, because of its autotrophic characteristics and its simple structure. The cells, moreover, are not surrounded by tissue, so observation and manipulation of the cells are very easy. Here I summarize a part of my knowledge of fern systems, which I have studied for nearly 40 years.

Widespread genome duplications throughout the history of flowering plants

Genomic comparisons provide evidence for ancient genome-wide duplications in a diverse array of animals and plants. We developed a birth-death model to identify evidence for genome duplication in EST data, and applied a mixture model to estimate the age distribution of paralogous pairs identified in EST sets for species representing the basal-most extant flowering plant lineages. We found evidence for episodes of ancient genome-wide duplications in the basal angiosperm lineages including Nuphar advena (yellow water lily: Nymphaeaceae) and the magnoliids Persea americana (avocado: Lauraceae), Liriodendron tulipifera (tulip poplar: Magnoliaceae), and Saruma henryi (Aristolochiaceae). In addition, we detected independent genome duplications in the basal eudicot Eschscholzia californica (California poppy: Papaveraceae) and the basal monocot Acorus americanus (Acoraceae), both of which were distinct from duplications documented for ancestral grass (Poaceae) and core eudicot lineages. Among gymnosperms, we found equivocal evidence for ancient polyploidy in Welwitschia mirabilis (Gnetales) and no evidence for polyploidy in pine, although gymnosperms generally have much larger genomes than the angiosperms investigated. Cross-species sequence divergence estimates suggest that synonymous substitution rates in the basal angiosperms are less than half those previously reported for core eudicots and members of Poaceae. These lower substitution rates permit inference of older duplication events. We hypothesize that evidence of an ancient duplication observed in the Nuphar data may represent a genome duplication in the common ancestor of all or most extant angiosperms, except Amborella.

Functional analysis and comparative genomics of expressed sequence tags from the lycophyte Selaginella moellendorffii

BACKGROUND

The lycophyte Selaginella moellendorffii is a member of one of the oldest lineages of vascular plants on Earth. Fossil records show that the lycophyte clade arose 400 million years ago, 150-200 million years earlier than angiosperms, a group of plants that includes the well-studied flowering plant Arabidopsis thaliana. S. moellendorffii has a genome size of approximately 100 Mbp, as small or smaller than that of A. thaliana. S. moellendorffii has the potential to provide significant comparative information to better understand the evolution of vascular plants.

RESULTS

We sequenced 2181 Expressed Sequence Tags (ESTs) from a S. moellendorffii cDNA library. One thousand three hundred and one non-redundant sequences were assembled, containing 291 contigs and 1010 singletons. Approximately 75% of the ESTs matched proteins in the non-redundant protein database. Among 1301 clusters, 343 were categorized according to Gene Ontology (GO) hierarchy and were compared to the GO mapping of A. thaliana tentative consensus sequences. We compared S. moellendorffii ESTs to the A. thaliana and Physcomitrella patens EST databases, using the tBLASTX algorithm. Approximately 60% of the ESTs exhibited similarity with both A. thaliana and P. patens ESTs; whereas, 13% and 1% of the ESTs had exclusive similarity with A. thaliana and P. patens ESTs, respectively. A substantial proportion of the ESTs (26%) had no match with A. thaliana or P. patens ESTs.

CONCLUSION

We discovered 1301 putative unigenes in S. moellendorffii. These results give an initial insight into its transcriptome that will aid in the study of the S. moellendorffii genome in the near future.