Sequence divergence at the putative flowering time locus COL1 in Brassicaceae

An insertion/deletion polymorphism (Ind2) in the Brassica nigra CONSTANS LIKE 1 (Bni COL1) gene was previously found to be associated with variation in flowering time. In the present study we examine the inter-specific divergence of COL1 in the family Brassicaceae. Analysis of codon substitution models did not reveal evidence of positive Darwinian selection, but comparisons of the COL1 gene in different species revealed a surprising number of indels. A total of 24 indels were found in the 650 bp of the middle variable region of the gene. This high number of indels could reflect a lack of constraint on length of this region of the protein, or the effect of positive selection. The number of indels was close to that expected in non-coding DNA, but the indels were longer in COL1 than those observed in non-coding regions. Reconstruction of indel evolution indicated that most indels resulted from deletions rather than insertions. The Ind2 indel that has shown association with flowering time in Brassica nigra exhibited a remarkable distribution in the Brassicaceae family, indicating that the polymorphism may have persisted more than ten million years. Considering presumed historic populations sizes of Brassicaceae species, such a long persistence time seems unlikely for a neutral polymorphism.     

Naturally occurring indel variation in the Brassica nigra COL1 gene is associated with variation in flowering time

Previous QTL mapping identified a Brassica nigra homolog to Arabidopsis thaliana CO as a candidate gene affecting flowering time in B. nigra. Transformation of an A. thaliana co mutant with two different alleles of the B. nigra CO (Bni COa) homolog, one from an early-flowering B. nigra plant and one from a late one, did not show any differential effect of the two alleles on flowering time. The DNA sequence of the coding region of the two alleles was also identical, showing that nucleotide variation influencing flowering time must reside outside the coding region of Bni COa. In contrast, the nucleotide sequence of the B. nigra COL1 (Bni COL1) gene located 3.5 kb upstream of Bni COa was highly diverged between the alleles from early and late plants. One indel polymorphism in the Bni COL1 coding region, present in several natural populations of B. nigra, displayed a significant association with flowering time within a majority of these populations. These data indicate that a quantitative trait nucleotide (QTN) affecting flowering time is located within or close to the Bni COL1 gene. The intergenic sequence between Bni COL1 and Bni COa displayed a prominent peak of divergence 1 kb downstream of the Bni COL1 coding region. This region could contain regulatory elements for the downstream Bni COa gene. Our data suggest that a naturally occurring QTN for flowering time affects the function or expression of either Bni COL1 or Bni COa.     

Sequence variation and haplotype structure at the putative flowering-time locus COL1 of Brassica nigra

Motivated by a previous study indicating that polymorphism at an indel, Ind2, within the Brassica nigra COL1 gene is significantly associated with flowering time, we searched for evidence of selection in a sample of 41 complete sequences of B. nigra COL1. The within-gene population recombination rate is moderate, and all neutrality tests used in the present study failed to detect departure from the standard neutral model or evidence of selection. The haplotype structure of the 5'-half of the gene is primarily associated with the demographic history of the species and more specifically with the split between European and Ethiopian populations, whereas the structure of the 3'-half reflects the polymorphism at Ind2. This could be the result of selection or a combination of recombination and migration during the history of the sample of sequences. Without additional information on polymorphism in flanking areas, these two alternatives are difficult to tell apart. If selection acted on the gene, we suggest that if the indel itself is not the target of selection, among the polymorphic sites cosegregating with the polymorphism at Ind2, replacement polymorphisms around sites 890 and 1260 are the most likely quantitative trait nucleotides within the gene.     

The Mhc class II of the Black grouse (Tetrao tetrix) consists of low numbers of B and Y genes with variable diversity and expression

We found that the Black grouse (Tetrao tetrix) possess low numbers of Mhc class II B (BLB) and Y (YLB) genes with variable diversity and expression. We have therefore shown, for the first time, that another bird species (in this case, a wild lek-breeding galliform) shares several features of the simple Mhc of the domestic chicken (Gallus gallus). The Black grouse BLB genes showed the same level of polymorphism that has been reported in chicken, and we also found indications of balancing selection in the peptide-binding regions. The YLB genes were less variable than the BLB genes, also in accordance with earlier studies in chicken, although their functional significance still remains obscure. We hypothesize that the YLB genes could have been under purifying selection, just as the mammal Mhc-E gene cluster.     

Use of the SSLP-based method for detection of rare apomictic events in a sexual AtSERK1 transgenic Arabidopsis population

Here we present a screening method to evaluate the potential of genes to transfer aspects of apomixis into sexual crop plants. Based on the assumption that an apomictic progeny is an exact genetic replica of the mother plant we employed a set of single sequence length polymorphism (SSLP) markers to identify individuals displaying heterozygosity fixation in segregating sexual populations as an indication of rare apomictic events. Here we present the results of such a study using the Arabidopsis thaliana SOMATIC EMBRYOGENESIS RECEPTOR KINASE 1 (AtSERK1) gene expressed under the control of the AtLTP1 promoter in sexual Arabidopsis plants. In one of the three tested F2 transgenic populations expressing the AtLTP1::AtSERK1 construct we observed two plants with heterozygosity maintenance for the full set of SSLP markers indicating a possible clonal inheritance. However, as their offspring revealed a close to binomial segregation for a number of heterozygous loci, it was concluded that these two putative apomictic plants either lost their clonal ability in the next generation or resulted from incidental recombination events displaying the genotype of the parent. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Local Patterns of Nucleotide Polymorphism Are Highly Variable in the Selfing Species Arabidopsis thaliana

Neighboring genes predictably share similar evolutionary histories to an extent delineated by recombination. This correlation should extend across multiple linked genes in a selfing species such as Arabidopsis thaliana due to its low effective recombination rate. To test this prediction, we performed a molecular population genetics analysis of nucleotide polymorphism and divergence in chromosomal regions surrounding four low-diversity loci. Three of these loci, At1g67140, At3g03700, and TERMINAL FLOWER1 (TFL1), have been previously implicated as targets of selection and we would predict stronger correlations in polymorphism between neighboring loci due to genetic hitchhiking around these loci. The remaining locus, At1g04300, was identified in a study of linkage disequilibrium surrounding the CRYPTOCHROME2 (CRY2) locus. Although we found broad valleys of reduced nucleotide variation around two of our focal genes, At1g67140 and At3g03700, all chromosomal regions exhibited extreme variation in the patterns of polymorphism and evolution between neighboring loci. Although three of our four regions contained potential targets of selection, application of the composite-likelihood-ratio test of selection in conjunction with a goodness-of-fit test supports the selection hypothesis only for the region containing At3g03700. The degree of discordance in evolutionary histories between linked loci within each region generally correlated with estimates of recombination and linkage disequilibrium for that region, with the exception of the region containing At1g04300. We discuss the implications of these data for future population genetics analyses and genomics studies in A. thaliana. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Comparative mapping between Arabidopsis thaliana and Brassica nigra indicates that Brassica genomes have evolved through extensive genome replication accompanied by chromosome fusions and frequent rearrangements.

Chromosome organization and evolution in the Brassicaceae family was studied using comparative linkage mapping. A total of 160 mapped Arabidopsis thaliana DNA fragments identified 284 homologous loci covering 751 cM in Brassica nigra. The data support that modern diploid Brassica species are descended from a hexaploid ancestor, and that the A. thaliana genome is similar in structure and complexity to those of each of the hypothetical diploid progenitors of the proposed hexaploid. Thus, the Brassica lineage probably went through a triplication after the divergence of the lineages leading to A. thaliana and B. nigra. These duplications were also accompanied by an exceptionally high rate of chromosomal rearrangements. The average length of conserved segments between A. thaliana and B. nigra was estimated at 8 cM. This estimate corresponds to approximately 90 rearrangements since the divergence of the two species. The estimated rate of chromosomal rearrangements is higher than any previously reported data based on comparative mapping. Despite the large number of rearrangements, fine-scale comparative mapping between model plant A. thaliana and Brassica crops is likely to result in the identification of a large number of genes that affect important traits in Brassica crops.     

Phylogenomic Analysis of the PEBP Gene Family in Cereals

The TFL1 and FT genes, which are key genes in the control of flowering time in Arabidopsis thaliana, belong to a small multigene family characterized by a specific phosphatidylethanolamine-binding protein domain, termed the PEBP gene family. Several PEBP genes are found in dicots and monocots, and act on the control of flowering time. We investigated the evolution of the PEBP gene family in cereals. First, taking advantage of the complete rice genome sequence and EST databases, we found 19 PEBP genes in this species, 6 of which were not previously described. Ten genes correspond to five pairs of paralogs mapped on known duplicated regions of the rice genome. Phylogenetic analysis of Arabidopsis and rice genes indicates that the PEBP gene family consists of three main homology classes (the so-called TFL1-LIKE, MFT-LIKE, and FT-LIKE subfamilies), in which gene duplication and/or loss occurred independently in Arabidopsis and rice. Second, phylogenetic analyses of genomic and EST sequences from five cereal species indicate that the three subfamilies of PEBP genes have been conserved in cereals. The tree structure suggests that the ancestral grass genome had at least two MFT-like genes, two TFL1-like genes, and eight FT-like genes. A phylogenomic approach leads to some hypotheses about conservation of gene function within the subfamilies. [Reviewing Editor: Dr. Yves Van de Peer]     

Origin and maintenance of a broad-spectrum disease resistance locus in Arabidopsis

The broad-spectrum mildew resistance genes RPW8.1 and RPW8.2 define a unique type of plant disease resistance (R) gene, and so far homologous sequences have been found in Arabidopsis thaliana only, which suggests a recent origin. In addition to RPW8.1 and RPW8.2, the RPW8 locus contains three homologs of RPW8, HR1, HR2, and HR3, which do not contribute to powdery mildew resistance. To investigate whether RPW8 has originated recently, and if so the processes involved, we have isolated and analyzed the syntenic RPW8 loci from Arabidopsis lyrata, and from Brassica rapa and B. oleracea. The A. lyrata locus contains four genes orthologous to HR1, HR2, HR3, and RPW8.2, respectively. Two syntenic loci have been characterized in Brassica; one locus contains three genes and is present in both B. oleracea and B. rapa, and the other locus contains a single gene and is detected in B. rapa only. The Brassica homologs have highest similarity to HR3. Sequence analyses suggested that the RPW8 gene family in Brassicaceae originated from an HR3-like ancestor gene through a series of duplications and that RPW8.1 and RPW8.2 evolved from functional diversification through positive selection several MYA. Examination of the sequence polymorphism of 32 A. thaliana accessions at the RPW8 locus and their disease reaction phenotypes revealed that the polymorphic RPW8 locus defines a major source of resistance to powdery mildew diseases. A possible evolutionary mechanism by which functional polymorphism at the AtRPW8 locus has been maintained in contemporary populations of A. thaliana is discussed.     

Asr genes belong to a gene family comprising at least three closely linked loci on chromosome 4 in tomato

Asr1, Asr2 andAsr3 are three homologous clones isolated from tomato whose expression is believed to be regulated by abscisic acid (ABA); the corresponding genes thus participate in physiological and developmental processes such as responses of leaf and root to water stress, and fruit ripening. In this report, results obtained with Near Isogenic Lines reveal thatAsr1, Asr2 andAsr3 represent three different loci. In addition, we map these genes on the restriction fragment length polymorphism (RFLP) map of the tomato genome by using an F₂ population derived from an interspecific hybrid crossL. esculentum × L. penelli. RFLP data allow us to map these genes on chromosome 4, suggesting that they belong to a gene family. The elucidation of the genomic organization of theAsr gene family may help in understanding the role of its members in the response to osmotic stress, as well as in fruit ripening, at the molecular level.     

Rapid genotyping of the osteoporosis-associated polymorphic transcription factor Sp1 binding site in the COL1A1 gene by pyrosequencing

We describe a novel polymerase chain reaction (PCR) and deoxyribonucleic acid (DNA) sequencingbased assay for rapid genotyping of the polymorphic Sp1 binding site in the COL1A1 gene (1). A single nucleotide G-->T substitution polymorphism at this GC-rich site has recently been reported to be a predictive genetic marker for low bone mineral density (BMD). To simplify screening for this marker, we optimized PCR conditions and subjected the amplicons to pyrosequencing, which is a convenient high-throughput sequence analysis technique, readily amenable to automation. The analysis of 200 deidentified convenience DNA samples extracted from blood revealed genotype frequences in Hardy-Weinberg equilibrium (SS 68.0%, Ss 28.5%, and ss 3.5%) in agreement with other studies of European populations. This study demonstrates for the first time that pyrosequencing can be used for rapid identification of the osteoporosis-associated single nucleotide polymorphism (SNP) in the COL1A1 gene.     

Recent degeneration of an old duplicated flowering time gene in Brassica nigra

Gene and genome duplications play a major role in the evolution of plant species. The Brassica nigra genome is highly replicated as a result of ancient polyploidization events. Two copies of the flowering time gene CONSTANS (COa and COb) have been identified in B. nigra, and previous studies showed that COa is functional. In the present study, the polymorphism of 92 COb alleles sampled in seven populations was analyzed. Both polymorphism and recombination levels were elevated and varied strongly among populations and 8% of COb alleles exhibit apparently disabling mutations. Sequence data, however, do not provide unambiguous support for the presence of relaxed selective constraint on COb as compared to known functional CO genes. On the one hand, some of the disabling mutations reached high-frequency arguing for a loss of function but, on the other hand, the ratio of nonsynonymous to synonymous nucleotide polymorphism and diversity is low and similar to that observed in other B. nigra CO and CO-like genes, supporting the conservation of some function. We also showed that COb is still transcribed. Finally, the flowering time of Arabidopsis thaliana co mutant plants transformed with COb alleles with and without apparent disabling mutations was similar. We propose that COb was retained for a long period after duplication, but a recent fixation of a detrimental mutation, possibly as an effect of a bottleneck, resulted in its nonfunctionalization. We also speculate as to the presence of subsequent selection for rapid degeneration of the gene.     

Comparative mapping in Arabidopsis and Brassica, fine scale genome collinearity and congruence of genes controlling flowering time

The model dicotyledonous plant, Arabidopsis thaliana , is closely related to Brassica crop species. It is intended that information concerning the genetic control of basic biological processes in Arabidopsis will be transferable to other species. Genome collinearity and its potential to facilitate the identification of candidate genes in Arabidopsis homologous to genes controlling important agronomic traits in Brassica was investigated. Genetic mapping in B. nigra identified two loci influencing flowering time (FT), with loci on linkage groups 2 and 8 explaining 53% and 12% of the total variation in FT, respectively. The CO gene exerts an important control over FT in A. thaliana , and B. nigra homologues of CO probably also play an important role in regulating FT. B. nigra homologues of CO were identified on linkage groups 2 and 8, the homologue on group 2 was coincident with the major locus controlling FT while the homologue on group 8 was within the 90% confidence interval of the weaker FT gene. The CO homologue on group 2 exhibits abundant allelic variation suggesting that it naturally controls a wide range of flowering times. Fine-scale A. thaliana/B. nigra comparative mapping demonstrated short-range collinearity between the genomes of Arabidopsis and Brassica . Eleven DNA fragments spaced over a 1.5 Mb contig in A. thaliana were used as RFLP probes in B. nigra . Three collinear representations of the A. thaliana contig were identified in B. nigra , with one interrupted by a large chromosomal inversion. Collinearity over this range will allow the resources generated by the Arabidopsis genome project to facilitate map-based cloning in Brassica crops.