Self-(in)compatibility in Tunisian apple accessions [Malus domestica. Borkh]: S-genotypes identification and pollen tube growth analysis.

MAIN CONCLUSION: Self-incompatibility studies have revealed a potential use of Tunisian apple resources for crop improvement and modern breeding programs and a likely correlation between the pollen tube growth and flowering period. Apples [Malus domestica. Borkh] exhibit an S-RNase-based gametophytic self-incompatibility (GSI) system. Four primer combinations were used to S-genotype eighteen Tunisian local apple accessions and twelve introduced accessions that served as references. Within the Tunisian local accessions, S2, S3, S7, and S28 S-alleles were the most frequent and were assigned to 14 S-genotypes; among them, S7S28, S3S7, S2S5, and S2S3 were the most abundant. PCA plot showed that population structuring was affected by the S-alleles frequencies and revealed a modern origin of the Tunisian varieties rather than being ancient ones. Nonetheless, the results obtained with 17 SSR markers showed a separate grouping of local Tunisian accessions that calls into question the hypothesis discussed. Pollination experiments showed that the pollen started to germinate within 24 h of pollination but 48 h after pollination in the "El Fessi" accession. The first pollen tubes arrived in the styles within 36 h of pollination in two early flowering accessions known as "Arbi" and "Bokri", and after 72 h of pollination in late flowering "El Fessi" and 48 h after pollination in remaining accessions. The first pollen tube arrests were observed in accessions "Arbi" and "Bokri" within 84 h of pollination, within 108 h of pollination in "El Fessi" and within 108 h of pollination in remaining accessions. In the apple accession called "Boutabgaya," the pollen tubes reached the base of the style within 120 h of pollination without being aborted. Nevertheless, the self-compatible nature of "Boutabgaya" needs more studies to be confirmed. However, our results revealed the malfunction of the female component of the GSI in this accession. To conclude, this work paved the path for further studies to enhance the insight (i) into the relation between the flowering period and the pollen tube growth, (ii) self-compatible nature of "Boutabgaya", and (iii) the origin of the Tunisian apple.

Isolation and Functional Analysis of a PISTILLATA-like MADS-Box Gene from Argan Tree (Argania spinosa).

Argan trees (Argania spinosa) belong to a species native to southwestern Morocco, playing an important role in the environment and local economy. Argan oil extracted from kernels has a unique composition and properties. Argan trees were introduced in Tunisia, where hundreds of trees can be found nowadays. In this study, we examined reproductive development in Argan trees from four sites in Tunisia and carried out the functional characterization of a floral homeotic gene in this non-model species. Despite the importance of reproductive development, nothing is known about the genetic network controlling flower development in Argania spinosa. Results obtained in several plant species established that floral organ development is mostly controlled by MADS-box genes and, in particular, APETALA3 (AP3) and PISTILLATA (PI) homologs are required for proper petal and stamen identity. Here, we describe the isolation and functional characterization of a MADS-box gene from Argania spinosa. Phylogenetic analyses showed strong homology with PI-like proteins, and the expression of the gene was found to be restricted to the second and third whorls. Functional homology with Arabidopsis PI was demonstrated by the ability of AsPI to confer petal and stamen identity when overexpressed in a pi-1 mutant background. The identification and characterization of this gene support the strong conservation of PI homologs among distant angiosperm plants.

Molecular and Evolutionary Characterization of Pollen S Determinant (SFB Alleles) in Four Diploid and Hexaploid Plum Species (Prunus spp.).

In more than 60 families of angiosperms, the self- and cross-fertilization is avoided through a complex widespread genetic system called self-incompatibility (SI). One of the major puzzling issues concerning the SI is the evolution of this system in species with complex polyploid genomes. Among plums, one of the first fruits species to attract human interest, polyploid species represent enormous genetic potential, which can be exploited in breeding programs. However, molecular studies in these species are very scarce due to the complexity of their genome. In order to study the SFB gene [the male component of gametophytic self-incompatibility system (GSI)] in plum species, 36 plum accessions belonging to diploid and hexaploid species were used. A total of 19 different alleles were identified; 1 of them was revealed after analyzing sequences. Peptide sequence analysis allowed identifying the five domains features of the SFB gene. Polymorphism analysis showed a subtle difference between domesticated and open pollinated Tunisian accessions and suggested a probable influence of the ploidy level. Divergence analysis between studied sequences showed that a new specificity may appear after 5.3% of divergence at synonymous sites between pairs of sequences in Prunus insititia, 6% in Prunus cerasifera, 8% and 9% in Prunus domestica and Prunus salicina respectively. Furthermore, sites under positive selection, the ones more likely to be responsible for specificity determination, were identified. A positive and significant Pearson correlation was found between the divergence between sequences, divergence time, fixed substitutions (MK test), and PSS number. These results supported the model assuming that functionally distinct proteins have arisen not as a result of chance fixation of neutral variants, but rather as a result of positive Darwinian selection. Further, the role that plays recombination can not be ruled out, since a rate of 0.08 recombination event per polymorphic sites was identified.

Self-compatibility in peach [Prunus persica (L.) Batsch]: patterns of diversity surrounding the S-locus and analysis of SFB alleles.

Self-incompatibility (SI) to self-compatibility (SC) transition is one of the most frequent and prevalent evolutionary shifts in flowering plants. Prunus L. (Rosaceae) is a genus of over 200 species most of which exhibit a Gametophytic SI system. Peach [Prunus persica (L.) Batsch; 2n = 16] is one of the few exceptions in the genus known to be a fully self-compatible species. However, the evolutionary process of the complete and irreversible loss of SI in peach is not well understood and, in order to fill that gap, in this study 24 peach accessions were analyzed. Pollen tube growth was controlled in self-pollinated flowers to verify their self-compatible phenotypes. The linkage disequilibrium association between alleles at the S-locus and linked markers at the end of the sixth linkage group was not significant (P > 0.05), except with the closest markers suggesting the absence of a signature of negative frequency dependent selection at the S-locus. Analysis of SFB1 and SFB2 protein sequences allowed identifying the absence of some variable and hypervariable domains and the presence of additional α-helices at the C-termini. Molecular and evolutionary analysis of SFB nucleotide sequences showed a signature of purifying selection in SFB2, while the SFB1 seemed to evolve neutrally. Thus, our results show that the SFB2 allele diversified after P. persica and P. dulcis (almond) divergence, a period which is characterized by an important bottleneck, while SFB1 diversified at a transition time between the bottleneck and population expansion.

Identification of conserved genes linked to responses to abiotic stresses in leaves among different plant species.

As a consequence of global climate change, certain stress factors that have a negative impact on crop productivity such as heat, cold, drought and salinity are becoming increasingly prevalent. We conducted a meta-analysis to identify genes conserved across plant species involved in (1) general abiotic stress conditions, and (2) specific and unique abiotic stress factors (drought, salinity, extreme temperature) in leaf tissues. We collected raw data and re-analysed eight RNA-Seq studies using our previously published bioinformatic pipeline. A total of 68 samples were analysed. Gene set enrichment analysis was performed using MapMan and PageMan whereas DAVID (Database for Annotation, Visualisation and Integrated Discovery) was used for metabolic process enrichment analysis. We identified of a total of 5122 differentially expressed genes when considering all abiotic stresses (3895 were upregulated and 1227 were downregulated). Jasmonate-related genes were more commonly upregulated by drought, whereas gibberellin downregulation was a key signal for drought and heat. In contrast, cold stress clearly upregulated genes involved in ABA (abscisic acid), cytokinin and gibberellins. A gene (non-phototrophic hypocotyl) involved in IAA (indoleacetic acid) response was induced by heat. Regarding secondary metabolism, as expected, MVA pathway (mevalonate pathway), terpenoids and alkaloids were generally upregulated by all different stresses. However, flavonoids, lignin and lignans were more repressed by heat (cinnamoyl coA reductase 1 and isopentenyl pyrophosphatase). Cold stress drastically modulated genes involved in terpenoid and alkaloids. Relating to transcription factors, AP2-EREBP, MADS-box, WRKY22, MYB, homoebox genes members were significantly modulated by drought stress whereas cold stress enhanced AP2-EREBPs, bZIP members, MYB7, BELL 1 and one bHLH member. C2C2-CO-LIKE, MADS-box and a homeobox (HOMEOBOX3) were mostly repressed in response to heat. Gene set enrichment analysis showed that ubiquitin-mediated protein degradation was enhanced by heat, which unexpectedly repressed glutaredoxin genes. Cold stress mostly upregulated MAP kinases (mitogen-activated protein kinase). Findings of this work will allow the identification of new molecular markers conserved across crops linked to major genes involved in quantitative agronomic traits affected by different abiotic stress.

Conserved DNA-derived polymorphism, new markers for genetic diversity analysis of Tunisian Pistacia vera L.

Pistachio trees (Pistacia vera L.) have been cultivated in Tunisia for decades and the plantation was extended mostly in the center of the country contributing to the economic growth of marginalized areas. Herein we used conserved DNA derived polymorphism (CDDP) technique, which target specifically conserved sequences of plant functional genes, to assess the genetic diversity and construct genetic relationships among 65 Tunisian pistachio trees. A set of nine primers were used and 157 CDDP markers were revealed with an average of 17.44 showing a high degree of polymorphism (99.37%). The average of polymorphism information content of CDDP markers was of 0.86, which indicates the efficiency of CDDP primers in the estimation of genetic diversity between pistachios. UPGMA dendrogram and the principal component analysis showed four clusters of analyzed pistachios trees. Our results showed that the genetic structure depends on: (1) the gene exchanges between groups, (2) the geographical origin and (3) the sex of the tree. The same result was revealed by the Bayesian analysis implemented in STRUCTURE at K = 4, in which the pistachio genotypes of El Guettar, Kasserine and Sfax were assigned with more than 80% of probability. Our results prove polymorphism and efficiency of CDDP markers in the characterization and genetic diversity analysis of P. vera L. genotypes to define conservation strategy.

Analysis of Self-Incompatibility and Genetic Diversity in Diploid and Hexaploid Plum Genotypes.

During the last decade, S-genotyping has been extensively investigated in fruit tree crops such as those belonging to the Prunus genus, including plums. In plums, S-allele typing has been largely studied in diploid species but works are scarcer in polyploid species due to the complexity of the polyploid genome. This study was conducted in order to analyze the S-genotypes of 30 diploid P. salicina, 17 of them reported here for the first time, and 29 hexaploid plums (24 of P. domestica and 5 of P. insititia). PCR analysis allowed identifying nine S-alleles in the P. salicina samples allocating the 30 accessions in 16 incompatibility groups, two of them identified here for the first time. In addition, pollen tube growth was studied in self-pollinated flowers of 17 Tunisian P. salicina under the microscope. In 16 samples, including one carrying the Se allele, which has been correlated with self-compatibility, the pollen tubes were arrested in the style. Only in one cultivar ("Bedri"), the pollen tubes reached the base of the style. Twelve S-alleles were identified in the 24 P. domestica and 5 P. insititia accessions, assigning accessions in 16 S-genotypes. S-genotyping results were combined with nine SSR loci to analyze genetic diversity. Results showed a close genetic relationship between P. domestica and P. salicina and between P. domestica and P. insititia corroborating that S-locus genotyping is suitable for molecular fingerprinting in diploid and polyploid Prunus species.

Combination of Simple Sequence Repeat, S-Locus Polymorphism and Phenotypic Data for Identification of Tunisian Plum Species (Prunus spp.).

Plums (Prunus spp.) are among the first fruit tree species that attracted human interest. Artificial crosses between wild and domesticated species of plums are still paving the way for creation of new phenotypic variability. In Tunisia, despite a considerable varietal richness of plum as well as a high economic value, the plum sector is experiencing a significant regression. The main reason of this regression is the absence of a national program of plum conservation. Hence, this work was aimed to phenotypically and genetically characterize 23 Tunisian plum accessions to preserve this patrimony. Closely related Prunus species from the same subgenus may be differing at two characteristics: ploidy level and phenotypic traits. In this study, single sequence repeat (SSR) markers allowed distinguishing between eighteen diploid accessions and five polyploid accessions, but SSR data alone precluded unambiguous ploidy estimation due to homozygosity. In contrast, S-allele markers were useful to identify the ploidy level between polyploid species, but they did not distinguish species with the same ploidy level. Seven out of 12 phenotypic traits were shown to be discriminant traits for plum species identification. Molecular and phenotypic traits were significantly correlated and revealed a powerful tool to draw taxonomic and genotypic keys. The results obtained in this work are of great importance for local Tunisian plum germplasm management.

Genotyping by sequencing reveals the interspecific C. maxima / C. reticulata admixture along the genomes of modern citrus varieties of mandarins, tangors, tangelos, orangelos and grapefruits.

The mandarin horticultural group is an important component of world citrus production for the fresh fruit market. This group formerly classified as C. reticulata is highly polymorphic and recent molecular studies have suggested that numerous cultivated mandarins were introgressed by C. maxima (the pummelos). C. maxima and C. reticulata are also the ancestors of sweet and sour oranges, grapefruit, and therefore of all the "small citrus" modern varieties (mandarins, tangors, tangelos) derived from sexual hybridization between these horticultural groups. Recently, NGS technologies have greatly modified how plant evolution and genomic structure are analyzed, moving from phylogenetics to phylogenomics. The objective of this work was to develop a workflow for phylogenomic inference from Genotyping By Sequencing (GBS) data and to analyze the interspecific admixture along the nine citrus chromosomes for horticultural groups and recent varieties resulting from the combination of the C. reticulata and C. maxima gene pools. A GBS library was established from 55 citrus varieties, using the ApekI restriction enzyme and selective PCR to improve the read depth. Diagnostic polymorphisms (DPs) of C. reticulata/C. maxima differentiation were identified and used to decipher the phylogenomic structure of the 55 varieties. The GBS approach was powerful and revealed 30,289 SNPs and 8,794 Indels with 12.6% of missing data. 11,133 DPs were selected covering the nine chromosomes with a higher density in genic regions. GBS combined with the detection of DPs was powerful for deciphering the "phylogenomic karyotypes" of cultivars derived from admixture of the two ancestral species after a limited number of interspecific recombinations. All the mandarins, mandarin hybrids, tangelos and tangors analyzed displayed introgression of C. maxima in different parts of the genome. C. reticulata/C. maxima admixture should be a major component of the high phenotypic variability of this germplasm opening up the way for association studies based on phylogenomics.

Towards a molecular taxonomic key of the Aurantioideae subfamily using chloroplastic SNP diagnostic markers of the main clades genotyped by competitive allele-specific PCR.

BACKGROUND: Chloroplast DNA is a primary source of molecular variations for phylogenetic analysis of photosynthetic eukaryotes. However, the sequencing and analysis of multiple chloroplastic regions is difficult to apply to large collections or large samples of natural populations. The objective of our work was to demonstrate that a molecular taxonomic key based on easy, scalable and low-cost genotyping method should be developed from a set of Single Nucleotide Polymorphisms (SNPs) diagnostic of well-established clades. It was applied to the Aurantioideae subfamily, the largest group of the Rutaceae family that includes the cultivated citrus species. RESULTS: The publicly available nucleotide sequences of eight plastid genomic regions were compared for 79 accessions of the Aurantioideae subfamily to search for SNPs revealing taxonomic differentiation at the inter-tribe, inter-subtribe, inter-genus and interspecific levels. Diagnostic SNPs (DSNPs) were found for 46 of the 54 clade levels analysed. Forty DSNPs were selected to develop KASPar markers and their taxonomic value was tested by genotyping 108 accessions of the Aurantioideae subfamily. Twenty-seven markers diagnostic of 24 clades were validated and they displayed a very high rate of transferability in the Aurantioideae subfamily (only 1.2 % of missing data on average). The UPGMA from the validated markers produced a cladistic organisation that was highly coherent with the previous phylogenetic analysis based on the sequence data of the eight plasmid regions. In particular, the monophyletic origin of the "true citrus" genera plus Oxanthera was validated. However, some clarification remains necessary regarding the organisation of the other wild species of the Citreae tribe. CONCLUSIONS: We validated the concept that with well-established clades, DSNPs can be selected and efficiently transformed into competitive allele-specific PCR markers (KASPar method) allowing cost-effective highly efficient cladistic analysis in large collections at subfamily level. The robustness of this genotyping method is an additional decisive advantage for network collaborative research. The availability of WGS data for the main "true citrus" species should soon make it possible to develop a set of DSNP markers allowing very fine resolution of this very important horticultural group.

SSR marker-assisted screening of commercial tomato genotypes under salt stress.

Salt stress was applied to tomato commercial genotypes to study adverse effects on their phenotypic traits. Three were saline tolerant (San Miguel, Romelia and Llanero), two were mildly tolerant (Perfect peel HF1 and Heinz 1350) whereas the remaining were sensitive. Genotyping cultivars using 19 polymorphic SSRs out of 25 tested produced a total of 70 alleles with an average of 3.68 alleles per locus and PIC values ranging from 0.22 (SSR 26, 92, 66 and TG35) to 0.82 (SSR 356). Principal component analysis (PCA) showed two contrasting panels discriminating tolerant and sensitive groups and one panel with scattered genotypes. STRUCTURE analysis clustered genotypes within three groups in accordance with their salt stress behavior. The success of tomato salt-tolerance breeding programs can be enhanced through molecular characterization of diversity within commercial cultivars that adapt differently to stress conditions. To this end, we combined phenotypes and SSR marker-genotypes to seek sources of salt tolerance that might be tomato species-specific. We integrated and represented genotype-phenotype associations from multiple loci into a multi-layer network representation. It is a systemic view linking discriminating genotypes to salt stress phenotypes, which may guide strategies for the introgression of valuable traits in target tomato varieties to overcome salinity.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations.

Since the publication of this paper, it has become apparent that an error was made in the legend to Fig. 3 and the colours referring to occidental and oriental are the wrong way round. The authors apologise for this error, and a correct version of the legend to Fig. 3 is given below.

Genetic structure of the date palm (Phoenix dactylifera) in the Old World reveals a strong differentiation between eastern and western populations.

BACKGROUND AND AIMS: Date palms (Phoenix dactylifera, Arecaceae) are of great economic and ecological value to the oasis agriculture of arid and semi-arid areas. However, despite the availability of a large date palm germplasm spreading from the Atlantic shores to Southern Asia, improvement of the species is being hampered by a lack of information on global genetic diversity and population structure. In order to contribute to the varietal improvement of date palms and to provide new insights on the influence of geographic origins and human activity on the genetic structure of the date palm, this study analysed the diversity of the species. METHODS: Genetic diversity levels and population genetic structure were investigated through the genotyping of a collection of 295 date palm accessions ranging from Mauritania to Pakistan using a set of 18 simple sequence repeat (SSR) markers and a plastid minisatellite. KEY RESULTS: Using a Bayesian clustering approach, the date palm genotypes can be structured into two different gene pools: the first, termed the Eastern pool, consists of accessions from Asia and Djibouti, whilst the second, termed the Western pool, consists of accessions from Africa. These results confirm the existence of two ancient gene pools that have contributed to the current date palm diversity. The presence of admixed genotypes is also noted, which points at gene flows between eastern and western origins, mostly from east to west, following a human-mediated diffusion of the species. CONCLUSIONS: This study assesses the distribution and level of genetic diversity of accessible date palm resources, provides new insights on the geographic origins and genetic history of the cultivated component of this species, and confirms the existence of at least two domestication origins. Furthermore, the strong genetic structure clearly established here is a prerequisite for any breeding programme exploiting the effective polymorphism related to each gene pool.

Male-specific DNA markers provide genetic evidence of an XY chromosome system, a recombination arrest and allow the tracing of paternal lineages in date palm.

Whether sex chromosomes are differentiated is an important aspect of our knowledge of dioecious plants, such as date palm (Phoenix dactylifera). In this crop plant, the female individuals produce dates, and are thus the more valuable sex. However, there is no way to identify the sex of date palm plants before reproductive age, and the sex-determining mechanism is still unclear. To identify sex-linked microsatellite markers, we surveyed a set of 52 male and 55 female genotypes representing the geographical diversity of the species. We found three genetically linked loci that are heterozygous only in males. Male-specific alleles allowed us to identify the gender in 100% of individuals. These results confirm the existence of an XY chromosomal system with a nonrecombining XY-like region in the date palm genome. The distribution of Y haplotypes in western and eastern haplogroups allowed us to trace two male ancestral paternal lineages that account for all known Y diversity in date palm. The very low diversity associated with Y haplotypes is consistent with clonal paternal transmission of a nonrecombining male-determining region. Our results establish the date palm as a biological model with one of the most ancient sex chromosomes in flowering plants.

Development of molecular tools for characterization and genetic diversity analysis in Tunisian fig (Ficus carica) cultivars.

Fig, Ficus carica L., is a useful genetic resource for commercial cultivation. In this study, RAPD (60), ISSR (48), RAMPO (63), and SSR (34) markers were compared to detect polymorphism and to establish genetic relationships among Tunisian fig tree cultivars. The statistical procedures conducted on the combined data show considerable genetic diversity, and the tested markers discriminated all fig genotypes studied. The identification key established on the basis of SSR permitted the unambiguous discrimination of cultivars and confirmed the reliability of SSR for fingerprinting fig genotypes. The study findings are discussed in relation to the establishment of a national reference collection that will aid in the conservation of Tunisian fig resources.

Genetic diversity of different Tunisian fig (Ficuscarica L.) collections revealed by RAPD fingerprints.

The genetic diversity in Tunisian fig (Ficus carica L.) was studied using RAPD markers. Thirty-five fig cultivars originating from diverse geographical areas and belonging to three collections were analysed. Random decamer primers were screened to assess their ability to detect polymorphisms in this crop. Forty-four RAPD markers were revealed and used to survey the genetic diversity and to detect cases of mislabelling. As a result, considerable genetic diversity was detected among the studied F. carica accessions. The relationships among the 35 varieties were studied by cluster analysis. The dendrogram showed two main groups composed of cultivars with similar geographic origin. Moreover, the male accessions (caprifigs) were clustered indistinctively within the female ones, suggesting a narrow genetic diversity among these accessions. Our data proved that RAPD markers are useful for germplasm discrimination as well as for investigation of patterns of variation in fig. Since this designed procedure has permitted to establish a molecular database of the reference collections, the opportunity of this study is discussed in relation to the improvement and rational management of fig germplasm.