Characterization of a deep-coverage carrot (Daucus carota L.) BAC library and initial analysis of BAC-end sequences

Impact of Genomic and Transcriptomic Resources on Apiaceae Crop Breeding Strategies

The Apiaceae taxon is one of the most important families of flowering plants and includes thousands of species used for food, flavoring, fragrance, medical and industrial purposes. This study had the specific intent of reviewing the main genomics and transcriptomic data available for this family and their use for the constitution of new varieties. This was achieved starting from the description of the main reproductive systems and barriers, with particular reference to cytoplasmic (CMS) and nuclear (NMS) male sterility. We found that CMS and NMS systems have been discovered and successfully exploited for the development of varieties only in Foeniculum vulgare, Daucus carota, Apium graveolens and Pastinaca sativa; whereas, strategies to limit self-pollination have been poorly considered. Since the constitution of new varieties benefits from the synergistic use of marker-assisted breeding in combination with conventional breeding schemes, we also analyzed and discussed the available SNP and SSR marker datasets (20 species) and genomes (8 species). Furthermore, the RNA-seq studies aimed at elucidating key pathways in stress tolerance or biosynthesis of the metabolites of interest were limited and proportional to the economic weight of each species. Finally, by aligning 53 plastid genomes from as many species as possible, we demonstrated the precision offered by the super barcoding approach to reconstruct the phylogenetic relationships of Apiaceae species. Overall, despite the impressive size of this family, we documented an evident lack of molecular data, especially because genomic and transcriptomic resources are circumscribed to a small number of species. We believe that our contribution can help future studies aimed at developing molecular tools for boosting breeding programs in crop plants of the Apiaceae family.

Using carrot centromeric repeats to study karyotype relationships in the genus Daucus (Apiaceae)

Background

In the course of evolution, chromosomes undergo evolutionary changes; thus, karyotypes may differ considerably among groups of organisms, even within closely related taxa. The genus Daucus seems to be a promising model for exploring the dynamics of karyotype evolution. It comprises some 40 wild species and the cultivated carrot, a crop of great economic significance. However, Daucus species are very diverse morphologically and genetically, and despite extensive research, the taxonomic and phylogenetic relationships between them have still not been fully resolved. Although several molecular cytogenetic studies have been conducted to investigate the chromosomal structure and karyotype evolution of carrot and other Daucus species, detailed karyomorphological research has been limited to carrot and only a few wild species. Therefore, to better understand the karyotype relationships within Daucus, we (1) explored the chromosomal distribution of carrot centromeric repeats (CentDc) in 34 accessions of Daucus and related species by means of fluorescence in situ hybridization (FISH) and (2) performed detailed karyomorphological analysis in 16 of them.

Results

We determined the genomic organization of CentDc in 26 accessions of Daucus (belonging to both Daucus I and II subclades) and one accession of closely related species. The CentDc repeats were present in the centromeric regions of all chromosomes of 20 accessions (representing 11 taxa). In the other Daucus taxa, the number of chromosome pairs with CentDc signals varied depending on the species, yet their centromeric localization was conserved. In addition, precise chromosome measurements performed in 16 accessions showed the inter- and intraspecific karyological relationships among them.

Conclusions

The presence of the CentDc repeats in the genomes of taxa belonging to both Daucus subclades and one outgroup species indicated the ancestral status of the repeat. The results of our study provide useful information for further evolutionary, cytotaxonomic, and phylogenetic research on the genus Daucus and may contribute to a better understanding of the dynamic evolution of centromeric satellites in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07853-2.

Insights into genetic diversity and population structure of Indian carrot (Daucus carota L.) accessions

Carrot (Daucus carota L.) is acknowledged as a highly valuable vegetable crop. Despite having high demand, limited breeding efforts have been made to develop the varieties and hybrids suitable to wider climatic conditions due to improper characterization of the available germplasm. An accession panel (AP) consisting of 144 accessions of five different root colors representing Asiatic and Western gene pools collected from different parts of India was utilized in the present study. This diverse AP was used to assess the population structure and genetic diversity from 80 polymorphic DNA markers distributed throughout the genome. Population structure, neighbor-joining (NJ) tree, and principal coordinate analysis (PCoA)-based diversity assessment divided the AP into three subpopulations/clusters. Greater than ninety percent polymorphism and the higher average polymorphic information content (͂> 0.50) coupled with higher gene diversity (He) indicating the broad genetic base of the population. Moderate to high F_st and gene flow (Nm) between the subpopulations revealed a moderate genetic differentiation among Indian carrot accessions owing to the highly outcrossing nature of carrot. Analysis of molecular variance (AMOVA) exhibited higher variation among individuals within the subpopulations (69.00%) or total populations (19.00%) than among the subpopulations (13%) as expected in the single Daucus species used here. The information obtained in the study would benefit the carrot breeders to explore the genetic diversity of the Indian carrots in the carrot breeding program for widening the genetic base and multi-color target trait improvement.

Dof transcription factors in carrot: genome-wide analysis and their response to abiotic stress

OBJECTIVES

The DNA-binding one zinc finger (Dof) family transcription factors (TF) are involved in stress response. Dof TFs in carrot were identified and the responses of DcDof genes to abiotic stresses were analyzed.

RESULTS

46 DcDofs in carrot were identified from carrot genome database. Based on the conserved domain in Dof TF family of Arabidopsis thaliana, the DcDof TFs were divided into four classes, named class A, B, C and D. Carrot and Arabidopsis shared most motifs in the same subgroup. Real-time quantification PCR analysis showed tissue-specific expression patterns in DcDofs. DcDofs from eight subgroups responded to four abiotic stress treatments.

CONCLUSIONS

The expression profiles were different with the abiotic stresses changed, indicating complicated regulatory mechanisms in Dof TF family in higher plant, and the response mechanisms of Dof genes may be influenced by different plant species.

New insights into domestication of carrot from root transcriptome analyses

BACKGROUND

Understanding the molecular basis of domestication can provide insights into the processes of rapid evolution and crop improvement. Here we demonstrated the processes of carrot domestication and identified genes under selection based on transcriptome analyses.

RESULTS

The root transcriptomes of widely differing cultivated and wild carrots were sequenced. A method accounting for sequencing errors was introduced to optimize SNP (single nucleotide polymorphism) discovery. 11,369 SNPs were identified. Of these, 622 (out of 1000 tested SNPs) were validated and used to genotype a large set of cultivated carrot, wild carrot and other wild Daucus carota subspecies, primarily of European origin. Phylogenetic analysis indicated that eastern carrot may originate from Western Asia and western carrot may be selected from eastern carrot. Different wild D. carota subspecies may have contributed to the domestication of cultivated carrot. Genetic diversity was significantly reduced in western cultivars, probably through bottlenecks and selection. However, a high proportion of genetic diversity (more than 85% of the genetic diversity in wild populations) is currently retained in western cultivars. Model simulation indicated high and asymmetric gene flow from wild to cultivated carrots, spontaneously and/or by introgression breeding. Nevertheless, high genetic differentiation exists between cultivated and wild carrots (Fst = 0.295) showing the strong effects of selection. Expression patterns differed radically for some genes between cultivated and wild carrot roots which may be related to changes in root traits. The up-regulation of water-channel-protein gene expression in cultivars might be involved in changing water content and transport in roots. The activated expression of carotenoid-binding-protein genes in cultivars could be related to the high carotenoid accumulation in roots. The silencing of allergen-protein-like genes in cultivated carrot roots suggested strong human selection to reduce allergy. These results suggest that regulatory changes of gene expressions may have played a predominant role in domestication.

CONCLUSIONS

Western carrots may originate from eastern carrots. The reduction in genetic diversity in western cultivars due to domestication bottleneck/selection may have been offset by introgression from wild carrot. Differential gene expression patterns between cultivated and wild carrot roots may be a signature of strong selection for favorable cultivation traits.

Mapping genes governing flower architecture and pollen development in a double mutant population of carrot

A linkage map of carrot (Daucus carota L.) was developed in order to study reproductive traits. The F2 mapping population derived from an initial cross between a yellow leaf (yel) chlorophyll mutant and a compressed lamina (cola) mutant with unique flower defects of the sporophytic parts of male and female organs. The genetic map has a total length of 781 cM and included 285 loci. The length of the nine linkage groups (LGs) ranged between 65 and 145 cM. All LGs have been anchored to the reference map. The objective of this study was the generation of a well-saturated linkage map of D. carota. Mapping of the cola-locus associated with flower development and fertility was successfully demonstrated. Two MADS-box genes (DcMADS3, DcMADS5) with prominent roles in flowering and reproduction as well as three additional genes (DcAOX2a, DcAOX2b, DcCHS2) with further importance for male reproduction were assigned to different loci that did not co-segregate with the cola-locus.

Begin at the beginning: A BAC-end view of the passion fruit (Passiflora) genome

Background

The passion fruit (Passiflora edulis) is a tropical crop of economic importance both for juice production and consumption as fresh fruit. The juice is also used in concentrate blends that are consumed worldwide. However, very little is known about the genome of the species. Therefore, improving our understanding of passion fruit genomics is essential and to some degree a pre-requisite if its genetic resources are to be used more efficiently. In this study, we have constructed a large-insert BAC library and provided the first view on the structure and content of the passion fruit genome, using BAC-end sequence (BES) data as a major resource.

Results

The library consisted of 82,944 clones and its levels of organellar DNA were very low. The library represents six haploid genome equivalents, and the average insert size was 108 kb. To check its utility for gene isolation, successful macroarray screening experiments were carried out with probes complementary to eight Passiflora gene sequences available in public databases. BACs harbouring those genes were used in fluorescent in situ hybridizations and unique signals were detected for four BACs in three chromosomes (n = 9). Then, we explored 10,000 BES and we identified reads likely to contain repetitive mobile elements (19.6% of all BES), simple sequence repeats and putative proteins, and to estimate the GC content (~42%) of the reads. Around 9.6% of all BES were found to have high levels of similarity to plant genes and ontological terms were assigned to more than half of the sequences analysed (940). The vast majority of the top-hits made by our sequences were to Populus trichocarpa (24.8% of the total occurrences), Theobroma cacao (21.6%), Ricinus communis (14.3%), Vitis vinifera (6.5%) and Prunus persica (3.8%).

Conclusions

We generated the first large-insert library for a member of Passifloraceae. This BAC library provides a new resource for genetic and genomic studies, as well as it represents a valuable tool for future whole genome study. Remarkably, a number of BAC-end pair sequences could be mapped to intervals of the sequenced Arabidopsis thaliana, V. vinifera and P. trichocarpa chromosomes, and putative collinear microsyntenic regions were identified.

Electronic supplementary material

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

Characterization of centromeric histone H3 (CENH3) variants in cultivated and wild carrots (Daucus sp.)

In eukaryotes, centromeres are the assembly sites for the kinetochore, a multi-protein complex to which spindle microtubules are attached at mitosis and meiosis, thereby ensuring segregation of chromosomes during cell division. They are specified by incorporation of CENH3, a centromere specific histone H3 variant which replaces canonical histone H3 in the nucleosomes of functional centromeres. To lay a first foundation of a putative alternative haploidization strategy based on centromere-mediated genome elimination in cultivated carrots, in the presented research we aimed at the identification and cloning of functional CENH3 genes in Daucus carota and three distantly related wild species of genus Daucus varying in basic chromosome numbers. Based on mining the carrot transcriptome followed by a subsequent PCR-based cloning, homologous coding sequences for CENH3s of the four Daucus species were identified. The ORFs of the CENH3 variants were very similar, and an amino acid sequence length of 146 aa was found in three out of the four species. Comparison of Daucus CENH3 amino acid sequences with those of other plant CENH3s as well as their phylogenetic arrangement among other dicot CENH3s suggest that the identified genes are authentic CENH3 homologs. To verify the location of the CENH3 protein in the kinetochore regions of the Daucus chromosomes, a polyclonal antibody based on a peptide corresponding to the N-terminus of DcCENH3 was developed and used for anti-CENH3 immunostaining of mitotic root cells. The chromosomal location of CENH3 proteins in the centromere regions of the chromosomes could be confirmed. For genetic localization of the CENH3 gene in the carrot genome, a previously constructed linkage map for carrot was used for mapping a CENH3-specific simple sequence repeat (SSR) marker, and the CENH3 locus was mapped on the carrot chromosome 9.

A BAC based physical map and genome survey of the rice false smut fungus Villosiclava virens

Background

Rice false smut caused by Villosiclava virens is a devastating fungal disease that spreads in major rice-growing regions throughout the world. However, the genomic information for this fungal pathogen is limited and the pathogenic mechanism of this disease is still not clear. To facilitate genetic, molecular and genomic studies of this fungal pathogen, we constructed the first BAC-based physical map and performed the first genome survey for this species.

Results

High molecular weight genomic DNA was isolated from young mycelia of the Villosiclava virens strain UV-8b and a high-quality, large-insert and deep-coverage Bacterial Artificial Chromosome (BAC) library was constructed with the restriction enzyme HindIII. The BAC library consisted of 5,760 clones, which covers 22.7-fold of the UV-8b genome, with an average insert size of 140 kb and an empty clone rate of lower than 1%. BAC fingerprinting generated successful fingerprints for 2,290 BAC clones. Using the fingerprints, a whole genome-wide BAC physical map was constructed that contained 194 contigs (2,035 clones) spanning 51.2 Mb in physical length. Bidirectional-end sequencing of 4,512 BAC clones generated 6,560 high quality BAC end sequences (BESs), with a total length of 3,030,658 bp, representing 8.54% of the genome sequence. Analysis of the BESs revealed general genome information, including 51.52% GC content, 22.51% repetitive sequences, 376.12/Mb simple sequence repeat (SSR) density and approximately 36.01% coding regions. Sequence comparisons to other available fungal genome sequences through BESs showed high similarities to Metarhizium anisopliae, Trichoderma reesei, Nectria haematococca and Cordyceps militaris, which were generally in agreement with the 18S rRNA gene analysis results.

Conclusion

This study provides the first BAC-based physical map and genome information for the important rice fungal pathogen Villosiclava virens. The BAC clones, physical map and genome information will serve as fundamental resources to accelerate the genetic, molecular and genomic studies of this pathogen, including positional cloning, comparative genomic analysis and whole genome sequencing. The BAC library and physical map have been opened to researchers as public genomic resources (http://gresource.hzau.edu.cn/resource/resource.html).

Exploring the genome of the salt-marsh Spartina maritima (Poaceae, Chloridoideae) through BAC end sequence analysis

Spartina species play an important ecological role on salt marshes. Spartina maritima is an Old-World species distributed along the European and North-African Atlantic coasts. This hexaploid species (2n = 6x = 60, 2C = 3,700 Mb) hybridized with different Spartina species introduced from the American coasts, which resulted in the formation of new invasive hybrids and allopolyploids. Thus, S. maritima raises evolutionary and ecological interests. However, genomic information is dramatically lacking in this genus. In an effort to develop genomic resources, we analysed 40,641 high-quality bacterial artificial chromosome-end sequences (BESs), representing 26.7 Mb of the S. maritima genome. BESs were searched for sequence homology against known databases. A fraction of 16.91% of the BESs represents known repeats including a majority of long terminal repeat (LTR) retrotransposons (13.67%). Non-LTR retrotransposons represent 0.75%, DNA transposons 0.99%, whereas small RNA, simple repeats and low-complexity sequences account for 1.38% of the analysed BESs. In addition, 4,285 simple sequence repeats were detected. Using the coding sequence database of Sorghum bicolor, 6,809 BESs found homology accounting for 17.1% of all BESs. Comparative genomics with related genera reveals that the microsynteny is better conserved with S. bicolor compared to other sequenced Poaceae, where 37.6% of the paired matching BESs are correctly orientated on the chromosomes. We did not observe large macrosyntenic rearrangements using the mapping strategy employed. However, some regions appeared to have experienced rearrangements when comparing Spartina to Sorghum and to Oryza. This work represents the first overview of S. maritima genome regarding the respective coding and repetitive components. The syntenic relationships with other grass genomes examined here help clarifying evolution in Poaceae, S. maritima being a part of the poorly-known Chloridoideae sub-family.

Diversity, genetic mapping, and signatures of domestication in the carrot (Daucus carota L.) genome, as revealed by Diversity Arrays Technology (DArT) markers

Carrot is one of the most economically important vegetables worldwide, but genetic and genomic resources supporting carrot breeding remain limited. We developed a Diversity Arrays Technology (DArT) platform for wild and cultivated carrot and used it to investigate genetic diversity and to develop a saturated genetic linkage map of carrot. We analyzed a set of 900 DArT markers in a collection of plant materials comprising 94 cultivated and 65 wild carrot accessions. The accessions were attributed to three separate groups: wild, Eastern cultivated and Western cultivated. Twenty-seven markers showing signatures for selection were identified. They showed a directional shift in frequency from the wild to the cultivated, likely reflecting diversifying selection imposed in the course of domestication. A genetic linkage map constructed using 188 F2 plants comprised 431 markers with an average distance of 1.1 cM, divided into nine linkage groups. Using previously anchored single nucleotide polymorphisms, the linkage groups were physically attributed to the nine carrot chromosomes. A cluster of markers mapping to chromosome 8 showed significant segregation distortion. Two of the 27 DArT markers with signatures for selection were segregating in the mapping population and were localized on chromosomes 2 and 6. Chromosome 2 was previously shown to carry the Vrn1 gene governing the biennial growth habit essential for cultivated carrot. The results reported here provide background for further research on the history of carrot domestication and identify genomic regions potentially important for modern carrot breeding.

BAC-end sequences analysis provides first insights into coffee (Coffea canephora P.) genome composition and evolution

Coffee is one of the world's most important agricultural commodities. Coffee belongs to the Rubiaceae family in the euasterid I clade of dicotyledonous plants, to which the Solanaceae family also belongs. Two bacterial artificial chromosome (BAC) libraries of a homozygous doubled haploid plant of Coffea canephora were constructed using two enzymes, HindIII and BstYI. A total of 134,827 high quality BAC-end sequences (BESs) were generated from the 73,728 clones of the two libraries, and 131,412 BESs were conserved for further analysis after elimination of chloroplast and mitochondrial sequences. This corresponded to almost 13 % of the estimated size of the C. canephora genome. 6.7 % of BESs contained simple sequence repeats, the most abundant (47.8 %) being mononucleotide motifs. These sequences allow the development of numerous useful marker sites. Potential transposable elements (TEs) represented 11.9 % of the full length BESs. A difference was observed between the BstYI and HindIII libraries (14.9 vs. 8.8 %). Analysis of BESs against known coding sequences of TEs indicated that 11.9 % of the genome corresponded to known repeat sequences, like for other flowering plants. The number of genes in the coffee genome was estimated at 41,973 which is probably overestimated. Comparative genome mapping revealed that microsynteny was higher between coffee and grapevine than between coffee and tomato or Arabidopsis. BESs constitute valuable resources for the first genome wide survey of coffee and provide new insights into the composition and evolution of the coffee genome.

First insights into the large genome of Epimedium sagittatum (Sieb. et Zucc) Maxim, a Chinese Ttaditional medicinal plant

Epimedium sagittatum (Sieb. et Zucc) Maxim is a member of the Berberidaceae family of basal eudicot plants, widely distributed and used as a traditional medicinal plant in China for therapeutic effects on many diseases with a long history. Recent data shows that E. sagittatum has a relatively large genome, with a haploid genome size of ~4496 Mbp, divided into a small number of only 12 diploid chromosomes (2n = 2x = 12). However, little is known about Epimedium genome structure and composition. Here we present the analysis of 691 kb of high-quality genomic sequence derived from 672 randomly selected plasmid clones of E. sagittatum genomic DNA, representing ~0.0154% of the genome. The sampled sequences comprised at least 78.41% repetitive DNA elements and 2.51% confirmed annotated gene sequences, with a total GC% content of 39%. Retrotransposons represented the major class of transposable element (TE) repeats identified (65.37% of all TE repeats), particularly LTR (Long Terminal Repeat) retrotransposons (52.27% of all TE repeats). Chromosome analysis and Fluorescence in situ Hybridization of Gypsy-Ty3 retrotransposons were performed to survey the E. sagittatum genome at the cytological level. Our data provide the first insights into the composition and structure of the E. sagittatum genome, and will facilitate the functional genomic analysis of this valuable medicinal plant.

DcSto: carrot Stowaway-like elements are abundant, diverse, and polymorphic

We investigated nine families of Stowaway-like miniature inverted-repeat transposable elements (MITEs) in the carrot genome, named DcSto1 to DcSto9. All of them were AT-rich and shared a highly conserved 6 bp-long TIR typical for Stowaways. The copy number of DcSto1 elements was estimated as ca. 5,000 per diploid genome. We observed preference for clustered insertions of DcSto and other MITEs. Distribution of DcSto1 hybridization signals revealed presence of DcSto1 clusters within euchromatic regions along all chromosomes. An arrangement of eight regions encompassing DcSto insertion sites, studied in detail, was highly variable among plants representing different populations of Daucus carota. All of these insertions were polymorphic which most likely suggests a very recent mobilization of those elements. Insertions of DcSto near carrot genes and presence of putative promoters, regulatory motifs, and polyA signals within their sequences might suggest a possible involvement of DcSto in the regulation of gene expression.

Transcriptome analysis of Thapsia laciniata Rouy provides insights into terpenoid biosynthesis and diversity in Apiaceae

Thapsia laciniata Rouy (Apiaceae) produces irregular and regular sesquiterpenoids with thapsane and guaiene carbon skeletons, as found in other Apiaceae species. A transcriptomic analysis utilizing Illumina next-generation sequencing enabled the identification of novel genes involved in the biosynthesis of terpenoids in Thapsia. From 66.78 million HQ paired-end reads obtained from T. laciniata roots, 64.58 million were assembled into 76,565 contigs (N50: 1261 bp). Seventeen contigs were annotated as terpene synthases and five of these were predicted to be sesquiterpene synthases. Of the 67 contigs annotated as cytochromes P450, 18 of these are part of the CYP71 clade that primarily performs hydroxylations of specialized metabolites. Three contigs annotated as aldehyde dehydrogenases grouped phylogenetically with the characterized ALDH1 from Artemisia annua and three contigs annotated as alcohol dehydrogenases grouped with the recently described ADH1 from A. annua. ALDH1 and ADH1 were characterized as part of the artemisinin biosynthesis. We have produced a comprehensive EST dataset for T. laciniata roots, which contains a large sample of the T. laciniata transcriptome. These transcriptome data provide the foundation for future research into the molecular basis for terpenoid biosynthesis in Thapsia and on the evolution of terpenoids in Apiaceae.

Molecular mapping of vernalization requirement and fertility restoration genes in carrot

Carrot (Daucus carota L.) is a cool-season vegetable normally classified as a biennial species, requiring vernalization to induce flowering. Nevertheless, some cultivars adapted to warmer climates require less vernalization and can be classified as annual. Most modern carrot cultivars are hybrids which rely upon cytoplasmic male-sterility for commercial production. One major gene controlling floral initiation and several genes restoring male fertility have been reported but none have been mapped. The objective of the present work was to develop the first linkage map of carrot locating the genomic regions that control vernalization response and fertility restoration. Using an F(2) progeny, derived from the intercross between the annual cultivar 'Criolla INTA' and a petaloid male sterile biennial carrot evaluated over 2 years, both early flowering habit, which we name Vrn1, and restoration of petaloid cytoplasmic male sterility, which we name Rf1, were found to be dominant traits conditioned by single genes. On a map of 355 markers covering all 9 chromosomes with a total map length of 669 cM and an average marker-to-marker distance of 1.88 cM, Vrn1 mapped to chromosome 2 with flanking markers at 0.70 and 0.46 cM, and Rf1 mapped to chromosome 9 with flanking markers at 4.38 and 1.12 cM. These are the first two reproductive traits mapped in the carrot genome, and their map location and flanking markers provide valuable tools for studying traits important for carrot domestication and reproductive biology, as well as facilitating carrot breeding.

Fluorescent in situ hybridization with arbitrarily amplified DNA fragments differentiates carrot (Daucus carota L.) chromosomes

Carrot (Daucus carota L.) chromosomes are small and poorly differentiated in size and morphology. Here we demonstrate that fluorescent in situ hybridization (FISH) signals derived from arbitrary PCR probes can be used for chromosome identification in carrot. To prepare probes, we searched for nonpolymorphic products abundantly amplified with arbitrary decamer primers in a group of accessions representing carrot genetic diversity. As a result, 13 fragments ranging in size from 517 to 1758 bp were selected, sequenced, and used as probes for fluorescent in situ hybridization. Four of these probes produced clear and reproducible hybridization signals. The sequences showed similarity to a number of carrot BAC-end sequences, indicating their repetitive character. Three of them were similar to internal portions of gypsy and copia LTR retrotransposons previously identified in plants. Hybridization signals for the four probes were observed as dotted tracks on chromosomes, differing in distribution and intensity. Generally, they were present in pericentromeric and (or) interstitial localizations on chromosome arms. The use of the four probes allowed discrimination of chromosome pairs and construction of more detailed karyotypes and idiograms of carrot.

A bacterial artificial chromosome library for 'Jefferson' hazelnut and identification of clones associated with eastern filbert blight resistance and pollen-stigma incompatibility

European hazelnut (Corylus avellana L.) is the only economically important nut crop in the family Betulaceae. Because of its small genome size (~385 Mb / 1C), relatively short life cycle, availability of a dense linkage map, and amenability to transformation by Agrobacterium, the European hazelnut could serve as a model plant for the Betulaceae. Here we report the construction of a bacterial artificial chromosome (BAC) library for 'Jefferson' hazelnut using the cloning enzyme MboI and the vector pECBAC1 (BamHI site). The library consists of 39,936 clones arrayed in 104,384-well microtitre plates with a mean insert size of 117 kb. The genomic coverage of the library is estimated to be about 12 genome equivalents. This library provides a valuable resource for the map-based cloning of two important genes, the resistance gene from 'Gasaway' that confers resistance to eastern filbert blight caused by the fungus Anisogramma anomala (Peck) E. Müller and the S locus that controls pollen-stigma incompatibility. Fine-resolution mapping near the two loci was carried out using random amplified polymorphic DNA (RAPD) and simple sequence repeat (SSR) markers. Fine mapping at the disease resistance locus showed that markers W07-375 and X01-825 flanked the resistance locus at distances of 0.06 and 0.05 cM, respectively. The S locus is flanked by markers 204-950 and KG819-200 at distances of 0.14 and 0.24 cM, respectively. Assuming that 1 cM corresponds to a physical distance of 430 kb, it will take approximately two to three chromosome walks to assemble BAC contigs that span both loci.

De novo assembly and characterization of the carrot transcriptome reveals novel genes, new markers, and genetic diversity

Background

Among next generation sequence technologies, platforms such as Illumina and SOLiD produce short reads but with higher coverage and lower cost per sequenced nucleotide than 454 or Sanger. A challenge now is to develop efficient strategies to use short-read length platforms for de novo assembly and marker development. The scope of this study was to develop a de novo assembly of carrot ESTs from multiple genotypes using the Illumina platform, and to identify polymorphisms.

Results

A de novo assembly of transcriptome sequence from four genetic backgrounds produced 58,751 contigs and singletons. Over 50% of these assembled sequences were annotated allowing detection of transposable elements and new carrot anthocyanin genes. Presence of multiple genetic backgrounds in our assembly allowed the identification of 114 computationally polymorphic SSRs, and 20,058 SNPs at a depth of coverage of 20× or more. Polymorphisms were predominantly between inbred lines except for the cultivated x wild RIL pool which had high intra-sample polymorphism. About 90% and 88% of tested SSR and SNP primers amplified a product, of which 70% and 46%, respectively, were of the expected size. Out of verified SSR and SNP markers 84% and 82% were polymorphic. About 25% of SNPs genotyped were polymorphic in two diverse mapping populations.

Conclusions

This study confirmed the potential of short read platforms for de novo EST assembly and identification of genetic polymorphisms in carrot. In addition we produced the first large-scale transcriptome of carrot, a species lacking genomic resources.

Microsatellite isolation and marker development in carrot - genomic distribution, linkage mapping, genetic diversity analysis and marker transferability across Apiaceae

BACKGROUND

The Apiaceae family includes several vegetable and spice crop species among which carrot is the most economically important member, with ~21 million tons produced yearly worldwide. Despite its importance, molecular resources in this species are relatively underdeveloped. The availability of informative, polymorphic, and robust PCR-based markers, such as microsatellites (or SSRs), will facilitate genetics and breeding of carrot and other Apiaceae, including integration of linkage maps, tagging of phenotypic traits and assisting positional gene cloning. Thus, with the purpose of isolating carrot microsatellites, two different strategies were used; a hybridization-based library enrichment for SSRs, and bioinformatic mining of SSRs in BAC-end sequence and EST sequence databases. This work reports on the development of 300 carrot SSR markers and their characterization at various levels.

RESULTS

Evaluation of microsatellites isolated from both DNA sources in subsets of 7 carrot F2 mapping populations revealed that SSRs from the hybridization-based method were longer, had more repeat units and were more polymorphic than SSRs isolated by sequence search. Overall, 196 SSRs (65.1%) were polymorphic in at least one mapping population, and the percentage of polymophic SSRs across F2 populations ranged from 17.8 to 24.7. Polymorphic markers in one family were evaluated in the entire F2, allowing the genetic mapping of 55 SSRs (38 codominant) onto the carrot reference map. The SSR loci were distributed throughout all 9 carrot linkage groups (LGs), with 2 to 9 SSRs/LG. In addition, SSR evaluations in carrot-related taxa indicated that a significant fraction of the carrot SSRs transfer successfully across Apiaceae, with heterologous amplification success rate decreasing with the target-species evolutionary distance from carrot. SSR diversity evaluated in a collection of 65 D. carota accessions revealed a high level of polymorphism for these selected loci, with an average of 19 alleles/locus and 0.84 expected heterozygosity.

CONCLUSIONS

The addition of 55 SSRs to the carrot map, together with marker characterizations in six other mapping populations, will facilitate future comparative mapping studies and integration of carrot maps. The markers developed herein will be a valuable resource for assisting breeding, genetic, diversity, and genomic studies of carrot and other Apiaceae.

Analysis of BAC end sequences in oak, a keystone forest tree species, providing insight into the composition of its genome

BACKGROUND

One of the key goals of oak genomics research is to identify genes of adaptive significance. This information may help to improve the conservation of adaptive genetic variation and the management of forests to increase their health and productivity. Deep-coverage large-insert genomic libraries are a crucial tool for attaining this objective. We report herein the construction of a BAC library for Quercus robur, its characterization and an analysis of BAC end sequences.

RESULTS

The EcoRI library generated consisted of 92,160 clones, 7% of which had no insert. Levels of chloroplast and mitochondrial contamination were below 3% and 1%, respectively. Mean clone insert size was estimated at 135 kb. The library represents 12 haploid genome equivalents and, the likelihood of finding a particular oak sequence of interest is greater than 99%. Genome coverage was confirmed by PCR screening of the library with 60 unique genetic loci sampled from the genetic linkage map. In total, about 20,000 high-quality BAC end sequences (BESs) were generated by sequencing 15,000 clones. Roughly 5.88% of the combined BAC end sequence length corresponded to known retroelements while ab initio repeat detection methods identified 41 additional repeats. Collectively, characterized and novel repeats account for roughly 8.94% of the genome. Further analysis of the BESs revealed 1,823 putative genes suggesting at least 29,340 genes in the oak genome. BESs were aligned with the genome sequences of Arabidopsis thaliana, Vitis vinifera and Populus trichocarpa. One putative collinear microsyntenic region encoding an alcohol acyl transferase protein was observed between oak and chromosome 2 of V. vinifera.

CONCLUSIONS

This BAC library provides a new resource for genomic studies, including SSR marker development, physical mapping, comparative genomics and genome sequencing. BES analysis provided insight into the structure of the oak genome. These sequences will be used in the assembly of a future genome sequence for oak.

Physical mapping and BAC-end sequence analysis provide initial insights into the flax (Linum usitatissimum L.) genome

BACKGROUND

Flax (Linum usitatissimum L.) is an important source of oil rich in omega-3 fatty acids, which have proven health benefits and utility as an industrial raw material. Flax seeds also contain lignans which are associated with reducing the risk of certain types of cancer. Its bast fibres have broad industrial applications. However, genomic tools needed for molecular breeding were non existent. Hence a project, Total Utilization Flax GENomics (TUFGEN) was initiated. We report here the first genome-wide physical map of flax and the generation and analysis of BAC-end sequences (BES) from 43,776 clones, providing initial insights into the genome.

RESULTS

The physical map consists of 416 contigs spanning ~368 Mb, assembled from 32,025 fingerprints, representing roughly 54.5% to 99.4% of the estimated haploid genome (370-675 Mb). The N50 size of the contigs was estimated to be ~1,494 kb. The longest contig was ~5,562 kb comprising 437 clones. There were 96 contigs containing more than 100 clones. Approximately 54.6 Mb representing 8-14.8% of the genome was obtained from 80,337 BES. Annotation revealed that a large part of the genome consists of ribosomal DNA (~13.8%), followed by known transposable elements at 6.1%. Furthermore, ~7.4% of sequence was identified to harbour novel repeat elements. Homology searches against flax-ESTs and NCBI-ESTs suggested that ~5.6% of the transcriptome is unique to flax. A total of 4064 putative genomic SSRs were identified and are being developed as novel markers for their use in molecular breeding.

CONCLUSION

The first genome-wide physical map of flax constructed with BAC clones provides a framework for accessing target loci with economic importance for marker development and positional cloning. Analysis of the BES has provided insights into the uniqueness of the flax genome. Compared to other plant genomes, the proportion of rDNA was found to be very high whereas the proportion of known transposable elements was low. The SSRs identified from BES will be valuable in saturating existing linkage maps and for anchoring physical and genetic maps. The physical map and paired-end reads from BAC clones will also serve as scaffolds to build and validate the whole genome shotgun assembly.

Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing

Background

Milkweeds (Asclepias L.) have been extensively investigated in diverse areas of evolutionary biology and ecology; however, there are few genetic resources available to facilitate and compliment these studies. This study explored how low coverage genome sequencing of the common milkweed (Asclepias syriaca L.) could be useful in characterizing the genome of a plant without prior genomic information and for development of genomic resources as a step toward further developing A. syriaca as a model in ecology and evolution.

Results

A 0.5× genome of A. syriaca was produced using Illumina sequencing. A virtually complete chloroplast genome of 158,598 bp was assembled, revealing few repeats and loss of three genes: accD, clpP, and ycf1. A nearly complete rDNA cistron (18S-5.8S-26S; 7,541 bp) and 5S rDNA (120 bp) sequence were obtained. Assessment of polymorphism revealed that the rDNA cistron and 5S rDNA had 0.3% and 26.7% polymorphic sites, respectively. A partial mitochondrial genome sequence (130,764 bp), with identical gene content to tobacco, was also assembled. An initial characterization of repeat content indicated that Ty1/copia-like retroelements are the most common repeat type in the milkweed genome. At least one A. syriaca microread hit 88% of Catharanthus roseus (Apocynaceae) unigenes (median coverage of 0.29×) and 66% of single copy orthologs (COSII) in asterids (median coverage of 0.14×). From this partial characterization of the A. syriaca genome, markers for population genetics (microsatellites) and phylogenetics (low-copy nuclear genes) studies were developed.

Conclusions

The results highlight the promise of next generation sequencing for development of genomic resources for any organism. Low coverage genome sequencing allows characterization of the high copy fraction of the genome and exploration of the low copy fraction of the genome, which facilitate the development of molecular tools for further study of a target species and its relatives. This study represents a first step in the development of a community resource for further study of plant-insect co-evolution, anti-herbivore defense, floral developmental genetics, reproductive biology, chemical evolution, population genetics, and comparative genomics using milkweeds, and A. syriaca in particular, as ecological and evolutionary models.

Advancing Eucalyptus genomics: identification and sequencing of lignin biosynthesis genes from deep-coverage BAC libraries

Background

Eucalyptus species are among the most planted hardwoods in the world because of their rapid growth, adaptability and valuable wood properties. The development and integration of genomic resources into breeding practice will be increasingly important in the decades to come. Bacterial artificial chromosome (BAC) libraries are key genomic tools that enable positional cloning of important traits, synteny evaluation, and the development of genome framework physical maps for genetic linkage and genome sequencing.

Results

We describe the construction and characterization of two deep-coverage BAC libraries EG_Ba and EG_Bb obtained from nuclear DNA fragments of E. grandis (clone BRASUZ1) digested with HindIII and BstYI, respectively. Genome coverages of 17 and 15 haploid genome equivalents were estimated for EG_Ba and EG_Bb, respectively. Both libraries contained large inserts, with average sizes ranging from 135 Kb (Eg_Bb) to 157 Kb (Eg_Ba), very low extra-nuclear genome contamination providing a probability of finding a single copy gene ≥ 99.99%. Libraries were screened for the presence of several genes of interest via hybridizations to high-density BAC filters followed by PCR validation. Five selected BAC clones were sequenced and assembled using the Roche GS FLX technology providing the whole sequence of the E. grandis chloroplast genome, and complete genomic sequences of important lignin biosynthesis genes.

Conclusions

The two E. grandis BAC libraries described in this study represent an important milestone for the advancement of Eucalyptus genomics and forest tree research. These BAC resources have a highly redundant genome coverage (> 15×), contain large average inserts and have a very low percentage of clones with organellar DNA or empty vectors. These publicly available BAC libraries are thus suitable for a broad range of applications in genetic and genomic research in Eucalyptus and possibly in related species of Myrtaceae, including genome sequencing, gene isolation, functional and comparative genomics. Because they have been constructed using the same tree (E. grandis BRASUZ1) whose full genome is being sequenced, they should prove instrumental for assembly and gap filling of the upcoming Eucalyptus reference genome sequence.

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).

An overview of the Phalaenopsis orchid genome through BAC end sequence analysis

BACKGROUND

Phalaenopsis orchids are popular floral crops, and development of new cultivars is economically important to floricultural industries worldwide. Analysis of orchid genes could facilitate orchid improvement. Bacterial artificial chromosome (BAC) end sequences (BESs) can provide the first glimpses into the sequence composition of a novel genome and can yield molecular markers for use in genetic mapping and breeding.

RESULTS

We used two BAC libraries (constructed using the BamHI and HindIII restriction enzymes) of Phalaenopsis equestris to generate pair-end sequences from 2,920 BAC clones (71.4% and 28.6% from the BamHI and HindIII libraries, respectively), at a success rate of 95.7%. A total of 5,535 BESs were generated, representing 4.5 Mb, or about 0.3% of the Phalaenopsis genome. The trimmed sequences ranged from 123 to 1,397 base pairs (bp) in size, with an average edited read length of 821 bp. When these BESs were subjected to sequence homology searches, it was found that 641 (11.6%) were predicted to represent protein-encoding regions, whereas 1,272 (23.0%) contained repetitive DNA. Most of the repetitive DNA sequences were gypsy- and copia-like retrotransposons (41.9% and 12.8%, respectively), whereas only 10.8% were DNA transposons. Further, 950 potential simple sequence repeats (SSRs) were discovered. Dinucleotides were the most abundant repeat motifs; AT/TA dimer repeats were the most frequent SSRs, representing 253 (26.6%) of all identified SSRs. Microsynteny analysis revealed that more BESs mapped to the whole-genome sequences of poplar than to those of grape or Arabidopsis, and even fewer mapped to the rice genome. This work will facilitate analysis of the Phalaenopsis genome, and will help clarify similarities and differences in genome composition between orchids and other plant species.

CONCLUSION

Using BES analysis, we obtained an overview of the Phalaenopsis genome in terms of gene abundance, the presence of repetitive DNA and SSR markers, and the extent of microsynteny with other plant species. This work provides a basis for future physical mapping of the Phalaenopsis genome and advances our knowledge thereof.

The first insight into the Taxus genome via fosmid library construction and end sequencing

Taxus mairei is a critically endangered and commercially important cultured medicinal gymnosperm in China and forms an important medicinal resource, but the research of its genome is absent. In this study, we constructed a T. mairei fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of one million clones with an average insert size of about 39 kb, amounting to 3.9 genome equivalents. Fosmid stability assays indicate that T. mairei DNA was stable during propagation in the fosmid system. End sequencing of both 5' and 3' ends of 968 individual clones generated 1,923 sequences after trimming, with an average sequence length of 839 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 560 (29.1%) significant hits (E < e(-5)). Repetitive sequences analysis revealed that 20.8% of end sequences are repetitive elements, which were composed of retroelements, DNA transposons, satellites, simple repeats, and low complexity sequences. The distribution pattern of various repeat types was found to be more similar to the gymnosperm Pinus and Picea than to the monocot and dicot. The satellites of T. mairei were significantly longer than those of P. taeda and P. glauca. The tetra-nucleotide repeats of T. mairei were much longer than those of P. glauca and P. taeda. The fosmid library and the fosmid end sequences, for the first time, will serve as a useful resource for large-scale genome sequencing, physical mapping, SSR marker development and positional cloning, and provide a better understanding of the Taxus genome.

Major cytogenetic landmarks and karyotype analysis inDaucus carota and other Apiaceae

Karyotype analysis provides insights into genome organization at the chromosome level and into chromosome evolution. Chromosomes were marked for comparative karyotype analysis using FISH localization of rDNA genes for the first time in Apioideae species including taxa of economic importance and several wild Daucus relatives. Interestingly, Daucus species did not vary in number of rDNA loci despite variation in chromosome number (2n = 18, 20, 22, and 44) and previous publications suggesting multiple loci. All had single loci for both 5S and 18S-25S (nucleolar organizing region) rDNA, located on two different chromosome pairs. The 5S rDNA was on the short arm of a metacentric chromosome pair in D. crinitus (2n = 22) and D. glochidiatus (2n = 44) and on the long arm of a metacentric pair in other Daucus species, suggesting possible rearrangement of this chromosome. For other Apiaceae, from two (Apium graveolens), to three (Orlaya grandiflora), to four (Cuminum cyminum) chromosomes had 18S-25S rDNA sites. Variability for number and position of the 5S rDNA was also observed. FISH signals enabled us to identify 20-40% of the chromosome complement among species examined. Comparative karyotype analysis provides insights into the fundamental aspects of chromosome evolution in Daucus.