High-Density Linkage Map Construction and QTL Identification in a Diploid Blueberry Mapping Population.

Genotyping by sequencing approaches have been widely applied in major crops and are now being used in horticultural crops like berries and fruit trees. As the original and largest producer of cultivated blueberry, the United States maintains the most diverse blueberry germplasm resources comprised of many species of different ploidy levels. We previously constructed an interspecific mapping population of diploid blueberry by crossing the parent F1#10 (Vaccinium darrowii Fla4B × diploid V. corymbosum W85-20) with the parent W85-23 (diploid V. corymbosum). Employing the Capture-Seq technology developed by RAPiD Genomics, with an emphasis on probes designed in predicted gene regions, 117 F1 progeny, the two parents, and two grandparents of this population were sequenced, yielding 131.7 Gbp clean sequenced reads. A total of 160,535 single nucleotide polymorphisms (SNPs), referenced to 4,522 blueberry genome sequence scaffolds, were identified and subjected to a parent-dependent sliding window approach to further genotype the population. Recombination breakpoints were determined and marker bins were deduced to construct a high density linkage map. Twelve blueberry linkage groups (LGs) consisting of 17,486 SNP markers were obtained, spanning a total genetic distance of 1,539.4 cM. Among 18 horticultural traits phenotyped in this population, quantitative trait loci (QTLs) that were significant over at least 2 years were identified for chilling requirement, cold hardiness, and fruit quality traits of color, scar size, and firmness. Interestingly, in 1 year, a QTL associated with timing of early bloom, full bloom, petal fall, and early green fruit was identified in the same region harboring the major QTL for chilling requirement. In summary, we report here the first high density bin map of a diploid blueberry mapping population and the identification of several horticulturally important QTLs.

The signal crayfish (Pacifastacus leniusculus) in Lake Tahoe (USA) hosts multiple Aphanomyces species.

The genus Aphanomyces (Oomycetes) comprises approximately 50 known species of water molds in three lineages. One of the most notorious is Aphanomyces astaci, the causative agent of crayfish plague. In this study, fresh isolates of Aphanomyces were collected from 20 live specimens of the signal crayfish Pacifastacus leniusculus (Dana, 1852) from Lake Tahoe, California, providing 35 axenic cultures of A. astaci as well as two apparently undescribed Aphanomyces spp. isolates. Based on the results of ITS-, chitinase-, mitochondrial rnnS- and rnnL-sequences and microsatellite markers combined, the Lake Tahoe A. astaci isolates were identical to isolates of A. astaci B-haplogroup commonly detected in Europe, and infection experiments confirmed their high virulence towards noble crayfish. One of the two undescribed Aphanomyces spp. isolates was highly similar to an Aphanomyces lineage detected previously in crustacean zooplankton (Daphnia) in Central Europe, while the other was distinct and most closely related (ITS sequence similarity of 93%) to either A. astaci or to Aphanomyces fennicus isolated recently from Astacus astacus in Finland. Neither of the two Aphanomyces spp. isolates caused crayfish mortality under experimental conditions. Our results indicate that the populations of North American signal crayfish can act as carriers of both pathogenic and non-pathogenic Aphanomyces at the same time. Furthermore, considering that a limited number of crayfish individuals from a single location yielded multiple distinct Aphanomyces isolates, our results suggest that substantial species diversity within this genus remains undescribed.

The genome sequence and transcriptome of Potentilla micrantha and their comparison to Fragaria vesca (the woodland strawberry).

Background: The genus Potentilla is closely related to that of Fragaria, the economically important strawberry genus. Potentilla micrantha is a species that does not develop berries but shares numerous morphological and ecological characteristics with Fragaria vesca. These similarities make P. micrantha an attractive choice for comparative genomics studies with F. vesca. Findings: In this study, the P. micrantha genome was sequenced and annotated, and RNA-Seq data from the different developmental stages of flowering and fruiting were used to develop a set of gene predictions. A 327 Mbp sequence and annotation of the genome of P. micrantha, spanning 2674 sequence contigs, with an N50 size of 335,712, estimated to cover 80% of the total genome size of the species was developed. The genus Potentilla has a characteristically larger genome size than Fragaria, but the recovered sequence scaffolds were remarkably collinear at the micro-syntenic level with the genome of F. vesca, its closest sequenced relative. A total of 33,602 genes were predicted, and 95.1% of bench-marking universal single-copy orthologous genes were complete within the presented sequence. Thus, we argue that the majority of the gene-rich regions of the genome have been sequenced. Conclusions: Comparisons of RNA-Seq data from the stages of floral and fruit development revealed genes differentially expressed between P. micrantha and F. vesca.The data presented are a valuable resource for future studies of berry development in Fragaria and the Rosaceae and they also shed light on the evolution of genome size and organization in this family.

Soil Microbiomes Associated with Verticillium Wilt-Suppressive Broccoli and Chitin Amendments are Enriched with Potential Biocontrol Agents.

Two naturally infested Verticillium wilt-conducive soils from the Salinas Valley of coastal California were amended with disease-suppressive broccoli residue or crab meal amendments, and changes to the soil prokaryote community were monitored using Illumina sequencing of a 16S ribosomal RNA gene library generated from 160 bulk soil samples. The experiment was run in a greenhouse, twice, with eggplant as the Verticillium wilt-susceptible host. Disease suppression, plant height, soil microsclerotia density, and soil chitinase activity were assessed at the conclusion of each experiment. In soil with high microsclerotia density, all amendments significantly reduced Verticillium wilt severity and microsclerotia density, and increased soil chitinase activity. Plant height was increased only in the broccoli-containing treatments. In total, 8,790 error-corrected sequence variants representing 1,917,893 different sequences were included in the analyses. The treatments had a significant impact on the soil microbiome community structure but measures of α diversity did not vary between treatments. Community structure correlated with disease score, plant height, microsclerotia density, and soil chitinase activity, suggesting that the prokaryote community may affect the disease-related response variables or vice versa. Similarly, the abundance of 107 sequence variants correlated with disease-related response variables, which included variants from genera with known antagonists of filamentous fungal plant pathogens, such as Pseudomonas and Streptomyces. Overall, genera with antifungal antagonists were more abundant in amended soils than unamended soils, and constituted up to 8.9% of all sequences in broccoli+crabmeal-amended soil. This study demonstrates that substrate-mediated shifts in soil prokaryote communities are associated with the transition of Verticillium wilt-conducive soils to Verticillium wilt-suppressive soils, and suggests that soils likely harbor numerous additional antagonists of fungal plant pathogens that contribute to the biological suppression of plant disease.

A High-Density Linkage Map of the Ancestral Diploid Strawberry, , Constructed with Single Nucleotide Polymorphism Markers from the IStraw90 Array and Genotyping by Sequencing.

Makino is recognized as an ancestor of the octoploid strawberry species, which includes the cultivated strawberry, × Duchesne ex Rozier. Here we report the construction of the first high-density linkage map for . The linkage map (Fii map) is based on two high-throughput techniques of single nucleotide polymorphism (SNP) genotyping: the IStraw90 Array (hereafter "Array"), and genotyping by sequencing (GBS). The F generation mapping population was derived by selfing hybrid F1D, the product of a cross between two divergent accessions collected from Hokkaido, Japan. The Fii map consists of seven linkage groups (LGs) and has an overall length of 451.7 cM as defined by 496 loci populated by 4173 markers: 3280 from the Array and 893 from GBS. Comparisons with two versions of the ssp. L. 'Hawaii 4' pseudo-chromosome (PC) assemblies reveal substantial conservation of synteny and colinearity, yet identified differences that point to possible genomic divergences between and , and/or to genomic assembly errors. The Fii map provides a basis for anchoring a genome assembly as a prerequisite for constructing a second diploid reference genome for .

Molecular genetics and genomics of the Rosoideae: state of the art and future perspectives.

The Rosoideae is a subfamily of the Rosaceae that contains a number of species of economic importance, including the soft fruit species strawberry (Fragaria ×ananassa), red (Rubus idaeus) and black (Rubus occidentalis) raspberries, blackberries (Rubus spp.) and one of the most economically important cut flower genera, the roses (Rosa spp.). Molecular genetics and genomics resources for the Rosoideae have developed rapidly over the past two decades, beginning with the development and application of a number of molecular marker types including restriction fragment length polymorphisms, amplified fragment length polymorphisms and microsatellites, and culminating in the recent publication of the genome sequence of the woodland strawberry, Fragaria vesca, and the development of high throughput single nucleotide polymorphism (SNP)-genotyping resources for Fragaria, Rosa and Rubus. These tools have been used to identify genes and other functional elements that control traits of economic importance, to study the evolution of plant genome structure within the subfamily, and are beginning to facilitate genomic-assisted breeding through the development and deployment of markers linked to traits such as aspects of fruit quality, disease resistance and the timing of flowering. In this review, we report on the developments that have been made over the last 20 years in the field of molecular genetics and structural genomics within the Rosoideae, comment on how the knowledge gained will improve the efficiency of cultivar development and discuss how these advances will enhance our understanding of the biological processes determining agronomically important traits in all Rosoideae species.

An evaluation of the PacBio RS platform for sequencing and de novo assembly of a chloroplast genome.

BACKGROUND: Second generation sequencing has permitted detailed sequence characterisation at the whole genome level of a growing number of non-model organisms, but the data produced have short read-lengths and biased genome coverage leading to fragmented genome assemblies. The PacBio RS long-read sequencing platform offers the promise of increased read length and unbiased genome coverage and thus the potential to produce genome sequence data of a finished quality containing fewer gaps and longer contigs. However, these advantages come at a much greater cost per nucleotide and with a perceived increase in error-rate. In this investigation, we evaluated the performance of the PacBio RS sequencing platform through the sequencing and de novo assembly of the Potentilla micrantha chloroplast genome. RESULTS: Following error-correction, a total of 28,638 PacBio RS reads were recovered with a mean read length of 1,902 bp totalling 54,492,250 nucleotides and representing an average depth of coverage of 320× the chloroplast genome. The dataset covered the entire 154,959 bp of the chloroplast genome in a single contig (100% coverage) compared to seven contigs (90.59% coverage) recovered from an Illumina data, and revealed no bias in coverage of GC rich regions. Post-assembly the data were largely concordant with the Illumina data generated and allowed 187 ambiguities in the Illumina data to be resolved. The additional read length also permitted small differences in the two inverted repeat regions to be assigned unambiguously. CONCLUSIONS: This is the first report to our knowledge of a chloroplast genome assembled de novo using PacBio sequence data. The PacBio RS data generated here were assembled into a single large contig spanning the P. micrantha chloroplast genome, with a higher degree of accuracy than an Illumina dataset generated at a much greater depth of coverage, due to longer read lengths and lower GC bias in the data. The results we present suggest PacBio data will be of immense utility for the development of genome sequence assemblies containing fewer unresolved gaps and ambiguities and a significantly smaller number of contigs than could be produced using short-read sequence data alone.

Saturated linkage map construction in Rubus idaeus using genotyping by sequencing and genome-independent imputation.

BACKGROUND: Rapid development of highly saturated genetic maps aids molecular breeding, which can accelerate gain per breeding cycle in woody perennial plants such as Rubus idaeus (red raspberry). Recently, robust genotyping methods based on high-throughput sequencing were developed, which provide high marker density, but result in some genotype errors and a large number of missing genotype values. Imputation can reduce the number of missing values and can correct genotyping errors, but current methods of imputation require a reference genome and thus are not an option for most species. RESULTS: Genotyping by Sequencing (GBS) was used to produce highly saturated maps for a R. idaeus pseudo-testcross progeny. While low coverage and high variance in sequencing resulted in a large number of missing values for some individuals, a novel method of imputation based on maximum likelihood marker ordering from initial marker segregation overcame the challenge of missing values, and made map construction computationally tractable. The two resulting parental maps contained 4521 and 2391 molecular markers spanning 462.7 and 376.6 cM respectively over seven linkage groups. Detection of precise genomic regions with segregation distortion was possible because of map saturation. Microsatellites (SSRs) linked these results to published maps for cross-validation and map comparison. CONCLUSIONS: GBS together with genome-independent imputation provides a rapid method for genetic map construction in any pseudo-testcross progeny. Our method of imputation estimates the correct genotype call of missing values and corrects genotyping errors that lead to inflated map size and reduced precision in marker placement. Comparison of SSRs to published R. idaeus maps showed that the linkage maps constructed with GBS and our method of imputation were robust, and marker positioning reliable. The high marker density allowed identification of genomic regions with segregation distortion in R. idaeus, which may help to identify deleterious alleles that are the basis of inbreeding depression in the species.

Strategies for transcriptome analysis in nonmodel plants.

Even with recent reductions in sequencing costs, most plants lack the genomic resources required for successful short-read transcriptome analyses as performed routinely in model species. Several approaches for the analysis of short-read transcriptome data are reviewed for nonmodel species for which the genome of a close relative is used as the reference genome. Two approaches using a data set from Phytophthora-challenged Rubus idaeus (red raspberry) are compared. Over 70000000 86-nt Illumina reads derived from R. idaeus roots were aligned to the Fragaria vesca genome using publicly available informatics tools (Bowtie/TopHat and Cufflinks). Alignment identified 16956 putatively expressed genes. De novo assembly was performed with the same data set and a publicly available transcriptome assembler (Trinity). A BLAST search with a maximum e-value threshold of 1.0 × 10(-3) revealed that over 36000 transcripts had matches to plants and over 500 to Phytophthora. Gene expression estimates from alignment to F. vesca and de novo assembly were compared for raspberry (Pearson's correlation = 0.730). Together, alignment to the genome of a close relative and de novo assembly constitute a powerful method of transcriptome analysis in nonmodel organisms. Alignment to the genome of a close relative provides a framework for differential expression testing if alignments are made to the predefined gene-space of a close relative and de novo assembly provides a more robust method of identifying unique sequences and sequences from other organisms in a system. These methods are considered experimental in nonmodel systems, but can be used to generate resources and specific testable hypotheses.