Duplication of floral regulatory genes in the Lamiales

Duplication of some floral regulatory genes has occurred repeatedly in angiosperms, whereas others are thought to be single-copy in most lineages. We selected three genes that interact in a pathway regulating floral development conserved among higher tricolpates (LFY/FLO, UFO/FIM, and AP3/DEF) and screened for copy number among families of Lamiales that are closely related to the model species Antirrhinum majus. We show that two of three genes have duplicated at least twice in the Lamiales. Phylogenetic analyses of paralogs suggest that an ancient whole genome duplication shared among many families of Lamiales occurred after the ancestor of these families diverged from the lineage leading to Veronicaceae (including the single-copy species A. majus). Duplication is consistent with previous patterns among angiosperm lineages for AP3/DEF, but this is the first report of functional duplicate copies of LFY/FLO outside of tetraploid species. We propose Lamiales taxa will be good models for understanding mechanisms of duplicate gene preservation and how floral regulatory genes may contribute to morphological diversity.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

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

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

Diversification of genes for carotenoid biosynthesis in aphids following an ancient transfer from a fungus

The pea aphid genome was recently found to harbor genes for carotenoid biosynthesis, reflecting an ancestral transfer from a fungus. To explore the evolution of the carotene desaturase gene family within aphids, sequences were retrieved from a set of 34 aphid species representing numerous deeply diverging lineages of aphids and analyzed together with fungal sequences retrieved from databases. All aphids have at least one copy of this gene and some aphid species have up to seven, whereas fungal genomes consistently have a single copy. The closest relatives of aphids, adelgids, also have carotene desaturase; these sequences are most closely related to those from aphids, supporting a shared origin from a fungal to insect transfer predating the divergence of adelgids and aphids. Likewise, all aphids, and adelgids, have carotenoid profiles that are consistent with their biosynthesis using the acquired genes of fungal origin rather than derivation from food plants. The carotene desaturase was acquired from a fungal species outside of Ascomycota or Basidiomycota and closest to Mucoromycotina among sequences available in databases. In aphids, an ongoing pattern of gene duplication is indicated by the presence of both anciently and recently diverged paralogs within genomes and by the presence of a high frequency of pseudogenes that appear to be recently inactivated. Recombination among paralogs is evident, making analyses of patterns of selection difficult, but tests of selection for a nonrecombining region indicates that duplications tend to be followed by bouts of positive selection. Species of Macrosiphini, which often show color polymorphisms, typically have a larger number of desaturase copies relative to other species sampled in the study. These results indicate that aphid evolution has been accompanied by ongoing evolution of carotenogenic genes, which have undergone duplication, recombination, and occasional positive selection to yield a wide variety of carotenoid profiles in different aphid species.     

Evolution of bract development and B-class MADS box gene expression in petaloid bracts of Cornus s. l. (Cornaceae)

? Despite increasing interest in the molecular mechanisms of floral diversity, few studies have investigated the developmental and genetic bases of petaloid bracts. This study examined morphological patterns of bract initiation and expression patterns of B-class MADS-box genes in bracts of several Cornus species. We suggest that petaloid bracts in this genus may not share a single evolutionary origin. ? Developmental pathways of bracts and spatiotemporal expression of B-class genes in bracts and flowers were examined for four closely related dogwood species. ? Divergent morphological progressions and gene expression patterns were found in the two sister lineages with petaloid bracts, represented by Cornus florida and Cornus canadensis. Phylogeny-based analysis identified developmental and gene expression changes that are correlated with the evolution of petaloid bracts in C.?florida and C.?canadensis. ? Our data support the existence of independent evolutionary origins of petaloid bracts in C.?canadensis and C.?florida. Additionally, we suggest that functional transference within B-class gene families may have contributed to the origin of bract petaloidy in C.?florida. However, the underlying mechanisms of petaloid bract development likely differ between C.?florida and C.?canadensis. In the future this hypothesis can be tested by functional analyses of Cornus B-class genes.     

Evolution of duplicate genes in a tetraploid animal, Xenopus laevis

To understand the evolution of duplicate genes, we compared rates of nucleotide substitution between 17 pairs of nonallelic duplicated genes in the tetraploid frog Xenopus laevis with rates between the orthologous loci of human and rodent. For all duplicated X. laevis genes, the number of synonymous substitutions per site (dS) was greater than the number of nonsynonymous substitutions per site (dN), indicating that these genes are subject to purifying selection. There was also a significant positive correlation (r = 0.915) between dN for the X. laevis genes and dN for the mammalian genes, suggesting that, at the amino acid level, the X. laevis genes and the mammalian genes are under similar constraints. Results of relative-rate tests showed nearly equal rates of nonsynonymous substitution in each copy of the X. laevis genes; apparently there are similar constraints on both copies. No correlation was found between dS for the X. laevis genes and dS for the mammalian genes. There was a significant positive correlation both between members of pairs of duplicated X. laevis genes (r = 0.951) and between human and rodent orthologues (r = 0.854) with respect to third- position G+C content but no such relationship between the X. laevis genes and either of their mammalian orthologues. The results indicate that both copies of a duplicate gene can be subject to purifying selection and thus support the hypothesis of selection against all genotypes containing a null allele at either of two duplicate loci.      

Cloning and characterization of a PI-like MADS-box gene in Phalaenopsis orchid

The highly evolved flowers of orchids have colorful sepals and fused columns that offer an opportunity to discover new genes involved in floral development in monocotyledon species. In this investigation, we cloned and characterized the homologous PISTALLATA-like (PI-like) gene PhPI15 (Phalaenopsis PI STILLATA # 15), from the Phalaenopsis hybrid cultivar. The protein sequence encoded by PhPI15 contains a typical PI-motif. Its sequence also formed a subclade with other monocot PI-type genes in phylogenetic analysis. Southern analysis showed that PhPI15 was present in the Phalaenopsis orchid genome as a single copy. Furthermore, it was expressed in all the whorls of the Phalaenopsis flower, while no expression was detected in vegetative organs. The flowers of transgenic tobacco plants ectopically expressing PhPI15 showed male-sterile phenotypes. Thus, as a Class-B MADS-box gene, PhPI15 specifies floral organ identity in orchids.     

The Population Genetics of dN/dS

Evolutionary pressures on proteins are often quantified by the ratio of substitution rates at non-synonymous and synonymous sites. The dN/dS ratio was originally developed for application to distantly diverged sequences, the differences among which represent substitutions that have fixed along independent lineages. Nevertheless, the dN/dS measure is often applied to sequences sampled from a single population, the differences among which represent segregating polymorphisms. Here, we study the expected dN/dS ratio for samples drawn from a single population under selection, and we find that in this context, dN/dS is relatively insensitive to the selection coefficient. Moreover, the hallmark signature of positive selection over divergent lineages, dN/dS>1, is violated within a population. For population samples, the relationship between selection and dN/dS does not follow a monotonic function, and so it may be impossible to infer selection pressures from dN/dS. These results have significant implications for the interpretation of dN/dS measurements among population-genetic samples.     

Expression pattern of two paralogs of the <Emphasis Type="Italic">PI</Emphasis>/<Emphasis Type="Italic">GLO</Emphasis>-like locus during <Emphasis Type="Italic">Orchis italica</Emphasis> (Orchidaceae,Orchidinae) flower development

The class B MADS-box genes belong to two distinct functional groups: the AP3/DEF-like and the PI/GLO-like sub-families. In orchids, AP3/DEF-like genes are present in four copies, each with a different role in floral organ formation, which is described in the “orchid code” model. Interestingly, the orchid PI/GLO-like genes are present in two copies in Orchidinae, whereas they are described as single copy in the other orchid lineages. The two PI/GLO-like paralogs have site-specific different selective constraints; in addition, they show relaxation of purifying selection when compared to the single-copy lineages. In this study, we present a comparative analysis of the expression patterns of the two PI/GLO-like paralogs, OrcPI and OrcPI2, in floral tissues of Orchis italica in different developmental stages using real-time PCR. The two genes show similar expression profiles in the tissue examined, with differences detectable between immature and mature inflorescence. In all cases, OrcPI2 is expressed at a higher level than OrcPI. Real-time PCR results reveal that the co-expression of the two duplicated loci could have a fully or partially redundant function. The possible evolutionary fate of OrcPI and OrcPI2 is discussed as well as their involvement in ovary development.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

Selection for functional uniformity of tuf duplicates in gamma-proteobacteria

Having an extra copy of a gene is thought to provide some functional redundancy, which results in a higher rate of evolution in duplicated genes. In this article, we estimate the impact of gene duplication on the selection of tuf paralogs, and we find that in the absence of gene conversion, tuf paralogs have evolved significantly slower than when gene conversion has been a factor in their evolution. Thus, tuf gene copies evolve under a selective pressure that ensures their functional uniformity, and gene conversion reduces selection against amino acid substitutions that affect the function of the encoded protein, EF-Tu.     

Phylogeny reconstruction in the Caesalpinieae grade (Leguminosae) based on duplicated copies of the sucrose synthase gene and plastid markers

The Caesalpinieae grade (Leguminosae) forms a morphologically and ecologically diverse group of mostly tropical tree species with a complex evolutionary history. This grade comprises several distinct lineages, but the exact delimitation of the group relative to subfamily Mimosoideae and other members of subfamily Caesalpinioideae, as well as phylogenetic relationships among the lineages are uncertain. With the aim of better resolving phylogenetic relationships within the Caesalpinieae grade, we investigated the utility of several nuclear markers developed from genomic studies in the Papilionoideae. We cloned and sequenced the low copy nuclear gene sucrose synthase (SUSY) and combined the data with plastid trnL and matK sequences. SUSY has two paralogs in the Caesalpinieae grade and in the Mimosoideae, but occurs as a single copy in all other legumes tested. Bayesian and maximum likelihood phylogenetic analyses suggest the two nuclear markers are congruent with plastid DNA data. The Caesalpinieae grade is divided into four well-supported clades (Cassia, Caesalpinia, Tachigali and Peltophorum clades), a poorly supported clade of Dimorphandra Group genera, and two paraphyletic groups, one with other Dimorphandra Group genera and the other comprising genera previously recognized as the Umtiza clade. A selection analysis of the paralogs, using selection models from PAML, suggests that SUSY genes are subjected to a purifying selection. One of the SUSY paralogs, under slightly stronger positive selection, may be undergoing subfunctionalization. The low copy SUSY gene is useful for phylogeny reconstruction in the Caesalpinieae despite the presence of duplicate copies. This study confirms that the Caesalpinieae grade is an artificial group, and highlights the need for further analyses of lineages at the base of the Mimosoideae.