Cultivated strawberry emerged from the hybridization of two wild octoploid species, both descendants from the merger of four diploid progenitor species into a single nucleus more than 1 million years ago. Here we report a near-complete chromosome-scale assembly for cultivated octoploid strawberry (Fragaria × ananassa) and uncovered the origin and evolutionary processes that shaped this complex allopolyploid. We identified the extant relatives of each diploid progenitor species and provide support for the North American origin of octoploid strawberry. We examined the dynamics among the four subgenomes in octoploid strawberry and uncovered the presence of a single dominant subgenome with significantly greater gene content, gene expression abundance, and biased exchanges between homoeologous chromosomes, as compared with the other subgenomes. Pathway analysis showed that certain metabolomic and disease-resistance traits are largely controlled by the dominant subgenome. These findings and the reference genome should serve as a powerful platform for future evolutionary studies and enable molecular breeding in strawberry.

BackgroundAlthough draft genomes are available for most agronomically important plant species, the majority are incomplete, highly fragmented, and often riddled with assembly and scaffolding errors. These assembly issues hinder advances in tool development for functional genomics and systems biology.FindingsHere we utilized a robust, cost-effective approach to produce high-quality reference genomes. We report a near-complete genome of diploid woodland strawberry (Fragaria vesca) using single-molecule real-time sequencing from Pacific Biosciences (PacBio). This assembly has a contig N50 length of ∼7.9 million base pairs (Mb), representing a ∼300-fold improvement of the previous version. The vast majority (>99.8%) of the assembly was anchored to 7 pseudomolecules using 2 sets of optical maps from Bionano Genomics. We obtained ∼24.96 Mb of sequence not present in the previous version of the F. vesca genome and produced an improved annotation that includes 1496 new genes. Comparative syntenic analyses uncovered numerous, large-scale scaffolding errors present in each chromosome in the previously published version of the F. vesca genome.ConclusionsOur results highlight the need to improve existing short-read based reference genomes. Furthermore, we demonstrate how genome quality impacts commonly used analyses for addressing both fundamental and applied biological questions.

BackgroundHighbush blueberry (Vaccinium corymbosum) has long been consumed for its unique flavor and composition of health-promoting phytonutrients. However, breeding efforts to improve fruit quality in blueberry have been greatly hampered by the lack of adequate genomic resources and a limited understanding of the underlying genetics encoding key traits. The genome of highbush blueberry has been particularly challenging to assemble due, in large part, to its polyploid nature and genome size.FindingsHere, we present a chromosome-scale and haplotype-phased genome assembly of the cultivar “Draper,” which has the highest antioxidant levels among a diversity panel of 71 cultivars and 13 wild Vaccinium species. We leveraged this genome, combined with gene expression and metabolite data measured across fruit development, to identify candidate genes involved in the biosynthesis of important phytonutrients among other metabolites associated with superior fruit quality. Genome-wide analyses revealed that both polyploidy and tandem gene duplications modified various pathways involved in the biosynthesis of key phytonutrients. Furthermore, gene expression analyses hint at the presence of a spatial-temporal specific dominantly expressed subgenome including during fruit development.ConclusionsThese findings and the reference genome will serve as a valuable resource to guide future genome-enabled breeding of important agronomic traits in highbush blueberry.

QTL analysis in multi-development stages with different QTL models identified 12 consensus QTLs underlying fruit elongation and radial growth presenting a dynamic view of genetic control of cucumber fruit development. Fruit size is an important quality trait in cucumber (Cucumis sativus L.) of different market classes. However, the genetic and molecular basis of fruit size variations in cucumber is not well understood. In this study, we conducted QTL mapping of fruit size in cucumber using F2, F2-derived F3 families and recombinant inbred lines (RILs) from a cross between two inbred lines Gy14 (North American picking cucumber) and 9930 (North China fresh market cucumber). Phenotypic data of fruit length and diameter were collected at three development stages (anthesis, immature and mature fruits) in six environments over 4 years. QTL analysis was performed with three QTL models including composite interval mapping (CIM), Bayesian interval mapping (BIM), and multiple QTL mapping (MQM). Twenty-nine consistent and distinct QTLs were detected for nine traits from multiple mapping populations and QTL models. Synthesis of information from available fruit size QTLs allowed establishment of 12 consensus QTLs underlying fruit elongation and radial growth, which presented a dynamic view of genetic control of cucumber fruit development. Results from this study highlighted the benefits of QTL analysis with multiple QTL models and different mapping populations in improving the power of QTL detection. Discussion was presented in the context of domestication and diversifying selection of fruit length and diameter, marker-assisted selection of fruit size, as well as identification of candidate genes for fruit size QTLs in cucumber.

BackgroundThe fragmented nature of most draft plant genomes has hindered downstream gene discovery, trait mapping for breeding, and other functional genomics applications. There is a pressing need to improve or finish draft plant genome assemblies.FindingsHere, we present a chromosome-scale assembly of the black raspberry genome using single-molecule real-time Pacific Biosciences sequencing and high-throughput chromatin conformation capture (Hi-C) genome scaffolding. The updated V3 assembly has a contig N50 of 5.1 Mb, representing an ∼200-fold improvement over the previous Illumina-based version. Each of the 235 contigs was anchored and oriented into seven chromosomes, correcting several major misassemblies. Black raspberry V3 contains 47 Mb of new sequences including large pericentromeric regions and thousands of previously unannotated protein-coding genes. Among the new genes are hundreds of expanded tandem gene arrays that were collapsed in the Illumina-based assembly. Detailed comparative genomics with the high-quality V4 woodland strawberry genome (Fragaria vesca) revealed near-perfect 1:1 synteny with dramatic divergence in tandem gene array composition. Lineage-specific tandem gene arrays in black raspberry are related to agronomic traits such as disease resistance and secondary metabolite biosynthesis.ConclusionsThe improved resolution of tandem gene arrays highlights the need to reassemble these highly complex and biologically important regions in draft plant genomes. The updated, high-quality black raspberry reference genome will be useful for comparative genomics across the horticulturally important Rosaceae family and enable the development of marker assisted breeding in Rubus.

In humans, the clinical and molecular characterization of sporadic syndromes is often hindered by the small number of patients and the difficulty in developing animal models for severe dominant conditions. Here we show that the availability of large data sets of whole-genome sequences, high-density SNP chip genotypes and extensive recording of phenotype offers an unprecedented opportunity to quickly dissect the genetic architecture of severe dominant conditions in livestock. We report on the identification of seven dominant de novo mutations in CHD7, COL1A1, COL2A1, COPA, and MITF and exploit the structure of cattle populations to describe their clinical consequences and map modifier loci. Moreover, we demonstrate that the emergence of recessive genetic defects can be monitored by detecting de novo deleterious mutations in the genome of bulls used for artificial insemination. These results demonstrate the attractiveness of cattle as a model species in the post genomic era, particularly to confirm the genetic aetiology of isolated clinical case reports in humans.

The oomycete, Phytophthora capsici, infects cucumber (Cucumis sativus L.) fruit. An age-related resistance (ARR) to this pathogen was previously observed in fruit of cultivar ‘Vlaspik’ and shown to be associated with the peel. Young fruits are highly susceptible, but develop resistance at ~10–12 days post pollination (dpp). Peels from resistant (16 dpp) versus susceptible (8 dpp) age fruit are enriched with genes associated with defense, and methanolic extracts from resistant age peels inhibit pathogen growth. Here we compared developing fruits from ‘Vlaspik’ with those of ‘Gy14’, a line that does not exhibit ARR. Transcriptomic analysis of peels of the two lines at 8 and 16 dpp identified 80 genes that were developmentally upregulated in resistant ‘Vlaspik’ 16 dpp versus 8 dpp, but not in susceptible ‘Gy14’ at 16 dpp. A large number of these genes are annotated to be associated with defense and/or specialized metabolism, including four putative resistance (R) genes, and numerous genes involved in flavonoid and terpenoid synthesis and decoration. Untargeted metabolomic analysis was performed on extracts from 8 and 16 dpp ‘Vlaspik’ and ‘Gy14’ fruit peels using Ultra-Performance Liquid Chromatography and Quadrupole Time-of-Flight Mass Spectrometry. Multivariate analysis of the metabolomes identified 113 ions uniquely abundant in resistant ‘Vlaspik’ 16 dpp peel extracts. The most abundant compounds in this group had relative mass defects consistent with terpenoid glycosides. Two of the three most abundant ions were annotated as glycosylated nor-terpenoid esters. Together, these analyses reveal potential mechanisms by which ARR to P. capsici may be conferred.

Very young cucumber (Cucumis sativus) fruit are highly susceptible to infection by the oomycete pathogen, Phytophthora capsici. As the fruit complete exponential growth, at approximately 10–12 days post pollination (dpp), they transition to resistance. The development of age-related resistance (ARR) is increasingly recognized as an important defense against pathogens, however, underlying mechanisms are largely unknown. Peel sections from cucumber fruit harvested at 8 dpp (susceptible) and 16 dpp (resistant) showed equivalent responses to inoculation as did whole fruit, indicating that the fruit surface plays an important role in defense against P. capsici. Exocarp from 16 dpp fruit had thicker cuticles, and methanolic extracts of peel tissue inhibited growth of P. capsici in vitro, suggesting physical or chemical components to the ARR. Transcripts specifically expressed in the peel vs. pericarp showed functional differentiation. Transcripts predominantly expressed in the peel were consistent with fruit surface associated functions including photosynthesis, cuticle production, response to the environment, and defense. Peel-specific transcripts that exhibited increased expression in 16 dpp fruit relative to 8 dpp fruit, were highly enriched (P<0.0001) for response to stress, signal transduction, and extracellular and transport functions. Specific transcripts included genes associated with potential physical barriers (i.e., cuticle), chemical defenses (flavonoid biosynthesis), oxidative stress, penetration defense, and molecular pattern (MAMP)-triggered or effector-triggered (R-gene mediated) pathways. The developmentally regulated changes in gene expression between peels from susceptible- and resistant- age fruits suggest programming for increased defense as the organ reaches full size.