Transposons played a major role in the diversification between the closely related almond and peach genomes: results from the almond genome sequence.

We sequenced the genome of the highly heterozygous almond Prunus dulcis cv. Texas combining short- and long-read sequencing. We obtained a genome assembly totaling 227.6 Mb of the estimated almond genome size of 238 Mb, of which 91% is anchored to eight pseudomolecules corresponding to its haploid chromosome complement, and annotated 27 969 protein-coding genes and 6747 non-coding transcripts. By phylogenomic comparison with the genomes of 16 additional close and distant species we estimated that almond and peach (Prunus persica) diverged around 5.88 million years ago. These two genomes are highly syntenic and show a high degree of sequence conservation (20 nucleotide substitutions per kb). However, they also exhibit a high number of presence/absence variants, many attributable to the movement of transposable elements (TEs). Transposable elements have generated an important number of presence/absence variants between almond and peach, and we show that the recent history of TE movement seems markedly different between them. Transposable elements may also be at the origin of important phenotypic differences between both species, and in particular for the sweet kernel phenotype, a key agronomic and domestication character for almond. Here we show that in sweet almond cultivars, highly methylated TE insertions surround a gene involved in the biosynthesis of amygdalin, whose reduced expression has been correlated with the sweet almond phenotype. Altogether, our results suggest a key role of TEs in the recent history and diversification of almond and its close relative peach.

Genome editing reagent delivery in plants.

Genome editing holds the potential for rapid crop improvement to meet the challenge of feeding the planet in a changing climate. The delivery of gene editing reagents into the plant cells has been dominated by plasmid vectors delivered using agrobacterium or particle bombardment. This approach involves the production of genetically engineered plants, which need to undergo regulatory approvals. There are various reagent delivery approaches available that have enabled the delivery of DNA-free editing reagents. They invariably involve the use of ribonucleoproteins (RNPs), especially in the case of CRISPR/Cas9-mediated gene editing. The explant of choice for most of the non-DNA approaches utilizes protoplasts as the recipient explant. While the editing efficiency is high in protoplasts, the ability to regenerate individual plants from edited protoplasts remains a challenge. There are various innovative delivery approaches being utilized to perform in planta edits that can be incorporated in the germline cells or inherited via seed. With the modification and adoption of various novel approaches currently being used in animal systems, it seems likely that non-transgenic genome editing will become routine in higher plants.

Genome editing in fruit, ornamental, and industrial crops.

The advent of genome editing has opened new avenues for targeted trait enhancement in fruit, ornamental, industrial, and all specialty crops. In particular, CRISPR-based editing systems, derived from bacterial immune systems, have quickly become routinely used tools for research groups across the world seeking to edit plant genomes with a greater level of precision, higher efficiency, reduced off-target effects, and overall ease-of-use compared to ZFNs and TALENs. CRISPR systems have been applied successfully to a number of horticultural and industrial crops to enhance fruit ripening, increase stress tolerance, modify plant architecture, control the timing of flower development, and enhance the accumulation of desired metabolites, among other commercially-important traits. As editing technologies continue to advance, so too does the ability to generate improved crop varieties with non-transgenic modifications; in some crops, direct transgene-free edits have already been achieved, while in others, T-DNAs have successfully been segregated out through crossing. In addition to the potential to produce non-transgenic edited crops, and thereby circumvent regulatory impediments to the release of new, improved crop varieties, targeted gene editing can speed up trait improvement in crops with long juvenile phases, reducing inputs resulting in faster market introduction to the market. While many challenges remain regarding optimization of genome editing in ornamental, fruit, and industrial crops, the ongoing discovery of novel nucleases with niche specialties for engineering applications may form the basis for additional and potentially crop-specific editing strategies.

Fruit crops in the era of genome editing: closing the regulatory gap.

The conventional breeding of fruits and fruit trees has led to the improvement of consumer-driven traits such as fruit size, yield, nutritional properties, aroma and taste, as well as the introduction of agronomic properties such as disease resistance. However, even with the assistance of modern molecular approaches such as marker-assisted selection, the improvement of fruit varieties by conventional breeding takes considerable time and effort. The advent of genetic engineering led to the rapid development of new varieties by allowing the direct introduction of genes into elite lines. In this review article, we discuss three such case studies: the Arctic® apple, the Pinkglow pineapple and the SunUp/Rainbow papaya. We consider these events in the light of global regulations for the commercialization of genetically modified organisms (GMOs), focusing on the differences between product-related systems (the USA/Canada comparative safety assessment) and process-related systems (the EU "precautionary principle" model). More recently, genome editing has provided an efficient way to introduce precise mutations in plants, including fruits and fruit trees, replicating conventional breeding outcomes without the extensive backcrossing and selection typically necessary to introgress new traits. Some jurisdictions have reacted by amending the regulations governing GMOs to provide exemptions for crops that would be indistinguishable from conventional varieties based on product comparison. This has revealed the deficiencies of current process-related regulatory frameworks, particularly in the EU, which now stands against the rest of the world as a unique example of inflexible and dogmatic governance based on political expediency and activism rather than rigorous scientific evidence.

CRISPR-associated nucleases: the Dawn of a new age of efficient crop improvement.

The world stands at a new threshold today. As a planet, we face various challenges, and the key one is how to continue to produce enough food, feed, fiber, and fuel to support the burgeoning population. In the past, plant breeding and the ability to genetically engineer crops contributed to increasing food production. However, both approaches rely on random mixing or integration of genes, and the process can be unpredictable and time-consuming. Given the challenge of limited availability of natural resources and changing environmental conditions, the need to rapidly and precisely improve crops has become urgent. The discovery of CRISPR-associated endonucleases offers a precise yet versatile platform for rapid crop improvement. This review summarizes a brief history of the discovery of CRISPR-associated nucleases and their application in genome editing of various plant species. Also provided is an overview of several new endonucleases reported recently, which can be utilized for editing of specific genes in plants through various forms of DNA sequence alteration. Genome editing, with its ever-expanding toolset, increased efficiency, and its potential integration with the emerging synthetic biology approaches hold promise for efficient crop improvement to meet the challenge of supporting the needs of future generations.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization.

BACKGROUND: Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies. RESULTS: By applying a two-step design process, which consisted of the construction of a first 'draft' array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping. CONCLUSIONS: The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Neuropeptides and peptide hormones identified in codling moth, Cydia pomonella (Lepidoptera: Tortricidae).

The codling moth, Cydia pomonella, is a worldwide pest of pome fruits. Neuropeptides regulate most physiological functions in insects and represent new targets for the development of control agents. The only neuropeptides reported from the codling moth to date are the allatostatin A family peptides. To identify other neuropeptides and peptide hormones from codling moth, we analyzed head transcriptomes, identified 50 transcripts, and predicted 120 prepropeptides for the codling moth neuropeptides and peptide hormones. All transcripts were amplified, and these sequences were verified. One of the notable findings in this study is that diapause hormones (DHs) reported from Tortricid moths, including the codling moth, do not have the WFGPRL sequence in C-terminal ends in the pban genes. The C-terminal motif is critical to characterize insect DH peptides, and always conserved in pban/dh genes in Lepidoptera and many insect orders. Interestingly, the WFGPRL sequence was produced only from the capa gene in the codling moth. The allatostatin A-family encoding transcript predicted nine peptides, seven of which, as expected, are identical to those previously isolated from the moth. We also identified new codling moth orthologs of insect neuropeptides including CCHamides, allatostatin CC, RYamides, and natalisins. The information provided in this study will benefit future codling moth investigations using peptidoproteomics to determine peptide presence and functions.

Comparative ultrastructure of fruit plastids in three genetically diverse genotypes of apple (Malus × domestica Borkh.) during development.

KEY MESSAGE: Comparative ultrastructural developmental time-course analysis has identified discrete stages at which the fruit plastids undergo structural and consequently functional transitions to facilitate subsequent development-guided understanding of the complex plastid biology. Plastids are the defining organelle for a plant cell and are critical for myriad metabolic functions. The role of leaf plastid, chloroplast, is extensively documented; however, fruit plastids-chromoplasts-are poorly understood, especially in the context of the diverse metabolic processes operating in these diverse plant organs. Recently, in a comparative study of the predicted plastid-targeted proteomes across seven plant species, we reported that each plant species is predicted to harbor a unique set of plastid-targeted proteins. However, the temporal and developmental context of these processes remains unknown. In this study, an ultrastructural analysis approach was used to characterize fruit plastids in the epidermal and collenchymal cell layers at 11 developmental timepoints in three genotypes of apple (Malus × domestica Borkh.): chlorophyll-predominant 'Granny Smith', carotenoid-predominant 'Golden Delicious', and anthocyanin-predominant 'Top Red Delicious'. Plastids transitioned from a proplastid-like plastid to a chromoplast-like plastid in epidermis cells, while in the collenchyma cells, they transitioned from a chloroplast-like plastid to a chloro-chromo-amyloplast plastid. Plastids in the collenchyma cells of the three genotypes demonstrated a diverse array of structures and features. This study enabled the identification of discrete developmental stages during which specific functions are most likely being performed by the plastids as indicated by accumulation of plastoglobuli, starch granules, and other sub-organeller structures. Information regarding the metabolically active developmental stages is expected to facilitate biologically relevant omics studies to unravel the complex biochemistry of plastids in perennial non-model systems.

Thioredoxin-insensitive plastid ATP synthase that performs moonlighting functions.

The chloroplast ATP synthase catalyzes the light-driven synthesis of ATP and acts as a key feedback regulatory component of photosynthesis. Arabidopsis possesses two homologues of the regulatory γ subunit of the ATP synthase, encoded by the ATPC1 and ATPC2 genes. Using a series of mutants, we show that both these subunits can support photosynthetic ATP synthesis in vivo with similar specific activities, but that in wild-type plants, only γ(1) is involved in ATP synthesis in photosynthesis. The γ(1)-containing ATP synthase shows classical light-induced redox regulation, whereas the mutant expressing only γ(2)-ATP synthase (gamma exchange-revised ATP synthase, gamera) shows equally high ATP synthase activity in the light and dark. In situ redox titrations demonstrate that the regulatory thiol groups on γ(2)-ATP synthase remain reduced under physiological conditions but can be oxidized by the strong oxidant diamide, implying that the redox potential for the thiol/disulphide transition in γ(2) is substantially higher than that for γ(1). This regulatory difference may be attributed to alterations in the residues near the redox-active thiols. We propose that γ(2)-ATP synthase functions to catalyze ATP hydrolysis-driven proton translocation in nonphotosynthetic plastids, maintaining a sufficient transthylakoid proton gradient to drive protein translocation or other processes. Consistent with this interpretation, ATPC2 is predominantly expressed in the root, whereas modifying its expression results in alteration of root hair development. Phylogenetic analysis suggests that γ(2) originated from ancient gene duplication, resulting in divergent evolution of functionally distinct ATP synthase complexes in dicots and mosses.

Evaluation of Parameters Affecting Agrobacterium-Mediated Transient Gene Expression in Industrial Hemp (Cannabis sativa L.).

Industrial hemp Cannabis sativa L. is an economically important crop mostly grown for its fiber, oil, and seeds. Due to its increasing applications in the pharmaceutical industry and a lack of knowledge of gene functions in cannabinoid biosynthesis pathways, developing an efficient transformation platform for the genetic engineering of industrial hemp has become necessary to enable functional genomic and industrial application studies. A critical step in the development of Agrobacterium tumefaciens-mediated transformation in the hemp genus is the establishment of optimal conditions for T-DNA gene delivery into different explants from which whole plantlets can be regenerated. As a first step in the development of a successful Agrobacterium tumefaciens-mediated transformation method for hemp gene editing, the factors influencing the successful T-DNA integration and expression (as measured by transient ß-glucuronidase (GUS) and Green Florescent Protein (GFP) expression) were investigated. In this study, the parameters for an agroinfiltration system in hemp, which applies to the stable transformation method, were optimized. In the present study, we tested different explants, such as 1- to 3-week-old leaves, cotyledons, hypocotyls, root segments, nodal parts, and 2- to 3-week-old leaf-derived calli. We observed that the 3-week-old leaves were the best explant for transient gene expression. Fully expanded 2- to 3-week-old leaf explants, in combination with 30 min of immersion time, 60 µM silver nitrate, 0.5 µM calcium chloride, 150 µM natural phenolic compound acetosyringone, and a bacterial density of OD600nm = 0.4 resulted in the highest GUS and GFP expression. The improved method of genetic transformation established in the present study will be useful for the introduction of foreign genes of interest, using the latest technologies such as genome editing, and studying gene functions that regulate secondary metabolites in hemp.

Comparative genomics analysis in Prunoideae to identify biologically relevant polymorphisms.

Prunus is an economically important genus with a wide range of physiological and biological variability. Using the peach genome as a reference, sequencing reads from four almond accessions and one sweet cherry cultivar were used for comparative analysis of these three Prunus species. Reference mapping enabled the identification of many biological relevant polymorphisms within the individuals. Examining the depth of the polymorphisms and the overall scaffold coverage, we identified many potentially interesting regions including hundreds of small scaffolds with no coverage from any individual. Non-sense mutations account for about 70 000 of the 13 million identified single nucleotide polymorphisms (SNPs). Blast2GO analyses on these non-sense SNPs revealed several interesting results. First, non-sense SNPs were not evenly distributed across all gene ontology terms. Specifically, in comparison with peach, sweet cherry is found to have non-sense SNPs in two 1-aminocyclopropane-1-carboxylate synthase (ACS) genes and two 1-aminocyclopropane-1-carboxylate oxidase (ACO) genes. These polymorphisms may be at the root of the nonclimacteric ripening of sweet cherry. A set of candidate genes associated with bitterness in almond were identified by comparing sweet and bitter almond sequences. To the best of our knowledge, this is the first report in plants of non-sense SNP abundance in a genus being linked to specific GO terms.

Role of bacterial communities in the natural suppression of Rhizoctonia solani bare patch disease of wheat (Triticum aestivum L.).

Rhizoctonia bare patch and root rot disease of wheat, caused by Rhizoctonia solani AG-8, develops as distinct patches of stunted plants and limits the yield of direct-seeded (no-till) wheat in the Pacific Northwest of the United States. At the site of a long-term cropping systems study near Ritzville, WA, a decline in Rhizoctonia patch disease was observed over an 11-year period. Bacterial communities from bulk and rhizosphere soil of plants from inside the patches, outside the patches, and recovered patches were analyzed by using pyrosequencing with primers designed for 16S rRNA. Taxa in the class Acidobacteria and the genus Gemmatimonas were found at higher frequencies in the rhizosphere of healthy plants outside the patches than in that of diseased plants from inside the patches. Dyella and Acidobacteria subgroup Gp7 were found at higher frequencies in recovered patches. Chitinophaga, Pedobacter, Oxalobacteriaceae (Duganella and Massilia), and Chyseobacterium were found at higher frequencies in the rhizosphere of diseased plants from inside the patches. For selected taxa, trends were validated by quantitative PCR (qPCR), and observed shifts of frequencies in the rhizosphere over time were duplicated in cycling experiments in the greenhouse that involved successive plantings of wheat in Rhizoctonia-inoculated soil. Chryseobacterium soldanellicola was isolated from the rhizosphere inside the patches and exhibited significant antagonism against R. solani AG-8 in vitro and in greenhouse tests. In conclusion, we identified novel bacterial taxa that respond to conditions affecting bare patch disease symptoms and that may be involved in suppression of Rhizoctonia root rot and bare batch disease.

An Arabidopsis mutant with high cyclic electron flow around photosystem I (hcef) involving the NADPH dehydrogenase complex.

Cyclic electron flow (CEFI) has been proposed to balance the chloroplast energy budget, but the pathway, mechanism, and physiological role remain unclear. We isolated a new class of mutant in Arabidopsis thaliana, hcef for high CEF1, which shows constitutively elevated CEF1. The first of these, hcef1, was mapped to chloroplast fructose-1,6-bisphosphatase. Crossing hcef1 with pgr5, which is deficient in the antimycin A-sensitive pathway for plastoquinone reduction, resulted in a double mutant that maintained the high CEF1 phenotype, implying that the PGR5-dependent pathway is not involved. By contrast, crossing hcef1 with crr2-2, deficient in thylakoid NADPH dehydrogenase (NDH) complex, results in a double mutant that is highly light sensitive and lacks elevated CEF1, suggesting that NDH plays a direct role in catalyzing or regulating CEF1. Additionally, the NdhI component of the NDH complex was highly expressed in hcef1, whereas other photosynthetic complexes, as well as PGR5, decreased. We propose that (1) NDH is specifically upregulated in hcef1, allowing for increased CEF1; (2) the hcef1 mutation imposes an elevated ATP demand that may trigger CEF1; and (3) alternative mechanisms for augmenting ATP cannot compensate for the loss of CEF1 through NDH.

Antisense expression of peach mildew resistance locus O (PpMlo1) gene confers cross-species resistance to powdery mildew in Fragaria x ananassa.

Powdery mildew (PM) is one of the major plant pathogens. The conventional method of PM control includes frequent use of sulfur-based fungicides adding to production costs and potential harm to the environment. PM remains a major scourge for Rosaceae crops where breeding approaches mainly resort to gene-for-gene resistance. We have tested an alternate source of PM resistance in Rosaceae. Mildew resistance locus O (MLO) has been well studied in barley due to its role in imparting broad spectrum resistance to PM. We identified PpMlo1 (Prunus persica Mlo) in peach and characterized it further to test if a similar mechanism of resistance is conserved in Rosaceae. Due to its recalcitrance in tissue culture, reverse genetic studies involving PpMloI were not feasible in peach. Therefore, Fragaria x ananassa LF9 line, a taxonomic surrogate, was used for functional analysis of PpMlo1. Agrobacterium-mediated transformation yielded transgenic strawberry plants expressing PpMlo1 in sense and antisense orientation. Antisense expression of PpMlo1 in transgenic strawberry plants conferred resistance to Fragaria-specific powdery mildew, Podosphaera macularis. Phylogenetic analysis of 208 putative Mlo gene copies from 35 plant species suggests a large number of duplications of this gene family prior to the divergence of monocots and eudicots, early in eudicot diversification. Our results indicate that the Mlo-based resistance mechanism is functional in Rosaceae, and that Fragaria can be used as a host to test mechanistic function of genes derived from related tree species. To the best of our knowledge, this work is one of the first attempts at testing the potential of using a Mlo-based resistance strategy to combat powdery mildew in Rosaceae.

Rootstock scion somatogenetic interactions in perennial composite plants.

The ancient plant production practice of grafting which instantly imparts new physiological properties to the desirable scion still remains shrouded in mystery. Yet, grafting remains a widely used technique in the production of several horticultural species. In a composite grafted plant, rootstocks control many aspects of scion growth and physiology including yield and quality attributes as well as biotic and abiotic stress tolerance. Broadly, physical, physiological, biochemical and molecular mechanisms have been reviewed to develop an integrated understanding of this enigmatic process that challenges existing genetic paradigms. This review summarizes the reported mechanisms underlying some of the economically important traits and identifies several key points to consider when conducting rootstock scion interaction experiments. Study of the somatogenetic interactions between rootstock and scion is a field that is ripe for discovery and vast improvements in the coming decade. Further, utilization of rootstocks based on a better understanding of the somatogenetic interactions is highly relevant in the current agricultural environment where there is a need for sustainable production practices. Rootstocks may offer a non-transgenic approach to rapidly respond to the changing environment and expand agricultural production of annual and perennial crops where grafting is feasible in order to meet the global food, fiber and fuel demands of the future.

Impact of heat stress, water stress, and their combined effects on the metabolism and transcriptome of grape berries.

Recurring heat and drought episodes present challenges to the sustainability of grape production worldwide. We investigated the impacts of heat and drought stress on transcriptomic and metabolic responses of berries from two wine grape varieties. Cabernet Sauvignon and Riesling grapevines were subjected to one of four treatments during early fruit ripening: (1) drought stress only, (2) heat stress only, (3) simultaneous drought and heat stress, (4) no drought or heat stress (control). Berry metabolites, especially organic acids, were analyzed, and time-course transcriptome analysis was performed on samples before, during, and after the stress episode. Both alone and in conjunction with water stress, heat stress had a much more significant impact on berry organic acid content, pH, and titratable acidity than water stress. This observation contrasts with previous reports for leaves, which responded more strongly to water stress, indicating that grape berries display a distinct, organ-specific response to environmental stresses. Consistent with the metabolic changes, the global transcriptomic analysis revealed that heat stress had a more significant impact on gene expression in grape berries than water stress in both varieties. The differentially expressed genes were those associated with the tricarboxylic acid cycle and glyoxylate cycle, mitochondrial electron transport and alternative respiration, glycolysis and gluconeogenesis, carbohydrate allocation, ascorbate metabolism, and abiotic stress signaling pathways. Knowledge regarding how environmental stresses, alone and in combination, impact the berry metabolism of different grape varieties will form the basis for developing recommendations for climate change mitigation strategies and genetic improvement.

Understanding protein import in diverse non-green plastids.

The spectacular diversity of plastids in non-green organs such as flowers, fruits, roots, tubers, and senescing leaves represents a Universe of metabolic processes in higher plants that remain to be completely characterized. The endosymbiosis of the plastid and the subsequent export of the ancestral cyanobacterial genome to the nuclear genome, and adaptation of the plants to all types of environments has resulted in the emergence of diverse and a highly orchestrated metabolism across the plant kingdom that is entirely reliant on a complex protein import and translocation system. The TOC and TIC translocons, critical for importing nuclear-encoded proteins into the plastid stroma, remain poorly resolved, especially in the case of TIC. From the stroma, three core pathways (cpTat, cpSec, and cpSRP) may localize imported proteins to the thylakoid. Non-canonical routes only utilizing TOC also exist for the insertion of many inner and outer membrane proteins, or in the case of some modified proteins, a vesicular import route. Understanding this complex protein import system is further compounded by the highly heterogeneous nature of transit peptides, and the varying transit peptide specificity of plastids depending on species and the developmental and trophic stage of the plant organs. Computational tools provide an increasingly sophisticated means of predicting protein import into highly diverse non-green plastids across higher plants, which need to be validated using proteomics and metabolic approaches. The myriad plastid functions enable higher plants to interact and respond to all kinds of environments. Unraveling the diversity of non-green plastid functions across the higher plants has the potential to provide knowledge that will help in developing climate resilient crops.

Development and Evaluation of an AxiomTM 60K SNP Array for Almond (Prunus dulcis).

A high-density single nucleotide polymorphism (SNP) array is essential to enable faster progress in plant breeding for new cultivar development. In this regard, we have developed an Axiom 60K almond SNP array by resequencing 81 almond accessions. For the validation of the array, a set of 210 accessions were genotyped and 82.8% of the SNPs were classified in the best recommended SNPs. The rate of missing data was between 0.4% and 2.7% for the almond accessions and less than 15.5% for the few peach and wild accessions, suggesting that this array can be used for peach and interspecific peach × almond genetic studies. The values of the two SNPs linked to the RMja (nematode resistance) and SK (bitterness) genes were consistent. We also genotyped 49 hybrids from an almond F2 progeny and could build a genetic map with a set of 1159 SNPs. Error rates, less than 1%, were evaluated by comparing replicates and by detection of departures from Mendelian inheritance in the F2 progeny. This almond array is commercially available and should be a cost-effective genotyping tool useful in the search for new genes and quantitative traits loci (QTL) involved in the control of agronomic traits.

Rapid gene-based SNP and haplotype marker development in non-model eukaryotes using 3'UTR sequencing.

BACKGROUND: Sweet cherry (Prunus avium L.), a non-model crop with narrow genetic diversity, is an important member of sub-family Amygdoloideae within Rosaceae. Compared to other important members like peach and apple, sweet cherry lacks in genetic and genomic information, impeding understanding of important biological processes and development of efficient breeding approaches. Availability of single nucleotide polymorphism (SNP)-based molecular markers can greatly benefit breeding efforts in such non-model species. RNA-seq approaches employing second generation sequencing platforms offer a unique avenue to rapidly identify gene-based SNPs. Additionally, haplotype markers can be rapidly generated from transcript-based SNPs since they have been found to be extremely utile in identification of genetic variants related to health, disease and response to environment as highlighted by the human HapMap project. RESULTS: RNA-seq was performed on two sweet cherry cultivars, Bing and Rainier using a 3' untranslated region (UTR) sequencing method yielding 43,396 assembled contigs. In order to test our approach of rapid identification of SNPs without any reference genome information, over 25% (10,100) of the contigs were screened for the SNPs. A total of 207 contigs from this set were identified to contain high quality SNPs. A set of 223 primer pairs were designed to amplify SNP containing regions from these contigs and high resolution melting (HRM) analysis was performed with eight important parental sweet cherry cultivars. Six of the parent cultivars were distantly related to Bing and Rainier, the cultivars used for initial SNP discovery. Further, HRM analysis was also performed on 13 seedlings derived from a cross between two of the parents. Our analysis resulted in the identification of 84 (38.7%) primer sets that demonstrated variation among the tested germplasm. Reassembly of the raw 3'UTR sequences using upgraded transcriptome assembly software yielded 34,620 contigs containing 2243 putative SNPs in 887 contigs after stringent filtering. Contigs with multiple SNPs were visually parsed to identify 685 putative haplotypes at 335 loci in 301 contigs. CONCLUSIONS: This approach, which leverages the advantages of RNA-seq approaches, enabled rapid generation of gene-linked SNP and haplotype markers. The general approach presented in this study can be easily applied to other non-model eukaryotes irrespective of the ploidy level to identify gene-linked polymorphisms that are expected to facilitate efficient Gene Assisted Breeding (GAB), genotyping and population genetics studies. The identified SNP haplotypes reveal some of the allelic differences in the two sweet cherry cultivars analyzed. The identification of these SNP and haplotype markers is expected to significantly improve the genomic resources for sweet cherry and facilitate efficient GAB in this non-model crop.