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.     

Isolation of circadian-associated genes in Brassica rapa by comparative genomics with Arabidopsis thaliana

Elucidation of the roles of circadian associated factors requires a better understanding of the molecular mechanisms of circadian rhythms, control of flowering time through photoperiodic pathways, and photosensory signal transduction. In Arabidopsis, the APRR1 quintet, APRRs 1, 3, 5, 7, and 9, are known as central oscillator genes. Other plants may share the molecular mechanism underlying the circadian rhythm. To identify and characterize these circadian response genes in Brassica crops whose genome was triplicated after divergence from Arabidopsis, we identified B. rapa BAC clones containing these genes by BLAST analysis of B. rapa BAC end sequences against the five corresponding Arabidopsis regions. Subsequent fingerprinting, Southern hybridization, and PCR allowed identification of five BAC clones, one for each of the five circadian-related genes. By draft shotgun sequencing of the BAC clones, we identified the complete gene sequences and cloned the five expressed B. rapa circadian-associated gene members, BrPRRs 1, 3, 5, 7, and 9. Phylogenetic analysis revealed that each BrPRR was orthologous to the corresponding APRR at the sequence level. Northern hybridization revealed that the five genes were transcribed at distinct points in the 24 hour period, and Southern hybridization revealed that they are present in 2, 1, 2, 2, and 1 copies, respectively in the B. rapa genome, which was triplicated and then diploidized during the last 15 million years.     

Segmental structure of the Brassica napus genome based on comparative analysis with Arabidopsis thaliana

Over 1000 genetically linked RFLP loci in Brassica napus were mapped to homologous positions in the Arabidopsis genome on the basis of sequence similarity. Blocks of genetically linked loci in B. napus frequently corresponded to physically linked markers in Arabidopsis. This comparative analysis allowed the identification of a minimum of 21 conserved genomic units within the Arabidopsis genome, which can be duplicated and rearranged to generate the present-day B. napus genome. The conserved regions extended over lengths as great as 50 cM in the B. napus genetic map, equivalent to approximately 9 Mb of contiguous sequence in the Arabidopsis genome. There was also evidence for conservation of chromosome landmarks, particularly centromeric regions, between the two species. The observed segmental structure of the Brassica genome strongly suggests that the extant Brassica diploid species evolved from a hexaploid ancestor. The comparative map assists in exploiting the Arabidopsis genomic sequence for marker and candidate gene identification within the larger, intractable genomes of the Brassica polyploids.     

Mapping and analysis of simple sequence repeats in the Arabidopsis thaliana genome

Simple sequence repeats (SSRs) are becoming standard DNA markers for plant genome analysis and are being used as markers in marker assisted breeding. And hence because of its great significance we have initiated this study to analyze complete genome of Arabidopsis thaliana for the prevalence of mono-, di-, tri-, tetra-, penta- and hexa- mer repeats in the coding and non-coding regions of the chromosome and to map their exact position on the sequence. We have developed a program that can search a repeat of any length, its exact position on the chromosome and also its frequency of occurrence in the genome. Analysis of the results reveal that maximum number of repeats were found in chromosome 1 followed by chromosome 2 and 4 whereas, chromosome 3 and 5 contain relatively less number of these repeats. Among the SSRs, hexamers and dimers were more predominant in the chromosomes. Overall data showed that Chromosome 5 has minimum number of repeats. The abundance or rarity of various simple repeats in different chromosomes is not explained by nucleotide composition of sequence or potential repeated motifs to form alternative DNA structures. This suggests that in addition to nucleotide composition of repeat motifs, characteristic DNA replication / repair / recombination machinery might play an important role in genesis of repeats. The positional information is given at www.geocities.com/amubioinfo/ARD. This positional information can help Arabidopsis researchers to identify new polymorphisms in chromosomal regions of interest based on the SSRs that map in the area.     

Simple sequence repeat-based comparative genomics between Brassica rapa and Arabidopsis thaliana: the genetic origin of clubroot resistance

An SSR-based linkage map was constructed in Brassica rapa. It includes 113 SSR, 87 RFLP, and 62 RAPD markers. It consists of 10 linkage groups with a total distance of 1005.5 cM and an average distance of 3.7 cM. SSRs are distributed throughout the linkage groups at an average of 8.7 cM. Synteny between B. rapa and a model plant, Arabidopsis thaliana, was analyzed. A number of small genomic segments of A. thaliana were scattered throughout an entire B. rapa linkage map. This points out the complex genomic rearrangements during the course of evolution in Cruciferae. A 282.5-cM region in the B. rapa map was in synteny with A. thaliana. Of the three QTL (Crr1, Crr2, and Crr4) for clubroot resistance identified, synteny analysis revealed that two major QTL regions, Crr1 and Crr2, overlapped in a small region of Arabidopsis chromosome 4. This region belongs to one of the disease-resistance gene clusters (MRCs) in the A. thaliana genome. These results suggest that the resistance genes for clubroot originated from a member of the MRCs in a common ancestral genome and subsequently were distributed to the different regions they now inhabit in the process of evolution.     

Coding sequence divergence between two closely related plant species: Arabidopsis thaliana and Brassica rapa ssp. pekinensis

To characterize the coding-sequence divergence of closely related genomes, we compared DNA sequence divergence between sequences from a Brassica rapa ssp. pekinensis EST library isolated from flower buds and genomic sequences from Arabidopsis thaliana. The specific objectives were (i) to determine the distribution of and relationship between K _a and K _s, (ii) to identify genes with the lowest and highest K _a:K _s values, and (iii) to evaluate how codon usage has diverged between two closely related species. We found that the distribution of K _a:K _s was unimodal, and that substitution rates were more variable at nonsynonymous than synonymous sites, and detected no evidence that K _a and K _s were positively correlated. Several genes had K _a:K _s values equal to or near zero, as expected for genes that have evolved under strong selective constraint. In contrast, there were no genes with K _a:K _s >1 and thus we found no strong evidence that any of the 218 sequences we analyzed have evolved in response to positive selection. We detected a stronger codon bias but a lower frequency of GC at synonymous sites in A. thaliana than B. rapa. Moreover, there has been a shift in the profile of most commonly used synonymous codons since these two species diverged from one another. This shift in codon usage may have been caused by stronger selection acting on codon usage or by a shift in the direction of mutational bias in the B. rapa phylogenetic lineage.     

Comparison of a Brassica oleracea genetic map with the genome of Arabidopsis thaliana

Brassica oleracea is closely related to the model plant, Arabidopsis thaliana. Despite this relationship, it has been difficult to both identify the most closely related segments between the genomes and determine the degree of genome replication within B. oleracea relative to A. thaliana. These difficulties have arisen in part because both species have replicated genomes, and the criteria used to identify orthologous regions between the genomes are often ambiguous. In this report, we compare the positions of sequenced Brassica loci with a known position on a B. oleracea genetic map to the positions of their putative orthologs within the A. thaliana genome. We use explicit criteria to distinguish orthologous from paralogous loci. In addition, we develop a conservative algorithm to identify collinear loci between the genomes and a permutation test to evaluate the significance of these regions. The algorithm identified 34 significant A. thaliana regions that are collinear with >28% of the B. oleracea genetic map. These regions have a mean of 3.3 markers spanning 2.1 Mbp of the A. thaliana genome and 2.5 cM of the B. oleracea genetic map. Our findings are consistent with the hypothesis that the B. oleracea genome has been highly rearranged since divergence from A. thaliana, likely as a result of polyploidization.     

The Arabidopsis genome sequence as a tool for genome analysis in Brassicaceae. A comparison of the Arabidopsis and Capsella rubella genomes

The annotated Arabidopsis genome sequence was exploited as a tool for carrying out comparative analyses of the Arabidopsis and Capsella rubella genomes. Comparison of a set of random, short C. rubella sequences with the corresponding sequences in Arabidopsis revealed that aligned protein-coding exon sequences differ from aligned intron or intergenic sequences in respect to the degree of sequence identity and the frequency of small insertions/deletions. Molecular-mapped markers and expressed sequence tags derived from Arabidopsis were used for genetic mapping in a population derived from an interspecific cross between Capsella grandiflora and C. rubella. The resulting eight Capsella linkage groups were compared to the sequence maps of the five Arabidopsis chromosomes. Fourteen colinear segments spanning approximately 85% of the Arabidopsis chromosome sequence maps and 92% of the Capsella genetic linkage map were detected. Several fusions and fissions of chromosomal segments as well as large inversions account for the observed arrangement of the 14 colinear blocks in the analyzed genomes. In addition, evidence for small-scale deviations from genome colinearity was found. Colinearity between the Arabidopsis and Capsella genomes is more pronounced than has been previously reported for comparisons between Arabidopsis and different Brassica species.     

Contrasting genome organisation: two regions of the Brassica oleracea genome compared with collinear regions of the Arabidopsis thaliana genome.

Brassica crop species are of worldwide importance and are closely related to the model plant Arabidopsis thaliana for which the complete genome sequence has recently been established. We investigated collinearity of marker order by comparing two contrasting regions of the Brassica oleracea genome with homologous regions of A. thaliana. Although there is widespread replication of marker loci in both A. thaliana and B. oleracea, we found that a combination of genetic markers mapped in B. oleracea, including RFLPs, CAPS, and SSRs allowed comparison and interpretation of medium-scale chromosomal organisation and rearrangements. The interpretation of data was facilitated by hybridising probes onto the whole A. thaliana genome, as represented by BAC contigs. Twenty marker loci were sampled from the whole length of the shortest B. oleracea linkage group, 06, and 21 from a 30.4-cM section of the longest linkage group, 03. There is evidence of locus duplication on linkage group 06. Locus order is well conserved between a putative duplicated region of 10.5 cM and a discrete region comprising 25 cM of A. thaliana chromosome I. This was supported by evidence from seven paralogous loci, three of which were duplicated in a 30.6-cM region of linkage group 06. The pattern of locus order for the remainder of linkage group 06 and the sampled section of linkage group 03 was more complex when compared with the A. thaliana genome. Although there was some conservation of locus order between markers on linkage group 03 and approximately 9 cM of A. thaliana chromosome I, this was superimposed upon a complex pattern of additional loci that were replicated in both A. thaliana and B. oleracea. The results are discussed in the context of the ability to use collinear information to assist map-based cloning.     

An integrated map of Arabidopsis thaliana for functional analysis of its genome sequence.

The genome of the model plant species Arabidopsis thaliana has recently been sequenced. To accelerate its current genome research, we developed a whole-genome, BAC/BIBAC-based, integrated physical, genetic, and sequence map of the A. thaliana ecotype Columbia. This new map was constructed from the clones of a new plant-transformation-competent BIBAC library and is integrated with the existing sequence map. The clones were restriction fingerprinted by DNA sequencing gel-based electrophoresis, assembled into contigs, and anchored to an existing genetic map. The map consists of 194 BAC/BIBAC contigs, spanning 126 Mb of the 130-Mb Arabidopsis genome. A total of 120 contigs, spanning 114 Mb, were anchored to the chromosomes of Arabidopsis. Accuracy of the integrated map was verified using the existing physical and sequence maps and numerous DNA markers. Integration of the new map with the sequence map has enabled gap closure of the sequence map and will facilitate functional analysis of the genome sequence. The method used here has been demonstrated to be sufficient for whole-genome physical mapping from large-insert random bacterial clones and thus is applicable to rapid development of whole-genome physical maps for other species.     

Genetic similarity among Arabidopsis thaliana ecotypes estimated by DNA sequence comparison

DNA polymorphisms among Arabidopsis thaliana ecotypes are widely used as genetic markers in map-based cloning strategies. New PCR-based molecular markers do not only facilitate molecular mapping, but can also be used to obtain reliable sequence information for cladistic analyses. We have used CAPS (cleaved amplified polymorphic sequences) markers and a direct sequencing strategy to estimate genetic similarity among eighteen Arabidopsis ecotypes. Sequences at four loci, two from the nuclear and two from a non-nuclaar genome, were analysed. For each ecotype more than 1000pb of sequence information was obtained, and genetic similarity was calculated from a total of 35 polymorphic sites using a character-based approach. Divergence ranged from zero up to 50 discordant characters among the 72 characters defined by the polymorphisms. Separate calculations based on the nuclear and the non-nuclear sequences were performed and revealed a number of common features, including the existence of small clusters of very closely related ecotypes separated from each other by extensive sequence divergence. Our results provide information useful especially to investigators setting up crosses for chromosome landing strategies.     

Toward unraveling the structure of Brassica rapa genome

Genomic research in any organism encompasses understanding structure of the target genome and genes, their function, and evolution. Brassica rapa , which is phylogenetically related to Arabidopsis thaliana , is an important species with respect to its uses as vegetable, oil, and fodder. The availability of suitable genetic and genomic resources is a prerequisite to undertake genomic research in B. rapa . We have developed reference mapping populations of Chinese cabbage ( B. rapa ssp. pekinensis ) comprising 78 doubled haploid lines and over 250 recombinant inbred lines. Two Bacterial Artificial Chromosome (BAC) libraries, generated by restriction enzymes Hin dIII (KBrH) and Bam HI (KBrB), comprise 56 592 and 50 688 clones, respectively. We have also constructed 22 cDNA libraries from different plant tissues consisting of 104 914 clones with an average length of 575 bp. Initial BAC-end sequence analysis of 1473 clones of the KBrH library led us to understand the structure of B. rapa genome with respect to extent of genic sequences and their annotation, and relative abundance of different types of repetitive DNAs. Full-length sequence analysis of BAC clones revealed extensive triplication of B. rapa DNA segments coupled with variable gene losses within the segments. The formulation of the 'Multinational Brassica Genome Project' has laid the foundation to sequence the complete genome of B. rapa ssp. pekinensis by the international Brassica research community. It has been proposed to undertake BAC-to-BAC sequencing of genetically mapped seed BACs. In recent years, development of bioinformatics tools in Brassica has given a boost to structural genomics research in Brassica species. The research undertaken with the availability of various genomic resources in the public domain has added to our understanding of the structure of B. rapa .     

Evidence for homology of flowering-time genes VFR2 from Brassica rapa and FLC from Arabidopsis thaliana

Variation in flowering time is important for the adaptation of plant species to different natural and agricultural environments. We previously identified VFR2 as one of two major quantitative trait loci (QTLs) controlling vernalization-responsive flowering time in a segregating population derived from a cross of annual and biennial Brassica rapa. The region containing VFR2 is homologous to a region in Brassica napus that controls the same trait, and also to a region on chromosome 5 of Arabidopsis thaliana that contains several flowering-time loci. In order to determine precisely the allelic effects and map position of VFR2, we backcrossed the late allele into an early flowering line and obtained monogenic segregation for flowering time in a BC₃S₁ population. The two homozygous genotypic classes differed by 43 and 95 days to flowering in the field and growth chamber, respectively; and the effect of the late allele was almost completely additive. DNA probes that were previously shown to detect RFLP loci in the VFR2 region, or in the homologous regions of B. napus or A. thaliana (including two DNA clones of flowering-time genes) were used to construct a high-resolution map around VFR2. An RFLP detected by an A. thaliana cDNA clone of flowering locus C (FLC) co-segregated exactly with VFR2 in 414 gametes analyzed (<0.24 cM). FLC is a repressor of flowering and is required for the winter-annual habit of late-flowering ecotypes of A. thaliana. The regulation of FLC RNA in B. rapa was consistent with that seen in A. thaliana; RNA levels were up-regulated in the biennial parent and down-regulated by cold treatment. Thus, VFR2 appears to be homologous to FLC and may control flowering time though a similar mechanism as in A. thaliana .     

Pong-like elements in Arabidopsis and Brassica rapa: its regulation of F-box protein gene in different ecotypes of Arabidopsis thaliana

Pong elements are active class 2 transposons in rice. We found Pong-like elements in Arabidopsis thaliana and named them as AtPong (Arabidopsis thaliana Pong) elements. The AtPong elements carried 13 bp of TIRs and two ORFs in which NAM DNA binding domain and DDE catalytic domain are encoded. Ping and mPing (miniature Ping) elements are deficient Pong elements and we found AtPong derived deficient AtPing and AtmPing elements. We also found a homologous element of the AtPong element in Brassica rapa. This element was a deficient element by internal deletion, but showed high sequence similarity with the AtmPing so that it was named as BrmPing (Brassica rapa miniature Ping). The AtPong, AtPing, and AtmPing elements were present in intergenic regions except one element, AtPing1, which was present in an exon of a F-box protein gene. Among the different A. thaliana ecotypes, the AtPing1 showed polymorphisms of present and absent in the F-box protein gene. The excision of the AtPing1 restored the expression of the F-box protein gene, indicating that the expression of the F-box protein genes is regulated by the AtPing1.     

Analysis of the myosins encoded in the recently completed Arabidopsis thaliana genome sequence

Background

Three types of molecular motors play an important role in the organization, dynamics and transport processes associated with the cytoskeleton. The myosin family of molecular motors move cargo on actin filaments, whereas kinesin and dynein motors move cargo along microtubules. These motors have been highly characterized in non-plant systems and information is becoming available about plant motors. The actin cytoskeleton in plants has been shown to be involved in processes such as transportation, signaling, cell division, cytoplasmic streaming and morphogenesis. The role of myosin in these processes has been established in a few cases but many questions remain to be answered about the number, types and roles of myosins in plants.

Results

Using the motor domain of an Arabidopsis myosin we identified 17 myosin sequences in the Arabidopsis genome. Phylogenetic analysis of the Arabidopsis myosins with non-plant and plant myosins revealed that all the Arabidopsis myosins and other plant myosins fall into two groups - class VIII and class XI. These groups contain exclusively plant or algal myosins with no animal or fungal myosins. Exon/intron data suggest that the myosins are highly conserved and that some may be a result of gene duplication.

Conclusions

Plant myosins are unlike myosins from any other organisms except algae. As a percentage of the total gene number, the number of myosins is small overall in Arabidopsis compared with the other sequenced eukaryotic genomes. There are, however, a large number of class XI myosins. The function of each myosin has yet to be determined.     

Characterization of the genome of Arabidopsis thaliana

The small crucifer Arabidopsis thaliana has many useful features as an experimental organism for the study of plant molecular biology. It has a four-week life-cycle, only five chromosomes and a genome size less than half that of Drosophila. To characterize the DNA sequence organization of this plant, we have randomly selected 50 recombinant lambda clones containing inserts with an average length of 12,800 base-pairs and analyzed their content of repetitive and unique DNA by various genome blot, restriction digestion and RNA blot procedures. The following conclusions can be drawn. The DNA represented in this random sample is composed predominantly of single-copy sequences. This presumably reflects the organization of the Arabidopsis genome as a whole and supports prior conclusions reached on the basis of kinetics of DNA reassociation. The DNA that encodes the ribosomal RNAs constitutes the only major class of cloned nuclear repetitive DNA. It consists of approximately 570 tandem copies of a heterogeneous 9900-base-pair repeat unit. There is an average of approximately 660 copies of the chloroplast genome per cell. Therefore, the chloroplast genome constitutes the major component of the repetitive sequences found in A. thaliana DNA made from whole plants. The inner cytosine residue in the sequence C-C-G-G is methylated more often than the outer in the tandem ribosomal DNA units, whereas very few differences in the methylation state of these two cytosine residues are detected in unique sequences.     

Identification and structural analysis of SINE elements in the Arabidopsis thaliana genome

An insertion sequence was found in a Mu homologue in the genome of Arabidopsis thaliana. The insertion sequence had poly(A) at the 3' end, and promoter motifs (A- and B-boxes) recognized by RNA polymerase III. The sequence was flanked by direct repeats of a 15-bp sequence of the Mu homologue, which appears to be a target-site sequence duplicated upon insertion. These findings indicate that the insertion sequence is a retroposon SINE, and it was therefore named AtSN (A. thaliana SINE). Many members of the AtSN family were identified through a computer-aided homology search of databases and classified into two subfamilies, AtSN1 and AtSN2, having consensus sequences 159 and 149 bp in length, respectively. These had no homology to SINEs in other organisms. About half of AtSN members were truncated through loss of a region at either end of the element. Most of them were truncated at the 5' end, and had a duplication of the target-site sequence. This suggests that the ones with 5' truncation retroposed by the same mechanism as those without truncation. Members of the AtSN1 or AtSN2 subfamilies had many base substitutions when compared with the consensus sequence. All of the members examined were present in three different ecotypes of A. thaliana (Columbia, Landsberg erecta, and Wassilewskija). These findings suggest that AtSN members had proliferatedbefore the A. thaliana ecotype strains diverged.     

The shrunken genome of Arabidopsis thaliana

This paper examines macro and micro-level patterns of genome size evolution in the Brassicaceae. A phylogeny of 25 relatives of Arabidopsis thaliana was reconstructed using four molecular markers under both parsimony and Bayesian methods. Reconstruction of genome size (C value) evolution as a discrete character and as a continuous character was also performed. In addition, size dynamics in small chromosomal regions were assessed by comparing genomic clones generated for Arabidopsis lyrata and for Boechera stricta to the fully sequenced genome of A. thaliana. The results reveal a sevenfold variation in genome size among the taxa investigated and that the small genome size of A. thaliana is derived. Our results also indicate that the genome is free to increase or decrease in size across these evolutionary lineages without a directional bias. These changes are accomplished by insertions and deletions at both large and small-scales occurring mostly in intergenic regions, with repetitive sequences and transposable elements implicated in genome size increases. The focus upon taxa relatively closely related to the model organism A. thaliana, and the combination of complementary approaches, allows for unique insights into the processes driving genome size changes.