Genetic and physical mapping inBrassica diploid species of a gene cluster defined inArabidopsis thaliana

We report the genetic and physical analysis by pulse field gel electrophoresis (PFGE) in threeBrassica diploid genomes for a cluster of five genes characterized in a selected segment of 15 kb on chromosome 3 ofArabidopsis thaliana, encoding aBradyrhizobium CycJ homologue (At1), a rat p67 translation factor homologue (At2), an Em-like (early methionine) protein (At3), chlorophyll synthase (At4) and a yeast Sac1 homologue (A5). TheArabidopsis gene array was found to be conserved on a single linkage group in each of theBrassica genomes. However, partial complexes were found to be duplicated in other chromosome segments on the same or other linkage groups. Some of the At genes, which could not be genetically mapped because of lack of polymorphism, were assigned to their respective linkage groups by physical mapping. The presence of multiple copies of theA. thaliana gene cluster in the threeBrassica genomes further establishes their complex nature, which results from extensive duplication and chromosomal rearrangement. In general, genetic distances between the At genes agreed with values expected for the physical distances determined inBrassica.     

The physical map ofMus musculus chromosome 11 reveals evolutionary relationship with different syntenic groups of genes inHomo sapiens

Summary The physical localization of sequences homologous to three cloned genes was determined by in situ hybridization to metaphase chromosomes. Previous work had assigned the skeletal myosin heavy chain gene cluster (Myh), the functional locus for the cellular tumor antigen p53 (Trp53-1), and the cellular homologue of the viral erb-B oncogene (Erbb) toMus musculus chromosome 11 (MMU11). Our results provide regional assignments ofMyh andTrp53-1 to chromosome bands B2C, and ofErbb to bands A1A4. Taken together with in situ mapping of three other loci on MMU 11 (Hox-2 homeobox-containing gene cluster, theSparc protein, and theColla-1 collagen gene), which have been reported elsewhere, these data allowed us to construct a physical map of MMU11 and to compare it with the linkage map of this chromosome. The map positions of the homologous genes on human chromosomes suggest evolutionary relationships of distinct regions of MMU11 with six different human chromosome arms: 1p, 5q, 7p, 16p, 17p, and 17q. The delineation of conserved chromosome regions has important implications for the understanding of karyotype evolution in mammalian species and for the development of animal models of human genetic diseases.     

Organization of an Arabidopsis thaliana gene cluster on chromosome 4 including the RPS2 gene, in the Brassica nigra genome

 Genetic and physical maps, consisting of a large number of DNA markers for Arabidopsis thaliana chromosomes, represent excellent tools to determine the organization of related genomes such as those of Brassica. In this paper we report the chromosomal localization and physical analysis by pulsed-field gel electrophoresis (PFGE) of a well-defined gene complex of A. thaliana in the Brassica nigra genome (B genome n=8). This complex is approximately 30 kb in length in A. thaliana and contains a cluster of six genes including ABI1 (ABA-responsive), RPS2 (resistance against Pseudomonas syringae, a bacterial disease), CK1 (casein kinase I), NAP (nucleosome-assembly protein), X9 and X14 (both of unknown function). The Arabidopsis chromosomal complex was found to be duplicated and conserved in gene number at different levels in the Brassica genome. Linkage group B1 had the most-conserved arrangement carrying all six genes tightly linked. Group B4 had an almost complete complex except for the absence of RPS2. Other partial complexes of fewer members were found on three other chromosomes. Our studies demonstrate that by this approach it is possible to identify ancestrally related chromosome segments in a complex and duplicated genome, such as the genome of B. nigra, permitting one to draw conclusions as to its origin and evolution. Received: 11 July 1997 / Accepted: 9 October 1997     

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.     

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.     

Mapping in the realm of polyploidy: The wheat model

Wheat is an allopolyploid containing three distinct but genetically related (homoeologous) genomes, A, B and D. Because of polyploid inheritance and large genome size (16×10¹² bp), the wheat genome is thought to be intractable to map-based cloning of agronomic and other genes of interest. We propose a targeted geneti mapping strategy that combines linkage and physical mapping and may facilitate map-based cloning. High-density linkage maps are either generated in wheat or in diploid Triticum tauschii, the donor of the D genome to wheat. Molecular marker-based chromosome maps are constructed, using an array of deletion lines in wheat. The conventional genetic linkage maps are aligned with chromosome maps to construct cytogenetic ladder maps (CLMs). The CLMs allow region-specific mapping and convert genetic distances into physical distances. The information from CLMs suggests that many genes in wheat are present in clusters that are highly recombiogenic, small, and may be amenable to cloning by chromosome walking. Therefore, the effective genome size of wheat is relatively small in comparison to the whole genome. The utility of using the smaller genome of rice for mapping and homologous gene cloning is discussed.     

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.     

Extensive Chromosome Homoeology among Brassiceae Species Were Revealed by Comparative Genetic Mapping with High-Density EST-Based SNP Markers in Radish (Raphanus sativus L.)

A linkage map of expressed sequence tag (EST)-based markers in radish (Raphanus sativus L.) was constructed using a low-cost and high-efficiency single-nucleotide polymorphism (SNP) genotyping method named multiplex polymerase chain reaction–mixed probe dot-blot analysis developed in this study. Seven hundred and forty-six SNP markers derived from EST sequences of R. sativus were assigned to nine linkage groups with a total length of 806.7 cM. By BLASTN, 726 markers were found to have homologous genes in Arabidopsis thaliana, and 72 syntenic regions, which have great potential for utilizing genomic information of the model species A. thaliana in basic and applied genetics of R. sativus, were identified. By construction and analysis of the genome structures of R. sativus based on the 24 genomic blocks within the Brassicaceae ancestral karyotype, 23 of the 24 genomic blocks were detected in the genome of R. sativus, and half of them were found to be triplicated. Comparison of the genome structure of R. sativus with those of the A, B, and C genomes of Brassica species and that of Sinapis alba L. revealed extensive chromosome homoeology among Brassiceae species, which would facilitate transfer of the genomic information from one Brassiceae species to another.     

BAC-HAPPY Mapping (BAP Mapping): A New and Efficient Protocol for Physical Mapping

Physical and linkage mapping underpin efforts to sequence and characterize the genomes of eukaryotic organisms by providing a skeleton framework for whole genome assembly. Hitherto, linkage and physical “contig” maps were generated independently prior to merging. Here, we develop a new and easy method, BAC HAPPY MAPPING (BAP mapping), that utilizes BAC library pools as a HAPPY mapping panel together with an Mbp-sized DNA panel to integrate the linkage and physical mapping efforts into one pipeline. Using Arabidopsis thaliana as an exemplar, a set of 40 Sequence Tagged Site (STS) markers spanning ∼10% of chromosome 4 were simultaneously assembled onto a BAP map compiled using both a series of BAC pools each comprising 0.7x genome coverage and dilute (0.7x genome) samples of sheared genomic DNA. The resultant BAP map overcomes the need for polymorphic loci to separate genetic loci by recombination and allows physical mapping in segments of suppressed recombination that are difficult to analyze using traditional mapping techniques. Even virtual “BAC-HAPPY-mapping” to convert BAC landing data into BAC linkage contigs is possible.     

The Tarenaya hassleriana Genome Provides Insight into Reproductive Trait and Genome Evolution of Crucifers

The Brassicaceae, including Arabidopsis thaliana and Brassica crops, is unmatched among plants in its wealth of genomic and functional molecular data and has long served as a model for understanding gene, genome, and trait evolution. However, genome information from a phylogenetic outgroup that is essential for inferring directionality of evolutionary change has been lacking. We therefore sequenced the genome of the spider flower (Tarenaya hassleriana) from the Brassicaceae sister family, the Cleomaceae. By comparative analysis of the two lineages, we show that genome evolution following ancient polyploidy and gene duplication events affect reproductively important traits. We found an ancient genome triplication in Tarenaya (Th-α) that is independent of the Brassicaceae-specific duplication (At-α) and nested Brassica (Br-α) triplication. To showcase the potential of sister lineage genome analysis, we investigated the state of floral developmental genes and show Brassica retains twice as many floral MADS (for MINICHROMOSOME MAINTENANCE1, AGAMOUS, DEFICIENS and SERUM RESPONSE FACTOR) genes as Tarenaya that likely contribute to morphological diversity in Brassica. We also performed synteny analysis of gene families that confer self-incompatibility in Brassicaceae and found that the critical SERINE RECEPTOR KINASE receptor gene is derived from a lineage-specific tandem duplication. The T. hassleriana genome will facilitate future research toward elucidating the evolutionary history of Brassicaceae genomes.     

Peaks of linkage are localized by a BAC/PAC contig of the 6p reading disability locus.

A gene for reading disability has been localized by nonparametric linkage to 6p21.3-p22 in several published reports. However, the lack of an uninterrupted genomic clone contig has made it difficult to determine accurate intermarker distances, precise marker order, and genetic boundaries and hinders direct comparisons of linkage. The search and discovery of the hemochromatosis gene (HFE) led to the creation of a bacterial artificial chromosome (BAC) and P-1 derived artificial chromosome (PAC) contig that extended physical maps 4 Mb from the MHC toward pter and localized new markers in that region [10-12]. Using this contig, we localized 124 sequence tagged sites, expressed sequence tags, and short tandem repeats including most of the markers in linkage with reading disability phenotypes, succinic semialdehyde dehydrogenase, GPLD1, prolactin, and 18 uncharacterized genes. This new contig joins and extends previously published physical maps to span the entire chromosome 6 reading disability genetic locus. Physical mapping data from the complete contig show overlap of the published linkage peaks for reading disability, provide accurate intermarker distances and order, and offer resources for generating additional markers and candidate genes for high resolution genetic studies in this region.     

A comparative linkage map of oilseed rape and its use for QTL analysis of seed oil and erucic acid content

We have developed a new DH mapping population for oilseed rape, named TNDH, using genetically and phenotypically diverse parental lines. We used the population in the construction of a high stringency genetic linkage map, consisting of 277 loci, for use in quantitative genetic analysis. A proportion of the markers had been used previously in the construction of linkage maps for Brassica species, thus permitting the alignment of maps. The map includes 68 newly developed Sequence Tagged Site (STS) markers targeted to the homologues of defined genes of A. thaliana. The use of these markers permits the alignment of our linkage map with the A. thaliana genome sequence. An additional 74 loci (31 newly developed STS markers and 43 loci defined by SSR and RFLP markers that had previously been used in published linkage maps) were added to the map. These markers increased the resolution of alignment of the newly constructed linkage map with existing Brassica linkage maps and the A. thaliana genome sequence. We conducted field trials with the TNDH population at two sites, and over 2 years, and identified reproducible QTL for seed oil content and erucic acid content. The results provide new insights into the genetic control of seed oil and erucic acid content in oilseed rape, and demonstrate the utility of the linkage map and population.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.D. Qiu and C. Morgan authors contributed equally to the work.     

Functional innovations of three chronological mesohexaploid Brassica rapa genomes

Background

The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa.

Results

Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were ‘stress response’ and ‘development’ both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization.

Conclusion

We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Electronic supplementary material

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

Microdissection and microcloning of plant chromosomes

Cytogenetic and molecular tools play an increasingly important role in plant genome research. A number of interesting applications that involve chromosome painting, the relationship between specific chromosomes and specific linkage groups, the relationships between physical and genetic distances on linkage maps, and the isolation of genes of interest, have been the subjects of recently published research. The aim of this paper is to review the different techniques available for chromosome microdissection and microcloning, and their use for the study of plant genomes. The quality of chromosomal DNA obtained is considered, and some recent results from our laboratory are presented.     

Cloning of TTG1 gene and PCR identification of genomes A,B and C in Brassica species

Arabidopsis Transparent Testa Glabra 1 (TTG1) genes were cloned from three diploid Brassica species (B. rapa, B. nigra and B. oleracea) and two amphidiploids species (B. juncea and B. carinata) by homology cloning. TTG1 homologues identified in all the accessions of the investigated species had a coding sequence of 1,014 bp. One copy was obtained from each diploid species and two copies from each amphidiploid species. Combined analysis of the TTG1 sequences cloned in this study with those obtained from public databases demonstrated that three, forty-five and seven nucleotides were specific variations in TTG1 genes from genomes A, B and C, respectively. Primers designed with genome-specific nucleotide variations were able to distinguish among TTG1 genes originating from genomes A, B and C in Brassica. Therefore, the TTG1 gene could serve as a candidate marker gene to detect the pollen flow of Brassica and provide an alternative method for the detection of pollen drift and risk assessment of gene flow in Brassica species.     

Inter- and intra-genomic homology of the Brassica genomes: implications for their origin and evolution

In order to determine the homologous regions shared by the cultivated Brassica genomes, linkage maps of the diploid cultivated B. rapa (A genome, n = 10), B. nigra (B genome, n = 8) and B. oleracea (C genome, n = 9), were compared. We found intergenomic conserved regions but with extensitve reordering among the genomes. Eighteen linkage groups from all three species could be associated on the basis of homologous segments based on at least three common markers. Intragenomic homologous conservation was also observed for some of the chromosomes of the A, B and C genomes. A possible chromosome phylogenetic pathway based on an ancestral genome of at least five, and no more than seven chromosomes, was drawn from the chromosomal inter-relationships observed. These results demonstrate that extensive duplication and rearrangement have been involved in the formation of the Brassica genomes from a smaller ancestral genome.     

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.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.