Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana

Genetic factors explaining anthocyanin pigmentation differences

BACKGROUND

Anthocyanins are important contributors to coloration across a wide phylogenetic range of plants. Biological functions of anthocyanins span from reproduction to protection against biotic and abiotic stressors. Owing to a clearly visible phenotype of mutants, the anthocyanin biosynthesis and its sophisticated regulation have been studied in numerous plant species. Genes encoding the anthocyanin biosynthesis enzymes are regulated by a transcription factor complex comprising MYB, bHLH and WD40 proteins.

RESULTS

A systematic comparison of anthocyanin-pigmented vs. non-pigmented varieties was performed within numerous plant species covering the taxonomic diversity of flowering plants. The literature was screened for cases in which genetic factors causing anthocyanin loss were reported. Additionally, transcriptomic data sets from four previous studies were reanalyzed to determine the genes possibly responsible for color variation based on their expression pattern. The contribution of different structural and regulatory genes to the intraspecific pigmentation differences was quantified. Differences concerning transcription factors are by far the most frequent explanation for pigmentation differences observed between two varieties of the same species. Among the transcription factors in the analyzed cases, MYB genes are significantly more prone to account for pigmentation differences compared to bHLH or WD40 genes. Among the structural genes, DFR genes are most often associated with anthocyanin loss.

CONCLUSIONS

These findings support previous assumptions about the susceptibility of transcriptional regulation to evolutionary changes and its importance for the evolution of novel coloration phenotypes. Our findings underline the particular significance of MYBs and their apparent prevalent role in the specificity of the MBW complex.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

Pathway position constrains the evolution of an ecologically important pathway in aspens (Populus tremula L.)

Many ecological interactions of aspens and their relatives (Populus spp.) are affected by products of the phenylpropanoid pathway synthesizing condensed tannins (CTs), whose production involves trade-offs with other ecologically important compounds and with growth. Genes of this pathway are candidates for investigating the role of selection on ecologically important, polygenic traits. We analysed sequences from 25 genes representing 10 steps of the CT synthesis pathway, which produces CTs used in defence and lignins used for growth, in 12 individuals of European aspen (Populus tremula). We compared these to homologs from P. trichocarpa, to a control set of 77 P. tremula genes, to genome-wide resequencing data and to RNA-seq expression levels, in order to identify signatures of selection distinct from those of demography. In Populus, pathway position exerts a strong influence on the evolution of these genes. Nonsynonymous diversity, divergence and allele frequency shifts (Tajima's D) were much lower than for synonymous measures. Expression levels were higher, and the direction of selection more negative, for upstream genes than for those downstream. Selective constraints act with increasing intensity on upstream genes, despite the presence of multiple paralogs in most gene families. Pleiotropy, expression level, flux control and codon bias appear to interact in determining levels and patterns of variation in genes of this pathway, whose products mediate a wide array of ecological interactions for this widely distributed species.

The factors affecting the evolution of the anthocyanin biosynthesis pathway genes in monocot and dicot plant species

BACKGROUND

The available data demonstrate that even in universal metabolic pathways, some species-specific regulatory features of structural genes are present. For instance, in the anthocyanin biosynthesis pathway (ABP), genes may be regulated by ABP-specific regulatory factors, and their expression levels may be strongly associated with anthocyanin pigmentation, or they may be expressed independently of pigmentation. A dataset of orthologous ABP genes (Chs, Chi, F3h, F3'h, Dfr, Ans) from monocot and dicot plant species that have distinct gene regulation patterns and different types of pollination was constructed to test whether these factors affect the evolution of the genes.

RESULTS

Using a maximum likelihood approach, we demonstrated that although the whole set of the ABP genes is under purifying selection, with greater selection acting on the upstream genes than on the downstream genes, genes from distinct groups of plant species experienced different strengths of selective pressure. The selective pressure on the genes was higher in dicots than in monocots (F3h and further downstream genes) and in pollinator-dependent plants than in pollinator-independent species (Chi and further downstream genes), suggesting an important role of pollination type in the evolution of the anthocyanin biosynthesis gene network. Contrasting effects of the regulation patterns on evolution were detected for the F3h and Dfr genes, with greater selective pressure on the F3h gene in plant species where the gene expression was not strongly associated with pigmentation and greater selective pressure on Dfr in plant species where the gene expression was associated with pigmentation.

CONCLUSIONS

We demonstrated the effects of pollination type and patterns of regulation on the evolution of the ABP genes, but the evolution of some of the genes could not be explained in the framework of these factors, such as the weaker selective pressure acting on Chs in species that attract pollinators or the stronger selective pressure on F3h in plant species where the gene expression was not associated with pigmentation. The observations suggest that additional factors could affect the evolution of these genes. One such factor could be an effect of gene duplication with further division of functions among gene copies and relaxed selective pressure acting on them. Additional tests with an appropriate dataset combining data on duplicated gene sequences and their functions in the flavonoid biosynthesis pathway are required to test this hypothesis.

A neutrality test for detecting selection on DNA methylation using single methylation polymorphism frequency spectrum

Inheritable epigenetic mutations (epimutations) can contribute to transmittable phenotypic variation. Thus, epimutations can be subject to natural selection and impact the fitness and evolution of organisms. Based on the framework of the modified Tajima’s D test for DNA mutations, we developed a neutrality test with the statistic “D^m” to detect selection forces on DNA methylation mutations using single methylation polymorphisms. With computer simulation and empirical data analysis, we compared the D^m test with the original and modified Tajima’s D tests and demonstrated that the D^m test is suitable for detecting selection on epimutations and outperforms original/modified Tajima’s D tests. Due to the higher resetting rate of epimutations, the interpretation of D^m on epimutations and Tajima’s D test on DNA mutations could be different in inferring natural selection. Analyses using simulated and empirical genome-wide polymorphism data suggested that genes under genetic and epigenetic selections behaved differently. We applied the D^m test to recently originated Arabidopsis and human genes, and showed that newly evolved genes contain higher level of rare epialleles, suggesting that epimutation may play a role in origination and evolution of genes and genomes. Overall, we demonstrate the utility of the D^m test to detect whether the loci are under selection regarding DNA methylation. Our analytical metrics and methodology could contribute to our understanding of evolutionary processes of genes and genomes in the field of epigenetics. The Perl script for the “D^m” test is available at http://fanlab.wayne.edu/ (last accessed December 18, 2014).

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.

Transcriptome analysis of a petal anthocyanin polymorphism in the arctic mustard, Parrya nudicaulis

Angiosperms are renown for their diversity of flower colors. Often considered adaptations to pollinators, the most common underlying pigments, anthocyanins, are also involved in plants’ stress response. Although the anthocyanin biosynthetic pathway is well characterized across many angiosperms and is composed of a few candidate genes, the consequences of blocking this pathway and producing white flowers has not been investigated at the transcriptome scale. We take a transcriptome-wide approach to compare expression differences between purple and white petal buds in the arctic mustard, Parrya nudicaulis, to determine which genes’ expression are consistently correlated with flower color. Using mRNA-Seq and de novo transcriptome assembly, we assembled an average of 722 bp per gene (49.81% coding sequence based on the A. thaliana homolog) for 12,795 genes from the petal buds of a pair of purple and white samples. Our results correlate strongly with qRT-PCR analysis of nine candidate genes in the anthocyanin biosynthetic pathway where chalcone synthase has the greatest difference in expression between color morphs (P/W = ∼7×). Among the most consistently differentially expressed genes between purple and white samples, we found 3× more genes with higher expression in white petals than in purple petals. These include four unknown genes, two drought-response genes (CDSP32, ERD5), a cold-response gene (GR-RBP2), and a pathogen defense gene (DND1). Gene ontology analysis of the top 2% of genes with greater expression in white relative to purple petals revealed enrichment in genes associated with stress responses including cold, drought and pathogen defense. Unlike the uniform downregulation of chalcone synthase that may be directly involved in the loss of petal anthocyanins, the variable expression of several genes with greater expression in white petals suggest that the physiological and ecological consequences of having white petals may be microenvironment-dependent.

Evolutionary study of the isoflavonoid pathway based on multiple copies analysis in soybean

BACKGROUND

Previous studies suggest that the metabolic pathway structure influences the selection and evolution rates of involved genes. However, most of these studies have exclusively considered a single gene copy encoding each enzyme in the metabolic pathway. Considering multiple-copy encoding enzymes could provide direct evidence of gene evolution and duplication patterns in metabolic pathways. We conducted a detailed analysis of the phylogeny, synteny, evolutionary rate and selection pressure of the genes in the isoflavonoid metabolic pathway of soybeans.

RESULTS

The results revealed that 1) only the phenylalanine ammonia-lyase (PAL) gene family most upstream from the pathway preserved all of the ancient and recent segmental duplications and maintained a strongly conserved synteny among these duplicated segments; gene families encoding branch-point enzymes with higher pleiotropy tended to retain more types of duplication; and genes encoding chalcone reductase (CHR) and isoflavone synthase (IFS) specific for legumes retained only recent segmental duplications; 2) downstream genes evolved faster than upstream genes and were subject to positive selection or relaxed selection constraints; 3) gene members encoding enzymes with high pleiotropy at the branching points were more likely to have undergone evolutionary differentiation, which may correspond to their functional divergences.

CONCLUSIONS

We reconciled our results with existing controversies and proposed that gene copies at branch points with higher connectivity might be under stronger selective constraints and that the gene copies controlling metabolic flux allocation underwent positive selection. Our analyses demonstrated that the structure and function of a metabolic pathway shapes gene duplication and the evolutionary constraints of constituent enzymes.

Divergent evolutionary and expression patterns between lineage specific new duplicate genes and their parental paralogs in Arabidopsis thaliana

Gene duplication is an important mechanism for the origination of functional novelties in organisms. We performed a comparative genome analysis to systematically estimate recent lineage specific gene duplication events in Arabidopsis thaliana and further investigate whether and how these new duplicate genes (NDGs) play a functional role in the evolution and adaption of A. thaliana. We accomplished this using syntenic relationship among four closely related species, A. thaliana, A. lyrata, Capsella rubella and Brassica rapa. We identified 100 NDGs, showing clear origination patterns, whose parental genes are located in syntenic regions and/or have clear orthologs in at least one of three outgroup species. All 100 NDGs were transcribed and under functional constraints, while 24% of the NDGs have differential expression patterns compared to their parental genes. We explored the underlying evolutionary forces of these paralogous pairs through conducting neutrality tests with sequence divergence and polymorphism data. Evolution of about 15% of NDGs appeared to be driven by natural selection. Moreover, we found that 3 NDGs not only altered their expression patterns when compared with parental genes, but also evolved under positive selection. We investigated the underlying mechanisms driving the differential expression of NDGs and their parents, and found a number of NDGs had different cis-elements and methylation patterns from their parental genes. Overall, we demonstrated that NDGs acquired divergent cis-elements and methylation patterns and may experience sub-functionalization or neo-functionalization influencing the evolution and adaption of A. thaliana.

Intraspecific sequence variation and differential expression in starch synthase genes of Arabidopsis thaliana

BACKGROUND

Natural accessions of Arabidopsis thaliana are a well-known system to measure levels of intraspecific genetic variation. Leaf starch content correlates negatively with biomass. Starch is synthesized by the coordinated action of many (iso)enzymes. Quantitatively dominant is the repetitive transfer of glucosyl residues to the non-reducing ends of α-glucans as mediated by starch synthases. In the genome of A. thaliana, there are five classes of starch synthases, designated as soluble starch synthases (SSI, SSII, SSIII, and SSIV) and granule-bound synthase (GBSS). Each class is represented by a single gene. The five genes are homologous in functional domains due to their common origin, but have evolved individual features as well. Here, we analyze the extent of genetic variation in these fundamental protein classes as well as possible functional implications on transcript and protein levels.

FINDINGS

Intraspecific sequence variation of the five starch synthases was determined by sequencing the entire loci including promoter regions from 30 worldwide distributed accessions of A. thaliana. In all genes, a considerable number of nucleotide polymorphisms was observed, both in non-coding and coding regions, and several amino acid substitutions were identified in functional domains. Furthermore, promoters possess numerous polymorphisms in potentially regulatory cis-acting regions. By realtime experiments performed with selected accessions, we demonstrate that DNA sequence divergence correlates with significant differences in transcript levels.

CONCLUSIONS

Except for AtSSII, all starch synthase classes clustered into two or three groups of haplotypes, respectively. Significant difference in transcript levels among haplotype clusters in AtSSIV provides evidence for cis-regulation. By contrast, no such correlation was found for AtSSI, AtSSII, AtSSIII, and AtGBSS, suggesting trans-regulation. The expression data presented here point to a regulation by common trans-regulatory transcription factors which ensures a coordinated action of the products of these four genes during starch granule biosynthesis. The apparent cis-regulation of AtSSIV might be related to its role in the initiation of de novo biosynthesis of granules.

Dissecting structural and nucleotide genome-wide variation in inbred Iberian pigs

BACKGROUND

In contrast to international pig breeds, the Iberian breed has not been admixed with Asian germplasm. This makes it an important model to study both domestication and relevance of Asian genes in the pig. Besides, Iberian pigs exhibit high meat quality as well as appetite and propensity to obesity. Here we provide a genome wide analysis of nucleotide and structural diversity in a reduced representation library from a pool (n=9 sows) and shotgun genomic sequence from a single sow of the highly inbred Guadyerbas strain. In the pool, we applied newly developed tools to account for the peculiarities of these data.

RESULTS

A total of 254,106 SNPs in the pool (79.6 Mb covered) and 643,783 in the Guadyerbas sow (1.47 Gb covered) were called. The nucleotide diversity (1.31x10-3 per bp in autosomes) is very similar to that reported in wild boar. A much lower than expected diversity in the X chromosome was confirmed (1.79x10-4 per bp in the individual and 5.83x10-4 per bp in the pool). A strong (0.70) correlation between recombination and variability was observed, but not with gene density or GC content. Multicopy regions affected about 4% of annotated pig genes in their entirety, and 2% of the genes partially. Genes within the lowest variability windows comprised interferon genes and, in chromosome X, genes involved in behavior like HTR2C or MCEP2. A modified Hudson-Kreitman-Aguadé test for pools also indicated an accelerated evolution in genes involved in behavior, as well as in spermatogenesis and in lipid metabolism.

CONCLUSIONS

This work illustrates the strength of current sequencing technologies to picture a comprehensive landscape of variability in livestock species, and to pinpoint regions containing genes potentially under selection. Among those genes, we report genes involved in behavior, including feeding behavior, and lipid metabolism. The pig X chromosome is an outlier in terms of nucleotide diversity, which suggests selective constraints. Our data further confirm the importance of structural variation in the species, including Iberian pigs, and allowed us to identify new paralogs for known gene families.

Polymorphism at genes involved in salt tolerance in Arabidopsis thaliana (Brassicaceae)

PREMISE OF THE STUDY

Genes involved in relevant functions for environmental adaptation can be considered primary candidates for their variation having been shaped by natural selection. Detecting recent selective events through their footprint on nucleotide variation constitutes a challenging task in species with a complex demographic history such as Arabidopsis thaliana. We have surveyed nucleotide variation in this species at nine genes involved in salt tolerance. The available genomewide information for this species has allowed us to contrast the levels and patterns of variation detected at the candidate genes with empirical distributions obtained from noncandidate regions.

METHODS

We sequenced nine genes involved in salt tolerance (~32 kb) in 20 ecotypes of A. thaliana and analyzed polymorphism and divergence at the individual gene and multilocus levels.

KEY RESULTS

Variation at the nine genes studied was characterized by a generalized skew toward polymorphisms with low-frequency variants. Except for genes RCD1 and NHX8, this pattern was similar to that generally detected in the A. thaliana genome and could thus be primarily explained by the species demographic history. The more extreme deviation at the NHX8 gene and its excess of polymorphism relative to divergence points to the recent action of selection on this gene.

CONCLUSIONS

The analysis of nucleotide polymorphism and divergence at nine genes involved in salt tolerance provided little evidence for the recent action of positive selection. Only the signals detected at NHX8 from both polymorphism and divergence were suggestive of the putative contribution of this gene to local adaptation.

Characterization, fine mapping and expression profiling of Ragged leaves1 in maize

The Ragged leaves1 (Rg1) maize mutant frequently develops lesions on leaves, leaf sheaths, and ear bracts. Lesion formation is independent of biotic stress. High-level accumulation of H(2)O(2) revealed by staining Rg1 leaves, with 3',3'-diaminobenzidine and trypan blue, suggested that lesion formation appeared to be due to cell death. Rg1 was initially mapped to an interval around 70.5 Mb in bin 3.04 on the short arm of chromosome 3. Utilizing 15 newly developed markers, Rg1 was delimitated to an interval around 17 kb using 16,356 individuals of a BC1 segregating population. There was only one gene, rp3, predicted in this region according to the B73 genome. Analysis of transcriptome data revealed that 441 genes significantly up-regulated in Rg1 leaves were functionally over-represented. Among those genes, several were involved in the production of reactive oxygen species (ROS). Our results suggested that lesions of Rg1 maize arose probably due to an aberrant rust resistance allele of Rp3, which elicited the accumulation of ROS independent of biotic stress.

Differential selection on carotenoid biosynthesis genes as a function of gene position in the metabolic pathway: a study on the carrot and dicots

Background

Selection of genes involved in metabolic pathways could target them differently depending on the position of genes in the pathway and on their role in controlling metabolic fluxes. This hypothesis was tested in the carotenoid biosynthesis pathway using population genetics and phylogenetics.

Methodology/Principal Findings

Evolutionary rates of seven genes distributed along the carotenoid biosynthesis pathway, IPI, PDS, CRTISO, LCYB, LCYE, CHXE and ZEP, were compared in seven dicot taxa. A survey of deviations from neutrality expectations at these genes was also undertaken in cultivated carrot (Daucus carota subsp. sativus), a species that has been intensely bred for carotenoid pattern diversification in its root during its cultivation history. Parts of sequences of these genes were obtained from 46 individuals representing a wide diversity of cultivated carrots. Downstream genes exhibited higher deviations from neutral expectations than upstream genes. Comparisons of synonymous and nonsynonymous substitution rates between genes among dicots revealed greater constraints on upstream genes than on downstream genes. An excess of intermediate frequency polymorphisms, high nucleotide diversity and/or high differentiation of CRTISO, LCYB1 and LCYE in cultivated carrot suggest that balancing selection may have targeted genes acting centrally in the pathway.

Conclusions/Significance

Our results are consistent with relaxed constraints on downstream genes and selection targeting the central enzymes of the carotenoid biosynthesis pathway during carrot breeding history.

A role in immunity for Arabidopsis cysteine protease RD21, the ortholog of the tomato immune protease C14

Secreted papain-like Cys proteases are important players in plant immunity. We previously reported that the C14 protease of tomato is targeted by cystatin-like EPIC proteins that are secreted by the oomycete pathogen Phytophthora infestans (Pinf) during infection. C14 has been under diversifying selection in wild potato species coevolving with Pinf and reduced C14 levels result in enhanced susceptibility for Pinf. Here, we investigated the role C14-EPIC-like interactions in the natural pathosystem of Arabidopsis with the oomycete pathogen Hyaloperonospora arabidopsidis (Hpa). In contrast to the Pinf-solanaceae pathosystem, the C14 orthologous protease of Arabidopsis, RD21, does not evolve under diversifying selection in Arabidopsis, and rd21 null mutants do not show phenotypes upon compatible and incompatible Hpa interactions, despite the evident lack of a major leaf protease. Hpa isolates express highly conserved EPIC-like proteins during infections, but it is unknown if these HpaEPICs can inhibit RD21 and one of these HpaEPICs even lacks the canonical cystatin motifs. The rd21 mutants are unaffected in compatible and incompatible interactions with Pseudomonas syringae pv. tomato, but are significantly more susceptible for the necrotrophic fungal pathogen Botrytis cinerea, demonstrating that RD21 provides immunity to a necrotrophic pathogen.

Genetic adaptation of the antibacterial human innate immunity network

BACKGROUND

Pathogens have represented an important selective force during the adaptation of modern human populations to changing social and other environmental conditions. The evolution of the immune system has therefore been influenced by these pressures. Genomic scans have revealed that immune system is one of the functions enriched with genes under adaptive selection.

RESULTS

Here, we describe how the innate immune system has responded to these challenges, through the analysis of resequencing data for 132 innate immunity genes in two human populations. Results are interpreted in the context of the functional and interaction networks defined by these genes. Nucleotide diversity is lower in the adaptors and modulators functional classes, and is negatively correlated with the centrality of the proteins within the interaction network. We also produced a list of candidate genes under positive or balancing selection in each population detected by neutrality tests and showed that some functional classes are preferential targets for selection.

CONCLUSIONS

We found evidence that the role of each gene in the network conditions the capacity to evolve or their evolvability: genes at the core of the network are more constrained, while adaptation mostly occurred at particular positions at the network edges. Interestingly, the functional classes containing most of the genes with signatures of balancing selection are involved in autoinflammatory and autoimmune diseases, suggesting a counterbalance between the beneficial and deleterious effects of the immune response.

A gene duplication/loss event in the ribulose-1,5-bisphosphate-carboxylase/oxygenase (rubisco) small subunit gene family among accessions of Arabidopsis thaliana

Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase; EC 4.1.1.39), the most abundant protein in nature, catalyzes the assimilation of CO(2) (worldwide about 10(11) t each year) by carboxylation of ribulose-1,5-bisphosphate. It is a hexadecamer consisting of eight large and eight small subunits. Although the Rubisco large subunit (rbcL) is encoded by a single gene on the multicopy chloroplast genome, the Rubisco small subunits (rbcS) are encoded by a family of nuclear genes. In Arabidopsis thaliana, the rbcS gene family comprises four members, that is, rbcS-1a, rbcS-1b, rbcS-2b, and rbcS-3b. We sequenced all Rubisco genes in 26 worldwide distributed A. thaliana accessions. In three of these accessions, we detected a gene duplication/loss event, where rbcS-1b was lost and substituted by a duplicate of rbcS-2b (called rbcS-2b*). By screening 74 additional accessions using a specific polymerase chain reaction assay, we detected five additional accessions with this duplication/loss event. In summary, we found the gene duplication/loss in 8 of 100 A. thaliana accessions, namely, Bch, Bu, Bur, Cvi, Fei, Lm, Sha, and Sorbo. We sequenced an about 1-kb promoter region for all Rubisco genes as well. This analysis revealed that the gene duplication/loss event was associated with promoter alterations (two insertions of 450 and 850 bp, one deletion of 730 bp) in rbcS-2b and a promoter deletion (2.3 kb) in rbcS-2b* in all eight affected accessions. The substitution of rbcS-1b by a duplicate of rbcS-2b (i.e., rbcS-2b*) might be caused by gene conversion. All four Rubisco genes evolve under purifying selection, as expected for central genes of the highly conserved photosystem of green plants. We inferred a single positive selected site, a tyrosine to aspartic acid substitution at position 72 in rbcS-1b. Exactly the same substitution compromises carboxylase activity in the cyanobacterium Anacystis nidulans. In A. thaliana, this substitution is associated with an inferred recombination. Functional implications of the substitution remain to be evaluated.

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.

Molecular evolution of the endosperm starch synthesis pathway genes in rice (Oryza sativa L.) and its wild ancestor, O. rufipogon L

The evolution of metabolic pathways is a fundamental but poorly understood aspect of evolutionary change. One approach for understanding the complexity of pathway evolution is to examine the molecular evolution of genes that together comprise an integrated metabolic pathway. The rice endosperm starch biosynthetic pathway is one of the most thoroughly characterized metabolic pathways in plants, and starch is a trait that has evolved in response to strong selection during rice domestication. In this study, we have examined six key genes (AGPL2, AGPS2b, SSIIa, SBEIIb, GBSSI, ISA1) in the rice endosperm starch biosynthesis pathway to investigate the evolution of these genes before and after rice domestication. Genome-wide sequence tagged sites data were used as a neutral reference to overcome the problems of detecting selection in species with complex demographic histories such as rice. Five variety groups of Oryza sativa (aus, indica, tropical japonica, temperate japonica, aromatic) and its wild ancestor (O. rufipogon) were sampled. Our results showed evidence of purifying selection at AGPL2 in O. rufipogon and strong evidence of positive selection at GBSSI in temperate japonica and tropical japonica varieties and at GBSSI and SBEIIb in aromatic varieties. All the other genes showed a pattern consistent with neutral evolution in both cultivated rice and its wild ancestor. These results indicate the important role of positive selection in the evolution of starch genes during rice domestication. We discuss the role of SBEIIb and GBSSI in the evolution of starch quality during rice domestication and the power and limitation of detecting selection using genome-wide data as a neutral reference.

Population history in Arabidopsis halleri using multilocus analysis

A. halleri is a psuedometallophyte with a patchy distribution in Europe and is often spread by human activity. To determine the population history and whether this history is consistent with potential human effects, we surveyed nucleotide variation using 24 loci from 12 individuals in a large A. halleri population. The means of total and silent nucleotide variation (theta(W)) are within the range expected for the species. The population genetic neutrality tests Tajima's D and Wall's B had significant composite results rejecting panmixia, and Approximate Bayesian Computation analysis revealed that a subdivision model better explained the variation than the standard neutral model, refugia (or admixture), bottleneck or change of population size models. A categorical regression analysis further supports the subdivision model, and under the subdivision model, the neutrality tests are no longer significant. The best support was for two source populations, a situation consistent with the mixing of two populations possibly mediated by human activity. This scenario might limit the genetic diversity and adaptive potential of the population. The non-neutral population variation described here should be considered in bioinformatic searches for adaptation.

Evolutionary study of a potential selection target region in the pig

Domestication, modern breeding and artificial selection have shaped dramatically the genomic variability of domestic animals. In livestock, the so-called FAT1 quantitative trait locus (QTL) in porcine chromosome 4 was the first QTL uncovered although, to date, its precise molecular nature has remained elusive. Here, we characterize the nucleotide variability of 13 fragments of ∼500 bp equally spaced in a 2 Mb region in the vicinity of the FAT1 region in a wide-diversity panel of 32 pigs. Asian and European animals, including local Mediterranean and international pig breeds, were sequenced. Patterns of genetic variability were very complex and varied largely across loci and populations; they did not reveal overall a clear signal of a selective sweep in any breed, although FABP4 fragment showed a significantly higher diversity. We used an approximate Bayesian computation approach to infer the evolutionary history of this SSC4 region. Notably, we found that European pig populations have a much lower effective size than their Asian counterparts: in the order of hundreds vs hundreds of thousands. We show also an important part of extant European variability is actually due to introgression of Asian germplasm into Europe. This study shows how a potential loss in diversity caused by bottlenecks and possible selective sweeps associated with domestication and artificial selection can be counterbalanced by migration, making it much more difficult the identification of selection footprints based on naive demographic assumptions. Given the small fragment analyzed here, it remains to be studied how these conclusions apply to the rest of the genome.

Comparative SNP diversity among four Eucalyptus species for genes from secondary metabolite biosynthetic pathways

Background

There is little information about the DNA sequence variation within and between closely related plant species. The combination of re-sequencing technologies, large-scale DNA pools and availability of reference gene sequences allowed the extensive characterisation of single nucleotide polymorphisms (SNPs) in genes of four biosynthetic pathways leading to the formation of ecologically relevant secondary metabolites in Eucalyptus. With this approach the occurrence and patterns of SNP variation for a set of genes can be compared across different species from the same genus.

Results

In a single GS-FLX run, we sequenced over 103 Mbp and assembled them to approximately 50 kbp of reference sequences. An average sequencing depth of 315 reads per nucleotide site was achieved for all four eucalypt species, Eucalyptus globulus, E. nitens, E. camaldulensis and E. loxophleba. We sequenced 23 genes from 1,764 individuals and discovered 8,631 SNPs across the species, with about 1.5 times as many SNPs per kbp in the introns compared to exons. The exons of the two closely related species (E. globulus and E. nitens) had similar numbers of SNPs at synonymous and non-synonymous sites. These species also had similar levels of SNP diversity, whereas E. camaldulensis and E. loxophleba had much higher SNP diversity. Neither the pathway nor the position in the pathway influenced gene diversity. The four species share between 20 and 43% of the SNPs in these genes.

Conclusion

By using conservative statistical detection methods, we were confident about the validity of each SNP. With numerous individuals sampled over the geographical range of each species, we discovered one SNP in every 33 bp for E. nitens and one in every 31 bp in E. globulus. In contrast, the more distantly related species contained more SNPs: one in every 16 bp for E. camaldulensis and one in 17 bp for E. loxophleba, which is, to the best of our knowledge, the highest frequency of SNPs described in woody plant species.

Application of the isolation with migration model demonstrates the pleistocene origin of geographic differentiation in Cardamine nipponica (Brassicaceae), an endemic Japanese alpine plant

The Pleistocene was characterized by a cyclic pattern of cold and warm climatic periods, or climatic oscillations, which caused fluctuations in the distributions of organisms. This resulted in drastic changes in demography, thereby accelerating the genetic divergence of populations. Phylogeographic studies have elucidated the history of populations during the Pleistocene. However, given the lack of model-based analysis of population histories, previous phylogeographic studies could not adequately evaluate the effect of these Pleistocene climatic oscillations on the genetic divergence and migration events between populations. Populations of Japanese alpine plants in central and northern Japan are highly differentiated, and a history of isolation between regions during the Pleistocene was inferred. Using sequences of 10 nuclear genes (ca. approximately 7,000 bp in total) from Cardamine nipponica (Brassicaceae), we applied an isolation with migration (IM) model to test the significance of the isolation history between central and northern Japan and to assess whether range shifts during the Pleistocene climatic oscillations were involved in the genetic differentiation between regions. The estimated divergence time indicates that the two regions were separated about 100,000-110,000 years ago. The exclusive occurrence of closely related haplotypes within each region (parsimony network) and the high level of genetic differentiation between the regions (mean F(ST) = 0.417) indicate that genetic divergence occurred following the isolation of the two regions. Therefore, the genetic differentiation between regions was shaped during the Pleistocene, especially during the last glacial and inter and postglacial periods. In addition, our multilocus analysis showed that populations in central and northern Japan were completely isolated after they split. Geographic separation and subsequent restricted migration events among mountains could explain this isolation history between regions. Furthermore, genetic drift in the reduced populations would remove evidence of occasional migration, emphasizing the isolation history. Therefore, our application of a demographic model demonstrated the Pleistocene origin of geographic differentiation statistically and provided a plausible migration history for C. nipponica.

The correlation of evolutionary rate with pathway position in plant terpenoid biosynthesis

Genes are expected to face stronger selective constraint and to evolve more slowly if they encode enzymes upstream as opposed to downstream in metabolic pathways, because upstream genes are more pleiotropic, being required for a wider range of end products. However, few clear examples of this trend in evolutionary rate variation exist. We examined whether genes involved in plant terpenoid biosynthesis exhibit such a pattern, using data for 40 genes from five fully sequenced angiosperms, Oryza, Vitis, Arabidopsis, Populus, and Ricinus. Our results show that d(N)/d(S) does in fact correlate with pathway position along pathways converting glucose to the terpenoid phytohormones abscissic acid, gibberellic acid (GA), and brassinosteroids. Upstream versus downstream rate variation is particularly strong in the GA pathway. In contrast, we found no or little apparent variation in d(N)/d(S) with gene copy number. We also introduce a new measure of pathway position, the Pathway Pleiotropy Index (PPI), which counts groups of enzymes between pathway branch points. We found that this measure is superior to pathway position in explaining variation in d(N)/d(S) along each pathway. Although at least 8 of the 40 genes showed evidence of positive selection, correlations of d(N)/d(S) with PPI remain significant when these genes are removed. Therefore, our results are consistent with the prediction that selective constraint is progressively relaxed along metabolic pathways.

Relaxed selection on the CBF/DREB1 regulatory genes and reduced freezing tolerance in the southern range of Arabidopsis thaliana

Elucidating the molecular basis of adaptive phenotypic variation represents a central aim in evolutionary biology. Traits exhibiting patterns of clinal variation represent excellent models for studies of molecular adaptation, especially when variation in phenotype can be linked to organismal fitness in different environments. Natural accessions of the model plant species Arabidopsis thaliana exhibit clinal variation in freezing tolerance that follows a gradient of temperature variability across the species' native range (Zhen Y, Ungerer MC. 2008. Clinal variation in freezing tolerance among natural accessions of A. thaliana. New Phytol. 177:419-427). Here, we report that this pattern of variation is attributable, at least in part, to relaxed purifying selection on members of a small family of transcriptional activators (the CBF/DREB1s) in the species' southern range. These regulatory genes play a critical role in the ability of A. thaliana plants to undergo cold acclimation and thereby achieve maximum freezing tolerance. Relative to accessions from northern regions, accessions of A. thaliana from the southern part of their geographic range exhibit levels of nonsynonymous nucleotide polymorphism that are approximately 2.8-fold higher across this small gene subfamily. Relaxed selection on the CBF/DREB1s in southern accessions also has resulted in multiple mutations in regulatory regions resulting in abrogated expression of particular subfamily members in particular accessions. These coding-region and regulatory mutations compromise the ability of these genes to act as efficient transcriptional activators during the cold acclimation process, as determined by reductions in rates of induction and maximum levels of expression in the downstream genes they regulate. This study highlights the potential role of regulatory genes in underlying adaptive phenotypic variation in nature.