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 2 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 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.

Phylogenomics and the rise of the angiosperms

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

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.

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.

Plastid transit peptides-where do they come from and where do they all belong? Multi-genome and pan-genomic assessment of chloroplast transit peptide evolution

Subcellular relocalization of proteins determines an organism's metabolic repertoire and thereby its survival in unique evolutionary niches. In plants, the plastid and its various morphotypes import a large and varied number of nuclear-encoded proteins to orchestrate vital biochemical reactions in a spatiotemporal context. Recent comparative genomics analysis and high-throughput shotgun proteomics data indicate that there are a large number of plastid-targeted proteins that are either semi-conserved or non-conserved across different lineages. This implies that homologs are differentially targeted across different species, which is feasible only if proteins have gained or lost plastid targeting peptides during evolution. In this study, a broad, multi-genome analysis of 15 phylogenetically diverse genera and in-depth analyses of pangenomes from Arabidopsis and Brachypodium were performed to address the question of how proteins acquire or lose plastid targeting peptides. The analysis revealed that random insertions or deletions were the dominant mechanism by which novel transit peptides are gained by proteins. While gene duplication was not a strict requirement for the acquisition of novel subcellular targeting, 40% of novel plastid-targeted genes were found to be most closely related to a sequence within the same genome, and of these, 30.5% resulted from alternative transcription or translation initiation sites. Interestingly, analysis of the distribution of amino acids in the transit peptides of known and predicted chloroplast-targeted proteins revealed monocot and eudicot-specific preferences in residue distribution.

Genome-Scale Characterization of Predicted Plastid-Targeted Proteomes in Higher Plants

Plastids are morphologically and functionally diverse organelles that are dependent on nuclear-encoded, plastid-targeted proteins for all biochemical and regulatory functions. However, how plastid proteomes vary temporally, spatially, and taxonomically has been historically difficult to analyze at a genome-wide scale using experimental methods. A bioinformatics workflow was developed and evaluated using a combination of fast and user-friendly subcellular prediction programs to maximize performance and accuracy for chloroplast transit peptides and demonstrate this technique on the predicted proteomes of 15 sequenced plant genomes. Gene family grouping was then performed in parallel using modified approaches of reciprocal best BLAST hits (RBH) and UCLUST. A total of 628 protein families were found to have conserved plastid targeting across angiosperm species using RBH, and 828 using UCLUST. However, thousands of clusters were also detected where only one species had predicted plastid targeting, most notably in Panicum virgatum which had 1,458 proteins with species-unique targeting. An average of 45% overlap was found in plastid-targeted protein-coding gene families compared with Arabidopsis, but an additional 20% of proteins matched against the full Arabidopsis proteome, indicating a unique evolution of plastid targeting. Neofunctionalization through subcellular relocalization is known to impart novel biological functions but has not been described before on a genome-wide scale for the plastid proteome. Further work to correlate these predicted novel plastid-targeted proteins to transcript abundance and high-throughput proteomics will uncover unique aspects of plastid biology and shed light on how the plastid proteome has evolved to influence plastid morphology and biochemistry.

Co-expression clustering across flower development identifies modules for diverse floral forms in Achimenes (Gesneriaceae)

BACKGROUND

Genetic pathways involved with flower color and shape are thought to play an important role in the development of flowers associated with different pollination syndromes, such as those associated with bee, butterfly, or hummingbird pollination. Because pollination syndromes are complex traits that are orchestrated by multiple genes and pathways, the gene regulatory networks have not been explored. Gene co-expression networks provide a systems level approach to identify important contributors to floral diversification.

METHODS

RNA-sequencing was used to assay gene expression across two stages of flower development (an early bud and an intermediate stage) in 10 species of Achimenes (Gesneriaceae). Two stage-specific co-expression networks were created from 9,503 orthologs and analyzed to identify module hubs and the network periphery. Module association with bee, butterfly, and hummingbird pollination syndromes was tested using phylogenetic mixed models. The relationship between network connectivity and evolutionary rates (d N/d S) was tested using linear models.

RESULTS

Networks contained 65 and 62 modules that were largely preserved between developmental stages and contained few stage-specific modules. Over a third of the modules in both networks were associated with flower color, shape, and pollination syndrome. Within these modules, several hub nodes were identified that related to the production of anthocyanin and carotenoid pigments and the development of flower shape. Evolutionary rates were decreased in highly connected genes and elevated in peripheral genes.

DISCUSSION

This study aids in the understanding of the genetic architecture and network properties underlying the development of floral form and provides valuable candidate modules and genes for future studies.

Lessons from Cleomaceae, the Sister of Crucifers

Cleomaceae is a diverse group well-suited to addressing fundamental genomic and evolutionary questions as the sister group to Brassicaceae, facilitating transfer of knowledge from the model Arabidopsis thaliana. Phylogenetic and taxonomic revisions provide a framework for examining the evolution of substantive morphological and physiology diversity in Cleomaceae, but not necessarily in Brassicaceae. The investigation of both nested and contrasting whole-genome duplications (WGDs) between Cleomaceae and Brassicaceae allows comparisons of independently duplicated genes and investigation of whether they may be drivers of the observed innovations. Further, a wealth of outstanding genetic research has provided insight into how the important alternative carbon fixation pathway, C₄ photosynthesis, has evolved via differential expression of a suite of genes, of which the underlying mechanisms are being elucidated.

Genome-wide identification and analysis of the ALTERNATIVE OXIDASE gene family in diploid and hexaploid wheat

A comprehensive understanding of wheat responses to environmental stress will contribute to the long-term goal of feeding the planet. ALERNATIVE OXIDASE (AOX) genes encode proteins involved in a bypass of the electron transport chain and are also known to be involved in stress tolerance in multiple species. Here, we report the identification and characterization of the AOX gene family in diploid and hexaploid wheat. Four genes each were found in the diploid ancestors Triticum urartu, and Aegilops tauschii, and three in Aegilops speltoides. In hexaploid wheat (Triticum aestivum), 20 genes were identified, some with multiple splice variants, corresponding to a total of 24 proteins for those with observed transcription and translation. These proteins were classified as AOX1a, AOX1c, AOX1e or AOX1d via phylogenetic analysis. Proteins lacking most or all signature AOX motifs were assigned to putative regulatory roles. Analysis of protein-targeting sequences suggests mixed localization to the mitochondria and other organelles. In comparison to the most studied AOX from Trypanosoma brucei, there were amino acid substitutions at critical functional domains indicating possible role divergence in wheat or grasses in general. In hexaploid wheat, AOX genes were expressed at specific developmental stages as well as in response to both biotic and abiotic stresses such as fungal pathogens, heat and drought. These AOX expression patterns suggest a highly regulated and diverse transcription and expression system. The insights gained provide a framework for the continued and expanded study of AOX genes in wheat for stress tolerance through breeding new varieties, as well as resistance to AOX-targeted herbicides, all of which can ultimately be used synergistically to improve crop yield.

Non-dispersive phloem-protein bodies (NPBs) of Populus trichocarpa consist of a SEOR protein and do not respond to cell wounding and Ca2

Differentiating sieve elements in the phloem of angiosperms produce abundant phloem-specific proteins before their protein synthesis machinery is degraded. These P-proteins initially form dense bodies, which disperse into individual filaments when the sieve element matures. In some cases, however, the dense protein agglomerations remain intact and are visible in functional sieve tubes as non-dispersive P-protein bodies, or NPBs. Species exhibiting NPBs are distributed across the entire angiosperm clade. We found that NPBs in the model tree, Populus trichocarpa, resemble the protein bodies described from other species of the order Malpighiales as they all consist of coaligned tubular fibrils bundled in hexagonal symmetry. NPBs of all Malpighiales tested proved unresponsive to sieve tube wounding and Ca²⁺. The P. trichocarpa NPBs consisted of a protein encoded by a gene that in the genome database of this species had been annotated as a homolog of SEOR1 (sieve element occlusion-related 1) in Arabidopsis. Sequencing of the gene in our plants corroborated this interpretation, and we named the gene PtSEOR1. Previously characterized SEOR proteins form irregular masses of P-protein slime in functional sieve tubes. We conclude that a subgroup of these proteins is involved in the formation of NPBs at least in the Malpighiales, and that these protein bodies have no role in rapid wound responses of the sieve tube network.

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.

Notas taxonómicas y de distribución de Carex (Cyperaceae) en el Neotrópico

Notas taxonómicas y de distribución de Carex (Cyperaceae) en el Neotrópico. Presentamos datos relevantes acerca de la taxonomía y distribución de 20 especies del género Carex en el Neotrópico(de la region Caribe a la Patagonia). Aportamos un total de 19 registros nuevos, dos de ellos para especies hasta el momento no registradas para Sudamérica, diez especies citadas como novedad para diferentes países y siete especies aquí reportadas como novedades regionales en Argentina, Bolivia y Chile. Lectotipificamos de nueve nombres de especies (C. brehmeri, C. catharinensis, C. ecuadorica, C. firmicaulis, C. fuscula, C. mandoniana, C. phleoides, C. polysticha y C. subantarctica) y epitipificamos uno de ellos (C. fuscula). Finalmente, proponemos a Carex koyamae como nuevo sinónimo de Carex phleoides.

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.

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.

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.

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.

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.

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.

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.

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.

Diversity in forms of C4 in the genus Cleome (Cleomaceae)

BACKGROUND AND AIMS

Cleomaceae is one of 19 angiosperm families in which C(4) photosynthesis has been reported. The aim of the study was to determine the type, and diversity, of structural and functional forms of C(4) in genus Cleome. Methods Plants of Cleome species were grown from seeds, and leaves were subjected to carbon isotope analysis, light and scanning electron microscopy, western blot analysis of proteins, and in situ immunolocalization for ribulose bisphosphate carboxylase oxygenase (Rubisco) and phosphoenolpyruvate carboxylase (PEPC).

KEY RESULTS

Three species with C(4)-type carbon isotope values occurring in separate lineages in the genus (Cleome angustifolia, C. gynandra and C. oxalidea) were shown to have features of C(4) photosynthesis in leaves and cotyledons. Immunolocalization studies show that PEPC is localized in mesophyll (M) cells and Rubisco is selectively localized in bundle sheath (BS) cells in leaves and cotyledons, characteristic of species with Kranz anatomy. Analyses of leaves for key photosynthetic enzymes show they have high expression of markers for the C(4) cycle (compared with the C(3)-C(4) intermediate C. paradoxa and the C(3) species C. africana). All three are biochemically NAD-malic enzyme sub-type, with higher granal development in BS than in M chloroplasts, characteristic of this biochemical sub-type. Cleome gynandra and C. oxalidea have atriplicoid-type Kranz anatomy with multiple simple Kranz units around individual veins. However, C. angustifolia anatomy is represented by a double layer of concentric chlorenchyma forming a single compound Kranz unit by surrounding all the vascular bundles and water storage cells.

CONCLUSIONS

NAD-malic enzyme-type C(4) photosynthesis evolved multiple times in the family Cleomaceae, twice with atriplicoid-type anatomy in compound leaves having flat, broad leaflets in the pantropical species C. gynandra and the Australian species C. oxalidea, and once by forming a single Kranz unit in compound leaves with semi-terete leaflets in the African species C. angustifolia. The leaf morphology of C. angustifolia, which is similar to that of the sister, C(3)-C(4) intermediate African species C. paradoxa, suggests adaptation of this lineage to arid environments, which is supported by biogeographical information.

Basic helix-loop-helix transcription factor gene family phylogenetics and nomenclature

A phylogenetic analysis of the basic helix-loop-helix (bHLH) gene superfamily was performed using seven different species (human, mouse, rat, worm, fly, yeast, and plant Arabidopsis) and involving over 600 bHLH genes (Stevens et al., 2008). All bHLH genes were identified in the genomes of the various species, including expressed sequence tags, and the entire coding sequence was used in the analysis. Nearly 15% of the gene family has been updated or added since the original publication. A super-tree involving six clades and all structural relationships was established and is now presented for four of the species. The wealth of functional data available for members of the bHLH gene superfamily provides us with the opportunity to use this exhaustive phylogenetic tree to predict potential functions of uncharacterized members of the family. This phylogenetic and genomic analysis of the bHLH gene family has revealed unique elements of the evolution and functional relationships of the different genes in the bHLH gene family.

Range persistence during the last glacial maximum: Carex macrocephala was not restricted to glacial refugia

The distribution of many species inhabiting northwestern North America has been heavily influenced by the climatic changes during the late Pleistocene. Several studies have suggested that species were restricted to glacial refugia north and/or south of the continental ice sheet front. It is also hypothesized that the coast of northwestern North America could have been a prime location for glacial refugia because of the lowering of the eustatic sea level and the concomitant rise of the continental shelf because of tectonic rebound. Alternatively, some coastal species distributions and demographics may have been unaffected in the long-term by the last glacial maximum (LGM). We tested the glacial refugium hypothesis on an obligate coastal plant species, Carex macrocephala by sampling 600 individuals from 41 populations with 11 nuclear microsatellite loci and the rpL16 plastid intragenic spacer region. The microsatellite data sets suggest a low level of population differentiation with a standardized G'(ST) = 0.032 and inbreeding was high with an F = 0.969. The homogenization of the populations along the coast was supported by a principal coordinate analysis, amovas and samova analyses. Analyses using the rpL16 data set support the results of the microsatellite analyses, with a low F(ST) of 0.042. Coalescent and mismatch analyses using rpL16 suggest that C. macrocephala has not gone through a significant bottleneck within the past 100,000 years, although a much earlier population expansion was indicated by the mismatch analysis. Carex macrocephala exhibits the characteristics of metapopulation dynamics and on the basis of these results, we concluded that it was not restricted to glacial refugia during the LGM, but that it existed as a large metapopulation.

Stem architecture in Eleocharis subgenus Limnochloa (Cyperaceae): Evidence of dynamic morphological evolution in a group of pantropical sedges

We examined phylogenetic relationships and patterns of stem structural evolution in Eleocharis subgenus Limnochloa, an ecologically and economically important group of tropical to temperate-growing sedges, whose stems serve as the primary photosynthetic organs. We used maximum parsimony, likelihood, and Bayesian inference to develop phylogenetic trees and stochastic mapping and a Markov one-rate model to develop character history reconstructions of stem architecture. A complex history of stem shape evolution characterized by a high degree of homoplasy and rapid rates of change (an average of 13 transitions per character history for about 25 species) was identified across subgenus Limnochloa. Character states transition much more frequently in some lineages than others, but tend to follow a consistent directional pattern of evolutionary change. Our data also suggest that changes in stem shape and anatomy may be associated with speciation events in the subgenus (Pagel's κ = 0.3503, P = 0.04579) and may have some adaptive significance. The potential adaptive roles of stem structural traits are unclear, but may be elucidated by further studies. This work serves as a starting point for future evolutionary studies of stem shape and structure in monocots and provides important background knowledge for further studies of ecological adaptations of Eleocharis.

Babesia bovis: a comprehensive phylogenetic analysis of plastid-encoded genes supports green algal origin of apicoplasts

Apicomplexan parasites commonly contain a unique, non-photosynthetic plastid-like organelle termed the apicoplast. Previous analyses of other plastid-containing organisms suggest that apicoplasts were derived from a red algal ancestor. In this report, we present an extensive phylogenetic study of apicoplast origins using multiple previously reported apicoplast sequences as well as several sequences recently reported. Phylogenetic analysis of amino acid sequences was used to determine the evolutionary origin of the organelle. A total of nine plastid genes from 37 species were incorporated in our study. The data strongly support a green algal origin for apicoplasts and Euglenozoan plastids. Further, the nearest green algae lineage to the Apicomplexans is the parasite Helicosporidium, suggesting that apicoplasts may have originated by lateral transfer from green algal parasite lineages. The results also substantiate earlier findings that plastids found in Heterokonts such as Odontella and Thalassiosira were derived from a separate secondary endosymbiotic event likely originating from a red algal lineage.

A comparative study in ancestral range reconstruction methods: retracing the uncertain histories of insular lineages

Island systems have long been useful models for understanding lineage diversification in a geographic context, especially pertaining to the importance of dispersal in the origin of new clades. Here we use a well-resolved phylogeny of the flowering plant genus Cyrtandra (Gesneriaceae) from the Pacific Islands to compare four methods of inferring ancestral geographic ranges in islands: two developed for character-state reconstruction that allow only single-island ranges and do not explicitly associate speciation with range evolution (Fitch parsimony [FP; parsimony-based] and stochastic mapping [SM; likelihood-based]) and two methods developed specifically for ancestral range reconstruction, in which widespread ranges (spanning islands) are integral to inferences about speciation scenarios (dispersal-vicariance analysis [DIVA; parsimony-based] and dispersal-extinction-cladogenesis [DEC; likelihood-based]). The methods yield conflicting results, which we interpret in light of their respective assumptions. FP exhibits the least power to unequivocally reconstruct ranges, likely due to a combination of having flat (uninformative) transition costs and not using branch length information. SM reconstructions generally agree with a prior hypothesis about dispersal-driven speciation across the Pacific, despite the conceptual mismatch between its character-based model and this mode of range evolution. In contrast with narrow extant ranges for species of Cyrtandra, DIVA reconstructs broad ancestral ranges at many nodes. DIVA results also conflict with geological information on island ages; we attribute these conflicts to the parsimony criterion not considering branch lengths or time, as well as vicariance being the sole means of divergence for widespread ancestors. DEC analyses incorporated geological information on island ages and allowed prior hypotheses about range size and dispersal rates to be evaluated in a likelihood framework and gave more nuanced inferences about range evolution and the geography of speciation than other methods tested. However, ancestral ranges at several nodes could not be conclusively resolved, due possibly to uncertainty in the phylogeny or the relative complexity of the underlying model. Of the methods tested, SM and DEC both converge on plausible hypotheses for area range histories in Cyrtandra, due in part to the consideration of branch lengths and/or timing of events. We suggest that DEC model-based methods for ancestral range inference could be improved by adopting a Bayesian SM approach, in which stochastic sampling of complete geographic histories could be integrated over alternative phylogenetic topologies. Likelihood-based estimates of ancestral ranges for Cyrtandra suggest a major dispersal route into the Pacific through the islands of Fiji and Samoa, motivating future biogeographic investigation of this poorly known region.

Evolutionary switch and genetic convergence on rbcL following the evolution of C4 photosynthesis

Rubisco is responsible for the fixation of CO2 into organic compounds through photosynthesis and thus has a great agronomic importance. It is well established that this enzyme suffers from a slow catalysis, and its low specificity results into photorespiration, which is considered as an energy waste for the plant. However, natural variations exist, and some Rubisco lineages, such as in C4 plants, exhibit higher catalytic efficiencies coupled to lower specificities. These C4 kinetics could have evolved as an adaptation to the higher CO2 concentration present in C4 photosynthetic cells. In this study, using phylogenetic analyses on a large data set of C3 and C4 monocots, we showed that the rbcL gene, which encodes the large subunit of Rubisco, evolved under positive selection in independent C4 lineages. This confirms that selective pressures on Rubisco have been switched in C4 plants by the high CO2 environment prevailing in their photosynthetic cells. Eight rbcL codons evolving under positive selection in C4 clades were involved in parallel changes among the 23 independent monocot C4 lineages included in this study. These amino acids are potentially responsible for the C4 kinetics, and their identification opens new roads for human-directed Rubisco engineering. The introgression of C4-like high-efficiency Rubisco would strongly enhance C3 crop yields in the future CO2-enriched atmosphere.

Phylogenetic and expression analysis of the basic helix-loop-helix transcription factor gene family: genomic approach to cellular differentiation

A phylogenetic analysis of seven different species (human, mouse, rat, worm, fly, yeast, and plant) utilizing all (541) basic helix-loop-helix (bHLH) genes identified, including expressed sequence tags (EST), was performed. A super-tree involving six clades and a structural categorization involving the entire coding sequence was established. A nomenclature was developed based on clade distribution to discuss the functional and ancestral relationships of all the genes. The position/location of specific genes on the phylogenetic tree in relation to known bHLH factors allows for predictions of the potential functions of uncharacterized bHLH factors, including EST's. A genomic analysis using microarrays for four different mouse cell types (i.e. Sertoli, Schwann, thymic, and muscle) was performed and considered all known bHLH family members on the microarray for comparison. Cell-specific groups of bHLH genes helped clarify those bHLH genes potentially involved in cell specific differentiation. This phylogenetic and genomic analysis of the bHLH gene family has revealed unique aspects of the evolution and functional relationships of the different genes in the bHLH gene family.

Structural, biochemical, and physiological characterization of photosynthesis in two C4 subspecies of Tecticornia indica and the C3 species Tecticornia pergranulata (Chenopodiaceae)

Among dicotyledon families, Chenopodiaceae has the most C(4) species and the greatest diversity in structural forms of C(4). In subfamily Salicornioideae, C(4) photosynthesis has, so far, only been found in the genus Halosarcia which is now included in the broadly circumscribed Tecticornia. Comparative anatomical, cytochemical, and physiological studies on these taxa, which have near-aphyllous photosynthetic shoots, show that T. pergranulata is C(3), and that two subspecies of T. indica (bidens and indica) are C(4) (Kranz-tecticornoid type). In T. pergranulata, the stems have two layers of chlorenchyma cells surrounding the centrally located water storage tissue. The two subspecies of T. indica have Kranz anatomy in reduced leaves and in the fleshy stem cortex. They are NAD-malic enzyme-type C(4) species, with mesophyll chloroplasts having reduced grana, characteristic of this subtype. The Kranz-tecticornoid-type anatomy is unique among C(4) types in the family in having groups of chlorenchymatous cells separated by a network of large colourless cells (which may provide mechanical support or optimize the distribution of radiation in the tissue), and in having peripheral vascular bundles with the phloem side facing the bundle sheath cells. Also, the bundle sheath cells have chloroplasts in a centrifugal position, which is atypical for C(4) dicots. Fluorescence analyses in fresh sections indicate that all non-lignified cell walls have ferulic acid, a cell wall cross-linker. Structural-functional relationships of C(4) photosynthesis in T. indica are discussed. Recent molecular studies show that the C(4) taxa in Tecticornia form a monophyletic group, with incorporation of the Australian endemic genera of Salicornioideae, including Halosarcia, Pachycornia, Sclerostegia, and Tegicornia, into Tecticornia.

Phylogenetic analysis of the insulin-like growth factor binding protein (IGFBP) and IGFBP-related protein gene families

Insulin-like growth factor (IGF) activity is regulated by six high affinity binding proteins (IGFBPs) and possibly by some of the nine IGFBP-related proteins (IGFBP-rPs). To determine the phylogenetic relationship of this proposed gene superfamily, we conducted maximum likelihood (ML) and Bayesian inference analyses on a matrix of amino acid sequences from a diversity of vertebrate species. A single most likely phylogram, ML bootstrap, and Bayesian consensus tree of 10,000,000 generations revealed a monophyletic IGFBP lineage independent of the IGFBP-rPs. The IGFBPs segregated into three distinct clades: IGFBP-1, -3, and -6. Subsequent gene duplication events within the IGFBP-1 and -3 clades resulted in the production and divergence of IGFBP-2 and -4 within the IGFBP-1 clade and IGFBP-5 in the IGFBP-3 clade. By contrast, the IGFBP-rPs were distributed paraphyletically into two clades: IGFBP-rP1, 5, and 6 in one clade and the CCN family (IGFBP-rP2-4,7-9) in another. A recently identified IGFBP-3 homolog in rainbow trout localized to the IGFBP-2 subclade. Subsequence analysis identified a RGD motif common to IGFBP-2 orthologs, but did not identify the nuclear localization sequence present in IGFBP-3 and -5 homologs. The putative trout IGFBP-3 was 36-55% identical to different IGFBP-2 proteins, but only 24-27% identical to IGFBP-3 proteins. These results suggest that the IGFBPs and IGFBP-rPs are at best distantly related and that the limited similarities likely resulted from exon shuffling. They also suggest that rainbow trout, and possibly other salmonids, possess two IGFBP-2 paralogs as the putative trout IGFBP-3 is misannotated.

C4 Photosynthesis: Convergence upon Convergence upon…

C(4) photosynthesis has evolved independently numerous times in plants. New work suggests that these multiple origins are the result of recurring selection on a few amino acid positions in the key enzyme phosphoenolpyruvate carboxylase.

Diversity and plasticity of C4 photosynthesis in Eleocharis (Cyperaceae)

Eleocharis contains many amphibious species, and displays diversity of photosynthetic mechanism (C_3, C₄ or C₃-C₄ intermediates). A unique feature of Eleocharis is the plasticity in the photosynthetic mechanism of some species in response to the environment. In this study, we have examined the culm anatomy and photosynthetic property of several Eleocharis species grown terrestrially and the changes in the newly produced culms over a short period time frame after switching from terrestrial to submerged condition. Eleocharis baldwinii (Torrey) Chapman is C₄-like in terrestrial habitat, exhibiting O₂ inhibition of photosynthesis with Rubisco expressed in both mesophyll and bundle sheath cells and PEPC strictly in the mesophyll cells, but switches to C₃-C₄ intermediacy when submerged. In addition to Eleocharis vivipara Link type 1 (which switches from C₄-like to C₃), two other photosynthetic types examined in this study were shown to have different responses to growth in either terrestrial or submerged conditions. E. vivipara type 2 is a typical C₄ plant in the terrestrial habitat, but becomes a C₃-C₄ intermediate under submerged conditions. Further, terrestrially, E. vivipara type 3 is a C₃-C₄ intermediate, but when submerged the δ¹³C value increases to -6.7‰, indicating its use of bicarbonate as a major carbon source. The submerged form of this plant exhibited about three times higher photosynthetic O₂ evolution rate, compared to the C₃ species Eleocharis erythropoda Steudel. These Eleocharis species possess different molecular switches for regulating C₄ gene expression in response to environmental stimuli both between different species, and in E. vivipara among different populations. The apparent expression of a bicarbonate transport system by E. vivipara type 3 while submerged represents a unique adaptation to low CO₂ availability.

Characterization of rainbow trout myostatin-2 genes (rtMSTN-2a and -2b): genomic organization, differential expression, and pseudogenization

Myostatin is an extremely potent negative regulator of vertebrate skeletal muscle development. A phylogenetic analysis suggests that salmonids should possess four distinct genes, although only MSTN-1 orthologs have been characterized. Described herein are the rainbow trout (rt) MSTN-2a and -2b genes and subsequence analysis of their promoters and their quantitative expression profiles. Both genes are similarly organized, contain several putative myogenic response elements, and are legitimate MSTN-2 orthologs based on Bayesian analyses. However, rtMSTN-2b contains two in-frame stop codons within the first exon and unspliced variants of both transcripts were expressed in a tissue-specific manner. Complete splicing of rtMSTN-2a occurred only in brain, where expression is highest, whereas rtMSTN-2b transcripts were mostly present in unspliced forms. The presence of stop codons in the rtMSTN-2b open reading frame and the expression of mostly unspliced transcripts indicate that this particular homolog is a pseudogene. These results confirm our previous phylogenetic analysis and suggest that all salmonids likely possess four distinct myostatin genes. The tissue-specific expression and differential processing of both rtMSTN-2 transcripts as well the pseudogenization of rtMSTN-2b may reflect compensatory and adaptive responses to tetraploidization and may help limit rtMSTN-2a's influences primarily to neural tissue.

Ecological speciation in the East Maui-endemic Dubautia (Asteraceae) species

Adaptive diversification can be viewed as a balance between the conservative force of interpopulation gene flow and selection for differential environments. In this paper, we examine ecological, morphological, and genetic differentiation in a small clade consisting of four East Maui-endemic species of Dubautia: D. menziesii, D. platyphylla, D. reticulata, and D. waianapanapaensis, in the Hawaiian silversword alliance (Asteraceae). The East Maui clade is apparently recently derived (less than 1 million years ago) and is geographically restricted yet displays significant ecological and morphological differences. We used geographic data from historical herbarium specimens, measurements of plant architecture and leaf morphometrics, and measures of genetic differentiation in both microsatellite and nuclear coding loci to examine the correlation of different forms of divergence in this small species flock. We found overlap in large-scale geographic distributions, significant differentiation in most habitat factors, significant microsatellite differentiation, and many shared alleles at nuclear coding loci suggesting on-going lineage sorting. Despite the presence of apparent hybrids in some populations, microsatellite variation is consistent with isolation among species. Using Mantel tests, we compared the direction and extent of diversification among different datasets, to determine whether ecological/morphological divergence was correlated with genetic divergence. Correlations among different datasets showed that habitat was strongly correlated with plant architecture but not leaf morphology. Taken together, these results indicate that ecological and morphological diversification has driven genetic divergence at rapidly evolving microsatellite loci, whereas there is continuing lineage sorting at neutral sites in nuclear coding loci.

Genetic divergence and migration patterns in a North American passerine bird: implications for evolution and conservation

Like many other migratory birds, the black-throated blue warbler (Dendroica caerulescens) shows pronounced differences in migratory behaviour and other traits between populations: birds in the southern part of the breeding range have darker plumage and migrate to the eastern Caribbean during the winter, whereas those in the north have lighter plumage and migrate to the western Caribbean. We examined the phylogeography of this species, using samples collected from northern and southern populations, to determine whether differentiation between these populations dates to the Pleistocene or earlier, or whether differences in plumage and migratory behaviour have arisen more recently. We analysed variation at 369 bp of the mitochondrial control region domain I and also at seven nuclear microsatellites. Analyses revealed considerable genetic variation, but the vast majority of this variation was found within rather than between populations, and there was little differentiation between northern and southern populations. Phylogeographic analyses revealed a very shallow phylogenetic tree, a star-like haplotype network, and a unimodal mismatch distribution, all indicative of a recent range expansion from a single refugium. Coalescent modelling approaches also indicated a recent common ancestor for the entire group of birds analysed, no split between northern and southern populations, and high levels of gene flow. These results show that Pleistocene or earlier events have played little role in creating differences between northern and southern populations, suggesting that migratory and other differences between populations have arisen very recently. The implications of these results for the evolution of migration and defining taxonomic groups for conservation efforts are discussed.

Phylogenetic analysis of the myostatin gene sub-family and the differential expression of a novel member in zebrafish

The myostatin (MSTN)-null phenotype in mammals is characterized by extreme gains in skeletal muscle mass or "double muscling" as the cytokine negatively regulates skeletal muscle growth. Recent attempts, however, to reproduce a comparable phenotype in zebrafish have failed. Several aspects of MSTN biology in the fishes differ significantly from those in mammals and at least two distinct paralogs have been identified in some species, which possibly suggests functional divergence between the different vertebrate classes or between fish paralogs. We therefore conducted a phylogenetic analysis of the entire MSTN gene sub-family. Maximum likelihood, Bayesian inference, and bootstrap analyses indicated a monophyletic distribution of all MSTN genes with two distinct fish clades: MSTN-1 and -2. These analyses further indicated that all Salmonid genes described are actually MSTN-1 orthologs and that additional MSTN-2 paralogs may be present in most, if not all, teleosts. An additional zebrafish homolog was identified by BLAST searches of the zebrafish Hierarchical Tets Generation System database and was subsequently cloned. Comparative sequence analysis of both genes (zebrafish MSTN (zfMSTN)-1 and -2) revealed many differences, primarily within the latency-associated peptide regions, but also within the bioactive domains. The 2-kb promoter region of zfMSTN-2 contained many putative cis regulatory elements that are active during myogenesis, but are lacking in the zfMSTN-1 promoter. In fact, zfMSTN-2 expression was limited to the early stages of somitogenesis, whereas zfMSTN-1 was expressed throughout embryogenesis. These data suggest that zfMSTN-2 may be more closely associated with skeletal muscle growth and development. They also resolve the previous ambiguity in classification of fish MSTN genes.

Repeated evolution of net venation and fleshy fruits among monocots in shaded habitats confirms a priori predictions: evidence from an ndhF phylogeny

We present a well-resolved, highly inclusive phylogeny for monocots, based on ndhF sequence variation, and use it to test a priori hypotheses that net venation and vertebrate-dispersed fleshy fruits should undergo concerted convergence, representing independent but often concurrent adaptations to shaded conditions. Our data demonstrate that net venation arose at least 26 times and was lost eight times over the past 90 million years; fleshy fruits arose at least 21 times and disappeared 11 times. Both traits show a highly significant pattern of concerted convergence (p<10(-9)), arising 16 times and disappearing four times in tandem. This phenomenon appears driven by even stronger tendencies for both traits to evolve in shade and be lost in open habitats (p<10(-13)-10(-29)). These patterns are among the strongest ever demonstrated for evolutionary convergence in individual traits and the predictability of evolution, and the strongest evidence yet uncovered for concerted convergence. The rate of adaptive shifts per taxon has declined exponentially over the past 90 million years, as expected when large-scale radiations fill adaptive zones.

Phylogenetic analysis of the nuclear alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) and combined analyses of Adh and nuclear ribosomal ITS and ETS sequences for inferring species relationships

We analyzed sequence variation for the alcohol dehydrogenase (Adh) gene family in Carex section Acrocystis (Cyperaceae) to reconstruct Adh gene trees for Acrocystis species and to characterize the structure of the Adh gene family in Carex. Two Adh loci were included with ITS and ETS sequences in a combined Bayesian inference analysis of Carex section Acrocystis to gain a better understanding of species relationships in the section. In addition, we comment on how the results presented here contribute to our knowledge of the birth-death process of the Adh gene family in angiosperms. It appears that the structure of the Adh gene family in Carex is complex with possibly six loci present in the gene family. Additionally, variation among Acrocystis species within loci is quite low, and there is little phylogenetic resolution in the individual datasets. Bayesian inference analysis of the combined ITS, ETS, Adh1, and Adh2 datasets resulted in a moderately well-supported phylogenetic hypothesis of relationships in the section which is discussed in relation to previous hypotheses of relationships.