Complementing model species with model clades.

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant Tree of Life continues to improve. The intersection of these two research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade". These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis thaliana and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Genome-wide characterization and expression analysis of the CINNAMYL ALCOHOL DEHYDROGENASE gene family in Triticum aestivum.

BACKGROUND: CINNAMYL ALCOHOL DEHYDROGENASE (CAD) catalyzes the NADPH-dependent reduction of cinnamaldehydes into cinnamyl alcohols and is a key enzyme found at the final step of the monolignol pathway. Cinnamyl alcohols and their conjugates are subsequently polymerized in the secondary cell wall to form lignin. CAD genes are typically encoded by multi-gene families and thus traditionally organized into general classifications of functional relevance. RESULTS: In silico analysis of the hexaploid Triticum aestivum genome revealed 47 high confidence TaCAD copies, of which three were determined to be the most significant isoforms (class I) considered bone fide CADs. Class I CADs were expressed throughout development both in RNAseq data sets as well as via qRT-PCR analysis. Of the 37 class II TaCADs identified, two groups were observed to be significantly co-expressed with class I TaCADs in developing tissue and under chitin elicitation in RNAseq data sets. These co-expressed class II TaCADs were also found to be phylogenetically unrelated to a separate clade of class II TaCADs previously reported to be an influential resistance factor to pathogenic fungal infection. Lastly, two groups were phylogenetically identified as class III TaCADs, which possess distinct conserved gene structures. However, the lack of data supporting their catalytic activity for cinnamaldehydes and their bereft transcriptional presence in lignifying tissues challenges their designation and function as CADs. CONCLUSIONS: Taken together, our comprehensive transcriptomic analyses suggest that TaCAD genes contribute to overlapping but nonredundant functions during T. aestivum growth and development across a wide variety of agroecosystems and provide tolerance to various stressors.

Tangled webs and spider-flowers: Phylogenomics, biogeography, and seed morphology inform the evolutionary history of Cleomaceae.

PREMISE: Cleomaceae is an important model clade for studies of evolutionary processes including genome evolution, floral form diversification, and photosynthetic pathway evolution. Diversification and divergence patterns in Cleomaceae remain tangled as research has been restricted by its worldwide distribution, limited genetic sampling and species coverage, and a lack of definitive fossil calibration points. METHODS: We used target sequence capture and the Angiosperms353 probe set to perform a phylogenetic study of Cleomaceae. We estimated divergence times and biogeographic analyses to explore the origin and diversification of the family. Seed morphology across extant taxa was documented with multifocal image-stacking techniques and morphological characters were extracted, analyzed, and compared to fossil records. RESULTS: We recovered a well-supported and resolved phylogenetic tree of Cleomaceae generic relationships that includes 236 (~86%) species. We identified 11 principal clades and confidently placed Cleomella as sister to the rest of the family. Our analyses suggested that Cleomaceae and Brassicaceae diverged ~56 mya, and Cleomaceae began to diversify ~53 mya in the Palearctic and Africa. Multiple transatlantic disjunct distributions were identified. Seeds were imaged from 218 (~80%) species in the family and compared to all known fossil species. CONCLUSIONS: Our results represent the most comprehensive phylogenetic study of Cleomaceae to date. We identified transatlantic disjunctions and proposed explanations for these patterns, most likely either long-distance dispersals or contractions in latitudinal distributions caused by climate change over geological timescales. We found that seed morphology varied considerably but mostly mirrored generic relationships.

BUZZ: an essential gene for postinitiation root hair growth and a mediator of root architecture in Brachypodium distachyon.

Here, we discover a player in root development. Recovered from a forward-genetic screen in Brachypodium distachyon, the buzz mutant initiates root hairs but they fail to elongate. In addition, buzz roots grow twice as fast as wild-type roots. Also, lateral roots show increased sensitivity to nitrate, whereas primary roots are less sensitive to nitrate. Using whole-genome resequencing, we identified the causal single nucleotide polymorphism as occurring in a conserved but previously uncharacterized cyclin-dependent kinase (CDK)-like gene. The buzz mutant phenotypes are rescued by the wild-type B. distachyon BUZZ coding sequence and by an apparent homolog in Arabidopsis thaliana. Moreover, T-DNA mutants in A. thaliana BUZZ have shorter root hairs. BUZZ mRNA localizes to epidermal cells and develops root hairs and, in the latter, partially colocalizes with the NRT1.1A nitrate transporter. Based on qPCR and RNA-Seq, buzz overexpresses ROOT HAIRLESS LIKE SIX-1 and -2 and misregulates genes related to hormone signaling, RNA processing, cytoskeletal, and cell wall organization, and to the assimilation of nitrate. Overall, these data demonstrate that BUZZ is required for tip growth after root hair initiation and root architectural responses to nitrate.

Phylogenomic analyses reveal extensive gene flow within the magic flowers (Achimenes).

PREMISE OF THE STUDY: The Neotropical Gesneriaceae is a lineage known for its colorful and diverse flowers, as well as an extensive history of intra- and intergeneric hybridization, particularly among Achimenes (the magic flowers) and other members of subtribe Gloxiniinae. Despite numerous studies seeking to elucidate the evolutionary relationships of these lineages, relatively few have sought to infer specific patterns of gene flow despite evidence of widespread hybridization. METHODS: To explore the utility of phylogenomic data for reassessing phylogenetic relationships and inferring patterns of gene flow among species of Achimenes, we sequenced 12 transcriptomes. We used a variety of methods to infer the species tree, examine gene tree discordance, and infer patterns of gene flow. KEY RESULTS: Phylogenomic analyses resolve clade relationships at the crown of the lineage with strong support. In contrast to previous analyses, we recovered strong support for several new relationships despite a significant amount of gene tree discordance. We present evidence for at least two introgression events between two species pairs that share pollinators, and suggest that the species status of Achimenes admirabilis be reexamined. CONCLUSIONS: Our study demonstrates the utility of transcriptome data for phylogenomic analyses, and inferring patterns of gene flow despite gene tree discordance. Moreover, these data provide another example of prevalent interspecific gene flow among Neotropical plants that share pollinators.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae).

BACKGROUND: Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae). RESULTS: The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species. CONCLUSIONS: The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Distinct Processes Drive Diversification in Different Clades of Gesneriaceae.

Using a time-calibrated phylogenetic hypothesis including 768 Gesneriaceae species (out of [Formula: see text]3300 species) and more than 29,000 aligned bases from 26 gene regions, we test Gesneriaceae for diversification rate shifts and the possible proximal drivers of these shifts: geographic distributions, growth forms, and pollination syndromes. Bayesian Analysis of Macroevolutionary Mixtures analyses found five significant rate shifts in Beslerieae, core Nematanthus, core Columneinae, core Streptocarpus, and Pacific Cyrtandra These rate shifts correspond with shifts in diversification rates, as inferred by Binary State Speciation and Extinction Model and Geographic State Speciation and Extinction model, associated with hummingbird pollination, epiphytism, unifoliate growth, and geographic area. Our results suggest that diversification processes are extremely variable across Gesneriaceae clades with different combinations of characters influencing diversification rates in different clades. Diversification patterns between New and Old World lineages show dramatic differences, suggesting that the processes of diversification in Gesneriaceae are very different in these two geographic regions.

Kranz and single-cell forms of C4 plants in the subfamily Suaedoideae show kinetic C4 convergence for PEPC and Rubisco with divergent amino acid substitutions.

The two carboxylation reactions performed by phosphoenolpyruvate carboxylase (PEPC) and ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) are vital in the fixation of inorganic carbon for C4 plants. The abundance of PEPC is substantially elevated in C4 leaves, while the location of Rubisco is restricted to one of two chloroplast types. These differences compared with C3 leaves have been shown to result in convergent enzyme optimization in some C4 species. Investigation into the kinetic properties of PEPC and Rubisco from Kranz C4, single cell C4, and C3 species in Chenopodiaceae s. s. subfamily Suaedoideae showed that these major carboxylases in C4 Suaedoideae species lack the same mutations found in other C4 systems which have been examined; but still have similar convergent kinetic properties. Positive selection analysis on the N-terminus of PEPC identified residues 364 and 368 to be under positive selection with a posterior probability >0.99 using Bayes empirical Bayes. Compared with previous analyses on other C4 species, PEPC from C4 Suaedoideae species have different convergent amino acids that result in a higher K m for PEP and malate tolerance compared with C3 species. Kinetic analysis of Rubisco showed that C4 species have a higher catalytic efficiency of Rubisco (k catc in mol CO2 mol(-1) Rubisco active sites s(-1)), despite lacking convergent substitutions in the rbcL gene. The importance of kinetic changes to the two-carboxylation reactions in C4 leaves related to amino acid selection is discussed.

Environmental and genetic correlates of allocation to sexual reproduction in the circumpolar plant Bistorta vivipara.

UNLABELLED: • PREMISE OF THE STUDY: Sexual reproduction often requires more energy and time than clonal reproduction. In marginal arctic conditions, species that can reproduce both sexually and clonally dominate. Plants with this capacity may thrive because they can alter reproduction depending on environmental conditions. Bistorta vivipara is a circumpolar herb that predominately reproduces clonally, but certain environmental conditions promote higher investment in flowers (and possible sexual reproduction). Despite largely reproducing clonally, the herb has high levels of genetic variation, and the processes underlying this paradoxical pattern of variation remain unclear. Here we identified environmental factors associated with sexual investment and examined whether sexual reproduction is associated with higher levels of genetic variation.• METHODS: We sampled 20 populations of B. vivipara across the high Arctic archipelago of Svalbard. In each population, we measured reproductive traits, environmental variables, and collected samples for genetic analyses. These samples permitted hypotheses to be tested regarding sexual investment and ecological and genetic correlates.• KEY RESULTS: Increased soil nitrogen and organic matter content and decreased elevation were positively associated with investment in flowers. Increased investment in flowers significantly correlated with more genotypes per population. Linkage disequilibrium was consistent with predominant clonality, but several populations showed higher genetic variation and lower differentiation than expected. There was no geographical genetic structure.• CONCLUSIONS: In B. vivipara, sexual investment is positively associated with habitat quality. Bistorta vivipara predominantly reproduces clonally, but occasional outcrossing, efficient clonal reproduction, and dispersal by bulbils can explain the considerable genetic variation and weak genetic structure in B. vivipara.

Insights into the evolution and diversification of the AT-hook Motif Nuclear Localized gene family in land plants.

BACKGROUND: Members of the ancient land-plant-specific transcription factor AT-Hook Motif Nuclear Localized (AHL) gene family regulate various biological processes. However, the relationships among the AHL genes, as well as their evolutionary history, still remain unexplored. RESULTS: We analyzed over 500 AHL genes from 19 land plant species, ranging from the early diverging Physcomitrella patens and Selaginella to a variety of monocot and dicot flowering plants. We classified the AHL proteins into three types (Type-I/-II/-III) based on the number and composition of their functional domains, the AT-hook motif(s) and PPC domain. We further inferred their phylogenies via Bayesian inference analysis and predicted gene gain/loss events throughout their diversification. Our analyses suggested that the AHL gene family emerged in embryophytes and further evolved into two distinct clades, with Type-I AHLs forming one clade (Clade-A), and the other two types together diversifying in another (Clade-B). The two AHL clades likely diverged before the separation of Physcomitrella patens from the vascular plant lineage. In angiosperms, Clade-A AHLs expanded into 5 subfamilies; while, the ones in Clade-B expanded into 4 subfamilies. Examination of their expression patterns suggests that the AHLs within each clade share similar expression patterns with each other; however, AHLs in one monophyletic clade exhibit distinct expression patterns from the ones in the other clade. Over-expression of a Glycine max AHL PPC domain in Arabidopsis thaliana recapitulates the phenotype observed when over-expressing its Arabidopsis thaliana counterpart. This result suggests that the AHL genes from different land plant species may share conserved functions in regulating plant growth and development. Our study further suggests that such functional conservation may be due to conserved physical interactions among the PPC domains of AHL proteins. CONCLUSIONS: Our analyses reveal a possible evolutionary scenario for the AHL gene family in land plants, which will facilitate the design of new studies probing their biological functions. Manipulating the AHL genes has been suggested to have tremendous effects in agriculture through increased seedling establishment, enhanced plant biomass and improved plant immunity. The information gleaned from this study, in turn, has the potential to be utilized to further improve crop production.

Positive selection of Kranz and non-Kranz C4 phosphoenolpyruvate carboxylase amino acids in Suaedoideae (Chenopodiaceae).

In subfamily Suaedoideae, four independent gains of C4 photosynthesis are proposed, which includes two parallel origins of Kranz anatomy (sections Salsina and Schoberia) and two independent origins of single-cell C4 anatomy (Bienertia and Suaeda aralocaspica). Additional phylogenetic support for this hypothesis was generated from sequence data of the C-terminal portion of the phosphoenolpyruvate carboxylase (PEPC) gene used in C4 photosynthesis (ppc-1) in combination with previous sequence data. ppc-1 sequence was generated for 20 species in Suaedoideae and two outgroup Salsola species that included all types of C4 anatomies as well as two types of C3 anatomies. A branch-site test for positively selected codons was performed using the software package PAML. From labelling of the four branches where C4 is hypothesized to have developed (foreground branches), residue 733 (maize numbering) was identified to be under positive selection with a posterior probability >0.99 and residue 868 at the >0.95 interval using Bayes empirical Bayes (BEB). When labelling all the branches within C4 clades, the branch-site test identified 13 codons to be under selection with a posterior probability >0.95 by BEB; this is discussed considering current information on functional residues. The signature C4 substitution of an alanine for a serine at position 780 in the C-terminal end (which is considered a major determinant of affinity for PEP) was only found in four of the C4 species sampled, while eight of the C4 species and all the C3 species have an alanine residue; indicating that this substitution is not a requirement for C4 function.

Using supermatrices for phylogenetic inquiry: an example using the sedges.

In this article, we use supermatrix data-mining methods to reconstruct a large, highly inclusive phylogeny of Cyperaceae from nucleotide data available on GenBank. We explore the properties of these trees and their utility for phylogenetic inference, and show that even the highly incomplete alignments characteristic of supermatrix approaches may yield very good estimates of phylogeny. We present a novel pipeline for filtering sparse alignments to improve their phylogenetic utility by maximizing the partial decisiveness of the matrices themselves through a technique we call "phylogenetic scaffolding," and we present a new method of scoring tip instability (i.e. "rogue taxa") based on the I statistic implemented in the software Mesquite. The modified statistic, which we call I(S), is somewhat more straightforward to interpret than similar statistics, and our implementation of it may be applied to large sets of large trees. The largest sedge trees presented here contain more than 1500 tips (about one quarter of all sedge species) and are based on multigene alignments with more than 20 000 sites and more than 90% missing data. These trees match well with previously supported phylogenetic hypotheses, but have lower overall support values and less resolution than more heavily filtered trees. Our best-resolved trees are characterized by stronger support values than any previously published sedge phylogenies, and show some relationships that are incongruous with previous studies. Overall, we show that supermatrix methods offer powerful means of pursuing phylogenetic study and these tools have high potential value for many systematic biologists.

Potential use of low-copy nuclear genes in DNA barcoding: a comparison with plastid genes in two Hawaiian plant radiations.

BACKGROUND: DNA barcoding of land plants has relied traditionally on a small number of markers from the plastid genome. In contrast, low-copy nuclear genes have received little attention as DNA barcodes because of the absence of universal primers for PCR amplification. RESULTS: From pooled-species 454 transcriptome data we identified two variable intron-less nuclear loci for each of two species-rich genera of the Hawaiian flora: Clermontia (Campanulaceae) and Cyrtandra (Gesneriaceae) and compared their utility as DNA barcodes with that of plastid genes. We found that nuclear genes showed an overall greater variability, but also displayed a high level of heterozygosity, intraspecific variation, and retention of ancient alleles. Thus, nuclear genes displayed fewer species-diagnostic haplotypes compared to plastid genes and no interspecies gaps. CONCLUSIONS: The apparently greater coalescence times of nuclear genes are likely to limit their utility as barcodes, as only a small proportion of their alleles were fixed and unique to individual species. In both groups, species-diagnostic markers from either genome were scarce on the youngest island; a minimum age of ca. two million years may be needed for a species flock to be barcoded. For young plant groups, nuclear genes may not be a superior alternative to slowly evolving plastid genes.

Gene discordance in phylogenomics of recent plant radiations, an example from Hawaiian Cyrtandra (Gesneriaceae).

Resolving species relationships within recent radiations requires analysis at the interface of phylogenetics and population genetics, where coalescence and hybridization may confound our understanding of relationships. We developed 18 new primer pairs for nuclear loci in Cyrtandra (Gesneriaceae), one of the largest plant radiations in the Pacific Islands, and tested the concordance of 14 loci in establishing the phylogenetic relationships of a small number of Hawaiian species. Four genes yielded tree topologies conflicting with the primary concordance tree, suggesting plastid capture and horizontal transfer via hybridization. Combining all concordant genes yielded a tree with stronger support and a different topology from the total-evidence tree. We conclude that a small number of genes may be insufficient for accurate reconstruction of the phylogenetic relationships among closely related species. Further, the combination of genes for phylogenetic analysis without preliminary concordance tests can yield an erroneous tree topology. It seems that the number of genes needed for phylogenetic analysis of closely related species is significantly greater than the small numbers commonly used, which fail to isolate coalescence, introgression and hybridization.

Structural and physiological analyses in Salsoleae (Chenopodiaceae) indicate multiple transitions among C3, intermediate, and C4 photosynthesis.

In subfamily Salsoloideae (family Chenopodiaceae) most species are C4 plants having terete leaves with Salsoloid Kranz anatomy characterized by a continuous dual chlorenchyma layer of Kranz cells (KCs) and mesophyll (M) cells, surrounding water storage and vascular tissue. From section Coccosalsola sensu Botschantzev, leaf structural and photosynthetic features were analysed on selected species of Salsola which are not performing C4 based on leaf carbon isotope composition. The results infer the following progression in distinct functional and structural forms from C3 to intermediate to C4 photosynthesis with increased leaf succulence without changes in vein density: From species performing C3 photosynthesis with Sympegmoid anatomy with two equivalent layers of elongated M cells, with few organelles in a discontinuous layer of bundle sheath (BS) cells (S. genistoides, S. masenderanica, S. webbii) > development of proto-Kranz BS cells having mitochondria in a centripetal position and increased chloroplast number (S. montana) > functional C3-C4 intermediates having intermediate CO2 compensation points with refixation of photorespired CO2, development of Kranz-like anatomy with reduction in the outer M cell layer to hypodermal-like cells, and increased specialization (but not size) of a Kranz-like inner layer of cells with increased cell wall thickness, organelle number, and selective expression of mitochondrial glycine decarboxylase (Kranz-like Sympegmoid, S. arbusculiformis; and Kranz-like Salsoloid, S. divaricata) > selective expression of enzymes between the two cell types for performing C4 with Salsoloid-type anatomy. Phylogenetic analysis of tribe Salsoleae shows the occurrence of C3 and intermediates in several clades, and lineages of interest for studying different forms of anatomy.

Antisense expression of peach mildew resistance locus O (PpMlo1) gene confers cross-species resistance to powdery mildew in Fragaria x ananassa.

Powdery mildew (PM) is one of the major plant pathogens. The conventional method of PM control includes frequent use of sulfur-based fungicides adding to production costs and potential harm to the environment. PM remains a major scourge for Rosaceae crops where breeding approaches mainly resort to gene-for-gene resistance. We have tested an alternate source of PM resistance in Rosaceae. Mildew resistance locus O (MLO) has been well studied in barley due to its role in imparting broad spectrum resistance to PM. We identified PpMlo1 (Prunus persica Mlo) in peach and characterized it further to test if a similar mechanism of resistance is conserved in Rosaceae. Due to its recalcitrance in tissue culture, reverse genetic studies involving PpMloI were not feasible in peach. Therefore, Fragaria x ananassa LF9 line, a taxonomic surrogate, was used for functional analysis of PpMlo1. Agrobacterium-mediated transformation yielded transgenic strawberry plants expressing PpMlo1 in sense and antisense orientation. Antisense expression of PpMlo1 in transgenic strawberry plants conferred resistance to Fragaria-specific powdery mildew, Podosphaera macularis. Phylogenetic analysis of 208 putative Mlo gene copies from 35 plant species suggests a large number of duplications of this gene family prior to the divergence of monocots and eudicots, early in eudicot diversification. Our results indicate that the Mlo-based resistance mechanism is functional in Rosaceae, and that Fragaria can be used as a host to test mechanistic function of genes derived from related tree species. To the best of our knowledge, this work is one of the first attempts at testing the potential of using a Mlo-based resistance strategy to combat powdery mildew in Rosaceae.

Phylogenetic relationships among the North American cleomoids (Cleomaceae): a test of Iltis's reduction series.

PREMISE OF STUDY: A monophyletic group composed of five genera of the Cleomaceae represents an intriguing lineage with outstanding taxonomic and evolutionary questions. Generic boundaries are poorly defined, and historical hypotheses regarding the evolution of fruit type and phylogenetic relationships provide testable questions. This is the first detailed phylogenetic investigation of all 22 species in this group. We use this phylogenetic framework to assess generic monophyly and test Iltis's evolutionary "reduction series" hypothesis regarding phylogeny and fruit type/seed number. • METHODS: Maximum likelihood and Bayesian analyses of four plastid intergenic spacer region sequences (rpl32-trnL, trnQ-rps16, ycf1-rps15, and psbA-trnH) and one nuclear (ITS) region were used to reconstruct phylogenetic relationships among the NA cleomoid species. Stochastic mapping and ancestral-state reconstruction were used to study the evolution of fruit type. • KEY RESULTS: Both analyses recovered nearly identical phylogenies. Three of the currently recognized genera (Wislizenia, Carsonia, and Oxystylis) are monophyletic while two (Cleomella and Peritoma) are para- or polyphyletic. There was a single origin of the two-seeded schizocarp in the ancestor of the Oxystylis-Wislizenia clade and a secondary derivation of elongated capsule-type fruits in Peritoma from a truncated capsule state in Cleomella. • CONCLUSIONS: Our well-resolved phylogeny supports most of the current species circumscriptions but not current generic circumscriptions. Additionally, our results are inconsistent with Iltis's hypothesis of species with elongated many-seed fruits giving rise to species with truncated few-seeded fruits. Instead, we find support for the reversion to elongated multiseeded fruits from a truncate few-seeded ancestor in Peritoma.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner.

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Comparative genomic analysis and phylogenetic position of Theileria equi.

BACKGROUND: Transmission of arthropod-borne apicomplexan parasites that cause disease and result in death or persistent infection represents a major challenge to global human and animal health. First described in 1901 as Piroplasma equi, this re-emergent apicomplexan parasite was renamed Babesia equi and subsequently Theileria equi, reflecting an uncertain taxonomy. Understanding mechanisms by which apicomplexan parasites evade immune or chemotherapeutic elimination is required for development of effective vaccines or chemotherapeutics. The continued risk of transmission of T. equi from clinically silent, persistently infected equids impedes the goal of returning the U. S. to non-endemic status. Therefore comparative genomic analysis of T. equi was undertaken to: 1) identify genes contributing to immune evasion and persistence in equid hosts, 2) identify genes involved in PBMC infection biology and 3) define the phylogenetic position of T. equi relative to sequenced apicomplexan parasites. RESULTS: The known immunodominant proteins, EMA1, 2 and 3 were discovered to belong to a ten member gene family with a mean amino acid identity, in pairwise comparisons, of 39%. Importantly, the amino acid diversity of EMAs is distributed throughout the length of the proteins. Eight of the EMA genes were simultaneously transcribed. As the agents that cause bovine theileriosis infect and transform host cell PBMCs, we confirmed that T. equi infects equine PBMCs, however, there is no evidence of host cell transformation. Indeed, a number of genes identified as potential manipulators of the host cell phenotype are absent from the T. equi genome. Comparative genomic analysis of T. equi revealed the phylogenetic positioning relative to seven apicomplexan parasites using deduced amino acid sequences from 150 genes placed it as a sister taxon to Theileria spp. CONCLUSIONS: The EMA family does not fit the paradigm for classical antigenic variation, and we propose a novel model describing the role of the EMA family in persistence. T. equi has lost the putative genes for host cell transformation, or the genes were acquired by T. parva and T. annulata after divergence from T. equi. Our analysis identified 50 genes that will be useful for definitive phylogenetic classification of T. equi and closely related organisms.