Most fruits in our daily diet are the products of domestication and breeding. Here we report a map of genome variation for a major fruit that encompasses ~3.6 million variants, generated by deep resequencing of 115 cucumber lines sampled from 3,342 accessions worldwide. Comparative analysis suggests that fruit crops underwent narrower bottlenecks during domestication than grain crops. We identified 112 putative domestication sweeps; 1 of these regions contains a gene involved in the loss of bitterness in fruits, an essential domestication trait of cucumber. We also investigated the genomic basis of divergence among the cultivated populations and discovered a natural genetic variant in a β-carotene hydroxylase gene that could be used to breed cucumbers with enhanced nutritional value. The genomic history of cucumber evolution uncovered here provides the basis for future genomics-enabled breeding.

BackgroundTargeting specificity has been a barrier to applying genome editing systems in functional genomics, precise medicine and plant breeding. In plants, only limited studies have used whole-genome sequencing (WGS) to test off-target effects of Cas9. The cause of numerous discovered mutations is still controversial. Furthermore, WGS-based off-target analysis of Cpf1 (Cas12a) has not been reported in any higher organism to date.ResultsWe conduct a WGS analysis of 34 plants edited by Cas9 and 15 plants edited by Cpf1 in T0 and T1 generations along with 20 diverse control plants in rice. The sequencing depths range from 45× to 105× with read mapping rates above 96%. Our results clearly show that most mutations in edited plants are created by the tissue culture process, which causes approximately 102 to 148 single nucleotide variations (SNVs) and approximately 32 to 83 insertions/deletions (indels) per plant. Among 12 Cas9 single guide RNAs (sgRNAs) and three Cpf1 CRISPR RNAs (crRNAs) assessed by WGS, only one Cas9 sgRNA resulted in off-target mutations in T0 lines at sites predicted by computer programs. Moreover, we cannot find evidence for bona fide off-target mutations due to continued expression of Cas9 or Cpf1 with guide RNAs in T1 generation.ConclusionsOur comprehensive and rigorous analysis of WGS data across multiple sample types suggests both Cas9 and Cpf1 nucleases are very specific in generating targeted DNA modifications and off-targeting can be avoided by designing guide RNAs with high specificity.Electronic supplementary materialThe online version of this article (10.1186/s13059-018-1458-5) contains supplementary material, which is available to authorized users.

To gain insight into the evolution of xanthophyll synthesis in Arabidopsis thaliana, we analyzed two pairs of duplicated carotenoid hydroxylase enzymes in Arabidopsis thaliana: the cytochrome P450 enzymes CYP97A3 and CYP97C1, and non-heme di-iron enzymes, BCH1 and BCH2. Hydroxylated carotenes did not accumulate in a quadruple mutant for these four genes, demonstrating that they encode the full complement of carotenoid hydroxylases in A. thaliana. We were thus able to infer definitively the activity of each enzyme in vivo based on the phenotypes of selected double and triple mutant genotypes. The CYP97 and BCH gene pairs are primarily responsible for hydroxylation of alpha- and beta-carotenes, respectively, but exhibit some overlapping activities, most notably in hydroxylation of the beta-ring of alpha-carotene. Surprisingly, triple mutants containing only CYP97C1 or CYP97A3 activity produced 74 and 6% of the wild-type lutein level, indicating that CYP97C1 can efficiently hydroxylate both the beta- and epsilon-rings of alpha-carotene and that CYP97A3 also has low activity toward the epsilon-ring of alpha-carotene. The modes of functional divergence for the gene pairs appear distinct, with the CYP97 duplicates being strongly co-expressed but encoding enzymes with different in vivo substrates, while the BCH duplicates encode isozymes that show significant expression divergence in reproductive organs. By integrating the evolutionary history and substrate specificities of each extant enzyme with the phenotypic responses of various mutant genotypes to high light stress, we propose two likely scenarios for the evolution of alpha-xanthophyll biosynthesis in plants from ancestral organisms.

Hop (Humulus lupulus L. Cannabaceae) is an economically important crop for the brewing industry, where it is used to impart flavor and aroma to beer, and has also drawn attention in recent years due to its potential pharmaceutical applications. Essential oils (mono- and sesquiterpenes), bitter acids (prenylated polyketides), and prenylflavonoids are the primary phytochemical components that account for these traits, and all accumulate at high concentrations in glandular trichomes of hop cones. To understand the molecular basis for terpene accumulation in hop trichomes, a trichome cDNA library was constructed and 9,816 cleansed expressed sequence tag (EST) sequences were obtained from random sequencing of 16,152 cDNA clones. The ESTs were assembled into 3,619 unigenes (1,101 contigs and 2,518 singletons). Putative functions were assigned to the unigenes based on their homology to annotated sequences in the GenBank database. Two mono- and two sesquiterpene synthases identified from the EST collection were expressed in Escherichia coli. Hop MONOTERPENE SYNTHASE2 formed the linear monterpene myrcene from geranyl pyrophosphate, whereas hop SESQUITERPENE SYNTHASE1 (HlSTS1) formed both caryophyllene and humulene from farnesyl pyrophosphate. Together, these enzymes account for the production of the major terpene constituents of the hop trichomes. HlSTS2 formed the minor sesquiterpene constituent germacrene A, which was converted to beta-elemene on chromatography at elevated temperature. We discuss potential functions for other genes expressed at high levels in developing hop trichomes.