Arnica (Asteraceae) phylogeny revisited using RPB2: complex patterns and multiple d-paralogues

The region coding for the second largest subunit of RNA polymerase II (RPB2) was explored for resolving interspecific relationships in Arnica and lower level taxa in general. The region between exons 17 and 23 was cloned and sequenced for 33 accessions of Arnica and four outgroup taxa. Three paralogues of the RPB2-d copy (RPB2-dA, B and C) were detected in Arnica and outgroup taxa, indicating that the duplications must have occurred before the divergence of Arnica. Parsimony and Bayesian analyses of separate alignments of the three copies reveal complex patterns in Arnica, likely reflecting a history of lineage sorting in combination with apomixis, polyploidization, and possibly hybridization. Cloned sequences of some taxa do not form monophyletic clades within paralogues, but form multiple strongly supported clades with sequences of other taxa. Some well supported groups are present in more than one paralogue and many groups are in line with earlier hypotheses regarding interspecific relationships within the genus. Low levels of homoplasy in combination with relatively high sequence variation indicates that the introns of the RPB2 region could be suitable for phylogenetic studies in low level taxonomy.     

Phylogenetic relationships among early-diverging eudicots based on four genes: were the eudicots ancestrally woody?

Based on analyses of combined data sets of three genes (18S rDNA, rbcL, and atpB), phylogenetic relationships among the early-diverging eudicot lineages (Ranunculales, Proteales, Trochodendraceae, Sabiaceae, and Buxaceae) remain unclear, as are relationships within Ranunculales, especially the placement of Eupteleaceae. To clarify relationships among these early-diverging eudicot lineages, we added entire sequences of 26S rDNA to the existing three-gene data set. In the combined analyses of four genes based on parsimony, ML, and Bayesian analysis, Ranunculales are strongly supported as a clade and are sister to other eudicots. Proteales appear as sister to the remaining eudicots, which are weakly (59%) supported as a clade. Relationships among Trochodendraceae, Buxaceae (including Didymeles), Sabiaceae, and Proteales remain unclear. Within Ranunculales, Eupteleaceae are sister to all other Ranunculales, with bootstrap support of 70% in parsimony analysis and with posterior probability of 1.00 in Bayesian analysis. Our character reconstructions indicate that the woody habit is ancestral, not only for the basal angiosperms, but also for the eudicots. Furthermore, Ranunculales may not be ancestrally herbaceous, as long maintained. The woody habit appears to have been ancestral for several major clades of eudicots, including Caryophyllales, and asterids.     

Conservation and divergence in the AGAMOUS subfamily of MADS-box genes: evidence of independent sub- and neofunctionalization events

The MADS-box gene AGAMOUS (AG) plays a key role in determining floral meristem and organ identities. We identified three AG homologs, EScaAG1, EScaAG2, and EScaAGL11 from the basal eudicot Eschscholzia californica (California poppy). Phylogenetic analyses indicate that EScaAG1 and EScaAG2 are recent paralogs within the AG clade, independent of the duplication in ancestral core eudicots that gave rise to the euAG and PLENA (PLE) orthologs. EScaAGL11 is basal to core eudicot AGL11 orthologs in a clade representing an older duplication event after the divergence of the angiosperm and gymnosperm lineages. Detailed in situ hybridization experiments show that expression of EScaAG1 and EScaAG2 is similar to AG; however, both genes appear to be expressed earlier in floral development than described in the core eudicots. A thorough examination of available expression and functional data in a phylogenetic context for members of the AG and AGL11 clades reveals that gene expression has been quite variable throughout the evolutionary history of the AG subfamily and that ovule-specific expression might have evolved more than twice. Although sub- and neofunctionalization are inferred to have occurred following gene duplication, functional divergence among orthologs is evident, as is convergence, among paralogs sampled from different species. We propose that retention of multiple AG homologs in several paralogous lineages can be explained by the conservation of ancestral protein activity combined with evolutionarily labile regulation of expression in the AG and AGL11 clades such that the collective functions of the AG subfamily in stamen and carpel development are maintained following gene duplication.     

Comparative sequence analyses indicate that Coffea (Asterids) and Vitis (Rosids) derive from the same paleo-hexaploid ancestral genome

The complete sequence of Vitis vinifera revealed that the rosid clade derives from a hexaploid ancestor. At present, no analysis of complete genome sequence is available for an asterid, the other large eudicot clade, which includes the economically important species potato, tomato and coffee. To elucidate the genomic history of asterids, we compared the sequence of an 800 kb region of diploid Coffea genome to the orthologous regions of V. vinifera, Populus trichocarpa and Arabidopsis thaliana. We found a very high level of collinearity between around 80 genes of the three rosid species and Coffea. Collinearity comparisons between orthologous and paralogous regions indicates that (1) the Coffea (and consequently all asterids) and rosids share the same hexaploid ancestor; (2) the diploidization process (loss of duplicated and redundant copies from the whole genome duplication) was very advanced in the most recent common ancestor of rosids and asterids. Finally, no additional polyploidization events were detected in the Coffea lineage. Differences in gene loss rates were detected among the three rosid species and linked to the divergence in protein sequences.     

Evolution of the closely related, sex-related genes DM-W and DMRT1 in African clawed frogs (Xenopus)

DM-W is a dominant, female-specific, regulator of sex determination in the African clawed frog Xenopus laevis. This gene is derived from partial duplication of DMRT1, a male-related autosomal gene. We set out to better understand sex determination in Xenopus by studying this pair of genes. We found that DM-W evolved in Xenopus after divergence from the sister genus Silurana but before divergence of X. laevis and X. clivii, and that DM-W arose from partial duplication of DMRT1β, which is one of the two DMRT1 paralogs in the tetraploid ancestor of Xenopus. Using the rate ratio of nonsynonymous to synonymous substitutions per site and multilocus polymorphism data, we show that DM-W evolved non-neutrally. By cloning paralogs and using a pyrosequencing assay, we also demonstrate that DMRT1 underwent phylogenetically biased pseudogenization after polyploidization, and that expression of this gene is regulated by mechanisms that vary through development. One explanation for these observations is that the expression domain of DMRT1β was marginalized, which would explain why this paralog is dispensable in Xenopus polyploids and why DM-W has a narrow expression domain. These findings illustrate how evolution of the genetic control of stable phenotypes is facilitated by redundancy, degeneration, and compartmentalized regulation.     

Evolution and diversity of axon guidance Robo receptor family genes

The diversity of axon guidance (AG) receptors reflects gains in complexity of the animal nervous system during evolution. Members of the Roundabout (Robo) family of receptors interact with Slit proteins and play important roles in many developmental processes, including AG and neural crest cell migration. There are four members of the Robo gene family. However, the evolutionary history of Robo family genes remain obscure. We analyzed the distribution of Robo family members in metazoan species ranging in complexity from hydras to humans. We undertook a phylogenetic analysis in metazoans, synteny analysis, and ancestral chromosome mapping in vertebrates, and detected selection pressure and functional divergence among four mammalian Robo paralogs. Based on our analysis, we proposed that the ancestral Robo gene could have undergone a tandem duplication in the vertebrate ancestor; then one round of whole genome duplication events occurred before the divergence of ancestral lamprey and gnathostome, generating four paralogs in early vertebrates. Robo4 paralog underwent segmental loss in the following evolutionary process. Our results showed that Robo3 paralog is under more powerful purifying selection pressure compared with other three paralogs, which could correlate with its unique expression pattern and function. Furthermore, we found four sites under positive selection pressure on the Ig1‐2 domains of Robo4 that might interfere with its binding to Slits ligand. Diverge analysis at the amino acid level showed that Robo4 paralog have relatively greater functional diversifications than other Robo paralogs. This coincides with the fact that Robo4 predominantly functions in vascular endothelial cells but not the nervous system.     

Retention of protein complex membership by ancient duplicated gene products in budding yeast

To investigate functional divergence of gene duplicates, we examined the protein-protein interactions and coexistence in complexes of paralogs resulting from an ancient whole-genome duplication in yeast. Strikingly, half the surveyed paralog pairs were found to be co-clustered in protein complexes, and were more conserved and highly expressed than non-co-clustered paralogs; however, their discordant expression patterns and conservation rates indicate differential regulation of subfunctionalized paralogs. These results highlight the value of protein complex membership in studying functional divergence among gene duplicates.     

Genes and Processed Paralogs Co-exist in Plant Mitochondria

RNA-mediated gene duplication has been proposed to create processed paralogs in the plant mitochondrial genome. A processed paralog may retain signatures left by the maturation process of its RNA precursor, such as intron removal and no need of RNA editing. Whereas it is well documented that an RNA intermediary is involved in the transfer of mitochondrial genes to the nucleus, no direct evidence exists for insertion of processed paralogs in the mitochondria (i.e., processed and un-processed genes have never been found simultaneously in the mitochondrial genome). In this study, we sequenced a region of the mitochondrial gene nad1, and identified a number of taxa were two different copies of the region co-occur in the mitochondria. The two nad1 paralogs differed in their (a) presence or absence of a group II intron, and (b) number of edited sites. Thus, this work provides the first evidence of co-existence of processed paralogs and their precursors within the plant mitochondrial genome. In addition, mapping the presence/absence of the paralogs provides indirect evidence of RNA-mediated gene duplication as an essential process shaping the mitochondrial genome in plants.     

Extensive concerted evolution of rice paralogs and the road to regaining independence

Many genes duplicated by whole-genome duplications (WGDs) are more similar to one another than expected. We investigated whether concerted evolution through conversion and crossing over, well-known to affect tandem gene clusters, also affects dispersed paralogs. Genome sequences for two Oryza subspecies reveal appreciable gene conversion in the approximately 0.4 MY since their divergence, with a gradual progression toward independent evolution of older paralogs. Since divergence from subspecies indica, approximately 8% of japonica paralogs produced 5-7 MYA on chromosomes 11 and 12 have been affected by gene conversion and several reciprocal exchanges of chromosomal segments, while approximately 70-MY-old "paleologs" resulting from a genome duplication (GD) show much less conversion. Sequence similarity analysis in proximal gene clusters also suggests more conversion between younger paralogs. About 8% of paleologs may have been converted since rice-sorghum divergence approximately 41 MYA. Domain-encoding sequences are more frequently converted than nondomain sequences, suggesting a sort of circularity--that sequences conserved by selection may be further conserved by relatively frequent conversion. The higher level of concerted evolution in the 5-7 MY-old segmental duplication may reflect the behavior of many genomes within the first few million years after duplication or polyploidization.     

Evolution of plant RNA polymerase IV/V genes: evidence of subneofunctionalization of duplicated <Emphasis Type="Italic">NRPD2/NRPE2</Emphasis>-like paralogs in <Emphasis Type="Italic">Viola</Emphasis> (Violaceae)

Background  

DNA-dependent RNA polymerase IV and V (Pol IV and V) are multi-subunit enzymes occurring in plants. The origin of Pol V, specific to angiosperms, from Pol IV, which is present in all land plants, is linked to the duplication of the gene encoding the largest subunit and the subsequent subneofunctionalization of the two paralogs (NRPD1 and NRPE1). Additional duplication of the second-largest subunit, NRPD2/NRPE2, has happened independently in at least some eudicot lineages, but its paralogs are often subject to concerted evolution and gene death and little is known about their evolution nor their affinity with Pol IV and Pol V.     

Gunnerales are sister to other core eudicots: implications for the evolution of pentamery

Phylogenetic relationships among many lineages of angiosperms have been clarified via the analysis of large molecular data sets. However, with a data set of three genes (18S rDNA, rbcL, and atpB), relationships among lineages of core eudicots (Berberidopsidales, Caryophyllales, Gunnerales, Santalales, Saxifragales, asterids, rosids) remain essentially unresolved. We added 26S rDNA sequences to a three-gene matrix for 201 eudicots (8430 base pair aligned nucleotides per taxon). Parsimony analyses provided moderate (84%) jackknife support for Gunnerales, which comprise the two enigmatic families Gunneraceae and Myrothamnaceae, as sister to all other core eudicots. This position of Gunnerales has important implications for floral evolution. A dimerous or trimerous perianth is frequently encountered in early-diverging eudicots (e.g., Buxaceae, Proteales, Ranunculales, Trochodendraceae), whereas in core eudicots, pentamery predominates. Significantly, dimery is found in Gunneraceae and perhaps Myrothamnaceae (the merosity of the latter has also been interpreted as labile). Parsimony reconstructions of perianth merosity demonstrate lability among early-diverging eudicots and further indicate that a dimerous perianth could be the immediate precursor to the pentamerous condition characteristic of core eudicots. Thus, the developmental canalization that yielded the pentamerous condition of core eudicots occurred after the node leading to Gunnerales.     

Molecular evolution of PISTILLATA-like genes in the dogwood genus Cornus (Cornaceae)

The MADS-box gene family encodes critical regulators determining floral organ development. Understanding evolutionary patterns and processes of MADS-box genes is an important step toward unraveling the molecular basis of floral morphological evolution. In this study, we investigated the evolution of PI-like genes of the MADS-box family in the dogwood genus Cornus (Cornaceae). Cornus is a eudicot lineage in the asterids clade, and is intriguing in evolving petaloid bract morphology in two major lineages within the genus. The gene genealogy reconstructed using genomic DNA and cDNA sequences suggests multiple PI-like gene duplication events in Cornus. An ancient duplication event resulted in two ancient paralogs, CorPI-A and CorPI-B, which have highly diverged intron regions. Duplication of CorPI-A further resulted in two paralogs in one subgroup of Cornus, the BW group that does not produce modified bracts. Most species analyzed were found to contain more than one copy of the PI-like gene with most copies derived recently within species. Estimation and comparison of dN/dS ratios revealed relaxed selection in the PI-like gene in Cornus in comparison with the gene in the closely related outgroups Alangium and Davidia, and in other flowering plants. Selection also differed among major gene copies, CorPI-A and CorPI-B, and among different morphological subgroups of Cornus. Variation in selection pressures may indicate functional changes in PI-like genes after gene duplication and among different lineages. Strong positive selection at three amino acid sites of CorPI was also detected from a region critical for dimerization activity. Total substitution rates of the CorPI gene also differ among lineages of Cornus, showing a trend similar to that found in dN/dS ratios. We also found that the CorPI-A copy contains informative phylogenetic information when compared across species of Cornus.     

Improving phylogenetic inference of mushrooms with RPB1 and RPB2 nucleotide sequences (Inocybe; Agaricales)

Approximately 3000 bp across 84 taxa have been analyzed for variable regions of RPB1, RPB2, and nLSU-rDNA to infer phylogenetic relationships in the large ectomycorrhizal mushroom genus Inocybe (Agaricales; Basidiomycota). This study represents the first effort to combine variable regions of RPB1 and RPB2 with nLSU-rDNA for low-level phylogenetic studies in mushroom-forming fungi. Combination of the three loci increases non-parametric bootstrap support, Bayesian posterior probabilities, and resolution for numerous clades compared to separate gene analyses. These data suggest the evolution of at least five major lineages in Inocybe-the Inocybe clade, the Mallocybe clade, the Auritella clade, the Inosperma clade, and the Pseudosperma clade. Additionally, many clades nested within each major lineage are strongly supported. These results also suggest the family Crepiodataceae sensu stricto is sister to Inocybe. Recognition of Inocybe at the family level, the Inocybaceae, is recommended.