Genotyping by Sequencing in Almond: SNP Discovery, Linkage Mapping, and Marker Design

Advancing tree genomics to future proof next generation orchard production

The challenges facing tree orchard production in the coming years will be largely driven by changes in the climate affecting the sustainability of farming practices in specific geographical regions. Identifying key traits that enable tree crops to modify their growth to varying environmental conditions and taking advantage of new crop improvement opportunities and technologies will ensure the tree crop industry remains viable and profitable into the future. In this review article we 1) outline climate and sustainability challenges relevant to horticultural tree crop industries, 2) describe key tree crop traits targeted for improvement in agroecosystem productivity and resilience to environmental change, and 3) discuss existing and emerging genomic technologies that provide opportunities for industries to future proof the next generation of orchards.

Almond population genomics and non-additive GWAS reveal new insights into almond dissemination history and candidate genes for nut traits and blooming time

Domestication drastically changed crop genomes, fixing alleles of interest and creating different genetic populations. Genome-wide association studies (GWASs) are a powerful tool to detect these alleles of interest (and so QTLs). In this study, we explored the genetic structure as well as additive and non-additive genotype-phenotype associations in a collection of 243 almond accessions. Our genetic structure analysis strongly supported the subdivision of the accessions into five ancestral groups, all formed by accessions with a common origin. One of these groups was formed exclusively by Spanish accessions, while the rest were mainly formed by accessions from China, Italy, France, and the USA. These results agree with archaeological and historical evidence that separate modern almond dissemination into four phases: Asiatic, Mediterranean, Californian, and southern hemisphere. In total, we found 13 independent QTLs for nut weight, crack-out percentage, double kernels percentage, and blooming time. Of the 13 QTLs found, only one had an additive effect. Through candidate gene analysis, we proposed Prudul26A013473 as a candidate gene responsible for the main QTL found in crack-out percentage, Prudul26A012082 and Prudul26A017782 as candidate genes for the QTLs found in double kernels percentage, and Prudul26A000954 as a candidate gene for the QTL found in blooming time. Our study enhances our knowledge of almond dissemination history and will have a great impact on almond breeding.

Transcriptome analysis reveals the mechanism of different fruit appearance between apricot (Armeniaca vulgaris Lam.) and its seedling

BACKGROUND

Apricot fruit has great economic value. In the process of apricot breeding using traditional breeding methods, we obtained a larger seedling (named Us) from the original variety (named U). And Us fruit is larger than U, taste better. Therefore, revealing its mechanism is very important for Apricot breeding.

METHODS

In this study, de novo assembly and transcriptome sequencing (RNA-Seq) was used to screen the differently expressed genes (DEGs) between U and Us at three development stages, including young fruits stage, mid-ripening stage and mature fruit stage.

RESULTS

The results showed that there were 6,753 DEGs at different sampling time. "Cellulose synthase (UDP-forming) activity" and "cellulose synthase activity" were the key GO terms enriched in GO, of which CESA and CSL family played a key role. "Photosynthesis-antenna proteins" and "Plant hormone signal transduction" were the candidate pathways and lhca, lhcb, Aux/IAA and SAUR were the main regulators.

CONCLUSION

The auxin signaling pathway was active in Us, of which Aux/IAAs and SAUR were the key fruit size regulators. The low level of lhca and lhcb in Us could reveal the low demand for exogenous carbon, but they increased at mature stage, which might be due to the role of aux, who was keeping the fruit growing. Aux and photosynthesis maight be the main causes of appearance formation of Us fruits. Interestingly, the higher expression of CESA and CSL proved that Us entered the hardening process earlier than U. The advanced developmental progress might also be due to the role of Aux.

Exploring genetic diversity and population structure of a large grapevine (Vitis vinifera L.) germplasm collection in Türkiye

Grapevine (Vitis Vinifera L.) has been one of the significant perennial crops in widespread temperate climate regions since its domestication around 6000 years ago. Grapevine and its products, particularly wine, table grapes, and raisins, have significant economic importance not only in grapevine-growing countries but also worldwide. Grapevine cultivation in Türkiye dates back to ancient times, and Anatolia is considered one of the main grapevine migration routes around the Mediterranean basin. Turkish germplasm collection, conserved at the Turkish Viticulture Research Institutes, includes cultivars and wild relatives mainly collected in Türkiye, breeding lines, rootstock varieties, and mutants, but also cultivars of international origin. Genotyping with high-throughput markers enables the investigation of genetic diversity, population structure, and linkage disequilibrium, which are crucial for applying genomic-assisted breeding. Here, we present the results of a high-throughput genotyping-by-sequencing (GBS) study of 341 genotypes from grapevine germplasm collection at Manisa Viticulture Research Institute. A total of 272,962 high-quality single nucleotide polymorphisms (SNP) markers on the nineteen chromosomes were identified using genotyping-by-sequencing (GBS) technology. The high-density coverage of SNPs resulted in an average of 14,366 markers per chromosome, an average polymorphism information content (PIC) value of 0.23 and an expected heterozygosity (He) value of 0.28 indicating the genetic diversity within 341 genotypes. LD decayed very fast when r² was between 0.45 and 0.2 and became flat when r² was 0.05. The average LD decay for the entire genome was 30 kb when r^2 = 0.2. The PCA and structure analysis did not distinguish the grapevine genotypes based on different origins, highlighting the occurrence of gene flow and a high amount of admixture. Analysis of molecular variance (AMOVA) results indicated a high level of genetic differentiation within populations, while variation among populations was extremely low. This study provides comprehensive information on the genetic diversity and population structure of Turkish grapevine genotypes.

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.

Analysis of Genetic Diversity and Resistance to Foliar Pathogens Based on Genotyping-by-Sequencing of a Para Rubber Diversity Panel and Progeny of an Interspecific Cross

Para rubber trees (Hevea brasiliensis) are the largest major source of natural rubber in the world. Its major pathogens are Phytophthora spp., Corynespora cassiicola, and Colletotrichum spp. A rubber diversity panel of 116 clones using over 12,000 single nucleotide polymorphisms (SNPs) from DArTSeq genotyping revealed clear phylogenetic differences in clones that originated from different geographical regions of the world. An integrated linkage map constructed with an F₁ progeny of 86 from an interspecific cross between H. brasiliensis and H. benthamiana using 23,978 markers [10,323 SNPs and 13,655 SilicoDArTs] spanned 3947.83 cM with 0.83 cM average marker-interval. The genome scaffolds that were anchored to the linkage map, covering 1.44 Gb of H. brasiliensis reference genome, revealed a high level of collinearity between the genetic map and reference genome. Association analysis identified 12 SNPs significantly associated with the resistance against Phytophthora, Corynespora, and Colletotrichum in six linkage groups: 2, 6, 12, 14, 17, and 18. Kompetitive Allele-Specific PCR marker assays were developed for those 12 SNPs, screened with 178 individuals, and detected clear separation between two genotypes. Within the proximity to those SNPs, 41 potentially key genes that have previously been reported to associate with plant disease resistance were predicted with high confidence.

QTLs Identification for Iron Chlorosis in a Segregating Peach-Almond Progeny Through Double-Digest Sequence-Based Genotyping (SBG)

Linkage maps are highly appreciated tools for cultivar and rootstock breeding programs because they are suitable for genetic and genomic studies. In this study, we report on using sequence-based genotyping (SBG) approach to simultaneously discover and genotype SNPs from two peach-based rootstocks ("Adafuel" and "Flordaguard") and their progeny (n = 118): from a initial mapping population composed of 131 seedlings. The plant material was developed at the EEAD-CSIC Prunus rootstocks breeding program, aiming to obtain a segregating progeny for a range of characters of agronomical interest to rootstock breeding (iron-chlorosis and root-asphyxia tolerance, nematode resistance, vigor traits, and other effects on scion cultivars). Sequence reads obtained from double-digest SBG were aligned to the P. persica reference genome (Peach v2.0). While eight linkage groups were constructed for "Adafuel," only four linkage groups were constructed for "Flordaguard," given the low heterozygosity of this last genotype. High synteny and co-linearity were observed between obtained maps and Peach v2.0. On the other hand, this work aimed to elucidate the genetic basis of leaf chlorosis tolerance using the phenotypic segregation of the progeny to iron-chlorosis tolerance, along with the QTLs responsible for leaf chlorosis. The F1 mapping population, composed initially of 131 seedlings, was growing in four field trials established on calcareous soils at the experimental field of the EEAD-CSIC in Zaragoza, Spain. From the initial mapping population, 131 individuals were selected for their phenotypical characterization with SPAD measurements of plants grown in the field, exhibiting a great variability. Significant QTLs associated with tolerance to iron chlorosis were found in LG1, LG5, LG7, and LG8. The significant QTLs detected in LG5 and LG7 have not been associated with this abiotic stress before in Prunus. Several candidate genes such as Prupe.1G541100, predicted as glutamyl-tRNA reductase 1, Prupe.1G468200, encoding a 2-oxoglutarate (2OG), and Fe(II)-dependent oxygenase superfamily protein or Prupe.1G577000 (ppa011050.m), a NIFU-like protein 2 (NIFU2) were detected. The exact biological function of some of these genes should be verified for the future development of marker-assisted selection for peach iron chlorosis tolerance.

Whole-genome sequence and methylome profiling of the almond [Prunus dulcis (Mill.) D.A. Webb] cultivar 'Nonpareil'

Almond [Prunus dulcis (Mill.) D.A. Webb] is an economically important, specialty nut crop grown almost exclusively in the United States. Breeding and improvement efforts worldwide have led to the development of key, productive cultivars, including 'Nonpareil,' which is the most widely grown almond cultivar. Thus far, genomic resources for this species have been limited, and a whole-genome assembly for 'Nonpareil' is not currently available despite its economic importance and use in almond breeding worldwide. We generated a 571X coverage genome sequence using Illumina, PacBio, and optical mapping technologies. Gene prediction revealed 49,321 putative genes using MinION Oxford nanopore and Illumina RNA sequencing, and genome annotation found that 68% of predicted models are associated with at least one biological function. Furthermore, epigenetic signatures of almond, namely DNA cytosine methylation, have been implicated in a variety of phenotypes including self-compatibility, bud dormancy, and development of noninfectious bud failure. In addition to the genome sequence and annotation, this report also provides the complete methylome of several almond tissues, including leaf, flower, endocarp, mesocarp, exocarp, and seed coat. Comparisons between methylation profiles in these tissues revealed differences in genome-wide weighted % methylation and chromosome-level methylation enrichment.

Genotyping by sequencing-based linkage map construction and identification of quantitative trait loci for yield-related traits and oil content in Jatropha (Jatropha curcas L.)

BACKGROUND

Jatropha (Jatropha curcas L.) has been considered as a potential bioenergy crop and its genetic improvement is essential for higher seed yield and oil content which has been hampered due to lack of desirable molecular markers.

METHODS AND RESULTS

An F₂ population was created using an intraspecific cross involving a Central American line RJCA9 and an Asiatic species RJCS-9 to develop a dense genetic map and for Quantitative trait loci (QTL) identification. The genotyping-by-sequencing (GBS) approach was used to genotype the mapping population of 136 F₂ individuals along with the two parental lines for classification of the genotypes based on single nucleotide polymorphism (SNPs). NextSeq 2500 sequencing technology provided a total of 517.23 million clean reads, with an average of ~ 3.8 million reads per sample. We analysed 411 SNP markers and developed 11 linkage groups. The total length of the genetic map was 4092.3 cM with an average marker interval of 10.04 cM. We have identified a total of 83 QTLs for various yield and oil content governing traits. The percentage of phenotypic variation (PV) was found to be in the range of 8.81 to 65.31%, and a QTL showed the maximum PV of 65.3% for a total seed number on the 6th linkage group (LG).

CONCLUSIONS

The QTLs detected in this study for various phenotypic traits will lay down the path for marker-assisted breeding in the future and cloning of genes that are responsible for phenotypic variation.

Polymorphisms and gene expression in the almond IGT family are not correlated to variability in growth habit in major commercial almond cultivars

Almond breeding programs aimed at selecting cultivars adapted to intensive orchards have recently focused on the optimization of tree architecture. This multifactorial trait is defined by numerous components controlled by processes such as hormonal responses, gravitropism and light perception. Gravitropism sensing is crucial to control the branch angle and therefore, the tree habit. A gene family, denominated IGT family after a shared conserved domain, has been described as involved in the regulation of branch angle in several species, including rice and Arabidopsis, and even in fruit trees like peach. Here we identified six members of this family in almond: LAZY1, LAZY2, TAC1, DRO1, DRO2, IGT-like. After analyzing their protein sequences in forty-one almond cultivars and wild species, little variability was found, pointing a high degree of conservation in this family. To our knowledge, this is the first effort to analyze the diversity of IGT family proteins in members of the same tree species. Gene expression was analyzed in fourteen cultivars of agronomical interest comprising diverse tree habit phenotypes. Only LAZY1, LAZY2 and TAC1 were expressed in almond shoot tips during the growing season. No relation could be established between the expression profile of these genes and the variability observed in the tree habit. However, some insight has been gained in how LAZY1 and LAZY2 are regulated, identifying the IPA1 almond homologues and other transcription factors involved in hormonal responses as regulators of their expression. Besides, we have found various polymorphisms that could not be discarded as involved in a potential polygenic origin of regulation of architectural phenotypes. Therefore, we have established that neither the expression nor the genetic polymorphism of IGT family genes are correlated to diversity of tree habit in currently commercialized almond cultivars, with other gene families contributing to the variability of these traits.

Mutation of a bHLH transcription factor allowed almond domestication

Wild almond species accumulate the bitter and toxic cyanogenic diglucoside amygdalin. Almond domestication was enabled by the selection of genotypes harboring sweet kernels. We report the completion of the almond reference genome. Map-based cloning using an F₁ population segregating for kernel taste led to the identification of a 46-kilobase gene cluster encoding five basic helix-loop-helix transcription factors, bHLH1 to bHLH5. Functional characterization demonstrated that bHLH2 controls transcription of the P450 monooxygenase-encoding genes PdCYP79D16 and PdCYP71AN24, which are involved in the amygdalin biosynthetic pathway. A nonsynonymous point mutation (Leu to Phe) in the dimerization domain of bHLH2 prevents transcription of the two cytochrome P450 genes, resulting in the sweet kernel trait.

Update of Genetic Linkage Map and QTL Analysis for Growth Traits in Eucommia ulmoides Oliver

Eucommia ulmoides (Tu-chung) is an economically and ecologically important tree species which has attracted worldwide attention due to its application in pharmacology, landscaping, wind sheltering and sand fixation. Molecular marker technologies can elucidate the genetic mechanism and substantially improve the breeding efficiency of E. ulmoides. The current research updated the original linkage map, and quantitative trait loci (QTL) analysis was performed on tree growth traits measured over 10 consecutive years in an E. ulmoides F1 population (“Xiaoye” × “Qinzhong No.1”). In total, 452 polymorphic markers were scored from 365 simple sequence repeat (SSR) primers, with an average of 1.24 polymorphic markers per primer combination. The integrated map was 1913.29 cM (centimorgan) long, covering 94.10% of the estimated genome and with an average marker density of 2.20 cM. A total of 869 markers were mapped into 19 major independent linkage groups. Growth-related traits measured over 10 consecutive years showed a significant correlation, and 89 hypothetical QTLs were forecasted and divided into 27 distinct loci. Three traits for tree height, ground diameter and crown diameter detected 25 QTLs (13 loci), 32 QTLs (17 loci) and 15 QTLs (10 loci), respectively. Based on BLASTX search results in the NCBI database, six candidate genes were obtained. It is important to explore the growth-related genetic mechanism and lay the foundation for the genetic improvement of E. ulmoides at the molecular level.

Genome-wide SNP identification in Prunus rootstocks germplasm collections using Genotyping-by-Sequencing: phylogenetic analysis, distribution of SNPs and prediction of their effect on gene function

Genotyping-by-Sequencing (GBS) was applied in a set of 53 diploid Prunus rootstocks and five scion cultivars from three subgenera (Amygdalus, Prunus and Cerasus) for genome-wide SNP identification and to assess genetic diversity of both Chilean and Spanish germplasm collections. A group of 45,382 high quality SNPs (MAF >0.05; missing data <5%) were selected for analysis of this group of 58 accessions. These SNPs were distributed in genic and intergenic regions in the eight pseudomolecules of the peach genome (Peach v2.0), with an average of 53% located in exonic regions. The genetic diversity detected among the studied accessions divided them in three groups, which are in agreement with their current taxonomic classification. SNPs were classified based on their putative effect on annotated genes and KOG analysis was carried out to provide a deeper understanding of the function of 119 genes affected by high-impact SNPs. Results demonstrate the high utility for Prunus rootstocks identification and studies of diversity in Prunus species. Also, given the high number of SNPs identified in exonic regions, this strategy represents an important tool for finding candidate genes underlying traits of interest and potential functional markers for use in marker-assisted selection.

Construction of a highly saturated linkage map in Japanese plum (Prunus salicina L.) using GBS for SNP marker calling

This study reports the construction of high density linkage maps of Japanese plum (Prunus salicina Lindl.) using single nucleotide polymorphism markers (SNPs), obtained with a GBS strategy. The mapping population (An x Au) was obtained by crossing cv. "Angeleno" (An) as maternal line and cv. "Aurora" (Au) as the pollen donor. A total of 49,826 SNPs were identified using the peach genome V2.1 as a reference. Then a stringent filtering was carried out, which revealed 1,441 high quality SNPs in 137 An x Au offspring, which were mapped in eight linkage groups. Finally, the consensus map was built using 732 SNPs which spanned 617 cM with an average of 0.96 cM between adjacent markers. The majority of the SNPs were distributed in the intragenic region in all the linkage groups. Considering all linkage groups together, 85.6% of the SNPs were located in intragenic regions and only 14.4% were located in intergenic regions. The genetic linkage analysis was able to co-localize two to three SNPs over 37 putative orthologous genes in eight linkage groups in the Japanese plum map. These results indicate a high level of synteny and collinearity between Japanese plum and peach genomes.