Rapid evolution of a pollen-specific oleosin-like gene family from Arabidopsis thaliana and closely related species

Identification of key genes for triacylglycerol biosynthesis and storage in herbaceous peony (Paeonia lactifolra Pall.) seeds based on full-length transcriptome

BACKGROUND

The herbaceous peony (Paeonia lactiflora Pall.) is extensively cultivated in China due to its root being used as a traditional Chinese medicine known as 'Radix Paeoniae Alba'. In recent years, it has been discovered that its seeds incorporate abundant unsaturated fatty acids, thereby presenting a potential new oilseed plant. Surprisingly, little is known about the full-length transcriptome sequencing of Paeonia lactiflora, limiting research into its gene function and molecular mechanisms.

RESULTS

A total of 484,931 Reads of Inserts (ROI) sequences and 1,455,771 full-Length non-chimeric reads (FLNC) sequences were obtained for CDS prediction, TF analysis, SSR analysis and lncRNA identification. In addition, gene function annotation and gene structure analysis were performed. A total of 4905 transcripts were related to lipid metabolism biosynthesis pathway, belonging to 28 enzymes. We use these data to identify 10 oleosin (OLE) and 5 diacylglycerol acyltransferase (DGAT) gene members after de-redundancy. The analysis of physicochemical properties and secondary structure showed them similarity in gene family respectively. The phylogenetic analysis showed that the distribution of OLE and DGAT family members was roughly the same as that of Arabidopsis. Quantitative real-time polymerase chain reaction (qRT-PCR) analyses revealed expression changes in different seed development stages, and showed a trend of increasing and then decreasing.

CONCLUSION

In summary, these results provide new insights into the molecular mechanism of triacylglycerol (TAG) biosynthesis and storage during the seedling stage in Paeonia lactiflora. It provides theoretical references for selecting and breeding oil varieties and understanding the functions of oil storage as well as lipid synthesis related genes in Paeonia lactiflora.

Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

Molecular characterization of oleosin genes in Cyperus esculentus, a Cyperaceae plant producing oil in underground tubers

KEY MESSAGE

CeOLE genes exhibit a tuber-predominant expression pattern and their mRNA/protein abundances are positively correlated with oil accumulation during tuber development. Overexpression could significantly increase the oil content of tobacco leaves. Oleosins (OLEs) are abundant structural proteins of lipid droplets (LDs) that function in LD formation and stabilization in seeds of oil crops. However, little information is available on their roles in vegetative tissues. In this study, we present the first genome-wide characterization of the oleosin family in tigernut (Cyperus esculentus L., Cyperaceae), a rare example accumulating high amounts of oil in underground tubers. Six members identified represent three previously defined clades (i.e. U, SL and SH) or six out of seven orthogroups (i.e. U, SL1, SL2, and SH1-3) proposed in this study. Comparative genomics analysis reveals that lineage-specific expansion of Clades SL and SH was contributed by whole-genome duplication and dispersed duplication, respectively. Moreover, presence of SL2 and SH3 in Juncus effuses implies their appearance sometime before Cyperaceae-Juncaceae divergence, whereas SH2 appears to be Cyperaceae specific. Expression analysis showed that CeOLE genes exhibit a tuber-predominant expression pattern and transcript levels are considerably more abundant than homologs in the close relative Cyperus rotundus. Moreover, CeOLE mRNA and protein abundances were shown to positively correlate with oil accumulation during tuber development. Additionally, two dominant isoforms (i.e. CeOLE2 and -5) were shown to locate in LDs as well as the endoplasmic reticulum of tobacco (Nicotiana benthamiana) leaves, and are more likely to function in homo and heteromultimers. Furthermore, overexpression of CeOLE2 and -5 in tobacco leaves could significantly increase the oil content, supporting their roles in oil accumulation. These findings provide insights into lineage-specific family evolution and putative roles of CeOLE genes in oil accumulation of vegetative tissues, which facilitate further genetic improvement for tigernut.

Metabolomics and genetics of reproductive bud development in Ficus carica var. sativa (edible fig) and in Ficus carica var. caprificus (caprifig): similarities and differences

In figs, reproductive biology comprises cultivars requiring or not pollination, with female trees (edible fig) and male trees (caprifig) bearing different types of fruits. Metabolomic and genetic studies may clarify bud differentiation mechanisms behind the different fruits. We used a targeted metabolomic analysis and genetic investigation through RNA sequence and candidate gene investigation to perform a deep analysis of buds of two fig cultivars, 'Petrelli' (San Pedro type) and 'Dottato' (Common type), and one caprifig. In this work, proton nuclear magnetic resonance (¹H NMR-based metabolomics) has been used to analyze and compare buds of the caprifig and the two fig cultivars collected at different times of the season. Metabolomic data of buds collected on the caprifig, 'Petrelli', and 'Dottato' were treated individually, building three separate orthogonal partial least squared (OPLS) models, using the "y" variable as the sampling time to allow the identification of the correlations among metabolomic profiles of buds. The sampling times revealed different patterns between caprifig and the two edible fig cultivars. A significant amount of glucose and fructose was found in 'Petrelli', differently from 'Dottato', in the buds in June, suggesting that these sugars not only are used by the ripening brebas of 'Petrelli' but also are directed toward the developing buds on the current year shoot for either a main crop (fruit in the current season) or a breba (fruit in the successive season). Genetic characterization through the RNA-seq of buds and comparison with the literature allowed the identification of 473 downregulated genes, with 22 only in profichi, and 391 upregulated genes, with 21 only in mammoni.

The lineage-specific evolution of the oleosin family in Theaceae

Oleosins play essential roles in stabilization of lipid droplets (LDs) and seed oil production. However, evolution of this gene family has not been reported in Theaceae, a large plant family that contains many important tea and oil tea species. In this study, a total of 65 oleosin genes were identified in nine genome-sequenced Theaceae species. Among these genomes, the gene number of oleosin showed significant difference, with Camellia sinensis var. sinensis cv. Shuchazao and Camellia lanceoleosa displayed more oleosin numbers than other species. Phylogenetic analyses revealed that Theaceae oleosin genes were classified into three clades (U, SL, SH) respectively. Proteins within the same clade had similar gene structure and motif composition. Segmental duplication was the primary driving force for the evolution of oleosin genes in Shuchazao (SCZ), Huangdan (HD), C.lanceoleosa (Cla), and wild tea (DASZ). Synteny analysis showed that most oleosin genes displayed inter-species synteny among tea and oil tea species. Expression analysis demonstrated that oleosin genes were specifically expressed in seed and kernel of Huangdan (HD) and C.lanceoleosa. Moreover, expression divergence was observed in paralogous pairs and ∼1-2 oleosin genes in each clade have become activate. This study leads to a comprehensive understanding of evolution of oleosin family in Theaceae, and provides a rich resource to further address the functions of oleosin in tea and oil tea species.

Pollen Coat Proteomes of Arabidopsis thaliana, Arabidopsis lyrata, and Brassica oleracea Reveal Remarkable Diversity of Small Cysteine-Rich Proteins at the Pollen-Stigma Interface

The pollen coat is the outermost domain of the pollen grain and is largely derived from the anther tapetum, which is a secretory tissue that degenerates late in pollen development. By being localised at the interface of the pollen-stigma interaction, the pollen coat plays a central role in mediating early pollination events, including molecular recognition. Amongst species of the Brassicaceae, a growing body of data has revealed that the pollen coat carries a range of proteins, with a number of small cysteine-rich proteins (CRPs) being identified as important regulators of the pollen-stigma interaction. By utilising a state-of-the-art liquid chromatography/tandem mass spectrometry (LC-MS/MS) approach, rich pollen coat proteomic profiles were obtained for Arabidopsis thaliana, Arabidopsis lyrata, and Brassica oleracea, which greatly extended previous datasets. All three proteomes revealed a strikingly large number of small CRPs that were not previously reported as pollen coat components. The profiling also uncovered a wide range of other protein families, many of which were enriched in the pollen coat proteomes and had functions associated with signal transduction, cell walls, lipid metabolism and defence. These proteomes provide an excellent source of molecular targets for future investigations into the pollen-stigma interaction and its potential evolutionary links to plant-pathogen interactions.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Identification of favorable SNP alleles and candidate genes responsible for inflorescence-related traits via GWAS in chrysanthemum

81 SNPs were identified for three inflorescence-related traits, in which 15 were highly favorable. Two dCAPS markers were developed for future MAS breeding, and six candidate genes were predicted. Chrysanthemum is a leading ornamental species worldwide and demonstrates a wealth of morphological variation. Knowledge about the genetic basis of its phenotypic variation for key horticultural traits can contribute to its effective management and genetic improvement. In this study, we conducted a genome-wide association study (GWAS) based on two years of phenotype data and a set of 92,617 single nucleotide polymorphisms (SNPs) using a panel of 107 diverse cut chrysanthemums to dissect the genetic control of three inflorescence-related traits. A total of 81 SNPs were significantly associated with the three inflorescence-related traits (capitulum diameter, number of ray florets and flowering time) in at least one environment, with an individual allele explaining 22.72-38.67% of the phenotypic variation. Fifteen highly favorable alleles were identified for the three target traits by computing the phenotypic effect values for the stable associations detected in 2 year-long trials at each locus. Dosage pyramiding effects of the highly favorable SNP alleles and significant linear correlations between highly favorable allele numbers and corresponding phenotypic performance were observed. Two highly favorable SNP alleles correlating to flowering time and capitulum diameter were converted to derived cleaved amplified polymorphic sequence (dCAPS) markers to facilitate future breeding. Finally, six putative candidate genes were identified that contribute to flowering time and capitulum diameter. These results serve as a foundation for analyzing the genetic mechanisms underlying important horticultural traits and provide valuable insights into molecular marker-assisted selection (MAS) in chrysanthemum breeding programs.

The Evolutionary Consequences of Selection at the Haploid Gametic Stage

As an immediate consequence of sexual reproduction, biphasic life cycles with alternating diploid and haploid phases are a common characteristic of sexually reproducing eukaryotes. Much of our focus in evolutionary biology has been directed toward dynamics in diploid or haploid populations, but we rarely consider selection occurring during both phases when studying evolutionary processes. One of the reasons for this apparent omission is the fact that many flowering plants and metazoans are predominantly diploid with a very short haploid gametic phase. While this gametic phase may be short, it can play a crucial role in fundamental processes including the rate of adaptation, the load of mutation, and the evolution of features such as recombination. In addition, if selection acts in different directions between the two phases, a genetic conflict will occur, impacting the maintenance of genetic variation. Here we provide an overview of theoretical and empirical studies investigating the importance of selection at the haploid gametic phase in predominantly diploid organisms and discuss future directions to improve our understanding of the underlying dynamics and the general implications of haploid selection.

An integrated proteomic and transcriptomic analysis of perivitelline fluid proteins in a freshwater gastropod laying aerial eggs

Proteins of the egg perivitelline fluid (PVF) that surrounds the embryo are critical for embryonic development in many animals, but little is known about their identities. Using an integrated proteomic and transcriptomic approach, we identified 64 proteins from the PVF of Pomacea maculata, a freshwater snail adopting aerial oviposition. Proteins were classified into eight functional groups: major multifunctional perivitellin subunits, immune response, energy metabolism, protein degradation, oxidation-reduction, signaling and binding, transcription and translation, and others. Comparison of gene expression levels between tissues showed that 22 PVF genes were exclusively expressed in albumen gland, the female organ that secretes PVF. Base substitution analysis of PVF and housekeeping genes between P. maculata and its closely related species Pomacea canaliculata showed that the reproductive proteins had a higher mean evolutionary rate. Predicted 3D structures of selected PVF proteins showed that some nonsynonymous substitutions are located at or near the binding regions that may affect protein function. The proteome and sequence divergence analysis revealed a substantial amount of maternal investment in embryonic nutrition and defense, and higher adaptive selective pressure on PVF protein-coding genes when compared with housekeeping genes, providing insight into the adaptations associated with the unusual reproductive strategy in these mollusks.

SIGNIFICANCE

There has been great interest in studying reproduction-related proteins as such studies may not only answer fundamental questions about speciation and evolution, but also solve practical problems of animal infertility and pest outbreak. Our study has demonstrated the effectiveness of an integrated proteomic and transcriptomic approach in understanding the heavy maternal investment of proteins in the eggs of a non-model snail, and how the reproductive proteins may have evolved during the transition from laying underwater eggs to aerial eggs.

Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes

Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.

The Pollen Coat Proteome: At the Cutting Edge of Plant Reproduction

The tapetum is a single layer of secretory cells which encloses the anther locule and sustains pollen development and maturation. Upon apoptosis, the remnants of the tapetal cells, consisting mostly of lipids and proteins, fill the pits of the sculpted exine to form the bulk of the pollen coat. This extracellular matrix forms an impermeable barrier that protects the male gametophyte from water loss and UV light. It also aids pollen adhesion and hydration and retains small signaling compounds involved in pollen-stigma communication. In this study, we have updated the list of the pollen coat's protein components and also discussed their functions in the context of sexual reproduction.

Analysis of Arabidopsis floral transcriptome: detection of new florally expressed genes and expansion of Brassicaceae-specific gene families

The flower is essential for sexual reproduction of flowering plants and has been extensively studied. However, it is still not clear how many genes are expressed in the flower. Here, we performed RNA-seq analysis as a highly sensitive approach to investigate the Arabidopsis floral transcriptome at three developmental stages. We provide evidence that at least 23, 961 genes are active in the Arabidopsis flower, including 8512 genes that have not been reported as florally expressed previously. We compared gene expression at different stages and found that many genes encoding transcription factors are preferentially expressed in early flower development. Other genes with expression at distinct developmental stages included DUF577 in meiotic cells and DUF220, DUF1216, and Oleosin in stage 12 flowers. DUF1216 and DUF577 are Brassicaceae specific, and together with other families experienced expansion within the Brassicaceae lineage, suggesting novel/greater roles in Brassicaceae floral development than other plants. The large dataset from this study can serve as a resource for expression analysis of genes involved in flower development in Arabidopsis and for comparison with other species. Together, this work provides clues regarding molecular networks underlying flower development.

Selection on plant male function genes identifies candidates for reproductive isolation of yellow monkeyflowers

Understanding the genetic basis of reproductive isolation promises insight into speciation and the origins of biological diversity. While progress has been made in identifying genes underlying barriers to reproduction that function after fertilization (post-zygotic isolation), we know much less about earlier acting pre-zygotic barriers. Of particular interest are barriers involved in mating and fertilization that can evolve extremely rapidly under sexual selection, suggesting they may play a prominent role in the initial stages of reproductive isolation. A significant challenge to the field of speciation genetics is developing new approaches for identification of candidate genes underlying these barriers, particularly among non-traditional model systems. We employ powerful proteomic and genomic strategies to study the genetic basis of conspecific pollen precedence, an important component of pre-zygotic reproductive isolation among yellow monkeyflowers (Mimulus spp.) resulting from male pollen competition. We use isotopic labeling in combination with shotgun proteomics to identify more than 2,000 male function (pollen tube) proteins within maternal reproductive structures (styles) of M. guttatus flowers where pollen competition occurs. We then sequence array-captured pollen tube exomes from a large outcrossing population of M. guttatus, and identify those genes with evidence of selective sweeps or balancing selection consistent with their role in pollen competition. We also test for evidence of positive selection on these genes more broadly across yellow monkeyflowers, because a signal of adaptive divergence is a common feature of genes causing reproductive isolation. Together the molecular evolution studies identify 159 pollen tube proteins that are candidate genes for conspecific pollen precedence. Our work demonstrates how powerful proteomic and genomic tools can be readily adapted to non-traditional model systems, allowing for genome-wide screens towards the goal of identifying the molecular basis of genetically complex traits.

Selection-driven evolution of sex-biased genes is consistent with sexual selection in Arabidopsis thaliana

Sex-biased genes are genes with a preferential or specific expression in one sex and tend to show an accelerated rate of evolution in animals. Various hypotheses--which are not mutually exclusive--have been put forth to explain observed patterns of rapid evolution. One possible explanation is positive selection, but this has been shown only in few animal species and mostly for male-specific genes. Here, we present a large-scale study that investigates evolutionary patterns of sex-biased genes in the predominantly self-fertilizing plant Arabidopsis thaliana. Unlike most animal species, A. thaliana does not possess sex chromosomes, its flowers develop both male and female sexual organs, and it is characterized by low outcrossing rates. Using cell-specific gene expression data, we identified genes whose expression is enriched in comparison with all other tissues in the male and female gametes (sperm, egg, and central cell), as well as in synergids, pollen, and pollen tubes, which also play an important role in reproduction. Genes specifically expressed in gametes and synergids show higher rates of protein evolution compared with the genome-wide average and no evidence for positive selection. In contrast, pollen- and pollen tube-specific genes not only have lower rates of protein evolution but also exhibit a higher proportion of adaptive amino acid substitutions. We show that this is the result of increased levels of purifying and positive selection among genes with pollen- and pollen tube-specific expression. The increased proportion of adaptive substitutions cannot be explained by the fact that pollen- and pollen tube-expressed genes are enriched in segmental duplications, are on average older, or have a larger effective population size. Our observations are consistent with prezygotic sexual selection as a result of interactions during pollination and pollen tube growth such as pollen tube competition.

Pollen-specific, but not sperm-specific, genes show stronger purifying selection and higher rates of positive selection than sporophytic genes in Capsella grandiflora

Selection on the gametophyte can be a major force shaping plant genomes as 7-11% of genes are expressed only in that phase and 60% of genes are expressed in both the gametophytic and sporophytic phases. The efficacy of selection on gametophytic tissues is likely to be influenced by sexual selection acting on male and female functions of hermaphroditic plants. Moreover, the haploid nature of the gametophytic phase allows selection to be efficient in removing recessive deleterious mutations and fixing recessive beneficial mutations. To assess the importance of gametophytic selection, we compared the strength of purifying selection and extent of positive selection on gametophyte- and sporophyte-specific genes in the highly outcrossing plant Capsella grandiflora. We found that pollen-exclusive genes had a larger fraction of sites under strong purifying selection, a greater proportion of adaptive substitutions, and faster protein evolution compared with seedling-exclusive genes. In contrast, sperm cell-exclusive genes had a smaller fraction of sites under strong purifying selection, a lower proportion of adaptive substitutions, and slower protein evolution compared with seedling-exclusive genes. Observations of strong selection acting on pollen-expressed genes are likely explained by sexual selection resulting from pollen competition aided by the haploid nature of that tissue. The relaxation of selection in sperm might be due to the reduced influence of intrasexual competition, but reduced gene expression may also be playing an important role.

Tandem oleosin genes in a cluster acquired in Brassicaceae created tapetosomes and conferred additive benefit of pollen vigor

During evolution, genomes expanded via whole-genome, segmental, tandem, and individual-gene duplications, and the emerged redundant paralogs would be eliminated or retained owing to selective neutrality or adaptive benefit and further functional divergence. Here we show that tandem paralogs can contribute adaptive quantitative benefit and thus have been retained in a lineage-specific manner. In Brassicaceae, a tandem oleosin gene cluster of five to nine paralogs encodes ample tapetum-specific oleosins located in abundant organelles called tapetosomes in flower anthers. Tapetosomes coordinate the storage of lipids and flavonoids and their transport to the adjacent maturing pollen as the coat to serve various functions. Transfer-DNA and siRNA mutants of Arabidopsis thaliana with knockout and knockdown of different tandem oleosin paralogs had quantitative and correlated loss of organized structures of the tapetosomes, pollen-coat materials, and pollen tolerance to dehydration. Complementation with the knockout paralog restored the losses. Cleomaceae is the family closest to Brassicaceae. Cleome species did not contain the tandem oleosin gene cluster, tapetum oleosin transcripts, tapetosomes, or pollen tolerant to dehydration. Cleome hassleriana transformed with an Arabidopsis oleosin gene for tapetum expression possessed primitive tapetosomes and pollen tolerant to dehydration. We propose that during early evolution of Brassicaceae, a duplicate oleosin gene mutated from expression in seed to the tapetum. The tapetum oleosin generated primitive tapetosomes that organized stored lipids and flavonoids for their effective transfer to the pollen surface for greater pollen vitality. The resulting adaptive benefit led to retention of tandem-duplicated oleosin genes for production of more oleosin and modern tapetosomes.

Selection is no more efficient in haploid than in diploid life stages of an angiosperm and a moss

The masking hypothesis predicts that selection is more efficient in haploids than in diploids, because dominant alleles can mask the deleterious effects of recessive alleles in diploids. However, gene expression breadth and noise can potentially counteract the effect of masking on the rate at which genes evolve. Land plants are ideal to ask whether masking, expression breadth, or expression noise dominate in their influence on the rate of molecular evolution, because they have a biphasic life cycle in which the duration and complexity of the haploid and diploid phase varies among organisms. Here, we generate and compile genome-wide gene expression, sequence divergence, and polymorphism data for Arabidopsis thaliana and for the moss Funaria hygrometrica to show that the evolutionary rates of haploid- and diploid-specific genes contradict the masking hypothesis. Haploid-specific genes do not evolve more slowly than diploid-specific genes in either organism. Our data suggest that gene expression breadth influence the evolutionary rate of phase-specific genes more strongly than masking. Our observations have implications for the role of haploid life stages in the purging of deleterious mutations, as well as for the evolution of ploidy.

Species-specific size expansion and molecular evolution of the oleosins in angiosperms

Oleosins are hydrophobic plant proteins thought to be important for the formation of oil bodies, which supply energy for seed germination and subsequent seedling growth. To better understand the evolutionary history and diversity of the oleosin gene family in plants, especially angiosperms, we systematically investigated the molecular evolution of this family using eight representative angiosperm species. A total of 73 oleosin members were identified, with six members in each of four monocot species and a greater but variable number in the four eudicots. A phylogenetic analysis revealed that the angiosperm oleosin genes belonged to three monophyletic lineages. Species-specific gene duplications, caused mainly by segmental duplication, led to the great expansion of oleosin genes and occurred frequently in eudicots after the monocot-eudicot divergence. Functional divergence analyses indicate that significant amino acid site-specific selective constraints acted on the different clades of oleosins. Adaptive evolution analyses demonstrate that oleosin genes were subject to strong purifying selection after their species-specific duplications and that rapid evolution occurred with a high degree of evolutionary dynamics in the pollen-specific oleosin genes. In conclusion, this study serves as a foundation for genome-wide analyses of the oleosins. These findings provide insight into the function and evolution of this gene family in angiosperms and pave the way for studies in other plants.

In Arabidopsis thaliana codon volatility scores reflect GC3 composition rather than selective pressure

BACKGROUND

Synonymous codon usage bias has typically been correlated with, and attributed to translational efficiency. However, there are other pressures on genomic sequence composition that can affect codon usage patterns such as mutational biases. This study provides an analysis of the codon usage patterns in Arabidopsis thaliana in relation to gene expression levels, codon volatility, mutational biases and selective pressures.

RESULTS

We have performed synonymous codon usage and codon volatility analyses for all genes in the A. thaliana genome. In contrast to reports for species from other kingdoms, we find that neither codon usage nor volatility are correlated with selection pressure (as measured by dN/dS), nor with gene expression levels on a genome wide level. Our results show that codon volatility and usage are not synonymous, rather that they are correlated with the abundance of G and C at the third codon position (GC3).

CONCLUSIONS

Our results indicate that while the A. thaliana genome shows evidence for synonymous codon usage bias, this is not related to the expression levels of its constituent genes. Neither codon volatility nor codon usage are correlated with expression levels or selective pressures but, because they are directly related to the composition of G and C at the third codon position, they are the result of mutational bias. Therefore, in A. thaliana codon volatility and usage do not result from selection for translation efficiency or protein functional shift as measured by positive selection.

The dynamic roles of intracellular lipid droplets: from archaea to mammals

During the past decade, there has been a paradigm shift in our understanding of the roles of intracellular lipid droplets (LDs). New genetic, biochemical and imaging technologies have underpinned these advances, which are revealing much new information about these dynamic organelles. This review takes a comparative approach by examining recent work on LDs across the whole range of biological organisms from archaea and bacteria, through yeast and Drosophila to mammals, including humans. LDs probably evolved originally in microorganisms as temporary stores of excess dietary lipid that was surplus to the immediate requirements of membrane formation/turnover. LDs then acquired roles as long-term carbon stores that enabled organisms to survive episodic lack of nutrients. In multicellular organisms, LDs went on to acquire numerous additional roles including cell- and organism-level lipid homeostasis, protein sequestration, membrane trafficking and signalling. Many pathogens of plants and animals subvert their host LD metabolism as part of their infection process. Finally, malfunctions in LDs and associated proteins are implicated in several degenerative diseases of modern humans, among the most serious of which is the increasingly prevalent constellation of pathologies, such as obesity and insulin resistance, which is associated with metabolic syndrome.

Organ and cell type-specific complementary expression patterns and regulatory neofunctionalization between duplicated genes in Arabidopsis thaliana

Duplicated genes can contribute to the evolution of new functions and they are common in eukaryotic genomes. After duplication, genes can show divergence in their sequence and/or expression patterns. Qualitative complementary expression, or reciprocal expression, is when only one copy is expressed in some organ or tissue types and only the other copy is expressed in others, indicative of regulatory subfunctionalization or neofunctionalization. From analyses of two microarray data sets with 83 different organ types, developmental stages, and cell types in Arabidopsis thaliana, we determined that 30% of whole-genome duplicate pairs and 38% of tandem duplicate pairs show reciprocal expression patterns. We reconstructed the ancestral state of expression patterns to infer that considerably more cases of reciprocal expression resulted from gain of a new expression pattern (regulatory neofunctionalization) than from partitioning of ancestral expression patterns (regulatory subfunctionalization). Pollen was an especially common organ type for expression gain, resulting in contrasting expression of some duplicates in pollen. Many of the gene pairs with reciprocal expression showed asymmetric sequence rate evolution, consistent with neofunctionalization, and the more rapidly evolving copy often showed a more restricted expression pattern. A gene with reciprocal expression in pollen, involved in brassinosteroid signal transduction, has evolved more rapidly than its paralog, and it shows evidence for a new function in pollen. This study indicates the evolutionary importance of reciprocal expression patterns between gene duplicates, showing that they are common, often associated with regulatory neofunctionalization, and may be a factor allowing for retention and divergence of duplicated genes.

Selection-driven divergence after gene duplication in Arabidopsis thaliana

Gene duplications are one of the most important mechanisms for the origin of evolutionary novelties. Even though various models of the fate of duplicated genes have been established, current knowledge about the role of divergent selection after gene duplication is rather limited. In this study, we analyzed sequence divergence in response to neo- and subfunctionalization of segmentally duplicated genes in the genome of Arabidopsis thaliana. We compared the genomes of A. thaliana and the poplar Populus trichocarpa to identify orthologous pairs of genes and their corresponding inparalogs. Maximum-likelihood analyses of the nonsynonymous and synonymous substitution rate ratio [Formula: see text] of pairs of A. thaliana inparalogs were used to detect differences in the evolutionary rates of protein coding sequences. We analyzed 1,924 A. thaliana paralogous pairs and our results indicate that around 6.9% show divergent ω values between the lineages for a fraction of sites. We observe an enrichment of regulatory sequences, a reduced level of co-expression and an increased number of substitutions that can be attributed to positive selection based on an McDonald-Kreitman type of analysis. Taken together, these results show that selection after duplication contributes substantially to gene novelties and hence functional divergence in plants.

Dynamic programming procedure for searching optimal models to estimate substitution rates based on the maximum-likelihood method

The substitution rate in a gene can provide valuable information for understanding its functionality and evolution. A widely used method to estimate substitution rates is the maximum-likelihood method implemented in the CODEML program in the PAML package. A limited number of branch models, chosen based on a priori information or an interest in a particular lineage(s), are tested, whereas a large number of potential models are neglected. A complementary approach is also needed to test all or a large number of possible models to search for the globally optional model(s) of maximum likelihood. However, the computational time for this search even in a small number of sequences becomes impractically long. Thus, it is desirable to explore the most probable spaces to search for the optimal models. Using dynamic programming techniques, we developed a simple computational method for searching the most probable optimal branch-specific models in a practically feasible computational time. We propose three search methods to find the optimal models, which explored O(n) (method 1) to O(n(2)) (method 2 and method 3) models when the given phylogeny has n branches. In addition, we derived a formula to calculate the number of all possible models, revealing the complexity of finding the optimal branch-specific model. We show that in a reanalysis of over 50 previously published studies, the vast majority obtained better models with significantly higher likelihoods than the conventional hypothesis model methods.

Arabidopsis and relatives as models for the study of genetic and genomic incompatibilities

The past few years have seen considerable advances in speciation research, but whether drift or adaptation is more likely to lead to genetic incompatibilities remains unknown. Some of the answers will probably come from not only studying incompatibilities between well-established species, but also from investigating incipient speciation events, to learn more about speciation as an evolutionary process. The genus Arabidopsis, which includes the widely used Arabidopsis thaliana, provides a useful set of model species for studying many aspects of population divergence. The genus contains both self-incompatible and incompatible species, providing a platform for studying the impact of mating system changes on genetic differentiation. Another important path to plant speciation is via formation of polyploids, and this can be investigated in the young allotetraploid species A. arenosa. Finally, there are many cases of intraspecific incompatibilities in A. thaliana, and recent progress has been made in discovering the genes underlying both F(1) and F(2) breakdown. In the near future, all these studies will be greatly empowered by complete genome sequences not only for all members of this relatively small genus, but also for many different individuals within each species.

Genome-wide evidence for efficient positive and purifying selection in Capsella grandiflora, a plant species with a large effective population size

Recent studies comparing genome-wide polymorphism and divergence in Drosophila have found evidence for a surprisingly high proportion of adaptive amino acid fixations, but results for other taxa are mixed. In particular, few studies have found convincing evidence for adaptive amino acid substitution in plants. To assess the generality of this finding, we have sequenced 257 loci in the outcrossing crucifer Capsella grandiflora, which has a large effective population size and low population structure. Using a new method that jointly infers selective and demographic effects, we estimate that 40% of amino acid substitutions were fixed by positive selection in this species, and we also infer a low proportion of slightly deleterious amino acid mutations. We contrast these estimates with those for a similar data set from the closely related Arabidopsis thaliana and find significantly higher rates of adaptive evolution and fewer nearly neutral mutations in C. grandiflora. In agreement with results for other taxa, genes involved in reproduction show the strongest evidence for positive selection in C. grandiflora. Taken together, these results imply that both positive and purifying selection are more effective in C. grandiflora than in A. thaliana, consistent with the contrasting demographic history and effective population sizes of these species.

Dynamic evolution of oryza genomes is revealed by comparative genomic analysis of a genus-wide vertical data set

Oryza (23 species; 10 genome types) contains the world's most important food crop - rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past approximately 15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.

Diverse cell signalling pathways regulate pollen-stigma interactions: the search for consensus

Siphonogamy, the delivery of nonmotile sperm to the egg via a pollen tube, was a key innovation that allowed flowering plants (angiosperms) to carry out sexual reproduction on land without the need for water. This process begins with a pollen grain (male gametophyte) alighting on and adhering to the stigma of a flower. If conditions are right, the pollen grain germinates to produce a pollen tube. The pollen tube invades the stigma and grows through the style towards the ovary, where it enters an ovule, penetrates the embryo sac (female gametophyte) and releases two sperm cells, one of which fertilizes the egg, while the other fuses with the two polar nuclei of the central cell to form the triploid endosperm. The events before fertilization (pollen-pistil interactions) comprise a series of complex cellular interactions involving a continuous exchange of signals between the haploid pollen and the diploid maternal tissue of the pistil (sporophyte). In recent years, significant progress has been made in elucidating the molecular identity of these signals and the cellular interactions that they regulate. Here we review our current understanding of the cellular and molecular interactions that mediate the earliest of these interactions between the pollen and the pistil that occur on or within the stigma - the 'pollen-stigma interaction'.

Rapid evolution of yeast centromeres in the absence of drive

To find the most rapidly evolving regions in the yeast genome we compared most of chromosome III from three closely related lineages of the wild yeast Saccharomyces paradoxus. Unexpectedly, the centromere appears to be the fastest-evolving part of the chromosome, evolving even faster than DNA sequences unlikely to be under selective constraint (i.e., synonymous sites after correcting for codon usage bias and remnant transposable elements). Centromeres on other chromosomes also show an elevated rate of nucleotide substitution. Rapid centromere evolution has also been reported for some plants and animals and has been attributed to selection for inclusion in the egg or the ovule at female meiosis. But Saccharomyces yeasts have symmetrical meioses with all four products surviving, thus providing no opportunity for meiotic drive. In addition, yeast centromeres show the high levels of polymorphism expected under a neutral model of molecular evolution. We suggest that yeast centromeres suffer an elevated rate of mutation relative to other chromosomal regions and they change through a process of "centromere drift," not drive.

Duplication and adaptive evolution of the COR15 genes within the highly cold-tolerant Draba lineage (Brassicaceae)

Plants have evolved diverse adaptive mechanisms that enable them to tolerate abiotic stresses, to varying degrees, and such stresses may have strongly influenced evolutionary changes at levels ranging from molecular to morphological. Previous studies on these phenomena have focused on the adaptive evolution of stress-related orthologous genes in specific lineages. However, heterogenetic evolution of the paralogous genes following duplication has only been examined in a very limited number of stress-response gene families. The COR15 gene encodes a low molecular weight protein that plays an important role in protecting plants from cold stresses. Although two different copies of this gene have been found in the model species, Arabidopsis thaliana, evolutionary patterns of this small gene family in plants have not been previously explored. In this study, we cloned COR15-like sequences and performed evolutionary analyses of these sequences (including those previously reported) in the highly cold-tolerant Draba lineage and related lineages of Brassicaceae. Our phylogenetic analyses indicate that all COR15-like sequences clustered into four clades that corresponded well to the morphological lineages. Gene conversions were found to have probably occurred before/during the divergence of Brassica and Draba lineage. However, repeated, independent duplications of this gene have occurred in different lineages of Brassicaceae. Further comparisons of all sequences suggest that there have been significant inter-lineage differences in evolutionary rates between the duplicated and original genes. We assessed the likelihood that the differences between two well-supported gene subfamilies that appear to have originated from a single duplication, COR15a and COR15b, within the Draba lineage have been driven by adaptive evolution. Comparisons of their non-synonymous/synonymous substitution ratios and rates of predicted amino acid changes indicate that these two gene groups are evolving under different selective pressures and may be functionally divergent. This functional divergence was confirmed by comparing site-specific shifts in evolution indexes of the two groups of predicted proteins. The evidence of differential selection and possible functional divergence suggests that the duplication may be of adaptive significance, with possible implications for the explosive diversification of the Draba lineage during the cooling Quaternary stages and the following worldwide colonization of arid alpine and artic regions.

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.

Comparative genetic mapping in Boechera stricta, a close relative of Arabidopsis

The angiosperm family Brassicaceae contains both the research model Arabidopsis (Arabidopsis thaliana) and the agricultural genus Brassica. Comparative genomics in the Brassicaceae has largely focused on direct comparisons between Arabidopsis and the species of interest. However, the reduced genome size and chromosome number (n = 5) of Arabidopsis complicates comparisons. Arabidopsis shows extensive genome and chromosome reshuffling compared to its close relatives Arabidopsis lyrata and Capsella rubella, both with n = 8. To facilitate comparative genomics across the Brassicaceae we recently outlined a system of 24 conserved chromosomal blocks based on their positions in an ancestral karyotype of n = 8, rather than by their position in Arabidopsis. In this report we use this system as a tool to understand genome structure and evolution in Boechera stricta (n = 7). B. stricta is a diploid, sexual, and highly self-fertilizing species occurring in mostly montane regions of western North America. We have created an F(2) genetic map of B. stricta based on 192 individuals scored at 196 microsatellite and candidate gene loci. Single-nucleotide polymorphism genotyping of 94 of the loci was done simultaneously using an Illumina bead array. The total map length is 725.8 cM, with an average marker spacing of 3.9 cM. There are no gaps greater than 19.3 cM. The chromosomal reduction from n = 8 to n = 7 and other genomic changes in B. stricta likely involved a pericentric inversion, a chromosomal fusion, and two reciprocal translocations that are easily visualized using the genomic blocks. Our genetic map will facilitate the analysis of ecologically relevant quantitative variation in Boechera. Sequence data from this article can be found in the GenBank/EMBL data libraries under accession numbers DU 667459 to DU 708532.

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).

The extracellular pollen coat in members of the Brassicaceae: composition, biosynthesis, and functions in pollination

I have used cellular and molecular genetic and bioinformatic approaches to characterise the components of the pollen coat in plants of the family Brassicaceae, including Arabidopsis thaliana and several brassicas including Brassica napus, B. oleracea, and B. rapa. The pollen coat in these species is mostly made up of a unique mixture of lipids that is highly enriched in acylated compounds, such as sterol esters and phospholipids. These acyl lipids are characterised by an unusually high degree of saturation. The fatty acids typically contain 70-90% saturated acyl residues such as myristate, palmitate, and stearate. The major sterol components of the pollen coat are saturated fatty acyl esters of stigmasterol, campesterol, and campestdienol. In addition to lipids, the second major component of the pollen coat is a specific group of proteins that is dominated by a family of proteins that we term pollenins. Although pollenins are by far the major protein components of the pollen coat of members of the Brassicaceae, proteomic analysis reveals that there are several additional protein components, including lipases, protein kinases, a pectin esterase, and a caleosin. The biosynthesis of these lipids and proteins and their significance for overall pollen function are reviewed and discussed.

Molecular Evolution of a Small Gene Family of Wound Inducible Kunitz Trypsin Inhibitors in Populus

Maximum likelihood models of codon substitutions were used to analyze the molecular evolution of a Kunitz trypsin inhibitor (KTI) gene family in Populus and Salix. The methods support previous assertions that the KTI genes comprise a rapidly evolving gene family. Models that allow for codon specific estimates of the ratio of nonsynonymous to synonymous substitutions (omega) among sites detect positive Darwinian selection at several sites in the KTI protein. In addition, branch-specific maximum likelihood models show that there is significant heterogeneity in omega among branches of the KTI phylogeny. In particular, omega is substantially higher following duplication than speciation. There is also evidence for significant rate heterogeneity following gene duplication, suggesting different evolutionary rates in newly arisen gene duplicates. The results indicate uneven evolutionary rates both between sites in the KTI protein and among different lineages in the KTI phylogeny, which is incompatible with a neutral model of sequence evolution.

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.

Evolution of reproductive proteins from animals and plants

Sexual reproduction is a fundamental biological process common among eukaryotes. Because of the significance of reproductive proteins to fitness, the diversity and rapid divergence of proteins acting at many stages of reproduction is surprising and suggests a role of adaptive diversification in reproductive protein evolution. Here we review the evolution of reproductive proteins acting at different stages of reproduction among animals and plants, emphasizing common patterns. Although we are just beginning to understand these patterns, by making comparisons among stages of reproduction for diverse organisms we can begin to understand the selective forces driving reproductive protein diversity and the functional consequences of reproductive protein evolution.

Sorghum bicolor's transcriptome response to dehydration, high salinity and ABA

Genome wide changes in gene expression were monitored in the drought tolerant C4 cereal Sorghum bicolor, following exposure of seedlings to high salinity (150 mM NaCl), osmotic stress (20% polyethylene glycol) or abscisic acid (125 microM ABA). A sorghum cDNA microarray providing data on 12,982 unique gene clusters was used to examine gene expression in roots and shoots at 3- and 27-h post-treatment. Expression of approximately 2200 genes, including 174 genes with currently unknown functions, of which a subset appear unique to monocots and/or sorghum, was altered in response to dehydration, high salinity or ABA. The modulated sorghum genes had homology to proteins involved in regulation, growth, transport, membrane/protein turnover/repair, metabolism, dehydration protection, reactive oxygen scavenging, and plant defense. Real-time PCR was used to quantify changes in relative mRNA abundance for 333 genes that responded to ABA, NaCl or osmotic stress. Osmotic stress inducible sorghum genes identified for the first time included a beta-expansin expressed in shoots, actin depolymerization factor, inositol-3-phosphate synthase, a non-C4 NADP-malic enzyme, oleosin, and three genes homologous to 9-cis-epoxycarotenoid dioxygenase that may be involved in ABA biosynthesis. Analysis of response profiles demonstrated the existence of a complex gene regulatory network that differentially modulates gene expression in a tissue- and kinetic-specific manner in response to ABA, high salinity and water deficit. Modulation of genes involved in signal transduction, chromatin structure, transcription, translation and RNA metabolism contributes to sorghum's overlapping but nonetheless distinct responses to ABA, high salinity, and osmotic stress. Overall, this study provides a foundation of information on sorghum's osmotic stress responsive gene complement that will accelerate follow up biochemical, QTL and comparative studies.

Species specificity in pollen-pistil interactions

For pollination to succeed, pollen must carry sperm through a variety of different floral tissues to access the ovules within the pistil. The pistil provides everything the pollen requires for success in this endeavor including distinct guidance cues and essential nutrients that allow the pollen tube to traverse enormous distances along a complex path to the unfertilized ovule. Although the pistil is a great facilitator of pollen function, it can also be viewed as an elaborate barrier that shields ovules from access from inappropriate pollen, such as pollen from other species. Each discrete step taken by pollen tubes en route to the ovules is a potential barrier point to ovule access and waste by inappropriate mates. In this review, we survey the current molecular understanding of how pollination proceeds, and ask to what extent is each step important for mate discrimination. As this field progresses, this synthesis of functional biology and evolutionary studies will provide insight into the molecular basis of the species barriers that maintain the enormous diversity seen in flowering plants.

Adaptive evolution in the SRZ chemoreceptor families of Caenorhabditis elegans and Caenorhabditis briggsae

We investigated the possibility of positive selection acting on members of the putative seven-pass chemoreceptor superfamily in Caenorhabditis elegans, which comprises approximately 1,300 genes encoding seven-pass G protein-coupled receptors (GPCRs). Using a maximum-likelihood approach, we conducted statistical tests for evidence of codon sites where the ratio of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (d(N)/d(S)) was >1. Evidence for positive selection was found only for the srz family, about which virtually nothing specific is known. We extended the annotation of the srz gene family, establishing gene models for 60 srz genes in C. elegans and 28 srz genes in Caenorhabditis briggsae. d(N)/d(S) ratios varied dramatically in different regions of the SRZ proteins, peaking in predicted extracellular regions. These regions included 23 sites where evidence of positive selection was highly significant, corresponding remarkably well with regions implicated in ligand binding in other GPCR family members. We interpret these results as indicating that the srz family is under positive selection, probably driven by ligand binding.