High-throughput protein characterization by complementation using DNA barcoded fragment libraries.

Our ability to predict, control, or design biological function is fundamentally limited by poorly annotated gene function. This can be particularly challenging in non-model systems. Accordingly, there is motivation for new high-throughput methods for accurate functional annotation. Here, we used complementation of auxotrophs and DNA barcode sequencing (Coaux-Seq) to enable high-throughput characterization of protein function. Fragment libraries from eleven genetically diverse bacteria were tested in twenty different auxotrophic strains of Escherichia coli to identify genes that complement missing biochemical activity. We recovered 41% of expected hits, with effectiveness ranging per source genome, and observed success even with distant E. coli relatives like Bacillus subtilis and Bacteroides thetaiotaomicron. Coaux-Seq provided the first experimental validation for 53 proteins, of which 11 are less than 40% identical to an experimentally characterized protein. Among the unexpected function identified was a sulfate uptake transporter, an O-succinylhomoserine sulfhydrylase for methionine synthesis, and an aminotransferase. We also identified instances of cross-feeding wherein protein overexpression and nearby non-auxotrophic strains enabled growth. Altogether, Coaux-Seq's utility is demonstrated, with future applications in ecology, health, and engineering.

Mutant phenotypes for thousands of bacterial genes of unknown function.

One-third of all protein-coding genes from bacterial genomes cannot be annotated with a function. Here, to investigate the functions of these genes, we present genome-wide mutant fitness data from 32 diverse bacteria across dozens of growth conditions. We identified mutant phenotypes for 11,779 protein-coding genes that had not been annotated with a specific function. Many genes could be associated with a specific condition because the gene affected fitness only in that condition, or with another gene in the same bacterium because they had similar mutant phenotypes. Of the poorly annotated genes, 2,316 had associations that have high confidence because they are conserved in other bacteria. By combining these conserved associations with comparative genomics, we identified putative DNA repair proteins; in addition, we propose specific functions for poorly annotated enzymes and transporters and for uncharacterized protein families. Our study demonstrates the scalability of microbial genetics and its utility for improving gene annotations.

Correction: Filling gaps in bacterial amino acid biosynthesis pathways with high-throughput genetics.

[This corrects the article DOI: 10.1371/journal.pgen.1007147.].

Filling gaps in bacterial amino acid biosynthesis pathways with high-throughput genetics.

For many bacteria with sequenced genomes, we do not understand how they synthesize some amino acids. This makes it challenging to reconstruct their metabolism, and has led to speculation that bacteria might be cross-feeding amino acids. We studied heterotrophic bacteria from 10 different genera that grow without added amino acids even though an automated tool predicts that the bacteria have gaps in their amino acid synthesis pathways. Across these bacteria, there were 11 gaps in their amino acid biosynthesis pathways that we could not fill using current knowledge. Using genome-wide mutant fitness data, we identified novel enzymes that fill 9 of the 11 gaps and hence explain the biosynthesis of methionine, threonine, serine, or histidine by bacteria from six genera. We also found that the sulfate-reducing bacterium Desulfovibrio vulgaris synthesizes homocysteine (which is a precursor to methionine) by using DUF39, NIL/ferredoxin, and COG2122 proteins, and that homoserine is not an intermediate in this pathway. Our results suggest that most free-living bacteria can likely make all 20 amino acids and illustrate how high-throughput genetics can uncover previously-unknown amino acid biosynthesis genes.

Monofluorophosphate is a selective inhibitor of respiratory sulfate-reducing microorganisms.

Despite the environmental and economic cost of microbial sulfidogenesis in industrial operations, few compounds are known as selective inhibitors of respiratory sulfate reducing microorganisms (SRM), and no study has systematically and quantitatively evaluated the selectivity and potency of SRM inhibitors. Using general, high-throughput assays to quantitatively evaluate inhibitor potency and selectivity in a model sulfate-reducing microbial ecosystem as well as inhibitor specificity for the sulfate reduction pathway in a model SRM, we screened a panel of inorganic oxyanions. We identified several SRM selective inhibitors including selenate, selenite, tellurate, tellurite, nitrate, nitrite, perchlorate, chlorate, monofluorophosphate, vanadate, molydate, and tungstate. Monofluorophosphate (MFP) was not known previously as a selective SRM inhibitor, but has promising characteristics including low toxicity to eukaryotic organisms, high stability at circumneutral pH, utility as an abiotic corrosion inhibitor, and low cost. MFP remains a potent inhibitor of SRM growing by fermentation, and MFP is tolerated by nitrate and perchlorate reducing microorganisms. For SRM inhibition, MFP is synergistic with nitrite and chlorite, and could enhance the efficacy of nitrate or perchlorate treatments. Finally, MFP inhibition is multifaceted. Both inhibition of the central sulfate reduction pathway and release of cytoplasmic fluoride ion are implicated in the mechanism of MFP toxicity.

Control of methionine metabolism by the SahR transcriptional regulator in Proteobacteria.

Sulphur is an essential element in the metabolism. The sulphur-containing amino acid methionine is a metabolic precursor for S-adenosylmethionine (SAM), which serves as a coenzyme for ubiquitous methyltrtansferases. Recycling of organic sulphur compounds, e.g. via the SAM cycle, is an important metabolic process that needs to be tightly regulated. Knowledge about transcriptional regulation of these processes is still limited for many free-living bacteria. We identified a novel transcription factor SahR from the ArsR family that controls the SAM cycle genes in diverse microorganisms from soil and aquatic ecosystems. By using comparative genomics, we predicted SahR-binding DNA motifs and reconstructed SahR regulons in the genomes of 62 Proteobacteria. The conserved core of SahR regulons includes all enzymes required for the SAM cycle: the SAH hydrolase AhcY, the methionine biosynthesis enzymes MetE/MetH and MetF, and the SAM synthetase MetK. By using a combination of experimental techniques, we validated the SahR regulon in the sulphate-reducing Deltaproteobacterium Desulfovibrio alaskensis. SahR functions as a negative regulator that responds to the S-adenosylhomocysteine (SAH). The elevated SAH level in the cell dissociates SahR from its DNA operators and induces the expression of SAM cycle genes. The effector-sensing domain in SahR is related to SAM-dependent methylases that are able to tightly bind SAH. SahR represents a novel type of transcriptional regulators for the control of sulphur amino acid metabolism.

The energy-conserving electron transfer system used by Desulfovibrio alaskensis strain G20 during pyruvate fermentation involves reduction of endogenously formed fumarate and cytoplasmic and membrane-bound complexes, Hdr-Flox and Rnf.

The adaptation capability of Desulfovibrio to natural fluctuations in electron acceptor availability was evaluated by studying Desulfovibrio alaskensis strain G20 under varying respiratory, fermentative and methanogenic coculture conditions in chemostats. Transition from lactate to pyruvate in coculture resulted in a dramatic shift in the population structure and closer interspecies cell-to-cell interactions. Lower methane production rates in coculture than predicted from pyruvate input was attributed to redirection of electron flow to fumarate reduction. Without a methanogenic partner, accumulation of H2and formate resulted in greater succinate production. Comparative transcript and gene fitness analysis in concert with physiological data of G20 wildtype and mutants demonstrated that pyruvate fermentation involves respiration of cytoplasmically formed fumarate using cytoplasmic and membrane-bound energy-conserving complexes, Rnf, Hdr-Flox-1 and Hmc. At the low H2/formate levels maintained in coculture, Rnf likely functions as proton-pumping ferredoxin (Fd): type-I cytochrome c oxidoreductase, which transitions to a proton-pumping Fd(red): nicotinamide adenine dinucleotide (NAD⁺) oxidoreductase at high H2/formate levels during fermentation in monoculture. Hdr-Flox-1 is postulated to recycle Fd(red) via a flavin-based electron bifurcation involving NADH, Fdox and the thiol/disulphide-containing DsrC. In a menaquinone (MQ)-based electron confurcation reaction, the high-molecular-weight cytochrome-c3complex, Hmc, is proposed to then couple DsrC(red) and periplasmic H2/formate oxidation using the MQ pool to fuel a membrane-bound fumarate reductase.

Indirect and suboptimal control of gene expression is widespread in bacteria.

Gene regulation in bacteria is usually described as an adaptive response to an environmental change so that genes are expressed when they are required. We instead propose that most genes are under indirect control: their expression responds to signal(s) that are not directly related to the genes' function. Indirect control should perform poorly in artificial conditions, and we show that gene regulation is often maladaptive in the laboratory. In Shewanella oneidensis MR-1, 24% of genes are detrimental to fitness in some conditions, and detrimental genes tend to be highly expressed instead of being repressed when not needed. In diverse bacteria, there is little correlation between when genes are important for optimal growth or fitness and when those genes are upregulated. Two common types of indirect control are constitutive expression and regulation by growth rate; these occur for genes with diverse functions and often seem to be suboptimal. Because genes that have closely related functions can have dissimilar expression patterns, regulation may be suboptimal in the wild as well as in the laboratory.

Genomic and environmental controls on Castellaniella biogeography in an anthropogenically disturbed subsurface.

Castellaniella species have been isolated from a variety of mixed-waste environments including the nitrate and multiple metal-contaminated subsurface at the Oak Ridge Reservation (ORR). Previous studies examining microbial community composition and nitrate removal at ORR during biostimulation efforts reported increased abundances of members of the Castellaniella genus concurrent with increased denitrification rates. Thus, we asked how genomic and abiotic factors control the Castellaniella biogeography at the site to understand how these factors may influence nitrate transformation in an anthropogenically impacted setting. We report the isolation and characterization of several Castellaniella strains from the ORR subsurface. Five of these isolates match at 100% identity (at the 16S rRNA gene V4 region) to two Castellaniella amplicon sequence variants (ASVs), ASV1 and ASV2, that have persisted in the ORR subsurface for at least 2 decades. However, ASV2 has consistently higher relative abundance in samples taken from the site and was also the dominant blooming denitrifier population during a prior biostimulation effort. We found that the ASV2 representative strain has greater resistance to mixed metal stress than the ASV1 representative strains. We attribute this resistance, in part, to the large number of unique heavy metal resistance genes identified on a genomic island in the ASV2 representative genome. Additionally, we suggest that the relatively lower fitness of ASV1 may be connected to the loss of the nitrous oxide reductase (nos) operon (and associated nitrous oxide reductase activity) due to the insertion at this genomic locus of a mobile genetic element carrying copper resistance genes. This study demonstrates the value of integrating genomic, environmental, and phenotypic data to characterize the biogeography of key microorganisms in contaminated sites.

Flexibility of syntrophic enzyme systems in Desulfovibrio species ensures their adaptation capability to environmental changes.

The mineralization of organic matter in anoxic environments relies on the cooperative activities of hydrogen producers and consumers obligately linked by interspecies metabolite exchange in syntrophic consortia that may include sulfate reducing species such as Desulfovibrio. To evaluate the metabolic flexibility of syntrophic Desulfovibrio to adapt to naturally fluctuating methanogenic environments, we studied Desulfovibrio alaskensis strain G20 grown in chemostats under respiratory and syntrophic conditions with alternative methanogenic partners, Methanococcus maripaludis and Methanospirillum hungatei, at different growth rates. Comparative whole-genome transcriptional analyses, complemented by G20 mutant strain growth experiments and physiological data, revealed a significant influence of both energy source availability (as controlled by dilution rate) and methanogen on the electron transfer systems, ratios of interspecies electron carriers, energy generating systems, and interspecies physical associations. A total of 68 genes were commonly differentially expressed under syntrophic versus respiratory lifestyle. Under low-energy (low-growth-rate) conditions, strain G20 further had the capacity to adapt to the metabolism of its methanogenic partners, as shown by its differing gene expression of enzymes involved in the direct metabolic interactions (e.g., periplasmic hydrogenases) and the ratio shift in electron carriers used for interspecies metabolite exchange (hydrogen/formate). A putative monomeric [Fe-Fe] hydrogenase and Hmc (high-molecular-weight-cytochrome c3) complex-linked reverse menaquinone (MQ) redox loop become increasingly important for the reoxidation of the lactate-/pyruvate oxidation-derived redox pair, DsrC(red) and Fd(red), relative to the Qmo-MQ-Qrc (quinone-interacting membrane-bound oxidoreductase; quinone-reducing complex) loop. Together, these data underscore the high enzymatic and metabolic adaptive flexibility that likely sustains Desulfovibrio in naturally fluctuating methanogenic environments.

Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough.

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

Extreme reconfiguration of plastid genomes in the angiosperm family Geraniaceae: rearrangements, repeats, and codon usage.

Geraniaceae plastid genomes (plastomes) have experienced a remarkable number of genomic changes. The plastomes of Erodium texanum, Geranium palmatum, and Monsonia speciosa were sequenced and compared with other rosids and the previously published Pelargonium hortorum plastome. Geraniaceae plastomes were found to be highly variable in size, gene content and order, repetitive DNA, and codon usage. Several unique plastome rearrangements include the disruption of two highly conserved operons (S10 and rps2-atpA), and the inverted repeat (IR) region in M. speciosa does not contain all genes in the ribosomal RNA operon. The sequence of M. speciosa is unusually small (128,787 bp); among angiosperm plastomes sequenced to date, only those of nonphotosynthetic species and those that have lost one IR copy are smaller. In contrast, the plastome of P. hortorum is the largest, at 217,942 bp. These genomes have experienced numerous gene and intron losses and partial and complete gene duplications. Some of the losses are shared throughout the family (e.g., trnT-GGU and the introns of rps16 and rpl16); however, other losses are homoplasious (e.g., trnG-UCC intron in G. palmatum and M. speciosa). IR length is also highly variable. The IR in P. hortorum was previously shown to be greatly expanded to 76 kb, and the IR is lost in E. texanum and reduced in G. palmatum (11 kb) and M. speciosa (7 kb). Geraniaceae plastomes contain a high frequency of large repeats (>100 bp) relative to other rosids. Within each plastome, repeats are often located at rearrangement end points and many repeats shared among the four Geraniaceae flank rearrangement end points. GC content is elevated in the genomes and also in coding regions relative to other rosids. Codon usage per amino acid and GC content at third position sites are significantly different for Geraniaceae protein-coding sequences relative to other rosids. Our findings suggest that relaxed selection and/or mutational biases lead to increased GC content, and this in turn altered codon usage. We propose that increases in genomic rearrangements, repetitive DNA, nucleotide substitutions, and GC content may be caused by relaxed selection resulting from improper DNA repair.

A Defined Medium for Cultivation and Exometabolite Profiling of Soil Bacteria.

Exometabolomics is an approach to assess how microorganisms alter, or react to their environments through the depletion and production of metabolites. It allows the examination of how soil microbes transform the small molecule metabolites within their environment, which can be used to study resource competition and cross-feeding. This approach is most powerful when used with defined media that enable tracking of all metabolites. However, microbial growth media have traditionally been developed for the isolation and growth of microorganisms but not metabolite utilization profiling through Liquid Chromatography Tandem Mass Spectrometry (LC-MS/MS). Here, we describe the construction of a defined medium, the Northen Lab Defined Medium (NLDM), that not only supports the growth of diverse soil bacteria but also is defined and therefore suited for exometabolomic experiments. Metabolites included in NLDM were selected based on their presence in R2A medium and soil, elemental stoichiometry requirements, as well as knowledge of metabolite usage by different bacteria. We found that NLDM supported the growth of 108 of the 110 phylogenetically diverse (spanning 36 different families) soil bacterial isolates tested and all of its metabolites were trackable through LC-MS/MS analysis. These results demonstrate the viability and utility of the constructed NLDM medium for growing and characterizing diverse microbial isolates and communities.

Evidence-based annotation of transcripts and proteins in the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough.

We used high-resolution tiling microarrays and 5' RNA sequencing to identify transcripts in Desulfovibrio vulgaris Hildenborough, a model sulfate-reducing bacterium. We identified the first nucleotide position for 1,124 transcripts, including 54 proteins with leaderless transcripts and another 72 genes for which a major transcript initiates within the upstream protein-coding gene, which confounds measurements of the upstream gene's expression. Sequence analysis of these promoters showed that D. vulgaris prefers -10 and -35 boxes different from those preferred by Escherichia coli. A total of 549 transcripts ended at intrinsic (rho-independent) terminators, but most of the other transcripts seemed to have variable ends. We found low-level antisense expression of most genes, and the 5' ends of these transcripts mapped to promoter-like sequences. Because antisense expression was reduced for highly expressed genes, we suspect that elongation of nonspecific antisense transcripts is suppressed by transcription of the sense strand. Finally, we combined the transcript results with comparative analysis and proteomics data to make 505 revisions to the original annotation of 3,531 proteins: we removed 255 (7.5%) proteins, changed 123 (3.6%) start codons, and added 127 (3.7%) proteins that had been missed. Tiling data had higher coverage than shotgun proteomics and hence led to most of the corrections, but many errors probably remain. Our data are available at http://genomics.lbl.gov/supplemental/DvHtranscripts2011/.

The Enhancer of split complex arose prior to the diversification of schizophoran flies and is strongly conserved between Drosophila and stalk-eyed flies (Diopsidae).

BACKGROUND: In Drosophila, the Enhancer of split complex (E(spl)-C) comprises 11 bHLH and Bearded genes that function during Notch signaling to repress proneural identity in the developing peripheral nervous system. Comparison with other insects indicates that the basal state for Diptera is a single bHLH and Bearded homolog and that the expansion of the gene complex occurred in the lineage leading to Drosophila. However, comparative genomic data from other fly species that would elucidate the origin and sequence of gene duplication for the complex is lacking. Therefore, in order to examine the evolutionary history of the complex within Diptera, we reconstructed, using several fosmid clones, the entire E(spl)-complex in the stalk-eyed fly, Teleopsis dalmanni and collected additional homologs of E(spl)-C genes from searches of dipteran EST databases and the Glossina morsitans genome assembly. RESULTS: Comparison of the Teleopsis E(spl)-C gene organization with Drosophila indicates complete conservation in gene number and orientation between the species except that T. dalmanni contains a duplicated copy of E(spl)m5 that is not present in Drosophila. Phylogenetic analysis of E(spl)-complex bHLH and Bearded genes for several dipteran species clearly demonstrates that all members of the complex were present prior to the diversification of schizophoran flies. Comparison of upstream regulatory elements and 3' UTR domains between the species also reveals strong conservation for many of the genes and identifies several novel characteristics of E(spl)-C regulatory evolution including the discovery of a previously unidentified, highly conserved SPS+A domain between E(spl)mγ and E(spl)mß. CONCLUSION: Identifying the phylogenetic origin of E(spl)-C genes and their associated regulatory DNA is essential to understanding the functional significance of this well-studied gene complex. Results from this study provide numerous insights into the evolutionary history of the complex and will help refine the focus of studies examining the adaptive consequences of this gene expansion.

Genome-wide analyses of Geraniaceae plastid DNA reveal unprecedented patterns of increased nucleotide substitutions.

Angiosperm plastid genomes are generally conserved in gene content and order with rates of nucleotide substitutions for protein-coding genes lower than for nuclear protein-coding genes. A few groups have experienced genomic change, and extreme changes in gene content and order are found within the flowering plant family Geraniaceae. The complete plastid genome sequence of Pelargonium X hortorum (Geraniaceae) reveals the largest and most rearranged plastid genome identified to date. Highly elevated rates of sequence evolution in Geraniaceae mitochondrial genomes have been reported, but rates in Geraniaceae plastid genomes have not been characterized. Analysis of nucleotide substitution rates for 72 plastid genes for 47 angiosperm taxa, including nine Geraniaceae, show that values of dN are highly accelerated in ribosomal protein and RNA polymerase genes throughout the family. Furthermore, dN/dS is significantly elevated in the same two classes of plastid genes as well as in ATPase genes. A relatively high dN/dS ratio could be interpreted as evidence of two phenomena, namely positive or relaxed selection, neither of which is consistent with our current understanding of plastid genome evolution in photosynthetic plants. These analyses are the first to use protein-coding sequences from complete plastid genomes to characterize rates and patterns of sequence evolution for a broad sampling of photosynthetic angiosperms, and they reveal unprecedented accumulation of nucleotide substitutions in Geraniaceae. To explain these remarkable substitution patterns in the highly rearranged Geraniaceae plastid genomes, we propose a model of aberrant DNA repair coupled with altered gene expression.

A Simple, Cost-Effective, and Automation-Friendly Direct PCR Approach for Bacterial Community Analysis.

Bacterial communities in water, soil, and humans play an essential role in environmental ecology and human health. PCR-based amplicon analysis, such as 16S rRNA sequencing, is a fundamental tool for quantifying and studying microbial composition, dynamics, and interactions. However, given the complexity of microbial communities, a substantial number of samples becomes necessary for analyses that parse the factors that determine microbial composition. A common bottleneck in performing these kinds of experiments is genomic DNA (gDNA) extraction, which is time-consuming, expensive, and often biased based on the types of species present. Direct PCR method is a potentially simpler and more accurate alternative to gDNA extraction methods that do not require the intervening purification step. In this study, we evaluated three variations of direct PCR methods using diverse heterogeneous bacterial cultures, including both Gram-positive and Gram-negative species, ZymoBIOMICS microbial community standards, and groundwater. By comparing direct PCR methods with DNeasy Blood and Tissue Kits for microbial isolates and DNeasy PowerSoil Kits for microbial communities, we found that a specific variant of the direct PCR method exhibits an overall efficiency comparable to that of the conventional DNeasy PowerSoil protocol in the circumstances we tested. We also found that the method showed higher efficiency for extracting gDNA from the Gram-negative strains compared to DNeasy Blood and Tissue protocol. This direct PCR method is 1,600 times less expensive ($0.34 for 96 samples) and 10 times simpler (15 min hands-on time for 96 samples) than the DNeasy PowerSoil protocol. The direct PCR method can also be fully automated and is compatible with small-volume samples, thereby permitting scaling of samples and replicates needed to support high-throughput large-scale bacterial community analysis. IMPORTANCE Understanding bacterial interactions and assembly in complex microbial communities using 16S rRNA sequencing normally requires a large experimental load. However, the current DNA extraction methods, including cell disruption and genomic DNA purification, are normally biased, costly, time-consuming, labor-intensive, and not amenable to miniaturization by droplets or 1,536-well plates due to the significant DNA loss during the purification step for tiny-volume and low-cell-density samples. A direct PCR method could potentially solve these problems. In this study, we developed a direct PCR method which exhibits similar efficiency as the widely used method, the DNeasy PowerSoil protocol, while being 1,600 times less expensive and 10 times faster to execute. This simple, cost-effective, and automation-friendly direct-PCR-based 16S rRNA sequencing method allows us to study the dynamics, microbial interaction, and assembly of various microbial communities in a high-throughput fashion.

Deletion Mutants, Archived Transposon Library, and Tagged Protein Constructs of the Model Sulfate-Reducing Bacterium Desulfovibrio vulgaris Hildenborough.

The dissimilatory sulfate-reducing deltaproteobacterium Desulfovibrio vulgaris Hildenborough (ATCC 29579) was chosen by the research collaboration ENIGMA to explore tools and protocols for bringing this anaerobe to model status. Here, we describe a collection of genetic constructs generated by ENIGMA that are available to the research community.

Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages.

BACKGROUND: Despite considerable progress in our understanding of land plant phylogeny, several nodes in the green tree of life remain poorly resolved. Furthermore, the bulk of currently available data come from only a subset of major land plant clades. Here we examine early land plant evolution using complete plastome sequences including two previously unexamined and phylogenetically critical lineages. To better understand the evolution of land plants and their plastomes, we examined aligned nucleotide sequences, indels, gene and nucleotide composition, inversions, and gene order at the boundaries of the inverted repeats. RESULTS: We present the plastome sequences of Equisetum arvense, a horsetail, and of Isoetes flaccida, a heterosporous lycophyte. Phylogenetic analysis of aligned nucleotides from 49 plastome genes from 43 taxa supported monophyly for the following clades: embryophytes (land plants), lycophytes, monilophytes (leptosporangiate ferns + Angiopteris evecta + Psilotum nudum + Equisetum arvense), and seed plants. Resolution among the four monilophyte lineages remained moderate, although nucleotide analyses suggested that P. nudum and E. arvense form a clade sister to A. evecta + leptosporangiate ferns. Results from phylogenetic analyses of nucleotides were consistent with the distribution of plastome gene rearrangements and with analysis of sequence gaps resulting from insertions and deletions (indels). We found one new indel and an inversion of a block of genes that unites the monilophytes. CONCLUSIONS: Monophyly of monilophytes has been disputed on the basis of morphological and fossil evidence. In the context of a broad sampling of land plant data we find several new pieces of evidence for monilophyte monophyly. Results from this study demonstrate resolution among the four monilophytes lineages, albeit with moderate support; we posit a clade consisting of Equisetaceae and Psilotaceae that is sister to the "true ferns," including Marattiaceae.

Chloroplast genome sequence of the moss Tortula ruralis: gene content, polymorphism, and structural arrangement relative to other green plant chloroplast genomes.

BACKGROUND: Tortula ruralis, a widely distributed species in the moss family Pottiaceae, is increasingly used as a model organism for the study of desiccation tolerance and mechanisms of cellular repair. In this paper, we present the chloroplast genome sequence of T. ruralis, only the second published chloroplast genome for a moss, and the first for a vegetatively desiccation-tolerant plant. RESULTS: The Tortula chloroplast genome is approximately 123,500 bp, and differs in a number of ways from that of Physcomitrella patens, the first published moss chloroplast genome. For example, Tortula lacks the approximately 71 kb inversion found in the large single copy region of the Physcomitrella genome and other members of the Funariales. Also, the Tortula chloroplast genome lacks petN, a gene found in all known land plant plastid genomes. In addition, an unusual case of nucleotide polymorphism was discovered. CONCLUSIONS: Although the chloroplast genome of Tortula ruralis differs from that of the only other sequenced moss, Physcomitrella patens, we have yet to determine the biological significance of the differences. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for mosses) of the generation of DNA markers for fine-level phylogenetic studies, or to investigate individual variation within populations.