DCAR_030401

An orange, doubled-haploid, Nantes-type carrot (DH1) was used for genome sequencing. We used BAC end sequences and a newly developed linkage map with 2,075 markers to correct 135 scaffolds with one or more chimeric regions. The resulting v2.0 assembly spans 421.5 Mb and contains 4,907 scaffolds (N50 of 12.7 Mb), accounting for ∼90% of the estimated genome size of 473 Mb. The scaftig N50 of 31.2 kb is similar to those of other high-quality genome assemblies such as potato and pepper. About 86% (362 Mb) of the assembled genome is included in only 60 superscaffolds anchored to the nine pseudomolecules. The longest superscaffold spans 30.2 Mb, 85% of chromosome 4.

There are a few different naming schemes for this assembly. First there is the
Authors' original naming scheme: Sequences with DCARv2 prefix are the original assembly as submitted to NCBI. These are labelled DCARv2_Chr1 through DCARv2_Chr9 for the chromosome pseudomolecules, DCARv2_MT and DCARv2_PT for the organellar assemblies, DCARv2_B1 and up for unincorporated superscaffolds, DCARv2_S26.1 and up for unincorporated scaffolds, and DCARv2_C10542132 and up for unincorporated contigs. A file with sequences using this naming scheme can be downloaded from the File: link below.
These sequences can be viewed in JBrowse here.

Phytozome genome ID 388: The authors' sequences and gene predictions were also submitted to Phytozome, and can be accessed at this address: https://phytozome-next.jgi.doe.gov/info/Dcarota_v2_0

LNRQ01: These sequences were then assigned GenBank accession numbers starting at LNRQ01000001.1 which corresponds to DCARv2_Chr1, up to LNRQ01004826.1 which corresponds to an unincorporated contig, DCARv2_C10750146. These reside in bioproject PRJNA268187, which is a subproject of umbrella project PRJNA285926.

Assembly GCA_001625215.1: The genome assembly was later defined an accession number GCA_001625215.1 for assembly ASM162521v1 which consists of only the 9 chromosome sequences and the plastid assembly, which have accession numbers from CM004278.1 to CM004286.1 for the chromosomes and CM004358.1 for the plastid. The mitochondrial genome was not included because it is classified as an incomplete sequence.

RefSeq: The assembly was then later added to RefSeq, and there another new set of identifiers was defined from NC_030381.1 to NC_030389.1 for the chromosomes, and from NW_016089425.1 to NW_016094239.1 for unincorporated scaffolds and contigs. These reside in bioproject PRJNA326436. Note that NCBI substituted different assembled organellar genomes from different genotypes for the RefSeq records.

The NCBI Sequence report lists the correspondences between the various naming methods

Link to the LNRQ01000000.1 master record at NCBI

Raw Reads: Link to SRA accessions used for the genome assembly

This genome is available in the CarrotOmics Blast Search

This analysis was part of the carrot genome assembly publication

In population 97837, root tissue was collected from plants with yellow (yyY₂Y₂) and white (YYY₂Y₂) genotypes, with two biological replications per genotype, at 80 days after planting (DAP). In population 70796, root tissue was collected from plants with dark orange (yyy₂y₂) and pale orange (YYy₂y₂) genotypes, with three biological replications per genotype, at 100 DAP. Total RNA was extracted from whole root tissue using the TRIzol® Plus RNA Purification Kit (Life Technologies, Carlsbad, CA) in accordance with the manufacturer’s protocol. RNA was treated for DNA contamination with the TurboDNA-free kit (Life Technologies, Carlsbad, CA). RNA quantity and integrity was confirmed with an Experion RNA StdSens Analysis kit (Bio-Rad, Hercules, CA). All samples had RQI values above 8.0.

For each biological replicate, a 133 nt insert size paired-end library was prepared at the Biotechnology Center, UW-Madison (WI, USA). Libraries were sequenced on Illumina HiSeq2000 lanes using 2 × 100 nt reads. Reads were filtered using Trimmomatic version 0.32 with adapter trimming and using a sliding window of length ≥50 and quality ≥28, i.e. “ILLUMINACLIP:adapterfna:2:40:15 LEADING:28 TRAILING:28 SLIDINGWINDOW:10:28 MINLEN:50”.

Filtered reads were aligned to the Daucus carota v2.0 genome assembly using the program TopHat v2.0.12 (ref. 30). Non-default parameters used were “--mate-inner-dist -67 --mate-std-dev 50 --min-intron-length 20 --max-intron-length 10000 --library-type fr-unstranded --num-threads 14”. The aligned read files were processed by Cufflinks v2.2.1 (ref. 61). Reads were assembled into transcripts with “cufflinks” using the carrot annotation v1.0 gene predictions as the reference gtf guide. Samples were combined with “cuffmerge”, and then differential expression analyzed with “cuffdiff”, using non-default parameters of “--multi-read-correct --min-alignment-count=5”. Using the abundance estimations, this performs tests for differential expression and regulation between the samples. Normalized counts of the mapped RNA sequences were used to calculate the relative abundances of transcripts expressed as Fragments Per Kilobase of exon per Million fragments mapped (FPKM). When testing for differential expression, biological replicates were included as a term in the mixed model analysis to account for experimental error. Testing for differential expression was done at the level of genes, isoforms, and promoters.

Some of the data from this analysis were published in
Sup. Table 44: Annotation of carrot genes upregulated in both yellow and dark orange (yy) storage roots of plants from mapping populations
and in
Sup. Table 45: Annotation of carrot genes downregulated in yellow or dark orange (yy) storage roots of plants from mapping populations
of the carrot genome assembly publication