Morphological, Biochemical, and Molecular Diversity of an Indian Ex Situ Collection of Pomegranate (Punica granatum L.).

Pomegranate (Punica granatum, L.) is a fruit tree that is increasingly popular worldwide due to the health-related properties of the fruit juice. While several studies highlighted the rich phytochemical diversity, few efforts have been devoted to an integrative understanding of the level of diversity of this species. This study investigated the diversity of 40 pomegranate accessions in an Indian ex situ collection by using twenty-nine morphological traits, six biochemical parameters, and twenty-nine Simple Sequence Repeats (SSR) markers. Among the evaluated traits, fruit volume (23.34% CV), fruit weight (21.12% CV), and fruit color (*a) (22.69 % CV) largely contributed to the morphological classification. Based on Mahalanobis D2 distance and Tocher's clustering, the 40 pomegranate accessions were grouped into eight clusters, partly consistent with their origin. Specifically, cultivars introduced from foreign countries were present in distinct clusters. The SSR marker analysis generated 66 alleles. The observed heterozygosity values ranged from 0.05 to 0.63, with a mean value of 0.30. Maximum molecular genetic dissimilarity was observed between 'IC-318720' and 'Gul-e-Shah Red' (0.30). The neighbor-joining dendrogram separated wild accessions from cultivated varieties. The combination of morphological, biochemical, and molecular characterization allowed for comprehensively characterizing the pomegranate diversity and provided information on the relationships between the different aspects of the diversity. This work also suggests that the origin of the accessions is an important factor of discrimination and that the level of admixture between local and foreign material is currently limited.

Corrigendum: Chloroplast genome sequencing, comparative analysis, and discovery of unique cytoplasmic variants in pomegranate (Punica granatum L.).

[This corrects the article DOI: 10.3389/fgene.2021.704075.].

Novel miRNA-SSRs for Improving Seed Hardness Trait of Pomegranate (Punica granatum L.).

Present research discovered novel miRNA-SSRs for seed type trait from 761 potential precursor miRNA sequences of pomegranate. SSR mining and BLASTx of the unique sequences identified 69 non-coding pre-miRNA sequences, which were then searched for BLASTn homology against Dabenzi genome. Sixty three true pri-miRNA contigs encoding 213 pre-miRNAs were predicted. Analysis of the resulting sequences enabled discovery of SSRs within pri-miRNA (227) and pre-miRNA sequences (79). A total of 132 miRNA-SSRs were developed for seed type trait from 63 true pri-miRNAs, of which 46 were specific to pre-miRNAs. Through ePCR, 123 primers were validated and mapped on eight Tunisia chromosomes. Further, 80 SSRs producing specific amplicons were ePCR-confirmed on multiple genomes i.e. Dabenzi, Taishanhong, AG2017 and Tunisia, yielding a set of 63 polymorphic SSRs (polymorphism information content ≥0.5). Of these, 32 miRNA-SSRs revealed higher polymorphism level (89.29%) when assayed on six pomegranate genotypes. Furthermore, target prediction and network analysis suggested a possible association of miRNA-SSRs i.e. miRNA_SH_SSR69, miRNA_SH_SSR36, miRNA_SH_SSR103, miRNA_SH_SSR35 and miRNA_SH_SSR53 with seed type trait. These miRNA-SSRs would serve as important genomic resource for rapid and targeted improvement of seed type trait of pomegranate.

In Vitro and In Planta Antagonistic Effect of Endophytic Bacteria on Blight Causing Xanthomonas axonopodis pv. punicae: A Destructive Pathogen of Pomegranate

Pomegranate bacterial blight caused by Xanthomonas axonopodis pv. punicae (Xap) is a highly destructive disease. In the absence of host resistance to the disease, we aimed to evaluate the biocontrol potential of endophytic bacteria against Xap. Thus, in this study, we isolated endophytes from pomegranate plants, identified them on the basis of 16S rDNA sequencing, tested them against Xap, and estimated the endophyte-mediated host defense response. The population of isolated endophytes ranged from 3 × 106 to 8 × 107 CFU/g tissue. Furthermore, 26 isolates were evaluated for their biocontrol activity against Xap, and all the tested isolates significantly reduced the in vitro growth of Xap (15.65% ± 1.25% to 56.35% ± 2.66%) as compared to control. These isolates could reduce fuscan, an uncharacterized factor of Xap involved in its aggressiveness. Lower blight incidence (11.6%) and severity (6.1%) were recorded in plants sprayed with endophytes 8 days ahead of Xap spray (Set-III) as compared to control plants which were not exposed to endophytes (77.33 and 50%, respectively%) during in vivo evaluation. Moreover, significantly high phenolic and chlorophyll contents were estimated in endophyte-treated plants as compared to control. The promising isolates mostly belonged to the genera Bacillus, Burkholderia, and Lysinibacillus, and they were deposited to the National Agriculturally Important Microbial Culture Collection, India.

Chloroplast Genome Sequencing, Comparative Analysis, and Discovery of Unique Cytoplasmic Variants in Pomegranate (Punica granatum L.).

Here we report on comprehensive chloroplast (cp) genome analysis of 16 pomegranate (Punica granatum L.) genotypes representing commercial cultivars, ornamental and wild types, through large-scale sequencing and assembling using next-generation sequencing (NGS) technology. Comparative genome analysis revealed that the size of cp genomes varied from 158,593 bp (in wild, "1201" and "1181") to 158,662 bp (cultivar, "Gul-e-Shah Red") among the genotypes, with characteristic quadripartite structures separated by a pair of inverted repeats (IRs). The higher conservation for the total number of coding and non-coding genes (rRNA and tRNA) and their sizes, and IRs (IR-A and IR-B) were observed across all the cp genomes. Interestingly, high variations were observed in sizes of large single copy (LSC, 88,976 to 89,044 bp) and small single copy (SSC, 18,682 to 18,684 bp) regions. Although, the structural organization of newly assembled cp genomes were comparable to that of previously reported cp genomes of pomegranate ("Helow," "Tunisia," and "Bhagawa"), the striking differences were observed with the Lagerstroemia lines, viz., Lagerstroemia intermedia (NC_0346620) and Lagerstroemia speciosa (NC_031414), which clearly confirmed previous findings. Furthermore, phylogenetic analysis also revealed that members outside the genus Punica were clubbed into a separate clade. The contraction and expansion analysis revealed that the structural variations in IRs, LSC, and SSC have significantly accounted for the evolution of cp genomes of Punica and L. intermedia over the periods. Microsatellite survey across cp genomes resulted in the identification of a total of 233 to 234 SSRs, with majority of them being mono- (A/T or C/G, 164-165 numbers), followed by di- (AT/AT or AG/CT, 54), tri- (6), tetra- (8), and pentanucleotides (1). Furthermore, the comparative structural variant analyses across cp genomes resulted in the identification of many varietal specific SNP/indel markers. In summary, our study has offered a successful development of large-scale cp genomics resources to leverage future genetic, taxonomical, and phylogenetic studies in pomegranate.

Comprehensive Characterization and Validation of Chromosome-Specific Highly Polymorphic SSR Markers From Pomegranate (Punica granatum L.) cv. Tunisia Genome.

The simple sequence repeat (SSR) survey of 'Tunisia' genome (296.85 Mb) identified a total of 365,279 perfect SSRs spanning eight chromosomes, with a mean marker density of 1,230.6 SSRs/Mb. We found a positive trend in chromosome length and the SSR abundance as marker density enhanced with a shorter chromosome length. The highest number of SSRs (60,708) was mined from chromosome 1 (55.56 Mb), whereas the highest marker density (1,294.62 SSRs/Mb) was recorded for the shortest chromosome 8 (27.99 Mb). Furthermore, we categorized all SSR motifs into three major classes based on their tract lengths. Across the eight chromosomes, the class III had maximum number of SSR motifs (301,684, 82.59%), followed by the class II (31,056, 8.50%) and the class I (5,003, 1.37%). Examination of the distribution of SSR motif types within a chromosome suggested the abundance of hexanucleotide repeats in each chromosome followed by dinucleotides, and these results are consistent with 'Tunisia' genome features as a whole. Concerning major repeat types, AT/AG was the most frequent (14.16%), followed by AAAAAT/AAAAAG (7.89%), A/C (7.54%), AAT/AAG (5.23%), AAAT/AAAG (4.37%), and AAAAT/AAAAG (1.2%) types. We designed and validated a total of 3,839 class I SSRs in the 'Tunisia' genome through electronic polymerase chain reaction (ePCR) and found 1,165 (30.34%) SSRs producing a single amplicon. Then, we selected 906 highly variable SSRs (> 40 nt) from the ePCR-verified class I SSRs and in silico validated across multiple draft genomes of pomegranate, which provided us a subset of 265 highly polymorphic SSRs. Of these, 235 primers were validated on six pomegranate genotypes through wet-lab experiment. We found 221 (94%) polymorphic SSRs on six genotypes, and 187 of these SSRs had ≥ 0.5 PIC values. The utility of the developed SSRs was demonstrated by analyzing genetic diversity of 30 pomegranate genotypes using 16 HvSSRs spanning eight pomegranate chromosomes. In summary, we developed a comprehensive set of highly polymorphic genome-wide SSRs. These chromosome-specific SSRs will serve as a powerful genomic tool to leverage future genetic studies, germplasm management, and genomics-assisted breeding in pomegranate.

Comparative transcriptome profiling of pomegranate genotypes having resistance and susceptible reaction to Xanthomonas axonopodis pv. punicae.

Pomegranate (Punica granatum L.) is an important fruit crop, rich in fiber, vitamins, antioxidants, minerals and source of different biologically active compounds. The bacterial blight caused by Xanthomonas axonopodispv. punicae is a serious threat to the crop leading to 60-80% yield loss under epiphytotic conditions. In this work, we have generated comparative transcriptome profile to mark the gene expression signatures during resistance and susceptible interactions. We analyzed leaf and fruits samples of moderately resistant genotype (IC 524207) and susceptible variety (Bhagawa) of pomegranate at three progressive infection stages upon inoculation with the pathogen. RNA-Seq with the Illumina HiSeq 2500 platform revealed 1,88,337 non-redundant (nr) transcript sequences from raw sequencing data, for a total of 34,626 unigenes with size >2 kb. Moreover, 85.3% unigenes were annotated in at least one of the seven databases examined. Comparative analysis of gene-expression signatures in resistant and susceptible varieties showed that the genes known to be involved in defense mechanism in plants were up-regulated in resistant variety. Gene Ontology (GO) analysis successfully annotated 90,485 pomegranate unigenes, of which 68,464 were assigned to biological, 78,107 unigenes molecular function and 44,414 to cellular components. Significantly enriched GO terms in DEGs were related to oxidations reduction biological process, protein binding and oxidoreductase activity. This transcriptome data on pomegranate could help in understanding resistance and susceptibility nature of cultivars and further detailed fine mapping and functional validation of identified candidate gene would provide scope for resistance breeding programme in pomegranate.

Genetic diversity and association mapping of bacterial blight and other horticulturally important traits with microsatellite markers in pomegranate from India.

This genetic diversity study aimed to estimate the population structure and explore the use of association mapping strategies to identify linked markers for bacterial resistance, growth and fruit quality in pomegranate collections from India. In total, 88 accessions including 37 cultivated types were investigated. A total of 112 alleles were amplified by use of 44 publicly available microsatellites for estimating molecular genetic diversity and population structure. Neighbor-joining analysis, model-based population structure and principal component analysis corroborated the genetic relationships among wild-type and cultivated pomegranate collections from India. Our study placed all 88 germplasm into four clusters. We identified a cultivated clade of pomegranates in close proximity to Daru types of wild-type pomegranates that grow naturally near the foothills of the Himalayas. Admixture analysis sorted various lineages of cultivated pomegranates to their respective ancestral forms. We identified four linked markers for fruit weight, titratable acidity and bacterial blight severity. PGCT001 was found associated with both fruit weight and bacterial blight, and the association with fruit weight during both seasons analyzed was significant after Bonferroni correction. This research demonstrates effectiveness of microsatellites to resolve population structure among the wild and cultivar collection of pomegranates and future use for association mapping studies.