Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes

BACKGROUND: Lateral gene transfer (LGT) is an important evolutionary mechanism among prokaryotes. The situation in eukaryotes is less clear; the human genome sequence failed to give strong support for any recent transfers from prokaryotes to vertebrates, yet a number of LGTs from prokaryotes to protists (unicellular eukaryotes) have been documented. Here, we perform a systematic analysis to investigate the impact of LGT on the evolution of diplomonads, a group of anaerobic protists.RESULTS: Phylogenetic analyses of 15 genes present in the genome of the Atlantic Salmon parasite Spironucleus barkhanus and/or the intestinal parasite Giardia lamblia show that most of these genes originated via LGT. Half of the genes are putatively involved in processes related to an anaerobic lifestyle, and this finding suggests that a common ancestor, which most probably was aerobic, of Spironucleus and Giardia adapted to an anaerobic environment in part by acquiring genes via LGT from prokaryotes. The sources of the transferred diplomonad genes are found among all three domains of life, including other eukaryotes. Many of the phylogenetic reconstructions show eukaryotes emerging in several distinct regions of the tree, strongly suggesting that LGT not only involved diplomonads, but also involved other eukaryotic groups.CONCLUSIONS: Our study shows that LGT is a significant evolutionary mechanism among diplomonads in particular and protists in general. These findings provide insights into the evolution of biochemical pathways in early eukaryote evolution and have important implications for studies of eukaryotic genome evolution and organismal relationships. Furthermore, "fusion" hypotheses for the origin of eukaryotes need to be rigorously reexamined in the light of these results.     

Recent events dominate interdomain lateral gene transfers between prokaryotes and eukaryotes and,with the exception of endosymbiotic gene transfers,few ancient transfer events persist

While there is compelling evidence for the impact of endosymbiotic gene transfer (EGT; transfer from either mitochondrion or chloroplast to the nucleus) on genome evolution in eukaryotes, the role of interdomain transfer from bacteria and/or archaea (i.e. prokaryotes) is less clear. Lateral gene transfers (LGTs) have been argued to be potential sources of phylogenetic information, particularly for reconstructing deep nodes that are difficult to recover with traditional phylogenetic methods. We sought to identify interdomain LGTs by using a phylogenomic pipeline that generated 13 465 single gene trees and included up to 487 eukaryotes, 303 bacteria and 118 archaea. Our goals include searching for LGTs that unite major eukaryotic clades, and describing the relative contributions of LGT and EGT across the eukaryotic tree of life. Given the difficulties in interpreting single gene trees that aim to capture the approximately 1.8 billion years of eukaryotic evolution, we focus on presence–absence data to identify interdomain transfer events. Specifically, we identify 1138 genes found only in prokaryotes and representatives of three or fewer major clades of eukaryotes (e.g. Amoebozoa, Archaeplastida, Excavata, Opisthokonta, SAR and orphan lineages). The majority of these genes have phylogenetic patterns that are consistent with recent interdomain LGTs and, with the notable exception of EGTs involving photosynthetic eukaryotes, we detect few ancient interdomain LGTs. These analyses suggest that LGTs have probably occurred throughout the history of eukaryotes, but that ancient events are not maintained unless they are associated with endosymbiotic gene transfer among photosynthetic lineages.     

Detection of Lateral Gene Transfer Events in the Prokaryotic tRNA Synthetases by the Ratios of Evolutionary Distances Method

The availability of large numbers of genomic sequences has demonstrated the importance of lateral gene transfer (LGT) in prokaryotic evolution. However, considerable uncertainty remains concerning the frequency of LGT compared to other evolutionary processes. To examine LGTs in ancient lineages of prokaryotes a method was developed that utilizes the ratios of evolutionary distances (RED) to distinguish between alternative evolutionary histories. The advantages of this approach are that the variability inherent in comparing protein sequences is transparent, the direction of LGT and the relative rates of evolution are readily identified, and it is possible to detect other types of evolutionary events. This method was standardized using 35 genes encoding ribosomal proteins that were believed to share a vertical evolution. Using RED-T, an original computer program designed to implement the RED method, the evolution of the genes encoding the 20 aminoacyl-tRNA synthetases was examined. Although LGTs were common in the evolution of the aminoacyl-tRNA synthetases, they were not sufficient to obscure the organismal phylogeny. Moreover, much of the apparent complexity of the gene tree was consistent with the formation of the paralogs in the ancestors to the modern lineages followed by more recent loss of one paralog or the other.     

Evolution of plant Ash1 SET genes: structural divergence and functional differentiation

Plant Ash1 SET proteins are involved in H3K36 methylation, and play a key role in plant reproductive development. Genes encoding Ash1 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 and functional differentiation of the Ash1 SET gene family, we made a comprehensive evolutionary analysis of this gene family from eleven major representatives of green plants. A novel deep sister relationship grouping previously resolved II-1 and II-2 orthologous groups was identified. The absence of AWS domain in the group II-2 suggests that the independent losses of AWS domain have occurred during evolution. A diversity of gene structures in plant Ash1 SET gene family have been presented since the divergence of Physcomitrella patens (moss) from the other land plants. A small proportion of codons in SET domain regions were detected to be under positive selection along the branches ancestral to land plant and angiosperms, which may have allowed changes of substrate specificity among different evolutionary groups while maintaining the primary function of SET domains. Our predictive subcellular localization and comparative anatomical meta-expression analyses can assort with the structural divergences of Ash1 SET proteins.     

Phylogeny and Molecular Evolution of the Green Algae

The green lineage (Viridiplantae) comprises the green algae and their descendants the land plants, and is one of the major groups of oxygenic photosynthetic eukaryotes. Current hypotheses posit the early divergence of two discrete clades from an ancestral green flagellate. One clade, the Chlorophyta, comprises the early diverging prasinophytes, which gave rise to the core chlorophytes. The other clade, the Streptophyta, includes the charophyte green algae from which the land plants evolved. Multi-marker and genome scale phylogenetic studies have greatly improved our understanding of broad-scale relationships of the green lineage, yet many questions persist, including the branching orders of the prasinophyte lineages, the relationships among core chlorophyte clades (Chlorodendrophyceae, Ulvophyceae, Trebouxiophyceae and Chlorophyceae), and the relationships among the streptophytes. Current phylogenetic hypotheses provide an evolutionary framework for molecular evolutionary studies and comparative genomics. This review summarizes our current understanding of organelle genome evolution in the green algae, genomic insights into the ecology of oceanic picoplanktonic prasinophytes, molecular mechanisms underlying the evolution of complexity in volvocine green algae, and the evolution of genetic codes and the translational apparatus in green seaweeds. Finally, we discuss molecular evolution in the streptophyte lineage, emphasizing the genetic facilitation of land plant origins.     

Assessing the effects of a sequestered germline on interdomain lateral gene transfer in Metazoa

A sequestered germline in Metazoa has been argued to be an obstacle to lateral gene transfer (LGT), though few studies have specifically assessed this claim. Here, we test the hypothesis that the origin of a sequestered germline reduced LGT events in Bilateria (i.e., triploblast lineages) as compared to early‐diverging Metazoa (i.e., Ctenophora, Cnidaria, Porifera, and Placozoa). We analyze single‐gene phylogenies generated with over 900 species sampled from among Bacteria, Archaea, and Eukaryota to identify well‐supported interdomain LGTs. We focus on ancient interdomain LGT (i.e., those between prokaryotes and multiple lineages of Metazoa) as systematic errors in single‐gene tree reconstruction create uncertainties for interpreting eukaryote‐to‐eukaryote transfer. The breadth of the sampled Metazoa enables us to estimate the timing of LGTs, and to examine the pattern before versus after the evolution of a sequestered germline. We identified 58 LGTs found only in Metazoa and prokaryotes (i.e., bacteria and/or archaea), and seven genes transferred from prokaryotes into Metazoa plus one other eukaryotic clade. Our analyses indicate that more interdomain transfers occurred before the development of a sequestered germline, consistent with the hypothesis that this feature is an obstacle to LGT.     

Complex distribution of EFL and EF-1α proteins in the green algal lineage

Background  

EFL (or elongation factor-like) is a member of the translation superfamily of GTPase proteins. It is restricted to eukaryotes, where it is found in a punctate distribution that is almost mutually exclusive with elongation factor-1 alpha (EF-1α). EF-1α is a core translation factor previously thought to be essential in eukaryotes, so its relationship to EFL has prompted the suggestion that EFL has spread by horizontal or lateral gene transfer (HGT or LGT) and replaced EF-1α multiple times. Among green algae, trebouxiophyceans and chlorophyceans have EFL, but the ulvophycean Acetabularia and the sister group to green algae, land plants, have EF-1α. This distribution singles out green algae as a particularly promising group to understand the origin of EFL and the effects of its presence on EF-1α.     

Green algae to land plants: An evolutionary transition

Studies focused upon the evolutionary transition from ancestral green algae to the earliest land plants are important from a range of ecological, molecular and evolutionary perspectives. A substantial suite of ultrastructural, biochemical and molecular data supports the concept that land plants (embryophytes) are monophyletically derived from an ancestral charophycean alga. However, the details of phylogenetic branching patterns linking extant charophytes and seedless embryophytes are currently unclear. Moreover, the fossil record has so far been mute regarding the algae-land plant transition. Nevertheless, an accurate reflection of major evolutionary events in the history of the earliest land plants can be obtained by comparative paleontological-neontological studies, and comparative molecular, cellular and developmental investigations of extant charophytes and bryophytes. This review focuses upon research progress toward understanding three clade-specific adaptations that were important in the successful colonization of land by plants: the histogenetic apical meristem, the matrotrophic embryo, and decay-resistant cell wall polymers.     

Evolution of fragmented mitochondrial ribosomal RNA genes inChlamydomonas

The fragmented mitochondrial ribosomal RNAs (rRNAs) of the green algaeChlamydomonas eugametos andChlamydomonas reinhardtii are discontinuously encoded in subgenic modules that are scrambled in order and interspersed with protein coding and tRNA genes. The mitochondrial rRNA genes of these two algae differ, however, in both the distribution and organization of rRNA coding information within their respective genomes. The objectives of this study were (1) to examine the phylogenetic relationships between the mitochondrial rRNA gene sequences ofC. eugametos andC. reinhardtii and those of the conventional mitochondrial rRNA genes of the green alga,Prototheca wickerhamii, and land plants and (2) to attempt to deduce the evolutionary pathways that gave rise to the unusual mitochondrial rRNA gene structures in the genusChlamydomonas. Although phylogenetic analysis revealed an affiliation between the mitochondrial rRNA gene sequences of the twoChlamydomonas taxa to the exclusion of all other mitochondrial rRNA gene sequences tested, no specific affiliation was noted between theChlamydomonas sequences andP. wickerhamii or land plants. Calculations of the minimal number of transpositions required to convert hypothetical ancestral rRNA gene organizations to the arrangements observed forC. eugametos andC. reinhardtii mitochondrial rRNA genes, as well as a limited survey of the size of mitochondrial rRNAs in other members of the genus, lead us to propose that the last common ancestor ofChlamydomonas algae contained fragmented mitochondrial rRNA genes that were nearly co-linear with conventional rRNA genes.     

Presence of a Bacterial-Like Citrate Synthase Gene in <Emphasis Type="Italic">Tetrahymena thermophila</Emphasis>: Recent Lateral Gene Transfers (LGT) or Multiple Gene Losses Subsequent to a Single Ancient LGT?

Citrate synthase is the initial enzyme in the tricarboxylic acid cycle of mitochondria. In plants and fungi, it is the second isozyme in the glyoxylate cycle of peroxisomes (or glyoxysomes), and it is also present in bacteria. Some of the biochemical reactions in the glyoxylate cycle of the ciliated protozoan Tetrahymena pyriformis depend upon mitochondrial enzymes, as T. pyriformis lacks some glyoxysome-specific enzymes. Here we demonstrate a new citrate synthase gene from Tetrahymena thermophila that is different from the mitochondrial counterpart. A potential peroxysome-targeted signal was detected in the N-terminus, suggesting the localization of the enzyme in peroxysomes. Phylogenetic analysis placed the Tetrahymena sequence in a clade consisting of a few sequences from eukaryotes such as cellular slime molds and two land plants, near a green sulfur bacterium and many proteobacteria as a sister group but not in a mitochondrial clade. Southern blot analysis revealed that this type of gene was absent from distantly related ciliates and other species of Tetrahymena except for the closest species, T. mallaccensis. The scattered presence of the bacterial-like genes among distantly related eukaryotes suggests three alternative interpretations of acquisition of the novel glyoxysomal citrate synthase gene via lateral gene transfer (LGT). (1) Some eukaryotes independently acquired the gene from a common bacterium or closely related bacteria via LGT. (2) A hypothetical eukaryote once acquired the gene, which was thereafter independently transferred from the eukaryote to other eukaryotes. (3) A single event of LGT (or duplication) occurred in a certain common ancestor of eukaryotes, followed by multiple losses in many eukaryotic lineages during the subsequent evolution. Considering the monophyly of the bacterial-like eukaryotic citrate synthase genes, the first model is somewhat unlikely, even though it is not impossible. The second and third models can rationally explain the present observation, so these models are discused in some detail.     

The Complete Mitochondrial Genome Sequence of the Hornwort <Emphasis Type="Italic">Megaceros</Emphasis><Emphasis Type="Italic">aenigmaticus</Emphasis> Shows a Mixed Mode of Conservative Yet Dynamic Evolution in Early Land Plant Mitochondrial Genomes

Land plants possess some of the most unusual mitochondrial genomes among eukaryotes. However, in early land plants these genomes resemble those of green and red algae or early eukaryotes. The question of when during land plant evolution the dramatic change in mtDNAs occurred remains unanswered. Here we report the first completely sequenced mitochondrial genome of the hornwort, Megaceros aenigmaticus, a member of the sister group of vascular plants. It is a circular molecule of 184,908 base pairs, with 32 protein genes, 3 rRNA genes, 17 tRNA genes, and 30 group II introns. The genome contains many genes arranged in the same order as in those of a liverwort, a moss, several green and red algae, and Reclinomonas americana, an early-branching eukaryote with the most ancestral form of mtDNA. In particular, the gene order between mtDNAs of the hornwort and Physcomitrella patens (moss) differs by only 8 inversions and translocations. However, the hornwort mtDNA possesses 4 derived features relative to green alga mtDNAs—increased genome size, RNA editing, intron gains, and gene losses—which were all likely acquired during the origin and early evolution of land plants. Overall, this genome and those of other 2 bryophytes show that mitochondrial genomes in early land plants, unlike their seed plant counterparts, exhibit a mixed mode of conservative yet dynamic evolution. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Libo Li and Bin Wang contributed equally to this work.     

Evolutionary origins of Hsp90 chaperones and a deep paralogy in their bacterial ancestors

The 82-90 kD family of molecular chaperone proteins has homologs in eukaryotes (Hsp90) and many eubacteria (HtpG) but not in Archaebacteria. We used representatives of all four different eukaryotic paralogs (cytosolic, endoplasmic reticulum (ER), chloroplast, mitochondrial) together with numerous eubacterial HtpG proteins for phylogenetic analyses to investigate their evolutionary origins. Our trees confirm that none of the organellar Hsp90s derives from the endosymbionts of early eukaryotes. Contrary to previous suggestions of distant origins through lateral gene transfer (LGT) all eukaryote Hsp90s are related to Gram-positive eubacterial HtpG proteins. The nucleocytosolic, ER and chloroplast Hsp90 paralogs are clearly mutually related. The origin of mitochondrial Hsp90 is more obscure, as these sequences are deeply nested within eubacteria. Our trees also reveal a deep split within eubacteria into a group of mainly long-branching sequences (including the eukaryote mitochondrial Hsp90s) and another group comprising exclusively short-branching HtpG proteins, from which the cytosolic/ER versions probably arose. Both versions are present in several eubacterial phyla, suggesting gene duplication very early in eubacterial evolution and multiple independent losses thereafter. We identified one probable case of LGT within eubacteria. However, multiple losses can simply explain the evolutionary pattern of the eubacterial HtpG paralogs and predominate over LGT. We suggest that the actinobacterial ancestor of eukaryotes harbored genes for both eubacterial HtpG paralogs, as the actinobacterium Streptomyces coelicolor still does; one could have given rise to the mitochondrial Hsp90 and the other, following another duplication event in the ancestral eukaryote, to the cytosolic and ER Hsp90 homologs.     

Ca2+ signalling in plants and green algae--changing channels

Eukaryotic cells generate cytosolic Ca(2+) signals via Ca(2+)-conducting channels in cellular membranes. Plants and animals exhibit substantial differences in their complement of Ca(2+) channels. In particular, the four-domain voltage-dependent Ca(2+) channels, transient receptor potential channels and inositol (1,4,5)-trisphosphate receptors, which have important roles in animal physiology, are all absent from land plants. Recent evidence from biochemical and genomic studies has indicated that representatives of these classes of Ca(2+) channels are present in members of the green plant lineage, the chlorophyte algae. This indicates that the Ca(2+)-signalling mechanisms absent from land plants were, in fact, present in ancestral eukaryotes and were lost by land plants after their divergence from the chlorophyte algae.     

The rate and pattern of cladogenesis in microbes

Theories of macroevolution rarely have been extended to include microbes; however, because microbes represent the most ancient and diverse assemblage of organismal diversity, such oversight limits our understanding of evolutionary history. Our analysis of phylogenetic trees for microbes suggests that macroevolution may differ between prokaryotes and both micro- and macroeukaryotes (mainly plants and animals). Phylogenetic trees inferred for prokaryotes and some microbial eukaryotes conformed to expectations assuming a constant rate of cladogenesis over time and among lineages: nevertheless, microbial eukaryote trees exhibited more variation in rates of cladogenesis than prokaryote trees. We hypothesize that the contrast of macroevolutionary dynamics between prokaryotes and many eukaryotes is due, at least in part, to differences in the prevalence of lateral gene transfer (LGT) between the two groups. Inheritance is predominantly, if not wholly, vertical within eukaryotes, a feature that allows for the emergence and maintenance of heritable variation among lineages. By contrast, frequent LGT in prokaryotes may ameliorate heritable variation in rate of cladogenesis resulting from the emergence of key innovations; thus, the inferred difference in macroevolution might reflect exclusivity of key innovations in eukaryotes and their promiscuous nature in prokaryotes.     

The evolutionary history of calreticulin and calnexin genes in green plants

Calreticulin and calnexin are Ca²⁺-binding chaperones localized in the endoplasmic reticulum of eukaryotes acting in glycoprotein folding quality control and Ca²⁺ homeostasis. The evolutionary histories of calreticulin and calnexin gene families were inferred by comprehensive phylogenetic analyses using 18 completed genomes and ESTs covering the major green plants groups, from green algae to angiosperms. Calreticulin and calnexin possibly share a common origin, and both proteins are present along all green plants lineages. The calreticulin founder gene within green plants duplicated in early tracheophytes leading to two possible groups of orthologs with specialized functions, followed by lineage-specific gene duplications in spermatophytes. Calnexin founder gene in land plants was inherited from basal green algae during evolution in a very conservative copy number. A comprehensive classification in possible groups of orthologs and a catalog of calreticulin and calnexin genes from green plants are provided.     

Transfer of energy pathway genes in microbial enhanced biological phosphorus removal communities

Background

Lateral gene transfer (LGT) is an important evolutionary process in microbial evolution. In sewage treatment plants, LGT of antibiotic resistance and xenobiotic degradation-related proteins has been suggested, but the role of LGT outside these processes is unknown. Microbial communities involved in Enhanced Biological Phosphorus Removal (EBPR) have been used to treat wastewater in the last 50 years and may provide insights into adaptation to an engineered environment. We introduce two different types of analysis to identify LGT in EBPR sewage communities, based on identifying assembled sequences with more than one strong taxonomic match, and on unusual phylogenetic patterns. We applied these methods to investigate the role of LGT in six energy-related metabolic pathways.

Results

The analyses identified overlapping but non-identical sets of transferred enzymes. All of these were homologous with sequences from known mobile genetic elements, and many were also in close proximity to transposases and integrases in the EBPR data set. The taxonomic method had higher sensitivity than the phylogenetic method, identifying more potential LGTs. Both analyses identified the putative transfer of five enzymes within an Australian community, two in a Danish community, and none in a US-derived culture.

Conclusions

Our methods were able to identify sequences with unusual phylogenetic or compositional properties as candidate LGT events. The association of these candidates with known mobile elements supports the hypothesis of transfer. The results of our analysis strongly suggest that LGT has influenced the development of functionally important energy-related pathways in EBPR systems, but transfers may be unique to each community due to different operating conditions or taxonomic composition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1752-5) contains supplementary material, which is available to authorized users.     

The complete mitochondrial DNA sequences of Nephroselmis olivacea and Pedinomonas minor. Two radically different evolutionary patterns within green algae.

Green plants appear to comprise two sister lineages, Chlorophyta (classes Chlorophyceae, Ulvophyceae, Trebouxiophyceae, and Prasinophyceae) and Streptophyta (Charophyceae and Embryophyta, or land plants). To gain insight into the nature of the ancestral green plant mitochondrial genome, we have sequenced the mitochondrial DNAs (mtDNAs) of Nephroselmis olivacea and Pedinomonas minor. These two green algae are presumptive members of the Prasinophyceae. This class is thought to include descendants of the earliest diverging green algae. We find that Nephroselmis and Pedinomonas mtDNAs differ markedly in size, gene content, and gene organization. Of the green algal mtDNAs sequenced so far, that of Nephroselmis (45,223 bp) is the most ancestral (minimally diverged) and occupies the phylogenetically most basal position within the Chlorophyta. Its repertoire of 69 genes closely resembles that in the mtDNA of Prototheca wickerhamii, a later diverging trebouxiophycean green alga. Three of the Nephroselmis genes (nad10, rpl14, and rnpB) have not been identified in previously sequenced mtDNAs of green algae and land plants. In contrast, the 25,137-bp Pedinomonas mtDNA contains only 22 genes and retains few recognizably ancestral features. In several respects, including gene content and rate of sequence divergence, Pedinomonas mtDNA resembles the reduced mtDNAs of chlamydomonad algae, with which it is robustly affiliated in phylogenetic analyses. Our results confirm the existence of two radically different patterns of mitochondrial genome evolution within the green algae.     

植物颗粒结合淀粉合酶GBSS基因家族的进化

为全面理解植物颗粒结合淀粉合酶(GBSS)基因在植物中的进化模式并重建其进化历史, 利用20种陆生植物和2种藻类植物的基因组数据, 通过生物信息学手段, 深入挖掘和分析植物类群基因组中GBSS基因家族的构成和基因特点, 推测其可能的扩增和丢失规律。结果共识别42条同源序列。系统发育和进化分析表明, GBSS基因起源古老, 可能在所有绿色植物的祖先中就已经出现, 之后在进化过程中不断发生谱系的特异扩张和拷贝丢失, 并最终通过功能分化的形式在植物类群中被固定。     

Characterization of an Ancient Lepidopteran Lateral Gene Transfer

Bacteria to eukaryote lateral gene transfers (LGT) are an important potential source of material for the evolution of novel genetic traits. The explosion in the number of newly sequenced genomes provides opportunities to identify and characterize examples of these lateral gene transfer events, and to assess their role in the evolution of new genes. In this paper, we describe an ancient lepidopteran LGT of a glycosyl hydrolase family 31 gene (GH31) from an Enterococcus bacteria. PCR amplification between the LGT and a flanking insect gene confirmed that the GH31 was integrated into the Bombyx mori genome and was not a result of an assembly error. Database searches in combination with degenerate PCR on a panel of 7 lepidopteran families confirmed that the GH31 LGT event occurred deep within the Order approximately 65–145 million years ago. The most basal species in which the LGT was found is Plutella xylostella (superfamily: Yponomeutoidea). Array data from Bombyx mori shows that GH31 is expressed, and low dN/dS ratios indicates the LGT coding sequence is under strong stabilizing selection. These findings provide further support for the proposition that bacterial LGTs are relatively common in insects and likely to be an underappreciated source of adaptive genetic material.     

A Multistep Process Gave Rise to RNA Polymerase IV of Land Plants

Since their discovery in Metazoa, the three nuclear RNA polymerases (RNAPs) have been found in fungi, plants, and diverse protists. In all eukaryotes studied to date, RNAPs I, II, and III collectively transcribe all major RNAs made in the nucleus. We have found genes for the largest subunit (RPD1/RPE1) of a new DNA-dependent RNAP, RNAP IV, in all major land plant taxa and in closely related green algae. Genes for the second-largest subunit (RPD2) of this enzyme were found in all land plants. Phylogenetic study indicates that RNAP IV genes are sister to the corresponding RNAP II genes. Our results show the genesis of RNAP IV to be a multistep process in which the largest and second-largest subunit genes evolved by independent duplication events in the ancestors of Charales and land plants. These findings provide insights into evolutionary mechanisms that can explain the origin of multiple RNAPs in the eukaryotic nucleus. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Patrick Keeling]