Homozygosity mapping identified loci and candidate genes responsible for freezing tolerance in Camelina sativa

Homozygosity mapping is an effective tool for detecting genomic regions responsible for a given trait when the phenotype is controlled by a limited number of dominant or co-dominant loci. Freezing tolerance is a major attribute in agricultural crops such as camelina. Previous studies indicated that freezing tolerance differences between a tolerant (Joelle) and susceptible (CO46) variety of camelina were controlled by a small number of dominant or co-dominant genes. We performed whole genome homozygosity mapping to identify markers and candidate genes responsible for freezing tolerance difference between these two genotypes. A total of 28 F3 RILs were sequenced to similar to 30x coverage, and parental lines were sequenced to >30-40x coverage with Pacific Biosciences high fidelity technology and 60x coverage using Illumina whole genome sequencing. Overall, about 126k homozygous single nucleotide polymorphism markers were identified that differentiate both parents. Moreover, 617 markers were also homozygous in F3 families fixed for freezing tolerance/susceptibility. All these markers mapped to two contigs forming a contiguous stretch of chromosome 11. The homozygosity mapping detected 9 homozygous blocks among the selected markers and 22 candidate genes with strong similarity to regions in or near the homozygous blocks. Two such genes were differentially expressed during cold acclimation in camelina. The largest block contained a cold-regulated plant thionin and a putative rotamase cyclophilin 2 gene previously associated with freezing resistance in arabidopsis (Arabidopsis thaliana). The second largest block contains several cysteine-rich RLK genes and a cold-regulated receptor serine/threonine kinase gene. We hypothesize that one or more of these genes may be primarily responsible for freezing tolerance differences in camelina varieties.