QTL analysis on a lemon population provides novel insights on the genetic regulation of the tolerance to the two-spotted spider mite attack.

BACKGROUND: Among the Citrus species, lemon (Citrus limon Burm f.) is one of the most affected by the two-spotted spider mite (Tetranychus urticae Koch). Moreover, chemical control is hampered by the mite's ability to develop genetic resistance against acaricides. In this context, the identification of the genetic basis of the host resistance could represent a sustainable strategy for spider mite control. In the present study, a marker-trait association analysis was performed on a lemon population employing an association mapping approach. An inter-specific full-sib population composed of 109 accessions was phenotyped through a detached-leaf assays performed in modified Huffaker cells. Those individuals, complemented with two inter-specific segregating populations, were genotyped using a target-sequencing approach called SPET (Single Primer Enrichment Technology), the resulting SNPs were employed for the generation of an integrated genetic map. RESULTS: The percentage of damaged area in the full-sib population showed a quantitative distribution with values ranging from 0.36 to 9.67%. A total of 47,298 SNPs were selected for an association mapping study and a significant marker linked with resistance to spider mite was detected on linkage group 5. In silico gene annotation of the QTL interval enabled the detection of 13 genes involved in immune response to biotic and abiotic stress. Gene expression analysis showed an over expression of the gene encoding for the ethylene-responsive transcription factor ERF098-like, already characterized in Arabidopsis and in rice for its involvement in defense response. CONCLUSION: The identification of a molecular marker linked to the resistance to spider mite attack can pave the way for the development of marker-assisted breeding plan for the development of novel selection coupling favorable agronomical traits (e.g. fruit quality, yield) with a higher resistance toward the mite.

The origin and the genetic regulation of the self-compatibility mechanism in clementine (Citrus clementina Hort. ex Tan.).

Self-incompatibility (SI) is a genetic mechanism common in flowering plants to prevent self-fertilization. Among citrus species, several pummelo, mandarin, and mandarin-like accessions show SI behavior. In these species, SI is coupled with a variable degree of parthenocarpy ensuring the production of seedless fruits, a trait that is highly appreciated by consumers. In Citrus, recent evidences have shown the presence of a gametophytic SI system based on S-ribonucleases (S-RNases) ability to impair self-pollen tube growth in the upper/middle part of the style. In the present study, we combined PCR analysis and next-generation sequencing technologies, to define the presence of S7- and S11-Rnases in the S-genotype of the Citrus clementina (Hort. ex Tan.), the self-incompatible 'Comune' clementine and its self-compatible natural mutant 'Monreal'. The reference genome of 'Monreal' clementine is presented for the first time, providing more robust results on the genetic sequence of the newly discovered S7-RNase. SNP discovery analysis coupled with the annotation of the variants detected enabled the identification of 7,781 SNPs effecting 5,661 genes in 'Monreal' compared to the reference genome of C. clementina. Transcriptome analysis of unpollinated pistils at the mature stage from both clementine genotypes revealed the lack of expression of S7-RNase in 'Monreal' suggesting its involvement in the loss of the SI response. RNA-seq analysis followed by gene ontology studies enabled the identification of 2,680 differentially expressed genes (DEGs), a significant number of those is involved in oxidoreductase and transmembrane transport activity. Merging of DNA sequencing and RNA data led to the identification of 164 DEGs characterized by the presence of at least one SNP predicted to induce mutations with a high effect on their amino acid sequence. Among them, four candidate genes referring to two Agamous-like MADS-box proteins, to MYB111 and to MLO-like protein 12 were validated. Moreover, the transcription factor MYB111 appeared to contain a binding site for the 2.0-kb upstream sequences of the S7- and S11-RNase genes. These results provide useful information about the genetic bases of SI indicating that SNPs present in their sequence could be responsible for the differential expression and the regulation of S7-RNase and consequently of the SI mechanism.

Genome-wide footprints in the carob tree (Ceratonia siliqua) unveil a new domestication pattern of a fruit tree in the Mediterranean.

Intense research efforts over the last two decades have renewed our understanding of plant phylogeography and domestication in the Mediterranean basin. Here we aim to investigate the evolutionary history and the origin of domestication of the carob tree (Ceratonia siliqua), which has been cultivated for millennia for food and fodder. We used >1000 microsatellite genotypes to delimit seven carob evolutionary units (CEUs). We investigated genome-wide diversity and evolutionary patterns of the CEUs with 3557 single nucleotide polymorphisms generated by restriction-site associated DNA sequencing (RADseq). To address the complex wild vs. cultivated status of sampled trees, we classified 56 sampled populations across the Mediterranean basin as wild, seminatural or cultivated. Nuclear and cytoplasmic loci were identified from RADseq data and separated for analyses. Phylogenetic analyses of these genomic-wide data allowed us to resolve west-to-east expansions from a single long-term refugium probably located in the foothills of the High Atlas Mountains near the Atlantic coast. Our findings support multiple origins of domestication with a low impact on the genetic diversity at range-wide level. The carob was mostly domesticated from locally selected wild genotypes and scattered long-distance westward dispersals of domesticated varieties by humans, concomitant with major historical migrations by Romans, Greeks and Arabs. Ex situ efforts to preserve carob genetic resources should prioritize accessions from both western and eastern populations, with emphasis on the most differentiated CEUs situated in southwest Morocco, south Spain and eastern Mediterranean. Our study highlights the relevance of wild and seminatural habitats in the conservation of genetic resources for cultivated trees.

Biotechnological Approaches for Genetic Improvement of Lemon (Citrus limon (L.) Burm. f.) against Mal Secco Disease.

Among Citrus species, lemon is one of the most susceptible to mal secco disease, a tracheomycosis caused by the mitosporic fungus Plenodomus tracheiphilus, which induces chlorosis followed by leaf drop and progressive desiccation of twigs and branches. Severe infection can cause the death of the plant. Since no effective control strategies are available to efficiently control the pathogen spread, host tolerance is the most desirable goal in the struggle against mal secco disease. To date, both traditional breeding programs and biotechnological techniques were not efficient in developing novel varieties coupling tolerance to mal secco with optimal fruit quality. Furthermore, the genetic basis of host resistance has not been fully deciphered yet, hampering the set-up of marker-assisted selection (MAS) schemes. This paper provides an overview of the biotechnological approaches adopted so far for the selection of mal secco tolerant lemon varieties and emphasizes the promising contribution of marker-trait association analysis techniques for both unraveling the genetic determinism of the resistance to mal secco and detecting molecular markers that can be readily used for MAS. Such an approach has already proved its efficiency in several crops and could represent a valuable tool to select novel lemon varieties coupling superior fruit quality traits and resistance to mal secco.

Bioprospecting of Beneficial Bacteria Traits Associated With Tomato Root in Greenhouse Environment Reveals That Sampling Sites Impact More Than the Root Compartment.

Tomato is subject to several diseases that affect both field- and greenhouse-grown crops. To select cost-effective potential biocontrol agents, we used laboratory throughput screening to identify bacterial strains with versatile characteristics suitable for multipurpose uses. The natural diversity of tomato root-associated bacterial communities was bioprospected under a real-world environment represented by an intensive tomato cultivation area characterized by extraseasonal productions in the greenhouse. Approximately 400 tomato root-associated bacterial isolates, in majority Gram-negative bacteria, were isolated from three compartments: the soil close to the root surface (rhizosphere, R), the root surface (rhizoplane, RP), and the root interior (endorhizosphere, E). A total of 33% of the isolates produced siderophores and were able to solubilize phosphates and grow on NA with 8% NaCl. A total of 30% of the root-associated bacteria showed antagonistic activity against all the tomato pathogens tested, i.e., Clavibacter michiganesis pv. michiganensis, Pseudomonas syringae pv. tomato, Pseudomonas corrugata and Xanthomonas euvesicatoria pv. perforans, and Fusarium oxysporum f. sp. lycopersici. We found that the sampling site rather than the root compartment of isolation influenced bacterial composition in terms of analyzed phenotype. This was demonstrated through a diversity analysis including general characteristics and PGPR traits, as well as biocontrol activity in vitro. Analysis of 16S rRNA gene (rDNA) sequencing of 77 culturable endophytic bacteria that shared multiple beneficial activity revealed a predominance of bacteria in Bacillales, Enterobacteriales, and Pseudomonadales. Their in vitro antagonistic activity showed that Bacillus species were significantly more active than the isolates in the other taxonomic group. In planta activity against phytopathogenic bacteria of a subset of Bacillus and Pseudomonas isolates was also assessed.

Genomic prediction of fruit texture and training population optimization towards the application of genomic selection in apple.

Texture is a complex trait and a major component of fruit quality in apple. While the major effect of MdPG1, a gene controlling firmness, has already been exploited in elite cultivars, the genetic basis of crispness remains poorly understood. To further improve fruit texture, harnessing loci with minor effects via genomic selection is therefore necessary. In this study, we measured acoustic and mechanical features in 537 genotypes to dissect the firmness and crispness components of fruit texture. Predictions of across-year phenotypic values for these components were calculated using a model calibrated with 8,294 SNP markers. The best prediction accuracies following cross-validations within the training set of 259 genotypes were obtained for the acoustic linear distance (0.64). Predictions for biparental families using the entire training set varied from low to high accuracy, depending on the family considered. While adding siblings or half-siblings into the training set did not clearly improve predictions, we performed an optimization of the training set size and composition for each validation set. This allowed us to increase prediction accuracies by 0.17 on average, with a maximal accuracy of 0.81 when predicting firmness in the 'Gala' × 'Pink Lady' family. Our results therefore identified key genetic parameters to consider when deploying genomic selection for texture in apple. In particular, we advise to rely on a large training population, with high phenotypic variability from which a 'tailored training population' can be extracted using a priori information on genetic relatedness, in order to predict a specific target population.

Platelet­to­lymphocyte ratio, neutrophil­to­lymphocyte ratio and monocyte­to­HDL cholesterol ratio as markers of peripheral artery disease in elderly patients.

Solid evidence underlines the pivotal role played by inflammation regarding atherosclerosis. Peripheral artery disease (PAD) is one of atherosclerotic cardiovascular diseases (CVDs), it is highly frequently diagnosed in older individuals. In the present study we carried out an investigation on the association between platelet­to­lymphocytes ratio (PLR), neutrophil­to­lymphocyte ratio (NLR), monocyte­to­HDL cholesterol ratio (MHR) with PAD as favourable markers. We identified 300 subjects aged over 70 years, without any concomitant CVDs. The PLR, NLR and MHR were assessed from peripheral venous blood routinely drawn in the ward during hospitalization. Patients were divided in groups according to ankle brachial index (ABI) value (>0.9; 0.9­0.99; 1­1.4; >1.4). Higher PLR (P=0.007), NLR (P=0.0001) and MHR (P=0.0001) were associated with <0.9 ABI. Patients with a >1.4 ABI showed NLR values higher compared to >0.9l ABI (P<0.01). Univariate linear regression analysis demonstrated the direct correlation between increase in PLR (P=0.0023)and MHR (P<0.0001) with the decrease in ABI value. In multivariate linear regression analysis including main cardiovascular risk factors we found that PLR, NLR and MHR were independently associated with lower ABI (P=0.0011). Results show and suggest that the elevated PLR, NLR and MHR are related to PAD evaluated with ABI measurement. PLR and MHR seem to be more reliable markers than NLR in PAD. NLR seems to be more related to incompressibility of arterial wall. It is hypothesized that these three indexes may play a role as simple and repetitive markers of PAD.

Substantial Equivalence of a Transgenic Lemon Fruit Showing Postharvest Fungal Pathogens Resistance.

The development of genetically modified (GM) crops speeds up the obtainment of novel varieties with improved agronomic characteristics. However, the risk evaluation of the use of GMs is mandatory before their release in the market. In this paper, an untargeted and comprehensive nuclear magnetic resonance-based metabolomic study was carried out on the peel and flesh of a transgenic lemon clone (E23) expressing the chit42 gene and exhibiting an increased tolerance to some pathogenic fungi and on its wild type. Results highlighted a substantial equivalence of the metabolomics profile of the transgenic clone compared to the wild type. In addition, an enhanced response of the E23 clone toward fungal pathogens affecting the postharvest management in lemon was evidenced. These results confirm the potential of genetic engineering for the punctual modification of specific agronomic traits without altering the whole pattern of metabolites and open new perspectives for a more sustainable and effective management of specific postharvest diseases in citrus.

Elucidating the contribution of wild related species on autochthonous pear germplasm: A case study from Mount Etna.

The pear (genus Pyrus) is one of the most ancient and widely cultivated tree fruit crops in temperate climates. The Mount Etna area claims a large number of pear varieties differentiated due to a long history of cultivation and environmental variability, making this area particularly suitable for genetic studies. Ninety-five pear individuals were genotyped using the simple sequence repeat (SSR) methodology interrogating both the nuclear (nDNA) and chloroplast DNA (cpDNA) to combine an investigation of maternal inheritance of chloroplast SSRs (cpSSRs) with the high informativity of nuclear SSRs (nSSRs). The germplasm was selected ad hoc to include wild genotypes, local varieties, and national and international cultivated varieties. The objectives of this study were as follows: (i) estimate the level of differentiation within local varieties; (ii) elucidate the phylogenetic relationships between the cultivated genotypes and wild accessions; and (iii) estimate the potential genetic flow and the relationship among the germplasms in our analysis. Eight nSSRs detected a total of 136 alleles with an average minor allelic frequency and observed heterozygosity of 0.29 and 0.65, respectively, whereas cpSSRs allowed identification of eight haplotypes (S4 Table). These results shed light on the genetic relatedness between Italian varieties and wild genotypes. Among the wild species, compared with P. amygdaliformis, few P. pyraster genotypes exhibited higher genetic similarity to local pear varieties. Our analysis revealed the presence of genetic stratification with a 'wild' subpopulation characterizing the genetic makeup of wild species and the international cultivated varieties exhibiting the predominance of the 'cultivated' subpopulation.

Genome-wide association study unravels the genetic control of the apple volatilome and its interplay with fruit texture.

Fruit quality represents a fundamental factor guiding consumers' preferences. Among apple quality traits, volatile organic compounds and texture features play a major role. Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-ToF-MS), coupled with an artificial chewing device, was used to profile the entire apple volatilome of 162 apple accessions, while the fruit texture was dissected with a TAXT-AED texture analyzer. The array of volatile compounds was classed into seven major groups and used in a genome-wide association analysis carried out with 9142 single nucleotide polymorphisms (SNPs). Marker-trait associations were identified on seven chromosomes co-locating with important candidate genes for aroma, such as MdAAT1 and MdIGS. The integration of volatilome and fruit texture data conducted with a multiple factor analysis unraveled contrasting behavior, underlying opposite regulation of the two fruit quality aspects. The association analysis using the first two principal components identified two QTLs located on chromosomes 10 and 2, respectively. The distinction of the apple accessions on the basis of the allelic configuration of two functional markers, MdPG1 and MdACO1, shed light on the type of interplay existing between fruit texture and the production of volatile organic compounds.

Deciphering the genetic control of fruit texture in apple by multiple family-based analysis and genome-wide association.

Fruit texture is a complex feature composed of mechanical and acoustic properties relying on the modifications occurring in the cell wall throughout fruit development and ripening. Apple is characterized by a large variation in fruit texture behavior that directly impacts both the consumer's appreciation and post-harvest performance. To decipher the genetic control of fruit texture comprehensively, two complementing quantitative trait locus (QTL) mapping approaches were employed. The first was represented by a pedigree-based analysis (PBA) carried out on six full-sib pedigreed families, while the second was a genome-wide association study (GWAS) performed on a collection of 233 apple accessions. Both plant materials were genotyped with a 20K single nucleotide polymorphism (SNP) array and phenotyped with a sophisticated high-resolution texture analyzer. The overall QTL results indicated the fundamental role of chromosome 10 in controlling the mechanical properties, while chromosomes 2 and 14 were more associated with the acoustic response. The latter QTL, moreover, showed a consistent relationship between the QTL-estimated genotypes and the acoustic performance assessed among seedlings. The in silico annotation of these intervals revealed interesting candidate genes potentially involved in fruit texture regulation, as suggested by the gene expression profile. The joint integration of these approaches sheds light on the specific control of fruit texture, enabling important genetic information to assist in the selection of valuable fruit quality apple varieties.

A high-density, multi-parental SNP genetic map on apple validates a new mapping approach for outcrossing species.

Quantitative trait loci (QTL) mapping approaches rely on the correct ordering of molecular markers along the chromosomes, which can be obtained from genetic linkage maps or a reference genome sequence. For apple (Malus domestica Borkh), the genome sequence v1 and v2 could not meet this need; therefore, a novel approach was devised to develop a dense genetic linkage map, providing the most reliable marker-loci order for the highest possible number of markers. The approach was based on four strategies: (i) the use of multiple full-sib families, (ii) the reduction of missing information through the use of HaploBlocks and alternative calling procedures for single-nucleotide polymorphism (SNP) markers, (iii) the construction of a single backcross-type data set including all families, and (iv) a two-step map generation procedure based on the sequential inclusion of markers. The map comprises 15 417 SNP markers, clustered in 3 K HaploBlock markers spanning 1 267 cM, with an average distance between adjacent markers of 0.37 cM and a maximum distance of 3.29 cM. Moreover, chromosome 5 was oriented according to its homoeologous chromosome 10. This map was useful to improve the apple genome sequence, design the Axiom Apple 480 K SNP array and perform multifamily-based QTL studies. Its collinearity with the genome sequences v1 and v3 are reported. To our knowledge, this is the shortest published SNP map in apple, while including the largest number of markers, families and individuals. This result validates our methodology, proving its value for the construction of integrated linkage maps for any outbreeding species.

Detecting QTLs and putative candidate genes involved in budbreak and flowering time in an apple multiparental population.

In temperate trees, growth resumption in spring time results from chilling and heat requirements, and is an adaptive trait under global warming. Here, the genetic determinism of budbreak and flowering time was deciphered using five related full-sib apple families. Both traits were observed over 3 years and two sites and expressed in calendar and degree-days. Best linear unbiased predictors of genotypic effect or interaction with climatic year were extracted from mixed linear models and used for quantitative trait locus (QTL) mapping, performed with an integrated genetic map containing 6849 single nucleotide polymorphisms (SNPs), grouped into haplotypes, and with a Bayesian pedigree-based analysis. Four major regions, on linkage group (LG) 7, LG10, LG12, and LG9, the latter being the most stable across families, sites, and years, explained 5.6-21.3% of trait variance. Co-localizations for traits in calendar days or growing degree hours (GDH) suggested common genetic determinism for chilling and heating requirements. Homologs of two major flowering genes, AGL24 and FT, were predicted close to LG9 and LG12 QTLs, respectively, whereas Dormancy Associated MADs-box (DAM) genes were near additional QTLs on LG8 and LG15. This suggests that chilling perception mechanisms could be common among perennial and annual plants. Progenitors with favorable alleles depending on trait and LG were identified and could benefit new breeding strategies for apple adaptation to temperature increase.

ASSIsT: an automatic SNP scoring tool for in- and outbreeding species.

UNLABELLED: ASSIsT (Automatic SNP ScorIng Tool) is a user-friendly customized pipeline for efficient calling and filtering of SNPs from Illumina Infinium arrays, specifically devised for custom genotyping arrays. Illumina has developed an integrated software for SNP data visualization and inspection called GenomeStudio (GS). ASSIsT builds on GS-derived data and identifies those markers that follow a bi-allelic genetic model and show reliable genotype calls. Moreover, ASSIsT re-edits SNP calls with null alleles or additional SNPs in the probe annealing site. ASSIsT can be employed in the analysis of different population types such as full-sib families and mating schemes used in the plant kingdom (backcross, F1, F2), and unrelated individuals. The final result can be directly exported in the format required by the most common software for genetic mapping and marker-trait association analysis. ASSIsT is developed in Python and runs in Windows and Linux. AVAILABILITY AND IMPLEMENTATION: The software, example data sets and tutorials are freely available at http://compbiotoolbox.fmach.it/assist/. CONTACT: eric.vandeweg@wur.nl.

Accuracy and responses of genomic selection on key traits in apple breeding.

The application of genomic selection in fruit tree crops is expected to enhance breeding efficiency by increasing prediction accuracy, increasing selection intensity and decreasing generation interval. The objectives of this study were to assess the accuracy of prediction and selection response in commercial apple breeding programmes for key traits. The training population comprised 977 individuals derived from 20 pedigreed full-sib families. Historic phenotypic data were available on 10 traits related to productivity and fruit external appearance and genotypic data for 7829 SNPs obtained with an Illumina 20K SNP array. From these data, a genome-wide prediction model was built and subsequently used to calculate genomic breeding values of five application full-sib families. The application families had genotypes at 364 SNPs from a dedicated 512 SNP array, and these genotypic data were extended to the high-density level by imputation. These five families were phenotyped for 1 year and their phenotypes were compared to the predicted breeding values. Accuracy of genomic prediction across the 10 traits reached a maximum value of 0.5 and had a median value of 0.19. The accuracies were strongly affected by the phenotypic distribution and heritability of traits. In the largest family, significant selection response was observed for traits with high heritability and symmetric phenotypic distribution. Traits that showed non-significant response often had reduced and skewed phenotypic variation or low heritability. Among the five application families the accuracies were uncorrelated to the degree of relatedness to the training population. The results underline the potential of genomic prediction to accelerate breeding progress in outbred fruit tree crops that still need to overcome long generation intervals and extensive phenotyping costs.

Development and validation of a 20K single nucleotide polymorphism (SNP) whole genome genotyping array for apple (Malus × domestica Borkh).

High-density SNP arrays for genome-wide assessment of allelic variation have made high resolution genetic characterization of crop germplasm feasible. A medium density array for apple, the IRSC 8K SNP array, has been successfully developed and used for screens of bi-parental populations. However, the number of robust and well-distributed markers contained on this array was not sufficient to perform genome-wide association analyses in wider germplasm sets, or Pedigree-Based Analysis at high precision, because of rapid decay of linkage disequilibrium. We describe the development of an Illumina Infinium array targeting 20K SNPs. The SNPs were predicted from re-sequencing data derived from the genomes of 13 Malus × domestica apple cultivars and one accession belonging to a crab apple species (M. micromalus). A pipeline for SNP selection was devised that avoided the pitfalls associated with the inclusion of paralogous sequence variants, supported the construction of robust multi-allelic SNP haploblocks and selected up to 11 entries within narrow genomic regions of ±5 kb, termed focal points (FPs). Broad genome coverage was attained by placing FPs at 1 cM intervals on a consensus genetic map, complementing them with FPs to enrich the ends of each of the chromosomes, and by bridging physical intervals greater than 400 Kbps. The selection also included ∼3.7K validated SNPs from the IRSC 8K array. The array has already been used in other studies where ∼15.8K SNP markers were mapped with an average of ∼6.8K SNPs per full-sib family. The newly developed array with its high density of polymorphic validated SNPs is expected to be of great utility for Pedigree-Based Analysis and Genomic Selection. It will also be a valuable tool to help dissect the genetic mechanisms controlling important fruit quality traits, and to aid the identification of marker-trait associations suitable for the application of Marker Assisted Selection in apple breeding programs.

A multidisciplinary approach providing new insight into fruit flesh browning physiology in apple (Malus x domestica Borkh.).

In terms of the quality of minimally processed fruit, flesh browning is fundamentally important in the development of an aesthetically unpleasant appearance, with consequent off-flavours. The development of browning depends on the enzymatic action of the polyphenol oxidase (PPO). In the 'Golden Delicious' apple genome ten PPO genes were initially identified and located on three main chromosomes (2, 5 and 10). Of these genes, one element in particular, here called Md-PPO, located on chromosome 10, was further investigated and genetically mapped in two apple progenies ('Fuji x Pink Lady' and 'Golden Delicious x Braeburn'). Both linkage maps, made up of 481 and 608 markers respectively, were then employed to find QTL regions associated with fruit flesh browning, allowing the detection of 25 QTLs related to several browning parameters. These were distributed over six linkage groups with LOD values spanning from 3.08 to 4.99 and showed a rate of phenotypic variance from 26.1 to 38.6%. Anchoring of these intervals to the apple genome led to the identification of several genes involved in polyphenol synthesis and cell wall metabolism. Finally, the expression profile of two specific candidate genes, up and downstream of the polyphenolic pathway, namely phenylalanine ammonia lyase (PAL) and polyphenol oxidase (PPO), provided insight into flesh browning physiology. Md-PPO was further analyzed and two haplotypes were characterised and associated with fruit flesh browning in apple.