The molecular evolution of development

Genetic and molecular analysis of trichome development in Arabis alpina

The genetic and molecular analysis of trichome development in Arabidopsis thaliana has generated a detailed knowledge about the underlying regulatory genes and networks. However, how rapidly these mechanisms diverge during evolution is unknown. To address this problem, we used an unbiased forward genetic approach to identify most genes involved in trichome development in the related crucifer species Arabis alpina In general, we found most trichome mutant classes known in A. thaliana We identified orthologous genes of the relevant A. thaliana genes by sequence similarity and synteny and sequenced candidate genes in the A. alpina mutants. While in most cases we found a highly similar gene-phenotype relationship as known from Arabidopsis, there were also striking differences in the regulation of trichome patterning, differentiation, and morphogenesis. Our analysis of trichome patterning suggests that the formation of two classes of trichomes is regulated differentially by the homeodomain transcription factor AaGL2 Moreover, we show that overexpression of the GL3 basic helix-loop-helix transcription factor in A. alpina leads to the opposite phenotype as described in A. thaliana Mathematical modeling helps to explain how this nonintuitive behavior can be explained by different ratios of GL3 and GL1 in the two species.

Differences in brain transcriptomes of closely related Baikal coregonid species

The aim of this work was to get deeper insight into genetic factors involved in the adaptive divergence of closely related species, specifically two representatives of Baikal coregonids—Baikal whitefish (Coregonus baicalensis Dybowski) and Baikal omul (Coregonus migratorius Georgi)—that diverged from a common ancestor as recently as 10–20 thousand years ago. Using the Serial Analysis of Gene Expression method, we obtained libraries of short representative cDNA sequences (tags) from the brains of Baikal whitefish and omul. A comparative analysis of the libraries revealed quantitative differences among ~4% tags of the fishes under study. Based on the similarity of these tags with cDNA of known organisms, we identified candidate genes taking part in adaptive divergence. The most important candidate genes related to the adaptation of Baikal whitefish and Baikal omul, identified in this work, belong to the genes of cell metabolism, nervous and immune systems, protein synthesis, and regulatory genes as well as to DTSsa4 Tc1-like transposons which are widespread among fishes.

Adaptive molecular evolution of a defence gene in sexual but not functionally asexual evening primroses

Theory predicts that sexual reproduction provides evolutionary advantages over asexual reproduction by reducing mutational load and increasing adaptive potential. Here, we test the latter prediction in the context of plant defences against pathogens because pathogens frequently reduce plant fitness and drive the evolution of plant defences. Specifically, we ask whether sexual evening primrose plant lineages (Onagraceae) have faster rates of adaptive molecular evolution and altered gene expression of a class I chitinase, a gene implicated in defence against pathogens, than functionally asexual evening primrose lineages. We found that the ratio of amino acid to silent substitutions (K(a) /K(s) = 0.19 vs. 0.11 for sexual and asexual lineages, respectively), the number of sites identified to be under positive selection (four vs. zero for sexual and asexual lineages, respectively) and the expression of chitinase were all higher in sexual than in asexual lineages. Our results are congruent with the conclusion that a loss of sexual recombination and segregation in the Onagraceae negatively affects adaptive structural and potentially regulatory evolution of a plant defence protein.

Nucleotide diversity and molecular evolution of the PSY1 gene in Zea mays compared to some other grass species

Phytoene synthase (PSY), which is encoded by the phytoene synthase 1 (PSY1) gene, is the first rate-limiting enzyme in the plant carotenoid biosynthetic pathway. In order to examine the genetic diversity and evolution pattern of PSY1 within the Andropogoneae, sequences of 76 accessions from 5 species (maize, teosinte, tripsacum, coix, and sorghum) of the Andropogoneae were tested, along with 4 accessions of rice (Oryza sativa L.) included as outliers. Both the number and the order of exons and introns were relatively conserved across the species tested. Three domains were identified in the coding sequence, including signal peptide (SP), PSY, and highly conserved squalene synthase (SQS) domain. Although no positive selection signal was detected at an overall coding level among all species tested, the SP domain and the region upstream of the SQS-PSY domain appear to have undergone rapid evolution, as evidenced by a high d (N)/d (S) ratio (>1.0). At the nucleotide level, positive selection and balancing selection were detected only among the yellow maize germplasm and the white maize germplasm, respectively. The phylogenetic tree based on full-length sequences of PSY1-like regions supported the monophyletic theory of the Andropogoneae and the closest relationship between Zea and Tripsacum among the Andropogoneae. Coix, which was theorized to have a closer relationship with maize due to similarities in morphology and chromosome number, has been shown in this study to have diverged relatively early from the other Andropogoneae, including maize.

Comparative analysis of vertebrate PEPT1 and PEPT2 genes

The plasma membrane transport proteins belong to SoLute Carrier 15 (SLC15) family and two members of this family have been characterized extensively in higher vertebrates, namely PEPT1 and PEPT2. Despite many efforts have made to define a pharmacophore model for efficient binding and transporting of substrates, there is not a comprehensive study performed to elucidate the evolutionary mechanisms among the SLC15 family members and to statistically evaluate sequence conservation and functional divergence between members. In this study, we compared and contrasted the rates and patterns of molecular evolution of 2 PEPT genes. Phylogenetic tree assembly with all available vertebrate PEPTs suggests that the PEPTs originated by duplications and diverged from a common protein at the base of the eukaryotic tree. Topological structure demonstrates both members share the similar hydrophobic domains (TMDs), which have been constrained by purifying selection. Although both genes show qualitatively similar patterns, their rates of evolution differ significantly due to an increased rate of synonymous substitutions in the structural domains in one copy, suggesting substantial differences in functional constraint on each gene. Site-specific profiles were established by posterior probability analysis revealing significantly divergent regions mainly locate at the hydrophobic region between predicted transmembrane domains 9 and 10 of the proteins. Thus, these results provide the evidence that several amino acid residues with reduced selective constraints are largely responsible for functional divergence between the paralogous PEPTs. These findings may provide a starting point for further experimental verifications.

[Genome similarity of Baikal omul and sig]

Two members of the Baikal sig family, a lake sig (Coregonus lavaretus baicalensis Dybovsky) and omul (C. autumnalis migratorius Georgi), are close relatives that diverged from the same ancestor 10-20 thousand years ago. In this work, we studied genomic polymorphism of these two fish species. The method of subtraction hybridization (SH) did not reveal the presence of extended sequences in the sig genome and their absence in the omul genome. All the fragments found by SH corresponded to polymorphous noncoding genome regions varying in mononucleotide substitutions and short deletions. Many of them are mapped close to genes of the immune system and have regions identical to the Tc-1-like transposons abundant among fish, whose transcription activity may affect the expression of adjacent genes. Thus, we showed for the first time that genetic differences between Baikal sig family members are extremely small and cannot be revealed by the SH method. This is another endorsement of the hypothesis on the close relationship between Baikal sig and omul and their evolutionarily recent divergence from a common ancestor.

Pervasive sex-linked effects on transcription regulation as revealed by expression quantitative trait loci mapping in lake whitefish species pairs (Coregonus sp., Salmonidae)

Mapping of expression quantitative trait loci (eQTL) is a powerful means for elucidating the genetic architecture of gene regulation. Yet, eQTL mapping has not been applied toward investigating the regulation architecture of genes involved in the process of population divergence, ultimately leading to speciation events. Here, we conducted an eQTL mapping experiment to compare the genetic architecture of transcript regulation in adaptive traits, differentiating the recently evolved limnetic (dwarf) and benthic (normal) species pairs of lake whitefish. The eQTL were mapped in three data sets derived from an F(1) hybrid-dwarf backcrossed family: the entire set of 66 genotyped individuals and the two sexes treated separately. We identified strikingly more eQTL in the female data set (174), compared to both male (54) and combined (33) data sets. The majority of these genes were not differentially expressed between male and female progeny of the backcross family, thus providing evidence for a strong pleiotropic sex-linked effect in transcriptomic regulation. The subtelomeric region of a linkage group segregating in females encompassed >50% of all eQTL, which exhibited the most pronounced additive effects. We also conducted a direct comparison of transcriptomic profiles between pure dwarf and normal progeny reared in controlled conditions. We detected 34 differentially expressed transcripts associated with eQTL segregating only in sex-specific data sets and mostly belonging to functional groups that differentiate dwarf and normal whitefish in natural populations. Therefore, these eQTL are not related to interindividual variation, but instead to the adaptive and historical genetic divergence between dwarf and normal whitefish. This study exemplifies how the integration of genetic and transcriptomic data offers a strong means for dissecting the functional genomic response to selection by separating mapping family-specific effects from genetic factors under selection, potentially involved in the phenotypic divergence of natural populations.

The role of regulatory genes during maize domestication: evidence from nucleotide polymorphism and gene expression

We investigated DNA sequence variation in 72 candidate genes in maize landraces and the wild ancestor of maize, teosinte. The candidate genes were chosen because they exhibit very low sequence diversity among maize inbreds and have sequence homology to known regulatory genes. We observed signatures of selection in 17 candidate genes, indicating that they were potential targets of artificial selection during domestication. In addition, 21 candidate genes were identified as potential targets of natural selection in teosinte. A comparison of the proportion of selected genes between our regulatory genes and genes unfiltered for their potential function (but also with very low sequence diversity among maize inbreds) provided some weak evidence that regulatory genes are overrepresented among selected genes. We detected no significant association between the positions of genes identified as potential targets of selection during domestication and quantitative trait loci (QTL) responsible for maize domestication traits. However, a subset of these genes, those identified by sequence homology as kinase/phosphatase genes, significantly cluster with the domestication QTL. We also analyzed expression profiles of genes in distinct maize tissues and observed that domestication genes are expressed on average at a significantly higher level than neutral genes in reproductive organs, including kernels.

Positive selection and expression divergence following gene duplication in the sunflower CYCLOIDEA gene family

Members of the CYCLOIDEA (CYC)/TEOSINTE-BRANCHED1 (TB1) group of transcription factors have been implicated in the evolution of zygomorphic (i.e., bilaterally symmetric) flowers in Antirrhinum and Lotus and the loss of branching phenotype during the domestication of maize. The composite inflorescences of sunflower (Helianthus annuus L. Asteraceae) contain both zygomorphic and actinomorphic (i.e., radially symmetric) florets (rays and disks, respectively), and the cultivated sunflower has evolved an unbranched phenotype in response to domestication from its highly branched wild progenitor; hence, genes related to CYC/TB1 are of great interest in this study system. We identified 10 members of the CYC/TB1 gene family in sunflower, which is more than found in any other species investigated to date. Phylogenetic analysis indicates that these genes occur in 3 distinct clades, consistent with previous research in other eudicot species. A combination of dating the duplication events and linkage mapping indicates that only some of the duplications were associated with polyploidization. Cosegregation between CYC-like genes and branching-related quantitative trait loci suggest a minor, if any, role for these genes in conferring differences in branching. However, the expression patterns of one gene suggest a possible role in the development of ray versus disk florets. Molecular evolutionary analyses reveal that residues in the conserved domains were the targets of positive selection following gene duplication. Taken together, these results indicate that gene duplication and functional divergence have played a major role in diversification of the sunflower CYC gene family.

Molecular evolution of the OrcPI locus in natural populations of Mediterranean orchids

The evolutionary analysis of OrcPI, the orchid homologue to the PISTILLATA/GLOBOSA gene, was conducted on some Mediterranean orchid species, measuring mean pairwise Ka/Ks ratios and nucleotide variability. Evidence for positive selection was tested using a maximum likelihood approach implemented in PAML, and neutrality tests were conducted to assess deviation from neutral evolution. Data were also examined partitioning the coding region into four regions, corresponding to different functional domains of the protein. The results show that OrcPI is subjected to different evolutionary forces: diffuse purifying selection, localized positive selection or selective sweep, and different partitions of selective constraints.

Evolutionary patterns in the antR-Cor gene in the dwarf dogwood complex (Cornus, Cornaceae)

The evolutionary pattern of the myc-like anthocyanin regulatory gene antR-Cor was examined in the dwarf dogwood species complex (Cornus Subgenus Arctocrania) that contains two diploid species (C. canadensis and C. suecica), their putative hybrids with intermediate phenotypes, and a tetraploid derivative (C. unalaschkensis). Full-length sequences of this gene ( approximately 4 kb) were sequenced and characterized for 47 dwarf dogwood samples representing all taxa categories from 43 sites in the Pacific Northwest. Analysis of nucleotide diversity indicated departures from neutral evolution, due most likely to local population structure. Neighbor-joining and haplotype network analyses show that sequences from the tetraploid and diploid intermediates are much more strongly diverged from C. suecica than from C. canadensis, and that the intermediate phenotypes may represent an ancestral group to C. canadensis rather than interspecific hybrids. Seven amino acid mutations that are potentially linked to myc-like anthocyanin regulatory gene function correlate with petal colors differences that characterize the divergence between two diploid species and the tetraploid species in this complex. The evidence provides a working hypothesis for testing the role of the gene in speciation and its link to the petal coloration. Sequencing and analysis of additional nuclear genes will be necessary to resolve questions about the evolution of the dwarf dogwood complex.

Parallelism in gene transcription among sympatric lake whitefish (Coregonus clupeaformis Mitchill) ecotypes

Abstract We tested the hypothesis that phenotypic parallelism between dwarf and normal whitefish ecotypes (Coregonus clupeaformis, Salmonidae) is accompanied by parallelism in gene transcription. The most striking phenotypic differences between these forms implied energetic metabolism and swimming activity. Therefore, we predicted that genes showing parallel expression should mainly belong to functional groups associated with these phenotypes. Transcriptome profiles were obtained from white muscle by using a 3557 cDNA gene microarray developed for the Atlantic salmon (Salmo salar). A total of 1181 genes expressed in both lake populations hybridized on the array. Significant differential expression between ecotypes was detected for 134 (11.3%) and 195 (16.5%) gene clones in Cliff Lake and Indian Pond, respectively. Fifty-one genes (4.3%) showed parallel differential expression between lakes, among which 35 were expressed in opposite directions. Sixteen genes (1.35%) showed true parallelism of transcription, which mainly belonged to energetic metabolism and regulation of muscle contraction functional groups. Variance in expression was significantly reduced for these genes compared to those not showing directionality in parallelism of expression. Candidate genes associated with parallelism in swimming activity and energetic metabolism based on their level and variance in expression were identified. These results add to the growing evidence that parallel phenotypic evolution also involves parallelism at both the genotypic and regulatory level, which may at least partly be associated with genetic constraints. It also provides further evidence for the determinant role of divergent natural selection in driving phenotypic divergence, and perhaps reproductive isolation, in the adaptive radiation of lake whitefish. This study adds to a nascent field employing microarrays as powerful tools for investigating the evolutionary processes of adaptive divergence among natural populations.

A new MADS-box gene (IbMADS10) from sweet potato (Ipomoea batatas (L.) Lam) is involved in the accumulation of anthocyanin

A new MADS-box gene designated as IbMADS10 was cloned and its expression was characterized from sweet potato (Ipomoea batatas (L.) Lam.) cv. Beniazuma. The deduced amino acid sequence of the gene indicated high homology with members of the MADS-box family of transcription factors. IbMADS10 shares high amino acid sequence similarity with the DEFH28 of Antirrhinum majus (64%) and with BpMADS4 of Betula pendula (61%) of the SQUA subfamily. Southern blot analysis revealed that the IbMADS10 is present in one or low copy number in the sweet potato genome. The gene is specifically expressed in the pigmented tissues such as in the flower bud, in the pink and in red roots, and hence, it was speculated that the IbMADS10 gene might be correlated with anthocyanin biosynthesis in sweet potato. RNA blot expression of the anthocyanin biosynthesis genes encoding for CHS, CHI, F3H, DFR, ANS and UFTG carried out in the tissues where the IbMADS10 gene was expressed revealed similar transcript levels in all tissues where the IbMADS10 gene is highly expressed, indicating that the IbMADS10 gene is highly correlated with the anthocyanin biosynthesis genes. Another important aspect is the pigmented phenotypes of transgenic calli that ectopically express the IbMADS10 gene, thereby supporting its involvement in the developmental regulation of pigment formation. Tissue printing result further strengthens the hypothesis that the IbMADS10 gene is indeed involved in anthocyanin pigmentation in sweet potato. As the purpose of the IbMADS10 gene is pigmentation, its function, therefore, resembles that of the transparent testa (tt) genes of Arabidopsis.

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.

The evolution of the SEPALLATA subfamily of MADS-box genes: a preangiosperm origin with multiple duplications throughout angiosperm history

Members of the SEPALLATA (SEP) MADS-box subfamily are required for specifying the "floral state" by contributing to floral organ and meristem identity. SEP genes have not been detected in gymnosperms and seem to have originated since the lineage leading to extant angiosperms diverged from extant gymnosperms. Therefore, both functional and evolutionary studies suggest that SEP genes may have been critical for the origin of the flower. To gain insights into the evolution of SEP genes, we isolated nine genes from plants that occupy phylogenetically important positions. Phylogenetic analyses of SEP sequences show that several gene duplications occurred during the evolution of this subfamily, providing potential opportunities for functional divergence. The first duplication occurred prior to the origin of the extant angiosperms, resulting in the AGL2/3/4 and AGL9 clades. Subsequent duplications occurred within these clades in the eudicots and monocots. The timing of the first SEP duplication approximately coincides with duplications in the DEFICIENS/GLOBOSA and AGAMOUS MADS-box subfamilies, which may have resulted from either a proposed genome-wide duplication in the ancestor of extant angiosperms or multiple independent duplication events. Regardless of the mechanism of gene duplication, these pairs of duplicate transcription factors provided new possibilities of genetic interactions that may have been important in the origin of the flower.

No evidence for an association between common nonsynonymous polymorphisms in delta and bristle number variation in natural and laboratory populations of Drosophila melanogaster

We test the hypothesis that naturally occurring nonsynonymous variants in the Delta ligand of the Notch signaling pathway contribute to standing variation in sternopleural and/or abdominal bristle number in Drosophila melanogaster, for both a large cohort of wild-caught flies and previously described laboratory lines. We sequenced the transcribed region of Delta for 16 naturally occurring chromosomes and 65 SNPs, including 7 nonsynonymous SNPs (nsSNPs), were observed. Identified nsSNPs and 6 additional common SNPs, all located in exon 6 and the 3' UTR, were genotyped in 2060 wild-caught flies using an OLA-based methodology and genotyped in 38 additional natural chromosomes via DNA sequencing. None of the genotyped nsSNPs were significantly associated with natural variation in bristle number as assessed by a permutation test. A 95% upper bound on the additive genetic variance attributable to each genotyped SNP in the large natural cohort is <2% of the total phenotypic variation. Results suggest that two previously detected genotype/phenotype associations between bristle number and variants in the introns of Delta cannot be explained by linkage disequilibrium between these variants and nearby nonsynonymous variants. Unidentified regulatory variants more parsimoniously explain previous observations.

Adaptive evolution in the Arabidopsis MADS-box gene family inferred from its complete resolved phylogeny

Gene duplication is a substrate of evolution. However, the relative importance of positive selection versus relaxation of constraints in the functional divergence of gene copies is still under debate. Plant MADS-box genes encode transcriptional regulators key in various aspects of development and have undergone extensive duplications to form a large family. We recovered 104 MADS sequences from the Arabidopsis genome. Bayesian phylogenetic trees recover type II lineage as a monophyletic group and resolve a branching sequence of monophyletic groups within this lineage. The type I lineage is comprised of several divergent groups. However, contrasting gene structure and patterns of chromosomal distribution between type I and II sequences suggest that they had different evolutionary histories and support the placement of the root of the gene family between these two groups. Site-specific and site-branch analyses of positive Darwinian selection (PDS) suggest that different selection regimes could have affected the evolution of these lineages. We found evidence for PDS along the branch leading to flowering time genes that have a direct impact on plant fitness. Sites with high probabilities of having been under PDS were found in the MADS and K domains, suggesting that these played important roles in the acquisition of novel functions during MADS-box diversification. Detected sites are targets for further experimental analyses. We argue that adaptive changes in MADS-domain protein sequences have been important for their functional divergence, suggesting that changes within coding regions of transcriptional regulators have influenced phenotypic evolution of plants.

Generation and analysis of an artificial gene dosage series in tomato to study the mechanisms by which the cloned quantitative trait locus fw2.2 controls fruit size

It has been proposed that fw2.2 encodes a negative fruit-growth regulator that underlies natural fruit-size variation in tomato (Lycopersicon spp.) via heterochronic allelic variation of fw2.2 expression, rather than by variation in the structural protein itself. To further test the negative regulator and the transcriptional control hypotheses, a gene dosage series was constructed, which produced a wider range of fw2.2 transcript accumulation than can be found in natural tomato populations. Fruit developmental analyses revealed that fw2.2 transcript levels were highly correlated (negatively) with fruit mass, supporting the negative regulator and transcriptional regulation hypotheses. Further, the effect of fw2.2 on fruit mass was mediated by repressing three- and two-dimensional cell division in placental and pericarp tissues, respectively. Finally, fw2.2 had little effect on fertility and seed size/number, indicating that fruit size effects of fw2.2 are due largely to expression in the maternal tissues of developing fruit and not mediated through fertility or seed-setting-related processes.

Developmental morphology of normal and atypical flowers ofPhilodendron insigne(Araceae): a new case of homeosis

The early stages of development of the inflorescence of Philodendron insigne were examined using scanning electron microscopy. Pistillate flowers are initiated on the lower portion of the inflorescence and staminate flowers are initiated on the distal portion. The male flowers have three to five stamens. The female flowers have a multilocular ovary consisting of three to five locules. A transition zone consisting of sterile male flowers and atypical bisexual flowers with fused or free carpels and staminodes is located between the male and female floral zones. Generally, the portion of the bisexual flower facing the male zone forms stamens, and the portion facing the female zone develops one or two carpels. In P. insigne, the incomplete separation of staminodes from the gynoecial portion of the whorl shows that the staminodes and carpels belong to the same whorl. The atypical bisexual flowers of P. insigne are believed to be a case of homeosis where carpels have been replaced by sterile stamens on the same whorl. However, there is no regularity in the number of organs involved in the homeotic transformation taking place in atypical bisexual flowers. The presence of atypical bisexual flowers may correspond to a morphogenetic gradient at the level of the inflorescence as a whole.Key words: flower, development, gradient, inflorescence.

Contrasting evolutionary forces in the Arabidopsis thaliana floral developmental pathway

The floral developmental pathway in Arabidopsis thaliana is composed of several interacting regulatory genes, including the inflorescence architecture gene TERMINAL FLOWER1 (TFL1), the floral meristem identity genes LEAFY (LFY), APETALA1 (AP1), and CAULIFLOWER (CAL), and the floral organ identity genes APETALA3 (AP3) and PISTILLATA (PI). Molecular population genetic analyses of these different genes indicate that the coding regions of AP3 and PI, as well as AP1 and CAL, share similar levels and patterns of nucleotide diversity. In contrast, the coding regions of TFL1 and LFY display a significant reduction in nucleotide variation, suggesting that these sequences have been subjected to a recent adaptive sweep. Moreover, the promoter of TFL1, unlike its coding region, displays high levels of diversity organized into two distinct haplogroups that appear to be maintained by selection. These results suggest that patterns of molecular evolution differ among regulatory genes in this developmental pathway, with the earlier acting genes exhibiting evidence of adaptive evolution.

Low rates of silent substitution in nuclear genes of two distantly related Scrophulariaceae (Antirrhinum and Verbascum)

Low levels of genetic diversity and divergence at nuclear loci have previously been observed for cycloidea and fil1-like genes within and between several Antirrhinum species, and divergence at these loci is also low between species in genera at different levels of relatedness in the former family Scrophulariaceae (Digitalis and Verbascum). The low divergence values are surprising, because (based on the sequences of chloroplast loci) the Scrophulariaceae are thought to be polyphyletic, with two anciently diverged clades, and the species we compared belonged to the two different clades. Here, we extend our studies of sequence divergence to more nuclear genes: fil2, far, globosa, and ADH: Detailed studies revealed that in Antirrhinum these genes belong to gene families. Low levels of divergence between Antirrhinum and Verbascum were observed for four of the loci studied, fil2-1, fil2-2, far-L, and globosa, similar to our previous observations. We discuss hypotheses to explain these low synonymous divergence values. For Adh, no cases of very similar sequences were found, but, rather, our sequences from the three different genera (Antirrhinum, Digitalis, and Verbascum) were all very diverged. Repeated gene duplication and loss of elements in the Adh gene family is likely in these lineages, making it impossible to determine orthology of the Adh genes.

Accelerated regulatory gene evolution in an adaptive radiation

The disparity between rates of morphological and molecular evolution remains a key paradox in evolutionary genetics. A proposed resolution to this paradox has been the conjecture that morphological evolution proceeds via diversification in regulatory loci, and that phenotypic evolution may correlate better with regulatory gene divergence. This conjecture can be tested by examining rates of regulatory gene evolution in species that display rapid morphological diversification within adaptive radiations. We have isolated homologues to the Arabidopsis APETALA3 (ASAP3/TM6) and APETALA1 (ASAP1) floral regulatory genes and the CHLOROPHYLL A/B BINDING PROTEIN9 (ASCAB9) photosynthetic structural gene from species in the Hawaiian silversword alliance, a premier example of plant adaptive radiation. We have compared rates of regulatory and structural gene evolution in the Hawaiian species to those in related species of North American tarweeds. Molecular evolutionary analyses indicate significant increases in nonsynonymous relative to synonymous nucleotide substitution rates in the ASAP3/TM6 and ASAP1 regulatory genes in the rapidly evolving Hawaiian species. By contrast, no general increase is evident in neutral mutation rates for these loci in the Hawaiian species. An increase in nonsynonymous relative to synonymous nucleotide substitution rate is also evident in the ASCAB9 structural gene in the Hawaiian species, but not to the extent displayed in the regulatory loci. The significantly accelerated rates of regulatory gene evolution in the Hawaiian species may reflect the influence of allopolyploidy or of selection and adaptive divergence. The analyses suggest that accelerated rates of regulatory gene evolution may accompany rapid morphological diversification in adaptive radiations.

Molecular evolution of the teosinte branched gene among maize and related grasses

Several authors have proposed that changes in a small number of regulatory genes may be sufficient for the evolution of novel morphologies. Recent analyses have indicated that teosinte branched1 (tb1), a putative bHLH transcription factor, played such a role during the morphological evolution of maize from its wild ancestor, teosinte. To address whether or not tb1 played a similar role during the evolution of the Andropogoneae, the tribe to which maize belongs, and to examine the rate and pattern of tb1 evolution within this tribe, we analyzed tb1-like sequences from 23 members of the Andropogoneae and five other grasses. Our analysis revealed that the TB1 protein evolves slowly within three conserved domains but rapidly outside these domains. The nonconserved regions of the gene are characterized by both a high nonsynonymous substitution rate and frequent indels. The ratio of nonsynonymous substitutions per nonsynonymous site (d(N)) to synonymous substitutions per synonymous site (d(S)) was not significantly greater than 1.0, providing no evidence for positive selection. However, the d(N)/d(S) ratio varied significantly among lineages and was high compared with those of other plant nuclear genes. Variation in the d(N)/d(S) ratio among the Andropogoneae could be explained by unequal levels of purifying selection among lineages. Consistent with this interpretation, the rate of nonsynonymous substitution differed along several lineages, while the synonymous substitution rate did not differ significantly. Finally, using tb1, we examined phylogenetic relationships within the Andropogoneae. The phylogeny suggests that the tribe underwent a rapid radiation during its early history and that the monoecious Andropogoneae are polyphyletic.

Memory and DNA

A model is presented for the storage of long-term memory. In our model consolidation takes place by specific DNA sequences. These DNA sequences are obtained by the recombination of DNA in a similar way to that during meiosis and the production of immunological antibodies. DNA has the potential of the production of large numbers of specific DNA sequences. These sequences can be attached to images of neural networks. The following considerations lead to the theory: (1)Most of the DNA is not used: approximately 3% of our DNA is used. (2)There are no cell divisions in the brain after adulthood is reached. Structural DNA arrangements will not be altered nor disrupted as a consequence of cell division and mitosis. (3)Chromosomal pairing is demonstrated in the brain, which could indicate the exchange of DNA. In addition, in our first survey experiments we found a positive reaction of components of the synaptonemal complex (SC) in the nuclei of brain cells. The SC is highly meiosis specific and plays a major role in genetic recombination.

Evolution of Hoxa-11 in lineages phylogenetically positioned along the fin-limb transition

HOXA11 is a transcription factor implicated in paired appendage development. To identify signatures of evolutionary change in the structural, and putative functional, domains of HOXA11, we studied its evolution in tetrapod and nontetrapod lineages that represent approximately 1.5 billion years of evolutionary time. Here, Hoxa-11 gene proper sequences were determined for frog (Xenopus tropicalis), coelacanth (Latimeria chalumnae), common zebrafish (Danio rerio; Hoxa-11a and Hoxa-11b paralogs), and giant zebrafish (D. aequipinnatus; Hoxa-11b) and aligned against previously published Hoxa-11 sequences of human, mouse, chick, and newt. Based on aligned Hoxa-11 amino acid sequences, the protein was demarcated into three segments: Domains I (N-terminal) and III (homeobox + C-terminal), which varied slightly in rates and patterns of evolution, and a variable, overall hydrophilic region (HyD), which partially overlaps with Domain I. As judged by character reconstructions of HOXA11 Domains I and III, no significant changes in rates of coding sequence evolution occurred in tetrapods (frog and chick), relative to coelacanth (a lobe-finned fish), i.e., across the fin-limb transition. Accelerated rates of Hoxa-11 coding sequence evolution were observed for the mammalian and newt lineages. This was shown to be a gene-specific phenomenon. The duplicated Hoxa-11a and Hoxa-11b genes of zebrafish exhibited accelerated rates of evolution and accumulated substitutions at sites that are conserved among coelacanth and all tetrapods examined. Amino acid sequence comparisons of the HyD of HOXA11 suggested that a putative repressor subdomain, containing stretches of consecutive alanine residues, emerged within the tetrapods. A high degree of nucleotide conservation in the 5' half of the Hoxa-11 intron was observed for tetrapod and nontetrapod lineages. Using electrophoretic mobility shift assays, a 35-bp intron sequence, which is 100% conserved in all Hoxa-11 loci except for the zebrafish Hoxa-11a paralog, was found to bind protein(s) in HeLa and chick whole-cell extracts.

MADS-box gene evolution beyond flowers: expression in pollen, endosperm, guard cells, roots and trichomes

MADS-box genes encode transcriptional regulators involved in diverse aspects of plant development. Here we describe the cloning and mRNA spatio-temporal expression patterns of five new MADS-box genes from Arabidopsis: AGL16, AGL18, AGL19, AGL27 and AGL31. These genes will probably become important molecular tools for both evolutionary and functional analyses of vegetative structures. We mapped our data and previous expression patterns onto a new MADS-box phylogeny. These analyses suggest that the evolution of the MADS-box family has involved a rapid and simultaneous functional diversification in vegetative as well as reproductive structures. The hypothetical ancestral genes had broader expression patterns than more derived ones, which have been co-opted for putative specialized functions as suggested by their expression patterns. AGL27 and AGL31, which are closely related to the recently described flowering-time gene FLC (previously AGL25), are expressed in most plant tissues. AGL19 is specifically expressed in the outer layers of the root meristem (lateral root cap and epidermis) and in the central cylinder cells of mature roots. AGL18, which is most similar in sequence to the embryo-expressed AGL15 gene, is expressed in the endosperm and in developing male and female gametophytes, suggesting a role for AGL18 that is distinct from previously characterized MADS-box genes. Finally, AGL16 RNA accumulates in leaf guard cells and trichomes. Our new phylogeny reveals seven new monophyletic clades of MADS-box sequences not specific to flowers, suggesting that complex regulatory networks involving several MADS-box genes, similar to those that control flower development, underlie development of vegetative structures.

Formation of regulator/target gene relationships during evolution

The Foxn1-like forkhead/winged-helix transcription factor genes have been maintained in single copy throughout chordate evolution. Among other functions, Foxn1 (formerly known as Whn) regulates the expression of hair keratin genes in the hair follicle, which represents an evolutionarily novel organ characteristic of mammals. We show here that fish and mouse Foxn1-like genes are functionally equivalent in hair keratin gene activation, suggesting the absence of functionally relevant changes over the course of several hundred million years of vertebrate evolution. In contrast, the Foxn1-like gene from the cephalochordate Branchiostoma lanceolatum is inactive in this assay because of changes in the region located N-terminal to DNA binding and transcriptional activation domains of the protein. Our results indicate that functionally relevant changes in cis-regulatory regions are not necessarily accompanied by corresponding changes in transcription factor proteins in the formation of evolutionarily novel regulator/target gene relationships.

The molecular population genetics of regulatory genes

Regulatory loci, which may encode both trans acting proteins as well as cis acting promoter regions, are crucial components of an organism's genetic architecture. Although evolution of these regulatory loci is believed to underlie the evolution of numerous adaptive traits, there is little information on natural variation of these genes. Recent molecular population genetic studies, however, have provided insights into the extent of natural variation at regulatory genes, the evolutionary forces that shape them and the phenotypic effects of molecular regulatory variants. These recent analyses suggest that it may be possible to study the molecular evolutionary ecology of regulatory diversification by examining both the extent and patterning of regulatory gene diversity, the phenotypic effects of molecular variation at these loci and their ecological consequences.

The LAMB1 gene from the clubmoss, Lycopodium annotinum, is a divergent MADS-box gene, expressed specifically in sporogenic structures

Transcription factors encoded by the large MADS-box gene family have important developmental functions in angiosperms, the flowering plants. Mutations in certain MADS-box genes are known to cause homeotic alterations in floral organ identity, and the establishment of floral organ identity is the most well-studied developmental process in which MADS-box genes are known to function. Our interest is in the potential connection between the duplication history of this gene family and the evolutionary origin of the structures that the different MADS-box genes developmentally regulate in plants. Previous studies have demonstrated that the origin of the MADS-box genes that control floral organ identity predate the evolutionary origin of the flower itself, since gymnosperms have genes that are orthologous to angiosperm floral homeotic MADS-box genes, whereas ferns appear to lack such genes. Here we report on the isolation of a MADS-box gene from Lycopodium annotinum, which belongs to the clubmosses, the phylogenetic sister group to other vascular plants. The gene, LAMB1, in the sporophyte is expressed exclusively in the reproductive structure, the strobilus, during sporogenesis. LAMB1 is similar to other plant MADS-box genes in that it contains a MADS-box as well as a second conserved element, a K-box. However, it differs in length and in exon/intron structure in the region between the MADS- and K-box, and also in the length and structure of the C-terminal region. A phylogenetic analysis indicates that LAMB1 is not closely related to other plant-type MADS-box genes, and may represent one of the basal branches in the phylogenetic tree of plant MADS-box genes.

Molecular evolution of flower development

Flowers, as reproductive structures of the most successful group of land plants, have been a central focus of study for both evolutionists and ecologists. Recent advances in unravelling the genetics of flower development have provided insight into the evolution of floral structures among angiosperms. The study of the evolution of genes that control floral morphogenesis permits us to draw inferences on the diversification of developmental systems, the origin of floral organs and the selective forces that drive evolutionary change among these plant reproductive structures.

Towards a new evolutionary synthesis

New concepts and information from molecular developmental biology, systematics, geology and the fossil record of all groups of organisms, need to be integrated into an expanded evolutionary synthesis. These fields of study show that large-scale evolutionary phenomena cannot be understood solely on the basis of extrapolation from processes observed at the level of modern populations and species. Patterns and rates of evolution are much more varied than had been conceived by Darwin or the evolutionary synthesis, and physical factors of the earth's history have had a significant, but extremely varied, impact on the evolution of life.

Genome evolution in polyploids

Polyploidy is a prominent process in plants and has been significant in the evolutionary history of vertebrates and other eukaryotes. In plants, interdisciplinary approaches combining phylogenetic and molecular genetic perspectives have enhanced our awareness of the myriad genetic interactions made possible by polyploidy. Here, processes and mechanisms of gene and genome evolution in polyploids are reviewed. Genes duplicated by polyploidy may retain their original or similar function, undergo diversification in protein function or regulation, or one copy may become silenced through mutational or epigenetic means. Duplicated genes also may interact through inter-locus recombination, gene conversion, or concerted evolution. Recent experiments have illuminated important processes in polyploids that operate above the organizational level of duplicated genes. These include inter-genomic chromosomal exchanges, saltational, non-Mendelian genomic evolution in nascent polyploids, inter-genomic invasion, and cytonuclear stabilization. Notwithstanding many recent insights, much remains to be learned about many aspects of polyploid evolution, including: the role of transposable elements in structural and regulatory gene evolution; processes and significance of epigenetic silencing; underlying controls of chromosome pairing; mechanisms and functional significance of rapid genome changes; cytonuclear accommodation; and coordination of regulatory factors contributed by two, sometimes divergent progenitor genomes. Continued application of molecular genetic approaches to questions of polyploid genome evolution holds promise for producing lasting insight into processes by which novel genotypes are generated and ultimately into how polyploidy facilitates evolution and adaptation.

Reconstructing the history of human limb development: lessons from birth defects

A major goal of biology has been to understand the developmental mechanisms behind evolutionary trends. This has led to a growing interest in studying the molecular basis of the evolution of developmental programs such as those mediating the diversification of tetrapod limbs. Over the last 10 y, it has become clear that the genes and general developmental programs used to build a limb are strongly conserved among widely disparate species. This finding suggests that altered regulation of the timing and locations of developmental events may be responsible for the morphologic variation observed among some species. However, genetic analyses of the regulatory regions of genes controlling vertebrate developmental programs are very limited. Characterization of the genetic basis of human birth defects of the limb provides an opportunity to dissect the developmental programs used to modify the architecture of the hominoid limb. This may allow us to assess the relative contributions of altered gene regulation to morphologic variation among species and reconstruct the evolutionary history of the hominid limb. Such insight is also important because morphologic differences in the hominid upper limb have been correlated with the use of tools, and tool making is often regarded as the milestone that marked the emergence of the genus Homo.

Homeobox genes, fossils, and the origin of species

Ever since Darwin there has been a history of debate on the tempo and mode of evolution. Is speciation a gradual process involving the accumulation of minute variations extant within a species, or is it rapid, the result of major organismal reorganization? Does one define a species on the basis of genes, morphology, or geographic or reproductive isolation? In this communication I present a model of evolutionary change that is based on the Mendelian inheritance of mutations in regulatory genes and the fact that most nonlethal mutations arise in the recessive state. Since the new recessive allele will spread through many generations without expression until there is a critical mass of heterozygotes capable of producing homozygotes for the mutation, the novel feature thus produced will appear abruptly in the population and in more than one individual. This picture of punctuation is consistent with the fossil record, which typically fails to provide evidence of smoothly transitional states of morphological change. Given that the first of their kind in the fossil record are organisms in which their novel characteristics are often more fully expressed or complex than in their descendants, it would seem that, after the mutation involving a regulatory gene is introduced, the general tendency is for its effects to become diminished. Among the implications for speciation is that this process does not depend on either reproductive isolation or genetic incompatibility. Rather, barring effects on reproductive organs or behavior, homozygotes for a novelty should be able to breed with heterozygotes and homozygotes for the wild state of the original population. This, in turn, suggests that the species barrier between individuals is probably a matter of mate recognition.

Molecular population genetics of floral homeotic loci. Departures from the equilibrium-neutral model at the APETALA3 and PISTILLATA genes of Arabidopsis thaliana

Molecular variation in genes that regulate development provides insights into the evolutionary processes that shape the diversification of morphogenetic pathways. Intraspecific sequence variation at the APETALA3 and PISTILLATA floral homeotic genes of Arabidopsis thaliana was analyzed to infer the extent and nature of diversity at these regulatory loci. Comparison of AP3 and PI diversity with three previously studied genes revealed several features in the patterning of nucleotide polymorphisms common between Arabidopsis nuclear loci, including an excess of low-frequency nucleotide polymorphisms and significantly elevated levels of intraspecific replacement variation. This pattern suggests that A. thaliana has undergone recent, rapid population expansion and now exists in small, inbred subpopulations. The elevated intraspecific replacement levels may thus represent slightly deleterious polymorphisms that differentiate distinct ecotypes. The distribution of replacement and synonymous changes in AP3 and PI core and noncore functional domains also indicates differences in the patterns of molecular evolution between these interacting floral regulatory genes.