Accelerated regulatory gene evolution in an adaptive radiation

Phylogenetic relationships and the evolution of regulatory gene sequences in the parrotfishes

Regulatory genes control the expression of other genes and are key components of developmental processes such as segmentation and embryonic construction of the skull in vertebrates. Here we examine the variability and evolution of three vertebrate regulatory genes, addressing issues of their utility for phylogenetics and comparing the rates of genetic change seen in regulatory loci to the rates seen in other genes in the parrotfishes. The parrotfishes are a diverse group of colorful fishes from coral reefs and seagrasses worldwide and have been placed phylogenetically within the family Labridae. We tested phylogenetic hypotheses among the parrotfishes, with a focus on the genera Chlorurus and Scarus, by analyzing eight gene fragments for 42 parrotfishes and eight outgroup species. We sequenced mitochondrial 12s rRNA (967 bp), 16s rRNA (577 bp), and cytochrome b (477 bp). From the nuclear genome, we sequenced part of the protein-coding genes rag2 (715 bp), tmo4c4 (485 bp), and the developmental regulatory genes otx1 (672 bp), bmp4 (488bp), and dlx2 (522 bp). Bayesian, likelihood, and parsimony analyses of the resulting 4903 bp of DNA sequence produced similar topologies that confirm the monophyly of the scarines and provide a phylogeny at the species level for portions of the genera Scarus and Chlorurus. Four major clades of Scarus were recovered, with three distributed in the Indo-Pacific and one containing Caribbean/Atlantic taxa. Molecular rates suggest a Miocene origin of the parrotfishes (22 mya) and a recent divergence of species within Scarus and Chlorurus, within the past 5 million years. Developmentally important genes made a significant contribution to phylogenetic structure, and rates of genetic evolution were high in bmp4, similar to other coding nuclear genes, but low in otx1 and the dlx2 exons. Synonymous and non-synonymous substitution patterns in developmental regulatory genes support the hypothesis of stabilizing selection during the history of these genes, with several phylogenetic regions of accelerated non-synonymous change detected in the phylogeny.

Natural variation in gene expression between wild and weedy populations of Helianthus annuus

The molecular genetic changes underlying the transformation of wild plants into agricultural weeds are poorly understood. Here we use a sunflower cDNA microarray to detect variation in gene expression between two wild (non-weedy) Helianthus annuus populations from Utah and Kansas and four weedy H. annuus populations collected from agricultural fields in Utah, Kansas, Indiana, and California. When grown in a common growth chamber environment, populations differed substantially in their gene expression patterns, indicating extensive genetic differentiation. Overall, 165 uni-genes, representing approximately 5% of total genes on the array, showed significant differential expression in one or more weedy populations when compared to both wild populations. This subset of genes is enriched for abiotic/biotic stimulus and stress response proteins, which may underlie niche transitions from the natural sites to agricultural fields for H. annuus. However, only a small proportion of the differentially expressed genes overlapped in multiple wild vs. weedy comparisons, indicating that most of the observed expression changes are due to local adaptation or neutral processes, as opposed to parallel genotypic adaptation to agricultural fields. These results are consistent with an earlier phylogeographic study suggesting that weedy sunflowers have evolved multiple times in different regions of the United States and further indicate that the evolution of weedy sunflowers has been accompanied by substantial gene expression divergence in different weedy populations.

Allohexaploidy, introgression, and the complex phylogenetic history of Elymus repens (Poaceae)

The phylogenetic position of hexaploid Elymus repens within the tribe Triticeae (Poaceae) was examined using cloned sequences from the low-copy nuclear genes encoding phosphoenolpyruvate carboxylase (pepC) and beta-amylase. A previous analysis of E. repens using data from the nuclear granule-bound starch synthase I (GBSSI) gene had yielded five phylogenetically distinct gene copies, two more than expected from hexaploidy alone. The three gene trees share three distinct E. repens clades, suggesting that E. repens contains three phylogenetically divergent genomes, contributed by Hordeum, Pseudoroegneria, and an unknown donor. The two additional GBSSI sequences, including one that was apparently derived from outside of the tribe, appear to reflect past introgression of GBSSI sequences into the E. repens genome. On all three trees, the Hordeum-like E. repens sequences are polyphyletic within Hordeum, and the trees are in conflict with regard to the placement of these sequences within Hordeum, highlighting multiple contributions from Hordeum to E. repens.

Diversity and divergence patterns in regulatory genes suggest differential gene flow in recently derived species of the Hawaiian silversword alliance adaptive radiation (Asteraceae)

The impact of gene flow and population size fluctuations in shaping genetic variation during adaptive radiation, at both the genome-wide and gene-specific levels, is very poorly understood. To examine how historical population size and gene flow patterns within and between loci have influenced lineage divergence in the Hawaiian silversword alliance, we have investigated the nucleotide sequence diversity and divergence patterns of four floral regulatory genes (ASAP1-A, ASAP1-B, ASAP3-A, ASAP3-B) and a structural gene (ASCAB9). Levels and patterns of molecular divergence across these five nuclear loci were estimated between two recently derived species (Dubautia ciliolata and Dubautia arborea) which are presumed to be sibling species. This multilocus analysis of genetic variation, haplotype divergence and historical demography indicates that population expansion and differential gene flow occurred subsequent to the divergence of these two lineages. Moreover, contrasting patterns of allele- sharing for regulatory loci vs. a structural locus between these two sibling species indicate alternative histories of genetic variation and partitioning among loci where alleles of the floral regulatory loci are shared primarily from D. arborea to D. ciliolata and alleles of the structural locus are shared in both directions. Taken together, these results suggest that adaptively radiating species can exhibit contrasting allele migration rates among loci such that allele movement at specific loci may supersede genetic divergence caused by drift and that lineage divergence during adaptive radiation can be associated with population expansion.

Molecular evolution and genome divergence at RPB2 gene of the St and H genome in Elymus species

Molecular evolution of the second largest subunit of low copy nuclear RNA polymerase II (RPB2) in allotetrploid StH genomic species of Elymus is characterized here. Our study first reported a 39-bp MITE stowaway element insertion in the genic region of RPB2 gene for all tetraploid Elymus St genome and diploid Pseudoroegneria spicata and P. stipifolia St genome. The sequences on 3'-end are highly conserved, with AGTA in all sequences but H10339 (E. fibrosis), in which the AGTA was replaced with AGAA. All 12 Stowaway-containing sequences encompassed a 9 bp conserved TIRs (GAGGGAGTA). Interestingly, the 5'-end sequence of GGTA which was changed to AGTA or deleted resulted in Stowaway excision in the H genome of Elymus sepcies, in which Stowaway excision did not leave footprint. Another two large insertions in all St genome sequences are also transposable-like elements detected in the genic region of RPB2 gene. Our results indicated that these three transposable element indels have occurred prior to polyploidization, and shaped the homoeologous RPB2 loci in St and H genome of Eymus species. Nucleotide diversity analysis suggested that the RPB2 sequence may evolve faster in the polyploid species than in the diploids. Higher level of polymorphism and genome-specific amplicons generated by this gene indicated that RPB2 is an excellent tool for investigating the phylogeny and evolutionary dynamics of speciation, and the mode of polyploidy formation in Elymus species.

Breakdown of distyly in a tetraploid variety of Ophiorrhiza japonica (Rubiaceae) and its phylogenetic analysis

We examined the floral morph of tetraploid Ophiorrhiza japonica Blume var. amamiana Hatus. and diploid O. japonica var. japonica to elucidate the association of distyly and ploidy levels. Chloroplast DNA phylogeny was reconstructed to determine the number of tetraploidization events and floral morph shifts in O. japonica. All individuals of O. japonica var. amamiana proved to be long-homostylous, whereas O. japonica var. japonica was distylous with typical long- and short-styled flowers. Distyly is related to the ploidy level. The bagging treatment of flowers indicated that O. japonica var. amamiana is self-compatible and potentially automatically self-pollinating. In cpDNA sequencing analysis, no haplotype was shared between the two varieties. The cpDNA haplotype network displayed the monophyly of O. japonica var. amamiana, suggesting a single origin of this variety. Hence, both tetraploidization and the breakdown of distyly to homostyly in O. japonica var. amamiana likely occurred just once. Because O. japonica var. amamiana having the morphological and cytological entity is recognized as a single lineage and clearly separated from O. japonica var. japonica, this variety can be considered to be a distinct species. We therefore propose to raise O. japonica var. amamiana to the rank of species.

What makes species unique? The contribution of proteins with obscure features

An analysis of proteins with obscure features in ten eukaryotic genomes revealed that the majority are species-specific.

Background

Proteins with obscure features (POFs), which lack currently defined motifs or domains, represent between 18% and 38% of a typical eukaryotic proteome. To evaluate the contribution of this class of proteins to the diversity of eukaryotes, we performed a comparative analysis of the predicted proteomes derived from 10 different sequenced genomes, including budding and fission yeast, worm, fly, mosquito, Arabidopsis, rice, mouse, rat, and human.

Results

Only 1,650 protein groups were found to be conserved among these proteomes (BLAST E-value threshold of 10^-6). Of these, only three were designated as POFs. Surprisingly, we found that, on average, 60% of the POFs identified in these 10 proteomes (44,236 in total) were species specific. In contrast, only 7.5% of the proteins with defined features (PDFs) were species specific (17,554 in total). As a group, POFs appear similar to PDFs in their relative contribution to biological functions, as indicated by their expression, participation in protein-protein interactions and association with mutant phenotypes. However, POF have more predicted disordered structure than PDFs, implying that they may exhibit preferential involvement in species-specific regulatory and signaling networks.

Conclusion

Because the majority of eukaryotic POFs are not well conserved, and by definition do not have defined domains or motifs upon which to formulate a functional working hypothesis, understanding their biochemical and biological functions will require species-specific investigations.

Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera

As a model system, Arabidopsis thaliana and its wild relatives have played an important role in the study of genomics and evolution in plants. In this study, we examined the genetic diversity of the chalcone synthase (Chs) gene, which encodes a key enzyme of the flavonoid pathway and is located on chromosome five, as well as two Chs-like genes on the first and fourth chromosomes of Arabidopsis. The objectives of the study are to determine if natural selection operates differentially on the paralogs of the Chs gene family in A. thaliana and Arabidopsis halleri ssp. gemmifera. The mode of selection was inferred from Tajima's D values from noncoding and coding regions, as well as from the ratio of nonsynonymous to synonymous substitutions. Both McDonald-Kreitman and HKA tests revealed the effects of selection on the allelic distribution, except for the chromosome 1 paralog in ssp. gemmifera. The Chs gene on chromosome 5 was under purifying selection in both species. Significant, negative Tajima's D values at synonymous sites and positive Fay and Wu's H values within coding region, plus reduced genetic variability in introns, indicated effects of background selection in shaping the evolution of this gene region in A. thaliana. The Chs paralog on chromosome 1 was under positive selection in A. thaliana, while interspecific introgression and balancing selection determined the fates of the paralog and resulted in high heterogeneity in ssp. gemmifera. Local adaptation differentiated populations of Japan and China at the locus. In contrast, the other Chs-paralog of chromosome 4 was shaped by purifying selection in A. thaliana, while under positive selection in ssp. gemmifera, as indicated by dn/ds>1. Moreover, these contrasting patterns of selection have likely resulted in functional divergence in Arabidopsis, as indicated by radical amino acid substitutions at the chalcone synthase/stilbene synthase motif of the Chs genes. Unlike previous studies of the evolutionary history of A. thaliana, the high levels of genetic diversity in most gene regions of Chs paralogs and nonsignificant Tajima's D in the intron sequences of the Chs gene family in A. thaliana did not reflect the effects of a recent demographic expansion.

Molecular adaptation during adaptive radiation in the Hawaiian endemic genus Schiedea

BACKGROUND

"Explosive" adaptive radiations on islands remain one of the most puzzling evolutionary phenomena. The rate of phenotypic and ecological adaptations is extremely fast during such events, suggesting that many genes may be under fairly strong selection. However, no evidence for adaptation at the level of protein coding genes was found, so it has been suggested that selection may work mainly on regulatory elements. Here we report the first evidence that positive selection does operate at the level of protein coding genes during rapid adaptive radiations. We studied molecular adaptation in Hawaiian endemic plant genus Schiedea (Caryophyllaceae), which includes closely related species with a striking range of morphological and ecological forms, varying from rainforest vines to woody shrubs growing in desert-like conditions on cliffs. Given the remarkable difference in photosynthetic performance between Schiedea species from different habitats, we focused on the "photosynthetic" Rubisco enzyme, the efficiency of which is known to be a limiting step in plant photosynthesis.

RESULTS

We demonstrate that the chloroplast rbcL gene, encoding the large subunit of Rubisco enzyme, evolved under strong positive selection in Schiedea. Adaptive amino acid changes occurred in functionally important regions of Rubisco that interact with Rubisco activase, a chaperone which promotes and maintains the catalytic activity of Rubisco. Interestingly, positive selection acting on the rbcL might have caused favorable cytotypes to spread across several Schiedea species.

SIGNIFICANCE

We report the first evidence for adaptive changes at the DNA and protein sequence level that may have been associated with the evolution of photosynthetic performance and colonization of new habitats during a recent adaptive radiation in an island plant genus. This illustrates how small changes at the molecular level may change ecological species performance and helps us to understand the molecular bases of extremely fast rate of adaptation during island adaptive radiations.

Convergence, constraint and the role of gene expression during adaptive radiation: floral anthocyanins in Aquilegia

Convergent phenotypes are testament to the role of natural selection in evolution. However, little is known about whether convergence in phenotype extends to convergence at the molecular level. We use the independent losses of floral anthocyanins in columbines (Aquilegia) to determine the degree of molecular convergence in gene expression across the anthocyanin biosynthetic pathway (ABP). Using a phylogeny of the North American Aquilegia clade, we inferred six independent losses of floral anthocyanins. Via semiquantitative reverse transcriptase-polymerase chain reaction (RT-PCR), we monitored developmental and tissue-specific variation in expression of the six major structural ABP loci in three Aquilegia species, two that produce anthocyanins (A+) and one that does not (A-). We then compared ABP expression in petals of old-bud and pre-anthesis flowers of 13 Aquilegia species, eight wild species and two horticultural lines representing seven independent A- lineages as well as three wild A+ species. We only found evidence of down-regulation of ABP loci in A- lineages and losses of expression were significantly more prevalent for genes late in the pathway. Independent contrast analysis indicates that changes in expression of dihydroflavonol reductase (DFR) and anthocyanidin synthase (ANS) are strongly phylogenetically correlated consistent with the multilocus targets of trans-regulatory elements in the ABP of other systems. Our findings strongly suggest that pleiotropy constrains the evolution of loss of floral anthocyanins to mutations affecting genes late in the ABP mostly through convergent changes in regulatory genes. These patterns support the hypothesis that rapid evolutionary change occurs largely through regulatory rather than structural mutations.

Adaptive radiation in microbial microcosms

It has often been argued that evolutionary diversification is the result of divergent natural selection for specialization on alternative resources. I provide a comprehensive review of experiments that examine the ecology and genetics of resource specialization and adaptive radiation in microbial microcosms. In these experiments, resource heterogeneity generates divergent selection for specialization on alternative resources. At a molecular level, the evolution of specialization is generally attributable to mutations that de-regulate the expression of existing biosynthetic and catabolic pathways. Trade-offs are associated with the evolution of resource specialization, but these trade-offs are often not the result of antagonistic pleiotropy. Replicate adaptive radiations result in the evolution of a similar assemblage of specialists, but the genetic basis of specialization differs in replicate radiations. The implications of microbial selection experiments for evolutionary theory are discussed and future directions of research are proposed.

Allopolyploid origin of Cardamine asarifolia (Brassicaceae): incongruence between plastid and nuclear ribosomal DNA sequences solved by a single-copy nuclear gene

Interspecific hybridization and polyploidization have played central roles in plant diversification. However, technical difficulties in the analyses of low-copy genes have limited the study of the origins of hybrid and polyploid plants. Here, we present a phylogenetic analysis of the hexaploid Cardamine asarifolia, distributed in the southern European Alps and northern Apennines. Our study included all relevant taxa of the genus found in Europe. A marked discrepancy was revealed between the trnL-trnF region of cpDNA and internal transcribed spacer (nrDNA ITS) sequences. To solve the incongruence, we sequenced a single-copy nuclear CHS gene (chalcone synthase) using a novel method to design homoeologue-specific PCR primers to bypass artefacts caused by artificial recombination of homoeologues during PCR and/or cloning. Three homoeologues were isolated from C. asarifolia, providing evidence for its allopolyploid origin. One homoeologue, showing the same phylogenetic position as the ITS sequences, most likely originated from an extinct parent. Furthermore, we documented recurrent polytopic hybridizations between C. asarifolia and diploid C. amara. The allohexaploidization and the following hybridization with a diploid species exemplify the ongoing dynamic processes of speciation in the genus Cardamine.

Sequence evolution and sex-specific expression patterns of the C class floral identity gene, SpAGAMOUS, in dioecious Spinacia oleracea L

Development in dioecious cultivated spinach, Spinacia oleracea, is distinguished by the absence of alternative reproductive organ primordia in male and female flowers. Given the highly derived floral developmental program in spinach, we wished to characterize a spinach C class floral identity gene and to determine the patterns of sequence evolution as well as compare the spatial and temporal expression patterns with those of AGAMOUS. The isolated cDNA sequence clusters phylogenetically within the AGAMOUS/FARINELLI C class clade. In comparison with the SLM1 sequence from the related Silene latifolia, amino acid replacements are highly conservative and non-randomly distributed, being predominantly found in hinge regions or on exposed surfaces of helices. The spinach gene (SpAGAMOUS) appears to be exclusively expressed in reproductive tissues and not in vegetative organs. Initial expression of SpAGAMOUS is similar in male and female floral primordia. However, upon initiation of the first whorl organs, SpAGAMOUS becomes restricted to meristemic regions from which the reproductive primordia will develop. This results in an early gender-specific pattern. Thus, the spinach C class gene is differentially expressed prior to reproductive organ development and is, at least, correlated with, if not directly involved in, the sexual dimorphism in spinach.

Growth form evolution and shifting habitat specialization in annual plants

Optimal plant growth form should vary across environments. We examined the potential for mutations causing large changes in growth form to produce new optimal phenotypes across light environments. We predicted that the upright growth form would be favoured in a light limiting environment as leaves were in a position to maximize light interception, while a rosette (leaves in a basal position) growth form would be favoured in a high light environment. Growth form genotypes of Brassica rapa (upright wild-type and rosette mutants) and Arabidopsis thaliana (large rosette wild-type and increasingly upright growth form mutants) were grown in a greenhouse in control (ambient) and filtered (low) light treatments. Compared to upright genotypes, rosette genotypes had relatively high fitness in control light but had a relatively large fitness reduction in filtered light. Our results demonstrate the potential importance of rapid growth form evolution in plant adaptation to new or changing environments.

Neutral evolution of the nonbinding region of the anthocyanin regulatory gene Ipmyb1 in Ipomoea

Plant transcription factors often contain domains that evolve very rapidly. Although it has been suggested that this rapid evolution may contribute substantially to phenotypic differentiation among species, this suggestion has seldom been tested explicitly. We tested the validity of this hypothesis by examining the rapidly evolving non-DNA-binding region of an R2R3-myb transcription factor that regulates anthocyanin expression in flowers of the genus Ipomoea. We first provide evidence that the W locus in Ipomoea purpurea, which determines whether flowers will be pigmented or white, corresponds to a myb gene segregating in southeastern U.S. populations for one functional allele and one nonfunctional allele. While the binding domain exhibits substantial selective constraint, the nonbinding region evolves at an average K(a)/K(s) ratio of 0.74. This elevated rate of evolution is due to relaxed constraint rather than to increased levels of positive selection. Despite this relaxed constraint, however, approximately 20-25% of the codons, randomly distributed throughout the nonbinding region, are highly constrained, with the remainder evolving neutrally, indicating that the entire region performs important function(s). Our results provide little indication that rapid evolution in this regulatory gene is driven by natural selection or that it is responsible for floral-color differences among Ipomoea species.

Heterogeneous evolution of the Myc-like Anthocyanin regulatory gene and its phylogenetic utility in Cornus L. (Cornaceae)

Anthocyanin is a major pigment in vegetative and floral organs of most plants and plays an important role in plant evolution. The anthocyanin regulatory genes are responsible for regulating transcription of genes in the anthocyanin synthetic pathway. To assess evolutionary significance of sequence variation and evaluate the phylogenetic utility of an anthocyanin regulatory gene, we compared nucleotide sequences of the myc-like anthocyanin regulatory gene in the genus of dogwoods (Cornus: Cornaceae). Phylogenetic analyses demonstrate that the myc-like anthocyanin regulatory gene has potential as an informative phylogenetic marker at different taxonomic levels, depending on the data set considered (DNA or protein sequences) and regions applied (exons or introns). Pairwise nonsynonymous and synonymous substitution rate tests and codon-based substitution models were applied to characterize variation and to identify sites under diversifying selection. Mosaic evolution and heterogeneous rates among different domains and sites were detected.

Do plants and animals differ in phenotypic plasticity?

This paper compares the flexibility in the nexus between phenotype and genotype in plants and animals. These taxa although considered to be fundamentally different are found to be surprisingly similar in the mechanisms used to achieve plasticity. Although non-cognitive behaviour occurs in plants, its range is limited, while morphological and developmental plasticity also occur to a considerable extent in animals. Yet both plants and animals are subject to unique constraints and thus need to find unique solutions to functional problems. A true comparison between the plant and animal phenotype would be a comparison between plants and sessile photosynthesizing colonial invertebrates. Such comparisons are lacking. However, they would provide important insights into the adaptive significance of plasticity in these groups. It is also suggested that a comparison of inflexible traits in these groups would provide an understanding of the constraints, as well as the costs and benefits,of a plastic versus non-plastic phenotype in plants and animals.

Adaptive evolution of HoxA-11 and HoxA-13 at the origin of the uterus in mammals

The evolution of morphological characters is mediated by the evolution of developmental genes. Evolutionary changes can either affect cis-regulatory elements, leading to differences in their temporal and spatial regulation, or affect the coding region. Although there is ample evidence for the importance of cis-regulatory evolution, it has only recently been shown that transcription factors do not remain functionally equivalent during evolution. These results suggest that the evolution of transcription factors may play an active role in the evolution of development. To test this idea we investigated the molecular evolution of two genes essential for the development and function of the mammalian female reproductive organs, HoxA-11 and HoxA-13. We predicted that if coding-region evolution plays an active role in developmental evolution, then these genes should have experienced adaptive evolution at the origin of the mammalian female reproductive system. We report the sequences of HoxA-11 from basal mammalian and amniote taxa and analyse HoxA-11 and HoxA-13 for signatures of adaptive molecular evolution. The data demonstrate that these genes were under strong positive (directional) selection in the stem lineage of therian and eutherian mammals, coincident with the evolution of the uterus and vagina. These results support the idea that adaptive evolution of transcription factors can be an integral part in the evolution of novel structures.

Coalescent processes and relaxation of selective constraints leading to contrasting genetic diversity at paralogs AtHVA22d and AtHVA22e in Arabidopsis thaliana

Duplicate loci offer a very powerful system for understanding the complicated genome structure and adaptive evolution of a gene family. In this study, the genetic variation at paralogs AtHVA22d and AtHVA22e, members of an ABA- and stress-inducible gene family, is examined in the selfing Arabidopsis thaliana. Population genetic analysis indicates contrasting levels of nucleotide diversity at overall exon sequence and nonsynonymous sites between AtHVA22d (pi = 0.00337, pi(rep) = 0.00158) and AtHVA22e (pi = 0.00054, pi(rep) = 0.00023). The fact of Ka/Ks ratios significantly less than 1 in all sequences indicates that both genes are functional and subjected to purifying selection. In addition, rooted at barley HVA22, accelerated evolution is detected at replacement changes in the AtHVA22d locus, indicating relaxation of purifying selection after gene duplication. However, relative rate tests reveal no deviation from the neutrality at synonymous sites between the two paralogs. Based on clock-like evolution, the rate of synonymous substitution is estimated at 1.83 x 10(-9) substitutions per site per year; and the divergence of the two paralogs is traced to 90 MYA, coinciding with a period of the diversification of angiosperms. Given no codon usage bias in both genes, natural selection alone cannot account for the 6.4-fold differences in the nucleotide variation at synonymous sites between the two paralogs. Random processes resulting in different coalescence times, 3.65 MYA at AtHVA22d vs. 1.20 MYA at AtHVA22e, may have predominantly contributed to the evident differences of the genetic diversity. Partially nonoverlapping modes of expression between the two functional paralogs suggest a subfunctionalization hypothesis for explaining the fates of duplicate loci.

Changes in gene expression as biochemical adaptations to environmental change: a tribute to Peter Hochachka

Changes in gene expression are likely to play a critical role in both acclimation and adaptation to a changing environment. There is a rapidly growing body of literature implicating quantitative changes in gene expression during acclimation to environmental change, but less is known about the role of qualitative changes in gene expression, such as switching between alternative isoforms. Alternative isoforms can arise via gene duplication, alternative splicing, or alternative promoter usage. Organisms that have undergone recent genome duplication events may make use of environment-specific isoforms coded by multiple genes, but their role in other organisms is less well known. However, recent data suggest that isoforms arising from alternative splicing may be an under-appreciated source of physiological variation. The role of changes in gene expression during evolutionary adaptation has received comparatively limited attention, but novel approaches to addressing the adaptive significance of changes in gene expression have been applied to a few cases of differences in gene expression among taxa. Recent advances in genomics, including microarray technology, knock-out and knock-down approaches, and the wealth of data coming from large-scale sequencing projects have provided (and will continue to provide at ever increasing rates) new insights into these classic questions in comparative biochemistry.

Molecular evidence on plant divergence times

Estimation of divergence times from sequence data has become increasingly feasible in recent years. Conflicts between fossil evidence and molecular dates have sparked the development of new methods for inferring divergence times, further encouraging these efforts. In this paper, available methods for estimating divergence times are reviewed, especially those geared toward handling the widespread variation in rates of molecular evolution observed among lineages. The assumptions, strengths, and weaknesses of local clock, Bayesian, and rate smoothing methods are described. The rapidly growing literature applying these methods to key divergence times in plant evolutionary history is also reviewed. These include the crown group ages of green plants, land plants, seed plants, angiosperms, and major subclades of angiosperms. Finally, attempts to infer divergence times are described in the context of two very different temporal settings: recent adaptive radiations and much more ancient biogeographic patterns.

RATES OF EVOLUTION OF HYBRID INVIABILITY IN BIRDS AND MAMMALS

Almost 30 years ago, A. C. Wilson and colleagues presented results indicating that hybrid inviability between species evolves 10 times faster in mammals than in birds and frogs. Here I revisit this question for birds and mammals using modern molecular data (mitochondrial cytochrome b DNA) and a more phylogenetically appropriate statistical approach. My analyses confirm that diverging mammals lose the ability to form viable hybrids faster than birds. To explain the difference in rates of evolutionary loss of hybridization potential, Wilson and coworkers proposed that mammals have higher rates of regulatory evolution, causing higher probabilities of developmental incompatibilities between mammal species. I briefly discuss this and other potential explanations.

Phylogeny and divergence of basal angiosperms inferred from APETALA3- and PISTILLATA-like MADS-box genes

The B-class MADS-box genes composed of APETALA3 ( AP3) and PISTILLATA ( PI) lineages play an important role in petal and stamen identity in previously studied flowering plants. We investigated the diversification of the AP3-like and PI-like MADS-box genes of eight species in five basal angiosperm families: Amborella trichopoda (Amborellaceae); Brasenia schreberi and Cabomba caroliniana (Cabombaceae); Euryale ferox, Nuphar japonicum, and Nymphaea tetragona (Nymphaeaceae); Illicium anisatum (Illiciaceae); and Kadsura japonica (Schisandraceae). Sequence analysis showed that a four amino acid deletion in the K domain, which was found in all previously reported angiosperm PI genes, exists in a PI homologue of Schisandraceae, but not in six PI homologues of the Amborellaceae, Cabombaceae, and Nymphaeaceae, suggesting that the Amborellaceae, Cabombaceae, and Nymphaeaceae are basalmost lineages in angiosperms. The results of molecular phylogenetic analyses were not inconsistent with this hypothesis. The AP3 and PI homologues from Amborella share a sequence of five amino acids in the 5' region of exon 7. Using the linearized tree and likelihood methods, the divergence time between the AP3 and PI lineages was estimated as somewhere between immediately after to several tens of millions of years after the split between angiosperms and extant gymnosperms. Estimates of the age of the most recent common ancestor of all extant angiosperms range from approximately 140-210 Ma, depending on the trees used and assumptions made.

Isolation of the LFY/FLO homologue in Orchis italica and evolutionary analysis in some European orchids

It has been suggested that the evolutionary analysis of floral development genes could explain the divergences between the rates of morphological and molecular evolution. The LEAFY (LFY) gene of Arabidopsis thaliana is one of the central regulatory genes in the control of flower development. We have identified the homologue of this gene (OrcLFY) in Orchis italica using 5'/3'RACE and primer walking, and compared the coding sequences of several orchid species. We analyzed nonsynonymous and synonymous substitution rates between the OrcLFY coding regions of 14 species, and performed a McDonald-Kreitman test on Orchis morio and Orchis laxiflora populations, showing that purifying selection is acting on this gene in these orchids. We have performed a phylogenetic analysis showing that OrcLFY is a new useful marker to reconstruct molecular phylogenies at low taxonomic levels.

Hybridization and adaptive radiation

Whether interspecific hybridization is important as a mechanism that generates biological diversity is a matter of controversy. Whereas some authors focus on the potential of hybridization as a source of genetic variation, functional novelty and new species, others argue against any important role, because reduced fitness would typically render hybrids an evolutionary dead end. By drawing on recent developments in the genetics and ecology of hybridization and on principles of ecological speciation theory, I develop a concept that reconciles these views and adds a new twist to this debate. Because hybridization is common when populations invade new environments and potentially elevates rates of response to selection, it predisposes colonizing populations to rapid adaptive diversification under disruptive or divergent selection. I discuss predictions and suggest tests of this hybrid swarm theory of adaptive radiation and review published molecular phylogenies of adaptive radiations in light of the theory.

Alternative splicing mechanisms for the modulation of protein function: conservation between human and other species

Alternative splicing (AS) is an important process in eukaryotic organisms by which a given gene may express a set of different protein isoforms depending on the tissue, or the developmental stage of the individual. In the present work, we have compared AS among species, focusing on the conservation of AS mechanisms for the modulation of protein function. For this purpose, we first analysed the frequency with which different species, human, mouse, rat and fruitfly, utilise them. Second, we focused more directly on the conservation among species of the mechanisms themselves. To this end, we compared biologically equivalent AS events between human and mouse, or rat. Our results indicate only minor differences in the frequency of use of these mechanisms, as well as a high degree of conservation among the species studied.

Functional mutants of the sequence-specific transcription factor p53 and implications for master genes of diversity

There are many sources of genetic diversity, ranging from programmed mutagenesis in antibody genes to random mutagenesis during species evolution or development of cancer. We propose that mutations in DNA sequence-specific transcription factors that target response elements (REs) in many genes can also provide for rapid and broad phenotypic diversity, if the mutations lead to altered binding affinities at individual REs. To test this concept, we examined the in vivo transactivation capacity of wild-type human and murine p53 and 25 partial function mutants. The p53s were expressed in yeast from a rheostatable promoter, and the transactivation capacities toward >15 promoter REs upstream of a reporter gene were measured. Surprisingly, there was wide variation in transactivation by the mutant p53s toward the various REs. This is the first study to address directly the impact of mutations in a sequence-specific transcription factor on transactivation from a wide array of REs. We propose a master gene hypothesis for phenotypic diversity where the master gene is a single transcriptional activator (or repressor) that regulates many genes through different REs. Mutations of the master gene can lead to a variety of simultaneous changes in both the selection of targets and the extent of transcriptional modulation at the individual targets, resulting in a vast number of potential phenotypes that can be created with minimal mutational changes without altering existing protein-protein interactions.

Patterns of nucleotide variation in homoeologous regulatory genes in the allotetraploid Hawaiian silversword alliance (Asteraceae)

Genome-wide duplication (polyploidization) is prevalent in a large number of eukaryotic organisms and is particularly widespread in flowering plants. Polyploid species appear to vary from their diploid progenitors in a variety of ecologically important traits, suggesting that genome duplications provide a mechanism for ecological diversification. Studies of nucleotide variation at duplicate genes that arise via polyploidization allow us to infer the evolutionary forces that act on these polyploid loci. In an effort to examine the evolutionary dynamics of homoeologous loci, molecular population genetic analyses were undertaken for duplicate regulatory genes in the allopolyploid Hawaiian silversword alliance, a premier example of adaptive radiation. The levels and patterns of nucleotide variation for the floral homeotic genes ASAPETALA1 (ASAP1) and ASAPETALA3/TM6 (ASAP3/TM6) were studied in two species representing different lineages within the Hawaiian silversword alliance: Argyroxiphium sandwicense ssp. macrocephalum and Dubautia ciliolata ssp. glutinosa. Homoeologueous copies of ASAP1 and ASAP3/TM6 show differing levels and patterns of nucleotide polymorphism. Duplicate ASAP1 copies have similar levels of nucleotide diversity and haplotype structure in both species; by contrast, duplicate ASAP3/TM6 genes display different levels and patterns of variation in D. ciliolata ssp. glutinosa. Additionally, D. ciliolata ssp. glutinosa appears to be segregating for a moderate frequency null allele in one ASAP3/TM6 homoeologue. These results suggest that differing evolutionary forces can affect duplicate loci arising from allopolyploidization.

Tracing floral adaptations from ecology to molecules

Flowers have long fascinated humans. The scientific study of floral biology unifies many diverse areas of research, ranging from systematics to ecology, and from genetics to molecular biology. Despite this unity, few plant species offer the experimental versatility to encompass all levels of biological investigation in a single system. An exception is the morning glory genus Ipomoea, in which a broad picture of floral evolution, ranging from ecology to molecular biology, is emerging.

Recently duplicated maize R2R3 Myb genes provide evidence for distinct mechanisms of evolutionary divergence after duplication

R2R3 Myb genes are widely distributed in the higher plants and comprise one of the largest known families of regulatory proteins. Here, we provide an evolutionary framework that helps explain the origin of the plant-specific R2R3 Myb genes from widely distributed R1R2R3 Myb genes, through a series of well-established steps. To understand the routes of sequence divergence that followed Myb gene duplication, we supplemented the information available on recently duplicated maize (Zea mays) R2R3 Myb genes (C1/Pl1 and P1/P2) by cloning and characterizing ZmMyb-IF35 and ZmMyb-IF25. These two genes correspond to the recently expanded P-to-A group of maize R2R3 Myb genes. Although the origins of C1/Pl1 and ZmMyb-IF35/ZmMyb-IF25 are associated with the segmental allotetraploid origin of the maize genome, other gene duplication events also shaped the P-to-A clade. Our analyses indicate that some recently duplicated Myb gene pairs display substantial differences in the numbers of synonymous substitutions that have accumulated in the conserved MYB domain and the divergent C-terminal regions. Thus, differences in the accumulation of substitutions during evolution can explain in part the rapid divergence of C-terminal regions for these proteins in some cases. Contrary to previous studies, we show that the divergent C termini of these R2R3 MYB proteins are subject to purifying selection. Our results provide an in-depth analysis of the sequence divergence for some recently duplicated R2R3 Myb genes, yielding important information on general patterns of evolution for this large family of plant regulatory genes.

GAI homologues in the Hawaiian silversword alliance (Asteraceae-Madiinae): molecular evolution of growth regulators in a rapidly diversifying plant lineage

Accelerated evolution of regulatory genes has been proposed as an explanation for decoupled rates of morphological and molecular evolution. The Hawaiian silversword alliance (Asteraceae-Madiinae) has evolved drastic differences in growth form, including rosette plants, cushion plants, shrubs, and trees, since its origin approximately 6 MYA. We have isolated genes in the DELLA subfamily of putative growth regulators from 13 taxa of Hawaiian and North American Madiinae. The Hawaiian taxa contain two copies of DaGAI that form separate clades within the Madiinae, consistent with an allotetraploid origin for the silversword alliance. DaGAI retains conserved features that have previously been identified in DELLA genes. Selective constraint in the Hawaiian DaGAI copies remains strong in spite of rapid growth form divergence in the silversword alliance, although the constraint was somewhat relaxed in the Hawaiian copies relative to the North American lineages. We failed to detect evidence for positive selection on individual codons. Notably, selective constraint remained especially strong in the gibberellin-responsive DELLA region for which the gene subfamily is named, which is truncated or deleted in all identified dwarf mutants in GAI homologues in different angiosperm species. In contrast with the coding region, however, approximately 900 bp of the upstream flanking region shows variable rates and patterns of evolution, which might reflect positive selection on regulatory regions.

Power analysis of tests for loss of selective constraint in cave crayfish and nonphotosynthetic plant lineages

Loss of selective constraint on a gene may be expected following changes in the environment or life history that render its function unnecessary. The long-term persistence of protein-coding genes after the loss of known functional necessity can occur by chance or because of selective maintenance of an unknown gene function. The selective maintenance of an alternative gene function is not demonstrated by the failure of statistical tests to reject the hypothesis that there has been no change in the degree of constraint on the evolution of coding genes. Maintenance may be inferred, however, when power analyses of such tests demonstrate that there has been a sufficient number of nucleotide substitutions to detect the loss of selective constraint. Here, we describe a power analysis for tests of loss of constraint on protein-coding genes. The power analysis was applied to loss-of-constraint tests for opsin gene evolution in cave-dwelling crayfish and rbcL evolution in nonphotosynthetic parasitic plants. The power of previously applied tests for loss of constraint on cave crayfish opsin genes was insufficient to distinguish between chance retention and selective maintenance of opsin genes. However, the power of codon-based likelihood ratio tests for change in d(N)/d(S) (=omega) (nonsynonymous to synonymous change) did have sufficient power to detect a loss of constraint on rbcL associated with a loss of photosynthesis in most examples but failed to detect such a change in three independent lineages. We conclude that rbcL has been selectively maintained in these holoparasitic plant lineages. This conclusion suggests that either these taxa are photosynthetic for at least a part of their life or rbcL may have an unknown function in these plants unrelated to photosynthesis.

Cytochrome b-based phylogeny of the Praomys group (Rodentia, Murinae): a new African radiation?

Complete cytochrome b gene sequences allows, for the first time, establishing a nearly complete phylogeny among the Praomys group sensu lato. The genera Praomys, Mastomys and Stenocephalemys appear paraphyletic. Myomys is polyphyletic and this genus name probably needs to be restricted to its type species, M. verreauxii. The genera Zelotomys and Colomys appear as sister groups. Mastomys pernanus and Malacomys verschureni nest within the Praomys group, but their generic assignation must be further clarified. The genus Heimyscus appears closest to Praomys than to Hylomyscus. The different lineages probably result from an adaptive radiation at the end of the Miocene.

Molecular phylogeny and biogeography of the Qinghai-Tibet Plateau endemic Nannoglottis (Asteraceae)

All taxa endemic to the Qinghai-Tibet Plateau are hypothesized to have originated in situ or from immediately adjacent areas because of the relatively recent formation of the plateau since the Pliocene, followed by the large-scaled biota extinction and recession caused by the Quaternary ice sheet. However, identification of specific progenitors remains difficult for some endemics, especially some endemic genera. Nannoglottis, with about eight species endemic to this region, is one such genus. Past taxonomic treatments have suggested its relationships with four different tribes of Asteraceae. We intend to identify the closest relatives of Nannoglottis by evaluating the level of monophyly, tribal delimitation, and systematic position of the genus by using molecular data from ndhF gene, trnL-F, and ITS region sequences. We find that all sampled species of Nannoglottis are a well-defined monophyly. This supports all recent taxonomic treatments of Nannoglottis, in which all sampled species were placed in one broadly re-circumscribed genus. Nannoglottis is most closely related to the Astereae, but stands as an isolated genus as the first diverging lineage of the tribe, without close relatives. A tentative relationship was suggested for Nannoglottis and the next lineage of the tribe was based on the ITS topology, the "basal group," which consists of seven genera from the Southern Hemisphere. Such a relationship is supported by some commonly shared plesiomorphic morphological characters. Despite the very early divergence of Nannoglottis in the Astereae, the tribe must be regarded to have its origin in Southern Hemisphere rather than in Asia, because based on all morphological, molecular, biogeographical, and fossil data, the Asteraceae and its major lineages (tribes) are supposed to have originated in the former area. Long-distance dispersal using Southeast Asia as a steppingstone from Southern Hemisphere to the Qinghai-Tibet Plateau is the most likely explanation for this unusual biogeographic link of Nannoglottis. The 23-32-million-year divergence time between Nannoglottis and the other Astereae estimated by DNA sequences predated the formation of the plateau. This estimation is further favored by the fossil record of the Asteraceae and the possible time of origin of the Astereae. Nannoglottis seems to have reached the Qinghai-Tibet area in the Oligocene-Eocene and then re-diversified with the uplift of the plateau. The molecular infragenetic phylogeny of the genus identifies two distinct clades, which reject the earlier infrageneric classification based on the arrangement of the involucral bracts and the length of the ligules, but agree well with the habits and ecological preferences of its current species. The "alpine shrub" vs. "coniferous forest" divergence within Nannoglottis was estimated at about 3.4 million years ago when the plateau began its first large-scale uplifting and the coniferous vegetation began to appear. Most of the current species at the "coniferous forest" clade of the genus are estimated to have originated from 1.02 to 1.94 million years ago, when the second and third uprisings of the plateau occurred, the climate oscillated and the habitats were strongly changed. The assumed evolution, speciation diversity, and radiation of Nannoglottis based on molecular phylogeny and divergence times agree well with the known geological and paleobotanical histories of the Qinghai-Tibet Plateau.