Characterization of a Genomic Region under Selection in Cultivated Carrot (Daucus carota subsp. sativus) Reveals a Candidate Domestication Gene

Linkage mapping of root shape traits in two carrot populations

This study investigated the genetic basis of carrot root shape traits using composite interval mapping in two biparental populations (n = 119 and n = 128). The roots of carrot F2:3 progenies were grown over 2 years and analyzed using a digital imaging pipeline to extract root phenotypes that compose market class. Broad-sense heritability on an entry-mean basis ranged from 0.46 to 0.80 for root traits. Reproducible quantitative trait loci (QTL) were identified on chromosomes 2 and 6 on both populations. Colocalization of QTLs for phenotypically correlated root traits was also observed and coincided with previously identified QTLs in published association and linkage mapping studies. Individual QTLs explained between 14 and 27% of total phenotypic variance across traits, while four QTLs for length-to-width ratio collectively accounted for up to 73% of variation. Predicted genes associated with the OFP-TRM (OVATE Family Proteins-TONNEAU1 Recruiting Motif) and IQD (IQ67 domain) pathway were identified within QTL support intervals. This observation raises the possibility of extending the current regulon model of fruit shape to include carrot storage roots. Nevertheless, the precise molecular mechanisms through which this pathway operates in roots characterized by secondary growth originating from cambium layers remain unknown.

Diverse and mobile: eccDNA-based identification of carrot low-copy-number LTR retrotransposons active in callus cultures

Long terminal repeat retrotransposons (LTR-RTs) are mobilized via an RNA intermediate using a 'copy and paste' mechanism, and account for the majority of repetitive DNA in plant genomes. As a side effect of mobilization, the formation of LTR-RT-derived extrachromosomal circular DNAs (eccDNAs) occurs. Thus, high-throughput sequencing of eccDNA can be used to identify active LTR-RTs in plant genomes. Despite the release of a reference genome assembly, carrot LTR-RTs have not yet been thoroughly characterized. LTR-RTs are abundant and diverse in the carrot genome. We identified 5976 carrot LTR-RTs, 2053 and 1660 of which were attributed to Copia and Gypsy superfamilies, respectively. They were further classified into lineages, families and subfamilies. More diverse LTR-RT lineages, i.e. lineages comprising many low-copy-number subfamilies, were more frequently associated with genic regions. Certain LTR-RT lineages have been recently active in Daucus carota. In particular, low-copy-number LTR-RT subfamilies, e.g. those belonging to the DcAle lineage, have significantly contributed to carrot genome diversity as a result of continuing activity. We utilized eccDNA sequencing to identify and characterize two DcAle subfamilies, Alex1 and Alex3, active in carrot callus. We documented 14 and 32 de novo insertions of Alex1 and Alex3, respectively, which were positioned in non-repetitive regions.

Detection of Chromosomal Segments Introgressed from Wild Species of Carrot into Cultivars: Quantitative Trait Loci Mapping for Morphological Features in Backcross Inbred Lines

Cultivated carrot is thought to have been domesticated from a wild species, and various phenotypes developed through human domestication and selection over the past several centuries. Little is known about the genomic contribution of wild species to the phenotypes of present-day cultivars, although several studies have focused on identifying genetic loci that contribute to the morphology of storage roots. A backcross inbred line (BIL) population derived from a cross between the wild species Daucus carota ssp. carota "Songzi" and the orange cultivar "Amsterdam forcing" was developed. The morphological features in the BIL population became more diverse after several generations of selfing BC₂F₁ plants. Only few lines retained features of wild parent. Genomic resequencing of the two parental lines and the BILs resulted in 3,223,651 single nucleotide polymorphisms (SNPs), and 13,445 bin markers were generated using a sliding window approach. We constructed a genetic map with 2027 bins containing 154,776 SNPs; the total genetic distance was 1436.43 cM and the average interval between the bins was 0.71 cm. Five stable QTLs related to root length, root shoulder width, dry material content of root, and ratio of root shoulder width to root middle width were consistently detected on chromosome 2 in both years and explained 23.4-66.9% of the phenotypic variance. The effects of introgressed genomic segments from the wild species on the storage root are reported and will enable the identification of functional genes that control root morphological traits in carrot.

The influence of the Or and Carotene Hydroxylase genes on carotenoid accumulation in orange carrots [Daucus carota (L.)]

The Or and CH genes are necessary for the accumulation of high amounts of β-carotene and other carotenoid pigments in carrot roots, in addition to the Y and Y₂ genes. Carrot taproot color results from the accumulation of various carotenoid and anthocyanin pigments. Recently, the Or gene was identified as a candidate gene associated with the accumulation of β-carotene and other provitamin A carotenoids in roots. The specific molecular mechanisms involved with this process, as well as the interactions between Or and the other genes involved in this process are not well understood. In order to better characterize the effect that Or alleles have on conditioning the accumulation of carotenoids in roots, we analyzed an F₃ family fixed homozygous recessive for y and y₂, derived from a cross between an orange carrot and a white wild carrot, segregating for the two known Or alleles, which we name Or_c and Or_w. QTL mapping across three different environments revealed that the accumulation of several carotenoids was associated with the Or_c allele, with consistent patterns across environments. A second QTL on chromosome 7, harboring a carotene hydroxylase gene homologous to Lut5 in Arabidopsis, was also associated with the accumulation of several carotenoids. Two alleles for this gene, which we name CH_c and CH_w, were discovered to be segregating in this population. Our study provides further evidence that Or and CH are likely involved with controlling the accumulation of β-carotene and may be involved with modulating carotenoid flux in carrot, demonstrating that both were important domestication genes in carrot.

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.

Characteristics of the AT-Hook Motif Containing Nuclear Localized (AHL) Genes in Carrot Provides Insight into Their Role in Plant Growth and Storage Root Development

The AT-hook motif containing nuclear localized (AHL) gene family, controlling various developmental processes, is conserved in land plants. They comprise Plant and Prokaryote Conserved (PPC) domain and one or two AT-hook motifs. DcAHLc1 has been proposed as a candidate gene governing the formation of the carrot storage root. We identified and in-silico characterized carrot AHL proteins, performed phylogenetic analyses, investigated their expression profiles and constructed gene coexpression networks. We found 47 AHL genes in carrot and grouped them into two clades, A and B, comprising 29 and 18 genes, respectively. Within Clade-A, we distinguished three subclades, one of them grouping noncanonical AHLs differing in their structure (two PPC domains) and/or cellular localization (not nucleus). Coexpression network analysis attributed AHLs expressed in carrot roots into four of the 72 clusters, some of them showing a large number of interactions. Determination of expression profiles of AHL genes in various tissues and samples provided basis to hypothesize on their possible roles in the development of the carrot storage root. We identified a group of rapidly evolving noncanonical AHLs, possibly differing functionally from typical AHLs, as suggested by their expression profiles and their predicted cellular localization. We pointed at several AHLs likely involved in the development of the carrot storage root.

Mining for Candidate Genes Controlling Secondary Growth of the Carrot Storage Root

Background: Diverse groups of carrot cultivars have been developed to meet consumer demands and industry needs. Varietal groups of the cultivated carrot are defined based on the shape of roots. However, little is known about the genetic basis of root shape determination. Methods: Here, we used 307 carrot plants from 103 open-pollinated cultivars for a genome wide association study to identify genomic regions associated with the storage root morphology. Results: A 180 kb-long region on carrot chromosome 1 explained 10% of the total observed phenotypic variance in the shoulder diameter. Within that region, DcDCAF1 and DcBTAF1 genes were proposed as candidates controlling secondary growth of the carrot storage root. Their expression profiles differed between the cultivated and the wild carrots, likely indicating that their elevated expression was required for the development of edible roots. They also showed higher expression at the secondary root growth stage in cultivars producing thick roots, as compared to those developing thin roots. Conclusions: We provided evidence for a likely involvement of DcDCAF1 and/or DcBTAF1 in the development of the carrot storage root and developed a genotyping assay facilitating the identification of variants in the region on carrot chromosome 1 associated with secondary growth of the carrot root.

Stowaway miniature inverted repeat transposable elements are important agents driving recent genomic diversity in wild and cultivated carrot

BACKGROUND

Miniature inverted repeat transposable elements (MITEs) are small non-autonomous DNA transposons that are ubiquitous in plant genomes, and are mobilised by their autonomous relatives. Stowaway MITEs are derived from and mobilised by elements from the mariner superfamily. Those elements constitute a significant portion of the carrot genome; however the variation caused by Daucus carota Stowaway MITEs (DcStos), their association with genes and their putative impact on genome evolution has not been comprehensively analysed.

RESULTS

Fourteen families of Stowaway elements DcStos occupy about 0.5% of the carrot genome. We systematically analysed 31 genomes of wild and cultivated Daucus carota, yielding 18.5 thousand copies of these elements, showing remarkable insertion site polymorphism. DcSto element demography differed based on the origin of the host populations, and corresponded with the four major groups of D. carota, wild European, wild Asian, eastern cultivated and western cultivated. The DcStos elements were associated with genes, and most frequently occurred in 5' and 3' untranslated regions (UTRs). Individual families differed in their propensity to reside in particular segments of genes. Most importantly, DcSto copies in the 2 kb regions up- and downstream of genes were more frequently associated with open reading frames encoding transcription factors, suggesting their possible functional impact. More than 1.5% of all DcSto insertion sites in different host genomes contained different copies in exactly the same position, indicating the existence of insertional hotspots. The DcSto7b family was much more polymorphic than the other families in cultivated carrot. A line of evidence pointed at its activity in the course of carrot domestication, and identified Dcmar1 as an active carrot mariner element and a possible source of the transposition machinery for DcSto7b.

CONCLUSION

Stowaway MITEs have made a substantial contribution to the structural and functional variability of the carrot genome.

An Automated Image Analysis Pipeline Enables Genetic Studies of Shoot and Root Morphology in Carrot (Daucus carota L.)

Carrot is a globally important crop, yet efficient and accurate methods for quantifying its most important agronomic traits are lacking. To address this problem, we developed an automated image analysis platform that extracts components of size and shape for carrot shoots and roots, which are necessary to advance carrot breeding and genetics. This method reliably measured variation in shoot size and shape, petiole number, petiole length, and petiole width as evidenced by high correlations with hundreds of manual measurements. Similarly, root length and biomass were accurately measured from the images. This platform also quantified shoot and root shapes in terms of principal components, which do not have traditional, manually measurable equivalents. We applied the pipeline in a study of a six-parent diallel population and an F2 mapping population consisting of 316 individuals. We found high levels of repeatability within a growing environment, with low to moderate repeatability across environments. We also observed co-localization of quantitative trait loci for shoot and root characteristics on chromosomes 1, 2, and 7, suggesting these traits are controlled by genetic linkage and/or pleiotropy. By increasing the number of individuals and phenotypes that can be reliably quantified, the development of a rapid, automated image analysis pipeline to measure carrot shoot and root morphology will expand the scope and scale of breeding and genetic studies.

Carotenoid Presence Is Associated with the Or Gene in Domesticated Carrot

Carrots are among the richest sources of provitamin A carotenes in the human diet, but genetic variation in the carotenoid pathway does not fully explain the high levels of carotenoids in carrot roots. Using a diverse collection of modern and historic domesticated varieties, and wild carrot accessions, an association analysis for orange pigmentation revealed a significant genomic region that contains the Or gene, advancing it as a candidate for carotenoid presence in carrot. Analysis of sequence variation at the Or locus revealed a nonsynonymous mutation cosegregating with carotenoid content. This mutation was absent in all wild carrot samples and nearly fixed in all orange domesticated samples. Or has been found to control carotenoid presence in other crops but has not previously been described in carrot. Our analysis also allowed us to more completely characterize the genetic structure of carrot, showing that the Western domesticated carrot largely forms one genetic group, despite dramatic phenotypic differences among market classes. Eastern domesticated and wild accessions form a second group, which reflects the recent cultivation history of carrots in Central Asia. Other wild accessions form distinct geographic groups, particularly on the Iberian peninsula and in Northern Africa. Using genome-wide F_st , nucleotide diversity, and the cross-population composite likelihood ratio, we analyzed the genome for regions putatively under selection during domestication and identified 12 regions that were significant for all three methods of detection, one of which includes the Or gene. The Or domestication allele appears to have been selected after the initial domestication of yellow carrots in the East, near the proposed center of domestication in Central Asia. The rapid fixation of the Or domestication allele in almost all orange and nonorange carrots in the West may explain why it has not been found with less genetically diverse mapping populations.

Comparative Transcriptomics of Root Development in Wild and Cultivated Carrots

The carrot is the most popular root vegetable worldwide. The genetic makeup underlying the development of the edible storage root are fragmentary. Here, we report the first comparative transcriptome analysis between wild and cultivated carrot roots at multiple developmental stages. Overall, 3285, 4637, and 570 genes were differentially expressed in the cultivated carrot in comparisons made for young plants versus developing roots, young plants versus mature roots, and developing roots versus mature roots, respectively. Of those, 1916, 2645, and 475, respectively, were retained after filtering out genes showing similar profiles of expression in the wild carrot. They were assumed to be of special interest with respect to the development of the storage root. Among them, transcription factors and genes encoding proteins involved in post-translational modifications (signal transduction and ubiquitination) were mostly upregulated, while those involved in redox signaling were mostly downregulated. Also, genes encoding proteins regulating cell cycle, involved in cell divisions, development of vascular tissue, water transport, and sugar metabolism were enriched in the upregulated clusters. Genes encoding components of photosystem I and II, together with genes involved in carotenoid biosynthesis, were upregulated in the cultivated roots, as opposed to the wild roots; however, they were largely downregulated in the mature storage root, as compared with the young and developing root. The experiment produced robust resources for future investigations on the regulation of storage root formation in carrot and Apiaceae.