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

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata . Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D , Fu's F_S and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's θ) and pattern of variation (Tajima's D ) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.     

Nucleotide variation at the CHALCONE ISOMERASE locus in Arabidopsis thaliana

An approximately 1.9-kb region encompassing the CHI gene, which encodes chalcone isomerase, was sequenced in 24 worldwide ecotypes of Arabidopsis thaliana (L.) Heynh. and in 1 ecotype of A. lyrata ssp. petraea. There was no evidence for dimorphism at the CHI region. A minimum of three recombination events was inferred in the history of the sampled ecotypes of the highly selfing A. thaliana. The estimated nucleotide diversity theta(TOTAL) = 0.004, theta(SIL) = 0. 005 was on the lower part of the range of the corresponding estimates for other gene regions. The skewness of the frequency spectrum toward an excess of low-frequency polymorphisms, together with the bell-shaped distribution of pairwise nucleotide differences at CHI, suggests that A. thaliana has recently experienced a rapid population growth. Although this pattern could also be explained by a recent selective sweep at the studied region, results from the other studied loci and from an AFLP survey seem to support the expansion hypothesis. Comparison of silent polymorphism and divergence at the CHI region and at the Adh1 and ChiA revealed in some cases a significant deviation of the direct relationship predicted by the neutral theory, which would be compatible with balancing selection acting at the latter regions.     

Nucleotide variation at the myrosinase-encoding locus, TGG1, and quantitative myrosinase enzyme activity variation in Arabidopsis thaliana

The Arabidopsis thaliana TGG1 gene encodes thioglucoside glucohydrolase (myrosinase), an enzyme catalysing the hydrolysis of glucosinolate compounds. The enzyme is involved in plant defence against some insect herbivores, and is present in species of the order Capparales (Brassicales). Nucleotide variation was surveyed by sequencing c. 2.4 kb of the TGG1 locus in a sample of 28 worldwide A. thaliana accessions, and one Arabidopsis lyrata ssp. lyrata individual. Myrosinase activity was quantified for 27 of these same A. thaliana accessions, plus five additional others. Overall, estimated nucleotide diversity in A. thaliana was low compared to other published A. thaliana surveys, and the frequency distribution was skewed toward an excess of low-frequency variants. Furthermore, comparison to the outgroup species A. lyrata demonstrated that A. thaliana exhibited an excess of high-frequency derived variants relative to a neutral equilibrium model, suggesting a selective sweep. A. thaliana accessions differed significantly in total myrosinase activity, but analysis of variance detected no statistical evidence for an association between quantitative enzyme activity and alleles at the TGG1 myrosinase-encoding locus. We thus conclude that other, unsurveyed factors primarily affect the observed myrosinase activity levels in this species. The pattern of nucleotide variation was consistent with a model of positive selection but might also be compatible with a completely neutral model that takes into account the metapopulation behaviour of this highly inbreeding species which experienced a relatively recent worldwide expansion.     

CYP98A3 from Arabidopsis thaliana is a 3'-hydroxylase of phenolic esters, a missing link in the phenylpropanoid pathway.

The 4- and 5-hydroxylations of phenolic compounds in plants are catalyzed by cytochrome P450 enzymes. The 3-hydroxylation step leading to the formation of caffeic acid from p-coumaric acid remained elusive, however, alternatively described as a phenol oxidase, a dioxygenase, or a P450 enzyme, with no decisive evidence for the involvement of any in the reaction in planta. In this study, we show that the gene encoding CYP98A3, which was the best possible P450 candidate for a 3-hydroxylase in the Arabidopsis genome, is highly expressed in inflorescence stems and wounded tissues. Recombinant CYP98A3 expressed in yeast did not metabolize free p-coumaric acid or its glucose or CoA esters, p-coumaraldehyde, or p-coumaryl alcohol, but very actively converted the 5-O-shikimate and 5-O-d-quinate esters of trans-p-coumaric acid into the corresponding caffeic acid conjugates. The shikimate ester was converted four times faster than the quinate derivative. Antibodies directed against recombinant CYP98A3 specifically revealed differentiating vascular tissues in stem and root. Taken together, these data show that CYP98A3 catalyzes the synthesis of chlorogenic acid and very likely also the 3-hydroxylation of lignin monomers. This hydroxylation occurs on depsides, the function of which was so far not understood, revealing an additional and unexpected level of networking in lignin biosynthesis.     

Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.     

Nucleotide sequences of two corn histone H3 genes. Genomic organization of the corn histone H3 and H4 genes

Summary Two histone H3 genes have been cloned from a gtWES.B corn genomic library. The nucleotide sequences show 96% homology and both encode the same protein, which differs from its counterpart in wheat and pea by one amino acid substitution. The 5-flanking regions of the two corn H3 genes contain the classical histone-gene-specific consensus sequences and possess several regions of extensive nucleotide homology. A conserved octanucleotide 5-CGCGGATC-3 occurs at approximately 200 nucleotides upstream from the initiation ATG codon. This octanucleotide was found to exist in all of the 7 plant histone genes sequenced so far. Codon usage is characterized by a very high frequency of C (67%) and G (28%) at the third position of the codons, those ending by A (1%) and T (4%) being practically excluded.Comparison of Southern blots of EcoRI, EcoRV and BamHI digested genomic DNA suggests that the corn H3 and H4 genes are not closely associated. The H3 genes exist as 60 to 80 copies and the H4 genes as 100 to 120 copies per diploid genome. re]19851002 rv]19851212 ac]19851216