Comparative analysis of methylthioalkylmalate synthase (MAM) gene family and flanking DNA sequences in Brassica oleracea and Arabidopsis thaliana

Gene BoGSL-PRO is associated with presence of 3-carbon side-chain glucosinolates (GSL). This gene is a member of the methylthioalkylmalate synthase (MAM) gene family. A BAC clone of Brassica oleracea, B21F5, containing this gene, was sequenced, annotated and compared to its corresponding region in Arabidopsis thaliana. Twelve protein-coding genes and 10 transposable elements were found in this clone. The corresponding region in A. thaliana chromosome I has 14 genes and no transposable elements. Analysis of MAM gene family in both species, which also include genes controlling 4-carbon side-chain GSL, separated the genes in two groups based on exon numbers and function. Phylogenetic analysis of the amino acid sequences encoded by these genes suggest that these two groups were produced by a duplication that must have occurred before the divergence of the Rosid and Asterid lineages of angiosperms. Comparison with putative orthologs from several prokaryotes further suggest that the members of the gene family with 10 exons, which encode proteins involved in 4-carbon side-chain GSL biosynthesis, were derived via truncation of the 3′ end from ancestral genes more similar in length to those with 12 exons, which encode proteins involved in 3-carbon side-chain GSL biosynthesis. Lower gene density in B. oleracea compared to A. thaliana is due in part to presence of transposable elements (TE) mostly in inter-genic regions.     

Expression profiles and genomic organisation of group A protein phosphatase 2C genes in Brassica oleracea

We investigated the expression profiles and genomic organisation of the ABA‐responsive genes encoding protein phosphatases 2C (PP2C, group A members) in Brassica oleracea to better understand their functional and genetic relations. Gene expression profiling of drought responsive genes in B. oleracea and Arabidopsis thaliana revealed significant differences in the gene expression pattern of a key regulator of ABA signalling—ABI1 PP2C. This finding prompted us to study genetic relations within the PP2Cs group A in the Brassica species. Twenty homologous B. oleracea sequences were identified and characterised as putative PP2C group A members. Phylogenetic analysis revealed that the B. oleracea homologues were closely related to the particular members of the A. thaliana PP2C. The genetic analysis corroborated the presence of two to three gene copies in B. oleracea in comparison to the nine unique PP2C genes in the A. thaliana genome. Gene expression analyses showed significant differences in PP2C gene expression pattern in B. oleracea. Our results indicate that PP2C‐based drought stress signalling in B. oleracea has evolved distinctly. Different reactions of particular B. oleracea PP2C genes to drought stress and ABA treatment indicate low conservation of gene expression patterns and functional divergence between B. oleracea and A. thaliana homologous genes.     

From Arabidopsis thaliana to Brassica napus: development of amplified consensus genetic markers (ACGM) for construction of a gene map

The evolution of genomes can be studied by comparing maps of homologous genes which show changes in nucleic acid sequences and chromosome rearrangements. In this study, we developed a set of 32 amplified consensus gene markers (ACGMs) that amplified gene sequences from Arabidopsis thaliana and Brassica napus. Our methodology, based on PCR, facilitated the rapid sequencing of homologous genes from various species of the same phylogenetic family and the detection of intragenic polymorphism. We found that such polymorphism principally concerned intron sequences and we used it to attribute a Brassica oleracea or Brassica rapa origin to the B. napus sequences and to map 43 rapeseed genes. We confirm that the genetic position of homologous genes varied between B. napus and A. thaliana. ACGMs are a useful tool for genome evolution studies and for the further development of single nucleotide polymorphism suitable for use in genetic mapping and genetic diversity analyses.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Comparative sequence analysis for Brassica oleracea with similar sequences in B. rapa and Arabidopsis thaliana

We sequenced five BAC clones of Brassica oleracea doubled haploid ‘Early Big' broccoli containing major genes in the aliphatic glucosinolate pathway, and comparatively analyzed them with similar sequences in A. thaliana and B. rapa. Additionally, we included in the analysis published sequences from three other B. oleracea BAC clones and a contig of this species corresponding to segments in A. thaliana chromosomes IV and V. A total of 2,946 kb of B. oleracea, 1,069 kb of B. rapa sequence and 2,607 kb of A. thaliana sequence were compared and analyzed. We found conserved collinearity for gene order and content restricted to specific chromosomal segments, but breaks in collinearity were frequent resulting in gene absence likely not due to gene loss but rearrangements. B. oleracea has the lowest gene density of the three species, followed by B. rapa. The genome expansion of the Brassica species, B. oleracea in particular, is due to larger introns and gene spacers resulting from frequent insertion of DNA transposons and retrotransposons. These findings are discussed in relation to the possible origin and evolution of the Brassica genomes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Sequence and structure of Brassica rapa chromosome A3

Background

The species Brassica rapa includes important vegetable and oil crops. It also serves as an excellent model system to study polyploidy-related genome evolution because of its paleohexaploid ancestry and its close evolutionary relationships with Arabidopsis thaliana and other Brassica species with larger genomes. Therefore, its genome sequence will be used to accelerate both basic research on genome evolution and applied research across the cultivated Brassica species.

Results

We have determined and analyzed the sequence of B. rapa chromosome A3. We obtained 31.9 Mb of sequences, organized into nine contigs, which incorporated 348 overlapping BAC clones. Annotation revealed 7,058 protein-coding genes, with an average gene density of 4.6 kb per gene. Analysis of chromosome collinearity with the A. thaliana genome identified conserved synteny blocks encompassing the whole of the B. rapa chromosome A3 and sections of four A. thaliana chromosomes. The frequency of tandem duplication of genes differed between the conserved genome segments in B. rapa and A. thaliana, indicating differential rates of occurrence/retention of such duplicate copies of genes. Analysis of 'ancestral karyotype' genome building blocks enabled the development of a hypothetical model for the derivation of the B. rapa chromosome A3.

Conclusions

We report the near-complete chromosome sequence from a dicotyledonous crop species. This provides an example of the complexity of genome evolution following polyploidy. The high degree of contiguity afforded by the clone-by-clone approach provides a benchmark for the performance of whole genome shotgun approaches presently being applied in B. rapa and other species with complex genomes.     

Comparative Analysis between Homoeologous Genome Segments of Brassica napus and Its Progenitor Species Reveals Extensive Sequence-Level Divergence

Homoeologous regions of Brassica genomes were analyzed at the sequence level. These represent segments of the Brassica A genome as found in Brassica rapa and Brassica napus and the corresponding segments of the Brassica C genome as found in Brassica oleracea and B. napus. Analysis of synonymous base substitution rates within modeled genes revealed a relatively broad range of times (0.12 to 1.37 million years ago) since the divergence of orthologous genome segments as represented in B. napus and the diploid species. Similar, and consistent, ranges were also identified for single nucleotide polymorphism and insertion-deletion variation. Genes conserved across the Brassica genomes and the homoeologous segments of the genome of Arabidopsis thaliana showed almost perfect collinearity. Numerous examples of apparent transduplication of gene fragments, as previously reported in B. oleracea, were observed in B. rapa and B. napus, indicating that this phenomenon is widespread in Brassica species. In the majority of the regions studied, the C genome segments were expanded in size relative to their A genome counterparts. The considerable variation that we observed, even between the different versions of the same Brassica genome, for gene fragments and annotated putative genes suggest that the concept of the pan-genome might be particularly appropriate when considering Brassica genomes.     

Nuclear organization in crucifer genomes: nucleolus-associated telomere clustering is not a universal interphase configuration in Brassicaceae

Arabidopsis thaliana has become a major plant research model, where interphase nuclear organization exhibits unique features, including nucleolus-associated telomere clustering. The chromocenter (CC)-loop model, or rosette-like configuration, describes intranuclear chromatin organization in Arabidopsis as megabase-long loops anchored in, and emanating from, peripherally positioned CCs, with those containing telomeres associating with the nucleolus. To investigate whether the CC-loop organization is universal across the mustard family (crucifers), the nuclear distributions of centromeres, telomeres and nucleoli were analyzed by fluorescence in situ hybridization in seven diploid species (2n = 10–16) representing major crucifer clades with an up to 26-fold variation in genome size (160–4260 Mb). Nucleolus-associated telomere clustering was confirmed in Arabidopsis (157 Mb) and was newly identified as the major nuclear phenotype in other species with a small genome (215–381 Mb). In large-genome species (2611–4264 Mb), centromeres and telomeres adopted a Rabl-like configuration or dispersed distribution in the nuclear interior; telomeres only rarely associated with the nucleolus. In Arabis cypria (381 Mb) and Bunias orientalis (2611 Mb), tissue-specific patterns deviating from the major nuclear phenotypes were observed in anther and stem tissues, respectively. The rosette-like configuration, including nucleolus-associated telomere clustering in small-genome species from different infrafamiliar clades, suggests that genomic properties rather than phylogenetic position determine the interphase nuclear organization. Our data suggest that nuclear genome size, average chromosome size and degree of longitudinal chromosome compartmentalization affect interphase chromosome organization in crucifer genomes.     

Fine mapping and targeted SNP survey using rice-wheat gene colinearity in the region of the <Emphasis Type="Italic">Bo1</Emphasis> boron toxicity tolerance locus of bread wheat

Toxicity due to high levels of soil boron (B) represents a significant limitation to cereal production in some regions, and the Bo1 gene provides a major source of B toxicity tolerance in bread wheat (Triticum aestivum L.). A novel approach was used to develop primers to amplify and sequence gene fragments specifically from the Bo1 region of the hexaploid wheat genome. Single-nucleotide polymorphisms (SNPs) identified were then used to generate markers close to Bo1 on the distal end of chromosome 7BL. In the 16 gene fragments totaling 19.6 kb, SNPs were observed between the two cultivars Cranbrook and Halberd at a low frequency (one every 613 bp). Furthermore, SNPs were distributed unevenly, being limited to only two genes. In contrast, RFLP provided a much greater number of genetic markers, with every tested gene identifying polymorphism. Bo1 previously known only as a QTL was located as a discrete Mendelian locus. In total, 28 new RFLP, PCR and SSR markers were added to the existing map. The 1.8 cM Bo1 interval of wheat corresponds to a 227 kb section of rice chromosome 6L encoding 21 predicted proteins with no homology to any known B transporters. The co-dominant PCR marker AWW5L7 co-segregated with Bo1 and was highly predictive of B tolerance status within a set of 94 Australian bread wheat cultivars and breeding lines. The markers and rice colinearity described here represent tools that will assist B tolerance breeding and the positional cloning of Bo1. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Combining classifiers to predict gene function in Arabidopsis thaliana using large-scale gene expression measurements

Background  

Arabidopsis thaliana is the model species of current plant genomic research with a genome size of 125 Mb and approximately 28,000 genes. The function of half of these genes is currently unknown. The purpose of this study is to infer gene function in Arabidopsis using machine-learning algorithms applied to large-scale gene expression data sets, with the goal of identifying genes that are potentially involved in plant response to abiotic stress.     

Biochemical genetics of glucosinolate modification in Arabidopsis and Brassica

Fine mapping of the glucosinolate biosynthesis gene OHP, which regulates the conversion of 3-methylsulphinylpropyl to 3-hydroxypropyl glucosinolate, in an Arabidopsis thaliana Columbia × Landsberg erecta RI line population positioned the gene within 54 kb of DNA on chromosome IV. Sequence data identified a family of genes encoding 2-oxoglutarate-dependent dioxygenases in this region. A probe based on these genes co-segregated with ALK in Brassica oleracea,a gene regulating the synthesis of alkenyl glucosinolates. The reactions catalysed by the OHP and ALK enzymes utilise similar substrates and may have a common mechanism. Thus, these dioxygenases are prime candidates for controlling the side chain modification of glucosinolates. Received: 12 May 2000 / Accepted: 29 May 2000     

Sucrose Phosphate Synthase Genes in Plants Belong to Three Different Families

We present phylogenetic analyses to demonstrate that there are three families of sucrose phosphate synthase (SPS) genes present in higher plants. Two data sets were examined, one consisting of full-length proteins and a second larger set that covered a highly conserved region including the 14-3-3 binding region and the UDPGlu active site. Analysis of both datasets showed a well supported separation of known genes into three families, designated A, B, and C. The genomic sequences of Arabidopsis thaliana include a member in each family: two genes on chromosome 5 belong to Family A, one gene on chromosome 1 to Family B, and one gene on chromosome 4 to Family C. Each of three Citrus genes belong to one of the three families. Intron/exon organization of the four Arabidopsis genes differed according to phylogenetic analysis, with members of the same family from different species having similar genomic organization of their SPS genes. The two Family A genes on Arabidopsis chromosome 5 appear to be due to a recent duplication. Analysis of published literature and ESTs indicated that functional differentiation of the families was not obvious, although B family members appear not to be expressed in roots. B family genes were cloned from two Actinidia species and southern analysis indicated the presence of a single gene family, which contrasts to the multiple members of Family A in Actinidia. Only two family C genes have been reported to date. Received: 17 April 2001 / Accepted: 27 August 2001     

Microcolinearity between a 2-cM region encompassing the grain protein content locus<Emphasis Type="Italic"> Gpc-6B1</Emphasis> on wheat chromosome 6B and a 350-kb region on rice chromosome 2

The conservation of the linear order (colinearity) of genetic markers along large chromosome segments in wheat and rice is well established, but less is known about the microcolinearity between both genomes at subcentimorgan distances. In this study we focused on the microcolinearity between a 2.6-cM interval flanked by markers Xcdo365 and Xucw65 on wheat chromosome 6B and rice chromosome 2. A previous study has shown that this wheat segment includes the Gpc-6B1 locus, which is responsible for large differences in grain protein content (GPC) and is the target of a positional cloning effort in our laboratories. Twenty-one recombination events between Xcdo365 and Xucw65 were found in a large segregating population (935 gametes) and used to map 17 genes selected from rice chromosome 2 in the wheat genetic map. We found a high level of colinearity between a 2.1-cM region flanked by loci Xucw75 and Xucw67 on wheat chromosome 6B and a 350-kb uninterrupted sequenced region in rice chromosome arm 2S. Colinearity between these two genomes was extended to the region proximal to Xucw67 (eight colinear RFLP markers), but was interrupted distal to Xucw75 (six non-colinear RFLP markers). Analysis of different comparative studies between rice and wheat suggests that microcolinearity is more frequently disrupted in the distal region of the wheat chromosomes. Fortunately, the region encompassing the Gpc-6B1 locus showed an excellent conservation between the two genomes, facilitating the saturation of the target region of the wheat genetic map with molecular markers. These markers were used to map the Gpc-6B1 locus into a 0.3-cM interval flanked by PCR markers Xucw79 and Xucw71, and to identify five candidate genes within the colinear 64-kb region in rice.     

Fertility restorer locus Rf1 of sorghum (Sorghum bicolor L.) encodes a pentatricopeptide repeat protein not present in the colinear region of rice chromosome 12

With an aim to clone the sorghum fertility restorer gene Rf1, a high-resolution genetic and physical map of the locus was constructed. The Rf1 locus was resolved to a 32-kb region spanning four open reading frames: a plasma membrane Ca²⁺-ATPase, a cyclin D-1, an unknown protein, and a pentatricopeptide repeat (PPR13) gene family member. An ~19-kb region spanning the cyclin D-1 and unknown protein genes was completely conserved between sterile and fertile plants as was the sequence spanning the coding region of the Ca²⁺-ATPase. In contrast, 19 sequence polymorphisms were located in an ~7-kb region spanning PPR13, and all markers cosegregated with the fertility restoration phenotype. PPR13 was predicted to encode a mitochondrial-targeted protein containing a single exon with 14 PPR repeats, and the protein is classified as an E-type PPR subfamily member. To permit sequence-based comparison of the sorghum and rice genomes in the Rf1 region, 0.53 Mb of sorghum chromosome 8 was sequenced and compared to the colinear region of rice chromosome 12. Genome comparison revealed a mosaic pattern of colinearity with an ~275-kb gene-poor region with little gene conservation and an adjacent, ~245-kb gene-rice region that is more highly conserved between rice and sorghum. Despite being located in a region of high gene conservation, sorghum PPR13 was not located in a colinear position on rice chromosome 12. The present results suggest that sorghum PPR13 represents a potential candidate for the sorghum Rf1 gene, and its presence in the sorghum genome indicates a single gene transposition event subsequent to the divergence of rice and sorghum ancestors.An erratum to this article can be found at     

OrthoParaMap: Distinguishing orthologs from paralogs by integrating comparative genome data and gene phylogenies

Background  

In eukaryotic genomes, most genes are members of gene families. When comparing genes from two species, therefore, most genes in one species will be homologous to multiple genes in the second. This often makes it difficult to distinguish orthologs (separated through speciation) from paralogs (separated by other types of gene duplication). Combining phylogenetic relationships and genomic position in both genomes helps to distinguish between these scenarios. This kind of comparison can also help to describe how gene families have evolved within a single genome that has undergone polyploidy or other large-scale duplications, as in the case of Arabidopsis thaliana – and probably most plant genomes.     

Identification of anther-specific genes in a cruciferous model plant, Arabidopsis thaliana, by using a combination of Arabidopsis macroarray and mRNA derived from Brassica oleracea

Arabidopsis thaliana, a member of the Brassicaceae, is a model plant whose genome was the first higher plant genome to be sequenced. Because of the small size of the flowers, it is difficult to dissect and separate reproductive organs (anthers and pistils) at different developmental stages in A. thaliana. In order to perform genome-wide identification of anther-specific genes in A. thaliana, an Arabidopsis cDNA macroarray was hybridized to cDNA derived from anthers and pistils of another crucifer, Brassica oleracea. After scanning the signal intensity for each clone, and cluster analysis, 52 anther-specific genes were identified. These clones contained several anther-specific genes that have already been characterized, as well as novel anther-specific genes. In RT-PCR analysis with mRNA of A. thaliana and B. oleracea, the expression pattern of one-third of the clones was similar to that determined by cDNA macroarray. This system of heterologous hybridization analysis (Arabidopsis cDNA macroarray vs Brassica tissue-specific mRNA) should be applicable to other model species and their close relatives.     

Characterisation of disease resistance gene-like sequences in Brassica oleracea L.

Several cloned disease resistance genes from a wide range of plant species are known to share conserved regions with similar structural motifs. Degenerate primers based on conserved sequences of the nucleotide binding site of the genes RPS2, N and L6 were used for polymerase chain reaction (PCR) amplification from genomic DNA of two doubled haploid lines of Brassica oleracea. Sequences of amplified products were highly variable, but most of them showed similarity to known disease resistance genes, including RPS5, RPS2 and N, and to disease resistance gene-like sequences (RGLs) from different species. Primers based on B. oleracea sequences amplified five groups of RGLs. Products were mapped through cleaved amplified polymorphic sequence assays onto four different linkage groups of B. oleracea. PCR amplification from cDNA and allele analysis indicated that four locus-specific RGL fragments are expressed in cauliflower. Screening of a B. oleracea bacterial artificial chromosome library (BAC) with four B. oleracea RGL probes identified a small number of clones, suggesting that the four RGLs may not be highly copied. Screening of a BAC library of A. thaliana with the same probes identified clones that mapped onto four different chromosomes. These map positions correspond to known disease resistance loci of A. thaliana. Received: 12 November 1999 / Accepted: 19 June 2000     

Retention of triplicated phytoene synthase (PSY) genes in Brassica napus L. and its diploid progenitors during the evolution of the Brassiceae

The extent of genome redundancy exhibited by Brassica species provides a model to study the evolutionary fate of multi-copy genes and the effects of polyploidy in economically important crops. Phytoene synthase (PSY) catalyzes the first committed reaction of the carotenoid biosynthetic pathway, which has been shown to be rate-limiting in Brassica napus seeds. In Arabidopsis thaliana, a single PSY gene (AtPSY) regulates phytoene synthesis in all tissues. Considering that diploid Brassica genomes contain three Arabidopsis-like subgenomes, the objectives of the present work were to determine whether PSY gene families exist in B. napus (AACC) and its diploid progenitor species, Brassica rapa (AA) and Brassica oleracea (CC); to establish the level of retention of Brassica PSY genes; to map PSY gene family members in the A and C genomes and to compare Brassica PSY gene expression patterns. A total of 12 PSY homologues were identified, 6 in B. napus (BnaX.PSY.a-f) and 3 in B. rapa (BraA.PSY.a-c) and B. oleracea (BolC.PSY.a-c). Indeed, with six members, B. napus has the largest PSY gene family described to date. Sequence comparison between AtPSY and Brassica PSY genes revealed a highly conserved gene structure and identity percentages above 85% at the coding sequence (CDS) level. Altogether, our data indicate that PSY gene family expansion preceded the speciation of B. rapa and B. oleracea, dating back to the paralogous subgenome triplication event. In these three Brassica species, all PSY homologues are expressed, exhibiting overlapping redundancy and signs of subfunctionalization among photosynthetic and non-photosynthetic tissues. This evidence supports the hypothesis that functional divergence of PSY gene expression facilitates the accumulation of high levels of carotenoids in chromoplast-rich tissues. Thus, functional retention of triplicated Brassica PSY genes could be at least partially explained by the selective advantage provided by increased levels of gene product in floral organs. A better understanding of carotenogenesis in Brassica will aid in the future development of transgenic and conventional cultivars with carotenoid-enriched oil.     

Isolation of single-copy-sequence clones from a yeast artificial chromosome library of randomly-sheared Arabidopsis thaliana DNA

We describe the construction of a yeast artificial chromosome (YAC) library from the Arabidopsis thaliana genome. Randomly sheared high molecular weight source DNA was extracted from frozen, ground leaf tissue and blunt-end-ligated to the vector pYAC3. By size-fractionating the ligation products, we achieved an average clone size of 150 kb. Approximately 6% of the YACs contained inserts from the chloroplast genome. We screened clones equivalent to greater than four A. thaliana haploid nuclear genomes and isolated YACs homologous to five single-copy-sequence probes. The library should be useful chromosome walking and genome mapping experiments. In addition, the approach used for its construction should be applicable to other higher plant species.