- Research article
- Open Access
Development of genomic resources for Citrus clementina: Characterization of three deep-coverage BAC libraries and analysis of 46,000 BAC end sequences
© Terol et al; licensee BioMed Central Ltd. 2008
- Received: 29 January 2008
- Accepted: 18 September 2008
- Published: 18 September 2008
Citrus species constitute one of the major tree fruit crops of the subtropical regions with great economic importance. However, their peculiar reproductive characteristics, low genetic diversity and the long-term nature of tree breeding mostly impair citrus variety improvement. In woody plants, genomic science holds promise of improvements and in the Citrus genera the development of genomic tools may be crucial for further crop improvements. In this work we report the characterization of three BAC libraries from Clementine (Citrus clementina), one of the most relevant citrus fresh fruit market cultivars, and the analyses of 46.000 BAC end sequences. Clementine is a diploid plant with an estimated haploid genome size of 367 Mb and 2n = 18 chromosomes, which makes feasible the use of genomics tools to boost genetic improvement.
Three genomic BAC libraries of Citrus clementina were constructed through Eco RI, Mbo I and Hind III digestions and 56,000 clones, representing an estimated genomic coverage of 19.5 haploid genome-equivalents, were picked. BAC end sequencing (BES) of 28,000 clones produced 28.1 Mb of genomic sequence that allowed the identification of the repetitive fraction (12.5% of the genome) and estimation of gene content (31,000 genes) of this species. BES analyses identified 3,800 SSRs and 6,617 putative SNPs. Comparative genomic studies showed that citrus gene homology and microsyntheny with Populus trichocarpa was rather higher than with Arabidopsis thaliana, a species phylogenetically closer to citrus.
In this work, we report the characterization of three BAC libraries from C. clementina, and a new set of genomic resources that may be useful for isolation of genes underlying economically important traits, physical mapping and eventually crop improvement in Citrus species. In addition, BAC end sequencing has provided a first insight on the basic structure and organization of the citrus genome and has yielded valuable molecular markers for genetic mapping and cloning of genes of agricultural interest. Paired end sequences also may be very helpful for whole-genome sequencing programs.
- Citrus Species
- Putative SNPs
- Pinot Noir
- Citrus Genome
- Putative Code Region
Citrus, one of the major fruit tree crops is widely cultivated throughout the globe and therefore has a tremendous economical, social and cultural impact in our society. Citrus improvement through traditional techniques, however, is highly impaired due to the unusual combination of biological characteristics of Citrus species, their low genetic diversity and the long-term nature of tree breeding. Citrus are diploid plants with an estimated haploid genome size of about 367 Mb and 2n = 18 chromosomes, which may facilitate the use of genomics tools for crop improvement. Expressed sequence tag (EST) analyses and molecular marker studies strongly suggest that the main commercial citrus cultivars (oranges, lemons and grapefruits) are mostly interspecific hybrids and therefore are heterozygous "species" . In addition, most of the cultivars in these groups, including Clementine varieties, may represent accumulated somatic mutations identified over centuries .
The development of citrus genomic resources is in its infancy although in recent years major efforts and goals mostly on functional genomics have certainly been undertaken . Critical functional and expression analyses through microarrays with several platforms have also been published and analyses of ESTs in public databases have been initiated [4, 5]. For instance, 401,692 citrus ESTs have been deposited at GenBank and are currently available. This collection constitutes a valuable source for the direct access to the genes of interest and for the development of molecular markers for map-based cloning purposes or marker-assisted selection programs [6–9]. Moreover, genetic linkage maps have been produced with increasing value and resolution, following the evolution of new marker systems [10, 11]. Genetic transformation in citrus is also available  and strategies based on genome-wide mutagenesis are being explored. Other innovative resources such as viral-induced gene silencing (VIGS) are being developed and work in citrus proteomics is in progress [13, 14]. Thus, current advances in citrus research include the rapid development of functional genomics and molecular biology resources  although, on the other hand, basic information on the organization and structure of citrus genome is lacking. The main challenge for a comprehensive and meaningful description of genomes is the integration of the DNA marker-based genetic maps with physical maps, and eventually with DNA sequence of the whole genome, the ultimate physical map. For the generation of high-resolution physical maps, the construction of BAC libraries containing clones with large DNA fragments appears to be indispensable. BAC end sequencing is indeed an important component of physical map development and can be considered a form of low coverage sequencing . Paired end sequences of BACs form an important part of scaffolding whole-genome shotgun programs  as well as in BAC based genome projects [18, 19]. In addition, BAC clone collections and BAC-based contig maps are powerful tools having multiple applications in genomics including positional cloning. The BAC end sequence provides a random survey of the information contents (genes, transposons, repeats) of unsequenced genomes [20–22], and yields molecular markers useful for genetic mapping [23–25], and cloning of genes of agricultural interest [26–28]. Furthermore, in many agriculturally important species BAC clones and physical maps are being rapidly developed since they are essential components in linking phenotypic traits to the responsible genetic variation, to integrate the genetic data, for the comparative analysis of genomes, and to speed up marker-assisted selection (MAS) for breeding. Thus, BAC libraries have become central for physical mapping, genome analysis, clone based sequencing and sequencing of complex genomes, for both model [18, 19] and main crop plants [29–32].
In citrus, two BAC libraries from Poncirus trifoliata  and a hybrid of Citrus × Poncirus  have been described in detail. These libraries were constructed as part of a map-based cloning strategy of genes conferring resistance to citrus tristeza virus that causes significant economic damage and losses to citrus worldwide. Poncirus is a non domesticated genus related to citrus species that produces inedible fruit. However, other efforts to generate BAC citrus resources, for instance in Satsuma or sweet orange have also been accomplished .
In this work, we report the characterization of three genomic BAC libraries from Clementine (Citrus clementina) mandarin, a cultivar of great economic importance that has been a main target of recent studies . The development of these new genomic tools also complements the functional genomics platform generated for this species including an extensive EST collection and a 20 k cDNA chip [4, 5], expanding further possibilities for isolation of genes of agronomical interest [33, 34]. This study provides the most comprehensive, large insert clone resource of any Citrus species and reports the analysis of 46,339 BAC end sequences offering a first detailed insight into the sequence composition of the Clementine genome. The analyses focused on protein coding regions, repeat element composition, microsatellite and single nucleotide polymorphism contents. Additionally, data on gene homology based comparative genomics with poplar and Arabidopsis are also presented. The annotated BAC-end sequences may well serve as useful resources for physical mapping, positional cloning, genetic marker development and genome sequencing of C. clementina.
BAC library characterization
Genomic Citrus clementina BAC Libraries
Partial digest enzyme
Average insert size
N° of clones
BAC end sequencing
Citrus clementina BAC end sequencing summary
BES mate pairsa
Total raw sequence (bp)
Average read length (bp)
Chloro/mito total sequence (bp)
Final genomic sequence (bp)d
Chloroplast and mitochondrion DNA analysis
In order to identify extranuclear sequences, BES were first compared against the Citrus sinensis chloroplastic  and the Arabidopsis thaliana  mitochondrial genome sequences with an stringent threshold of 1e-15. The comparison indicated that 736 (1.68%) and 46 (0.1%) sequences produced significant matches with the chloroplastic and mitochondrial genomes, respectively. Chloroplastic BESs were assembled and the consensus sequences obtained spanned 101,470 bp, approximately 70% of the chloroplast genome. Chloroplastic and mitochondrial DNA summed up to 480 kb, and therefore, the total genomic DNA obtained was 28,1 Mb. Considering the average size of the BAC clones, the coverage provided by the BAC end sequencing was higher than 8.4 genomic equivalents.
Repetitive DNA Analysis
Repetitive DNA in Citrus clementina BESs identified by Repeat Masker
Number of elements
Number of elements (%)
Total DNA transposons
Total interspersed repeats
Analysis of coding regions
After filtration of mitochondrial, chloroplastic and repetitive sequences, the remaining 30,787 BES were analyzed for coding region identification via homology search. Parallel searches with BLASTX and BLASTN were performed against the non-redundant database (e value cut off of 1e-4) and a database of Citrus ESTs from GenBank (e value cut off of 1e-15), respectively. The BLASTX search identified 14,030 BES (36% of the total BESs) with significant protein hits, while the BLASTN search revealed a similar number of reads, i.e. 14,023 reads that rendered significant homology with 40,536 citrus ESTs. Overall, 20,185 BESs produced BLASTX and/or BLASN hits, and the 3 BAC libraries rendered a similar number of clones carrying potential coding regions (see Additional File 1). The total number of BAC clones that produced protein and/or EST hits was 15.658, while 4,527 of them provided hits in both 5' and 3' ends, an observation that suggested high gene contents in these BAC clones and indeed their possible location at the euchromatin. It was also found that 7,868 BES produced hits for both proteins and ESTs, strongly indicating that they contained active transcription units.
Assuming that each significant BLAST hit corresponds to a different transcription unit, the total number of hypothetical genes described in the BAC ends was 20,185. Thus, considering the amount of analyzed sequence (28.1 Mb), the estimated genome size (367 Mb), and the number of genomic equivalents analyzed (8.4), the gene contents of the genome of Citrus clementina was assessed as 31,000. This gene number is comparable to the estimate reported for three species of similar genome size that have been sequenced to completion to date: rice (Oryza sativa), with a genome size of 430 Mb and 41,042 genes identified , Black cottonwood with 55,000 genes in 485 Mb , and grapevine (Vitis vinifera), with 30,434 genes in 487 Mb . Other estimates based on BES analysis obtained for plant species were also analogous. For instance, in papaya, the estimated gene contents was 35,526 with a 372 Mb genome size , and in Chinese cabbage, with a genome of 529 Mb, 43,000 genes were calculated .
Lastly, Blast2GO  was used to analyze the different functions associated with the putative coding regions, and GO terms were assigned to 10,598 sequences. Additional file 2 shows the results obtained for the gene ontology categories Molecular Function and Process. This classification may be useful to identify and locate genes of agronomic interest, such as those related to sugar and cellulose synthesis [GenBank:ET090832, GenBank:ET070671, GenBank:ET074358], ion transport [GenBank:ET110643, GenBank:ET110643, GenBank:ET074792], or calcium metabolism [GenBank:ET068865, GenBank:ET091143, GenBank:ET084632], for instance.
BAC end sequences have proved to be excellent sources to identify simple sequence repeats (SSRs or microsatellites) in many plant species such as cotton or soybean [23, 25, 49]. In this work, Sputnik  was used to identify a total of 3,814 SSRs longer than 15 bp in the BESs that did not carry repetitive sequences. The occurrence of SSRs in the Clementine genome had a frequency of 0.20 SSR per kb, a value almost identical to that reported in a study based on citrus ESTs (0.19 SSR per kb)  and in B. rapa (0.18 SSR per kb) . This frequency, however, was lower than the ratios found in grapevine (0.48 SSR per kb) , papaya (0.43 SSRs longer than 12 bp per kb) , and A. thaliana (0.33 SSR per kb).
Microsatellites are co-dominant, highly polymorphic, and simple to use markers that have been successfully used in studies of genetic diversity  and genetic mapping in citrus [8, 52] and in many other plants [25, 53–55]. The additional markers reported here will certainly contribute to improve the coverage of the Citrus genome for many purposes, including the development of accurate linkage and genetic maps.
Contig Assembly and SNP analysis
Assembly of BESs that did not contain repetitive sequences was performed with CAP3 , and a total of 6,461 contigs including 19,057 reads and covering 6.14 Mb of sequence were produced. It has been suggested that C. clementina is an offspring of a C. sinensis × C. reticulata cross and therefore has a heterozygous genome . In order to identify possible polymorphisms affecting the assembly of the readings and the construction of the physical map of this species, a BLASTN search of all the contigs and singlets was performed against themselves. One hundred thirty sequence pairs that presented a single reciprocal BLASTN hit and the same protein hit from the previous BLASTX search (identical accession number) were selected for further analysis. A total of 81 pairs displayed sequence identities higher than 90% and were considered as originated from the same genomic region. These sequences were not assembled in the same contig due to the presence of polymorphisms, similarly to what was reported in the sequencing of the highly heterozygous grapevine variety, Pinot Noir .
The contigs generated in the assembly were analyzed in order to identify single nucleotide polymorphisms (SNPs) in the diploid genome of C. Clementina. SNPs are the most abundant and powerful polymorphic markers, since they provide gene-based markers that can be used in the creation of dense genetic linkage maps  and, more important, in the identification of genes associated with specific trait loci. BAC end sequences of heterozygous genotypes have also been successfully utilized in SNP discovery and construction of linkage maps [59–61]. PolyBayes , a software designed to use genomic sequence as a template and base quality values to discern true allelic variations from sequencing errors, was used to reveal SNP polymorphisms. The number of putative SNPs identified in the Clementine sequences that showed P ≥ 0.9 and SNP depth lower than 10 was 6,617, corresponding to 1.08 SNPs per kb. PolyBayes has been successfully utilized in automated high throughput identification of SNPs from EST collections, and it has been shown that when using P ≥ 0.95, 85% of the predicted SNPs are generally validated experimentally .
SNPs in Citrus clementina BESs identified by PolyBayes
Total seq (kb)c
Detailed relation of SNPs
It should be noted that the SNP listing reported in this work constitutes the first set of putative SNPs identified in any Citrus species, and hence provides a completely new resource for genome analysis in this genera. It is also worth mentioning that 4,500 SNPs were located in or close to putative coding regions, and therefore these 'functional SNPs' may provide an inestimable resource for the identification of genes associated with specific trait loci in addition to their utility as molecular markers for genetic and comparative mapping, nucleotide diversity analysis and association studies.
Two heterozygous genomes have been sequenced to completion, the Nisqually-1 poplar and the Pinot Noir grapevine strains. In both cases the frequency of polymorphisms found within these heterozygous genomes, was 2.6  and 4  SNPs per kb, respectively. These rates are in contrast with the frequency of 1.08 SNPs per kb found in Clementine. The reproductive biology of the different species could explain these differences. Gametophytic self- and cross-incompatibility, and apomixy would produce low variability within Citrus species , while outcrossing by means of insect and wind pollination, which is the norm for poplar and vitis, would result in highly heterozygous cultivars [17, 57].
Genomic BLASTN results with non-repetitive Citrus clementina BESs
Total length (kb)a
Average e value
Average distance between hits (kb)
Average n° of hits per chromosome
Mapping of citrus BESs hits on poplar chromosomes
Total Length (bp)
n° BES mapped
Average distance between BES
Following the approach of Lai et al. , we used forward and reverse BES read pairs separated by the approximate length of BAC clone inserts (~120 kb), to analyze the microsynteny between C. clementina and Arabidopsis, rice and poplar. To be considered as potentially collinear with the target genome, the citrus mate pairs had to map in the heterologous genome into a region comprised between 10 and 300 kb and be also oriented properly. The analyses of the sequences identified 108 Clementine BAC end pairs that met these criteria in poplar, while no one was found in Arabidopsis or rice. Furthermore, the majority of these BES pairs mapped on the Populus genome at a distance similar to the insert size of the Clementine libraries, suggesting the microsynteny between citrus and poplar is higher than between citrus and Arabidopsis. These results are striking since C. clementina and A. thaliana belong to Sapindales and Brassicales orders (eurosids II clade) that probably split approximately 87 MYA, while P. trichocarpa belongs to Malpighiales order, (eurosids I clade) that diverged from eurosids II around 109 MYA . Moreover, similar results were obtained by Lai et al  with papaya, that also exhibited higher level of colinearity with the poplar than with the Arabidopsis genome despite that C. papaya is a basal member of the Brassicales. Although a definitive explanation has not been provided yet, it is currently believed that the genome of A. thaliana has undergone a recent whole genome duplication, followed by subsequent gene losses and extensive local gene duplications , which might be responsible of the lack of colinearity with other eurosid II species. Comparative genomics with the recently sequenced genome of the grapevine, provides additional evidence that the genome of Arabidopsis has been thoroughly rearranged as related to an ancient angiosperm genome . The fact that papaya, Clementine, grapevine and poplar are long lived, clonally propagated, woody plants, might apparently cause a deceleration of their molecular clocks, resulting in genomes with higher resemblance to the ancestral eurosids genome .
We report here the construction of three genomic BAC libraries of Citrus clementina, with three restriction enzymes (Eco RI, Mbo I and Hind III). The number of picked BAC clones (56,000) and the average length of the inserts provide coverage of 19.5 haploid genome equivalents, ensuring a wide representation of the genome of this species. These libraries are adequate for the construction of the physical map of C. clementina. The analysis of 28.1 Mb of genomic sequence produced by BAC end sequencing has provided a first insight of the genome organization of C. clementina. The repetitive fraction of this species corresponding with transposable elements comprised 12.5% of the genome, while the gene number was estimated to be 31,000. This work also describes a set of 3,814 SSRs and a collection of the first 6,617 putative SNPs described in citrus that may be very useful for positional cloning, genetic and comparative mapping, nucleotide diversity analysis, and association studies. Finally, comparative genomics through gene homology searches has shown that, in spite of their taxonomic classification, microsynteny between Citrus and Populus is higher than with Arabidopsis that is a phylogenetically closer species.
Citrus clementina (Clementine mandarin, var. clemenules) developing leaves were used for BAC library construction.
BAC library construction
BAC libraries were constructed from high molecular weight (HMW) genomic DNA processed at Amplicon Express, (Pullman, Washington) using the method described in . DNA digestion was performed with varying amounts of Mbo I, Eco RI, and Hind III to identify appropriate partial digestion conditions. pECBAC1 and contained two FRT and one oriV elements, thus resulting in the pBAC(FRT-oriV) vector . Ligations were transformed into DH10B Escherichia coli cells (Invitrogen) and plated on LB agar with chloramphenicol (30 μg/ml), X-gal (20 mg/ml) and IPTG (0.1 M). Clones were robotically picked with a Genomic Solution G3 into 384 well plates containing LB freezing media. Plates were incubated for 18 h, replicated and then frozen at -80°C. The replicated copy was used for BAC end sequencing.
Insert size estimation
To estimate insert sizes, 10 μl aliquots of BAC miniprep DNA were digested with 5 U of Not I enzyme for three h at 37°C. The digestion products were separated by pulsed-field Weld gel electrophoresis (CHEF-DRIII system, Bio-Rad) in a 1% agarose gel in TBE buffer 0,5×. Insert sizes were compared to those of the Lambda Ladder PFG Marker (New England Biolabs). Electrophoresis was carried out for 18 h at 14°C with an initial switch time of 5 s, a Wnal switch time of 15 s, in a voltage gradient of 6 V/cm.
BAC End Sequencing
BAC clones were inoculated into 96-deep well macroplates and grown for 20 hs at 37°C. Cells were harvested by centrifugation and BACs were purified in 96-well plates by a standard alkaline lysis protocol developed by Genoscope (Paris, France). BAC DNA was precipitated with isopropanol and washed with 70% ethanol. Sequencing was carried out on ABI3730 equipment with "Dye Terminator" process using ABI kit version 3.1. in the Genoscope facility.
The software phred was used for base calling, and Crossmatch for vector masking. Repetitive DNA was identified with the RepeatMasker software , using the viridiplantae section of the RepBase Update  as database. Assembly was performed with CAP3 , using read quality and default parameters. Similarity searches were performed with the standalone version of BLAST [51, 70], against the NCBI non redundant protein, nucleotide and EST databases available on November 2007 . Parsing of the BLAST results was performed with the Bio::SearchIO module from the Bioperl package . Coding sequences were annotated with GO terms using Blast2GO . SPUTNIK  was used to identify simple sequence repeats (SSRs), and POLYBAYES  to search for SNPs.
DNA extraction was done from leaf tissues of C. clementina cv Nules using the DNeasy® Plant Mini Kit (Qiagen).
PCR amplifications of the samples were performed using a Mastercycler epgradiend S thermocycler (Eppendorf) in 100 μL final volume containing 0.025 U/μL of Pfu DNA polymerase (Fermentas), 0.2 ng/μL of genomic DNA, 0.2 mM of each dNTP, 2 mM MgSO4, 75 mM Tris-HCl (pH 8.8), 20 mM (NH4)2SO4, 0.2 μM of each primers. The following PCR program was applied: denaturation at 94° C for 5 min and 35 repeats of the following cycle: 30 s at 94°C, 1 min at 55°C or 60°C (according to primers Tm), 45 s at 72°C; and final elongation step of 4 min at 72°C.
PCR product purification was done directly or after cutting single bands on agarose gel, using respectively QIAquick® PCR Purification Kit and QIAquick® Gel Extraction Kit (Qiagen)
Sequencing was carried out on ABI3730 equipment with "Dye Terminator" process using ABI kit version 3.1. SNPs were identified as double peaks by manual inspection of chromatograms, and subsequently validated
Work at Centro de Genómica was supported by INIA grant RTA04-013, INCO contract 015453 and Ministerio de Educación y Ciencia grant AGL2007-65437-C04-01/AGR.
- Nicolosi E, Deng ZN, Gentile A, La Malfa S, Continella G, Tribulato E: Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theor Appl Genet. 2000, 100: 1155-1166.View ArticleGoogle Scholar
- Gmitter FG: Origin, evolution and breeding of the grapefruit. Plant Breeding Reviews. 1995, 13: 345-363.Google Scholar
- Talon M, Gmitter FG: Citrus Genomics. Int J Plant Genomics. 2008, 2008: 528361-PubMedPubMed CentralView ArticleGoogle Scholar
- Terol J, Conesa A, Colmenero JM, Cercos M, Tadeo F, Agusti J: Analysis of 13000 unique Citrus clusters associated with fruit quality, production and salinity tolerance. BMC Genomics. 2007, 8: 31-PubMedPubMed CentralView ArticleGoogle Scholar
- Forment J, Gadea J, Huerta L, Abizanda L, Agusti J, Alamar S: Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies. Plant Mol Biol. 2005, 57: 375-391.PubMedView ArticleGoogle Scholar
- Chen C, Zhou P, Choi YA, Huang S, Gmitter FG: Mining and characterizing microsatellites from citrus ESTs. Theor Appl Genet. 2006, 112: 1248-1257.PubMedView ArticleGoogle Scholar
- Barkley NA, Roose ML, Krueger RR, Federici CT: Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs). Theor Appl Genet. 2006, 112: 1519-1531.PubMedView ArticleGoogle Scholar
- Jiang D, Zhong GY, Hong QB: Analysis of microsatellites in citrus unigenes. Yi Chuan Xue Bao. 2006, 33: 345-353.PubMedGoogle Scholar
- Ruiz C, Asins MJ: Comparison between Poncirus and Citrus genetic linkage maps. Theor Appl Genet. 2003, 106: 826-836.PubMedGoogle Scholar
- Deng Z, Tao Q, Chang L, Huang S, Ling P, Yu C: Construction of a bacterial artificial chromosome (BAC) library for citrus and identification of BAC contigs containing resistance gene candidates. Theor Appl Genet. 2001, 102: 1177-1184.View ArticleGoogle Scholar
- Yang ZN, Ye XR, Choi S, Molina J, Moonan F, Wing RA: Construction of a 1.2-Mb contig including the citrus tristeza virus resistance gene locus using a bacterial artificial chromosome library of Poncirus trifoliata (L.) Raf. Genome. 2001, 44: 382-393.PubMedView ArticleGoogle Scholar
- Fagoaga C, Tadeo FR, Iglesias DJ, Huerta L, Lliso I, Vidal AM: Engineering of gibberellin levels in citrus by sense and antisense overexpression of a GA 20-oxidase gene modifies plant architecture. J Exp Bot. 2007, 58: 1407-1420.PubMedView ArticleGoogle Scholar
- Katz E, Fon M, Lee YJ, Phinney BS, Sadka A, Blumwald E: The citrus fruit proteome: insights into citrus fruit metabolism. Planta. 2007, 226: 989-1005.PubMedView ArticleGoogle Scholar
- Lliso I, Tadeo FR, Phinney BS, Wilkerson CG, Talon M: Protein changes in the albedo of citrus fruits on postharvesting storage. J Agric Food Chem. 2007, 55: 9047-9053.PubMedView ArticleGoogle Scholar
- Tadeo F, Cercos M, Colmenero-Flores JM, Iglesias DJ, Naranjo MA, Rios G: Molecular physiology of development and quality of citrus. Advances in Botanical Research. 2008,Google Scholar
- Town CD: Large-scale DNA sequencing. The handbook of Plant Genome Mapping. Edited by: Meksem K, Kahl G. 2005, Wiley-VCH, 337-351.View ArticleGoogle Scholar
- Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U: The Genome of Black Cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313: 1596-1604.PubMedView ArticleGoogle Scholar
- The Arabidopsis Genome Initiative A: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815.View ArticleGoogle Scholar
- International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800.View ArticleGoogle Scholar
- Lai C, Yu Q, Hou S, Skelton R, Jones M, Lewis K: Analysis of papaya BAC end sequences reveals first insights into the organization of a fruit tree genome. Mol Genet Genomics. 2006, 276: 1-12.PubMedView ArticleGoogle Scholar
- Cheung F, Town CD: A BAC end view of the Musa Accuminata genome. BMC Plant Biol. 2007, 7: 29-PubMedPubMed CentralView ArticleGoogle Scholar
- Hong CP, Plaha P, Koo DH, Yang TJ, Choi SR, Lee YK: A Survey of the Brassica rapa Genome by BAC-End Sequence Analysis and Comparison with Arabidopsis thaliana. Mol Cells. 2006, 22: 300-307.PubMedGoogle Scholar
- Frelichowski JE, Palmer MB, Main D, Tomkins JP, Cantrell RG, Stelly DM: Cotton genome mapping with new microsatellites from Acala 'Maxxa' BAC-ends. Mol Genet Genomics. 2006, 275: 479-491.PubMedView ArticleGoogle Scholar
- Marek LF, Mudge J, Darnielle L, Grant D, Hanson N, Paz M: Soybean genomic survey: BAC-end sequences near RFLP and SSR markers. Genome. 2001, 44: 572-581.PubMedView ArticleGoogle Scholar
- Shultz JL, Samreen K, Rabia B, Jawaad AA, Lightfoot DA: The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean. Theor Appl Genet. 2007, 114: 1081-1090.PubMedView ArticleGoogle Scholar
- Coyne CJ, McClendon MT, Walling JG, Timmerman-Vaughan GM, Murray S, Meksem K: Construction and characterization of two bacterial artificial chromosome libraries of pea (Pisum sativum L.) for the isolation of economically important genes. Genome. 2007, 50: 871-875.PubMedView ArticleGoogle Scholar
- Nam W, Penmetsa RV, Endre G, Uribe P, Kim D, Cook DR: Construction of a bacterial artificial chromosome library of Medicago truncatula and identification of clones containing ethylene-response genes. Theor Appl Genet. 1999, 98: 638-646.View ArticleGoogle Scholar
- Liang H, Fang E, Tomkins J, Luo M, Kudrna D, Kim H: Development of a BAC library for yellow-poplar (Liriodendron tulipifera) and the identification of genes associated with flower development and lignin biosynthesis. Tree Genetics & Genomes. 2007, 3: 215-225.View ArticleGoogle Scholar
- Han Y, Gasic K, Marron B, Beever JE, Korban SS: A BAC-based physical map of the apple genome. Genomics. 2007, 89: 630-637.PubMedView ArticleGoogle Scholar
- Wu C, Sun S, Nimmakayala P, Santos FA, Meksem K, Springman R: A BAC- and BIBAC-Based Physical Map of the Soybean Genome. Genome Res. 2004 Feb;14(2):319-26. 2004, 14 (2): 319-326.Google Scholar
- Yim YS, Davis GL, Duru NA, Musket TA, Linton EW, Messing JW: Characterization of Three Maize Bacterial Artificial Chromosome Libraries toward Anchoring of the Physical Map to the Genetic Map Using High-Density Bacterial Artificial Chromosome Filter Hybridization. Plant Physiol. 2002, 130: 1686-1696.PubMedPubMed CentralView ArticleGoogle Scholar
- Cenci A, Chantret N, Kong X, Gu Y, Anderson OD, Fahima T: Construction and characterization of a half million clone BAC library of durum wheat (Triticum turgidum ssp. durum). Theor Appl Genet. 2003, 107: 931-939.PubMedView ArticleGoogle Scholar
- Cercos M, Soler G, Iglesias D, Gadea J, Forment J, Talon M: Global Analysis of Gene Expression During Development and Ripening of Citrus Fruit Flesh. A Proposed Mechanism for Citric Acid Utilization. Plant Mol Biol. 2006, 62: 513-527.PubMedView ArticleGoogle Scholar
- Alos E, Cercos M, Rodrigo MJ, Zacarias L, Talon M: Regulation of color break in citrus fruits. Changes in pigment profiling and gene expression induced by gibberellins and nitrate, two ripening retardants. J Agric Food Chem. 2006, 54: 4888-4895.PubMedView ArticleGoogle Scholar
- Sun S, Xu Z, Wu C, Ding K, Zhang HB: Genome properties and their influences on library construction and physical mapping. Proceedings of Plant and Animal Genome XI Conference. 2003, P77-Google Scholar
- Ren C, Xu Z, Sun S, Lee MK, Wu C, Scheuring C: Genomic DNA Libraries and Physical Mapping. The handbook of plant genome mapping. Edited by: Meksem K, Kahl G. 2005, Wiley-VCH, 173-213.View ArticleGoogle Scholar
- Chang YL, Tao Q, Scheuring C, Ding K, Meksem K, Zhang HB: An integrated map of Arabidopsis thaliana for functional analysis of its genome sequence. Genetics. 2001, 159: 1231-1242.PubMedPubMed CentralGoogle Scholar
- Chen M, Presting G, Barbazuk WB, Goicoechea JL, Blackmon B, Fang G: An integrated physical and genetic map of the rice genome. Plant Cell. 2002, 14: 537-545.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu Z, Sun S, Covaleda L, Ding K, Zhang A, Wu C: Genome physical mapping with large-insert bacterial clones by fingerprint analysis: methodologies, source clone genome coverage, and contig map quality. Genomics. 2004, 84: 941-951.PubMedView ArticleGoogle Scholar
- Bausher M, Singh N, Lee SB, Jansen R, Daniell H: The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var 'Ridge Pineapple': organization and phylogenetic relationships to other angiosperms. BMC Plant Biology. 2006, 6: 21-PubMedPubMed CentralView ArticleGoogle Scholar
- Unseld M, Marienfeld JR, Brandt P, Brennicke A: The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides. Nat Genet. 1997, 15: 57-61.PubMedView ArticleGoogle Scholar
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467.PubMedView ArticleGoogle Scholar
- Rico-Cabanas L, Martinez-Izquierdo JA: CIRE1, a novel transcriptionally active Ty1-copia retrotransposon from Citrus sinensis. Mol Genet Genomics. 2007, 277: 365-377.PubMedView ArticleGoogle Scholar
- Bernet GP, Asins MJ: Identification and genomic distribution of gypsy like retrotransposons in Citrus and Poncirus. Theor Appl Genet. 2003, 108: 121-130.PubMedView ArticleGoogle Scholar
- Yuan Q, Ouyang S, Liu J, Suh B, Cheung F, Sultana R: The TIGR rice genome annotation resource: annotating the rice genome and creating resources for plant biologists. Nucl Acids Res. 2003, 31: 229-233.PubMedPubMed CentralView ArticleGoogle Scholar
- Bogunic F, Muratovic E, Brown SC, Siljak-Yakovlev S: Genome size and base composition of five Pinus species from the Balkan region. Plant Cell Reports. 2003, 22: 59-63.PubMedView ArticleGoogle Scholar
- Carels N, Hatey P, Jabbari K, Bernardi G: Compositional Properties of Homologous Coding Sequences from Plants. J Mol Evol. 1998, 46: 45-53.PubMedView ArticleGoogle Scholar
- Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676.PubMedView ArticleGoogle Scholar
- Clement D, Lanaud C, Sabau X, Fouet O, Cunff LL, Ruiz E: Creation of BAC genomic resources for cocoa (Theobroma cacao L.) for physical mapping of RGA containing BAC clones. Theor Appl Genet. 2004, 108: 1627-1634.PubMedView ArticleGoogle Scholar
- Abajian C: SPUTNIK. Computer Program. 1994Google Scholar
- Chen C, Zhou P, Choi YA, Huang S, Gmitter FG: Mining and characterizing microsatellites from citrus ESTs. Theor Appl Genet. 2006, 112 (7): 1248-1257.PubMedView ArticleGoogle Scholar
- Kijas JMH, Thomas MR, Fowler JCS, Roose ML: Integration of trinucleotide microsatellites into a linkage map of Citrus. Theor Appl Genet. 1997, 94: 701-706.View ArticleGoogle Scholar
- Varshney RK, Grosse I, H+ U, Siefken R, Prasad M, Stein N: Genetic mapping and BAC assignment of EST-derived SSR markers shows non-uniform distribution of genes in the barley genome. Theor Appl Genet. 2006, 113: 239-250.PubMedView ArticleGoogle Scholar
- Mun JH, Kim DJ, Choi HK, Gish J, Debelle F, Mudge J: Distribution of microsatellites in the genome of Medicago truncatula: a resource of genetic markers that integrate genetic and physical maps. Genetics. 2006, 172: 2541-2555.PubMedPubMed CentralView ArticleGoogle Scholar
- Feingold S, Lloyd J, Norero N, Bonierbale M, Lorenzen J: Mapping and characterization of new EST-derived microsatellites for potato (Solanum tuberosum L.). Theor Appl Genet. 2005, 111: 456-466.PubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877.PubMedPubMed CentralView ArticleGoogle Scholar
- Velasco R, Zharkikh A, Troggio M, Cartwright DA, Cestaro A, Pruss D: A High Quality Draft Consensus Sequence of the Genome of a Heterozygous Grapevine Variety. PLoS ONE. 2007, 2: e1326-PubMedPubMed CentralView ArticleGoogle Scholar
- Rafalski A: Applications of single nucleotide polymorphisms in crop genetics. Curr Opin Plant Biol. 2002, 5: 94-100.PubMedView ArticleGoogle Scholar
- Weil MM, Pershad R, Wang R, Zhao S: Use of BAC end sequences for SNP discovery. Methods Mol Biol. 2004, 256: 1-6.PubMedGoogle Scholar
- Abe K, Noguchi H, Tagawa K, Yuzuriha M, Toyoda A, Kojima T: Contribution of Asian mouse subspecies Mus musculus molossinus to genomic constitution of strain C57BL/6J, as defined by BAC-end sequence-SNP analysis. Genome Res. 2004, 14: 2439-2447.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamamoto K, Narukawa J, Kadono-Okuda K, Nohata J, Sasanuma M, Suetsugu Y: Construction of a Single Nucleotide Polymorphism Linkage Map for the Silkworm, Bombyx mori, Based on Bacterial Artificial Chromosome End Sequences. Genetics. 2006, 173: 151-161.PubMedPubMed CentralView ArticleGoogle Scholar
- Marth GT, Korf I, Yandell MD, Yeh RT, Gu Z, Zakeri H: A general approach to single-nucleotide polymorphism discovery. Nat Genet. 1999, 23: 452-456.PubMedView ArticleGoogle Scholar
- Pavy N, Parsons LS, Paule C, MacKay J, Bousquet J: Automated SNP detection from a large collection of white spruce expressed sequences: contributing factors and approaches for the categorization of SNPs. BMC Genomics. 2006, 7 (174): 174-PubMedPubMed CentralView ArticleGoogle Scholar
- Wikström N, Savolainen V, Chase MW: Evolution of the angiosperms: calibrating the family tree. Proceedings of the Royal Society B: Biological Sciences. 2001, 268: 2211-2220.PubMedPubMed CentralView ArticleGoogle Scholar
- Tao Q, Wang A, Zhang HB: One large-insert plant-transformation-competent BIBAC library and three BAC libraries of Japonica rice for genome research in rice and other grasses. Theor Appl Genet. 2002, 105: 1058-1066.PubMedView ArticleGoogle Scholar
- Frijters CJ, Zhang Z, Damme Mv, Wang L, Ronald PC, Michelmore RW: Construction of a bacterial artificial chromosome library containing large Eco RI and Hin dIII genomic fragments of lettuce. Theor Appl Genet. 1997, 94: 390-399.View ArticleGoogle Scholar
- Ewing B, Green P: Base-Calling of Automated Sequencer Traces Using Phred. II Error Probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. Computer Program. 1996Google Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- National Center for Biotechnology Information. Electronic Citation.Google Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12: 1611-1618.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.