Development of novel g-SSR markers in guava (Psidium guajava L.) cv. Allahabad Safeda and their application in genetic diversity, population structure and cross species transferability studies

Genetic analyses and dispersal patterns unveil the Amazonian origin of guava domestication

Guava (Psidium guajava L.) is a semi-domesticated fruit tree of moderate importance in the Neotropics, utilized for millennia due to its nutritional and medicinal benefits, but its origin of domestication remains unknown. In this study, we examine genetic diversity and population structure in 215 plants from 11 countries in Mesoamerica, the Andes, and Amazonia using 25 nuclear microsatellite loci to propose an origin of domestication. Genetic analyses reveal one gene pool in Mesoamerica (Mexico) and four in South America (Brazilian Amazonia, Peruvian Amazonia and Andes, and Colombia), indicating greater differentiation among localities, possibly due to isolation between guava populations, particularly in the Amazonian and Andean regions. Moreover, Mesoamerican populations show high genetic diversity, with moderate genetic structure due to gene flow from northern South American populations. Dispersal scenarios suggest that Brazilian Amazonia is the probable origin of guava domestication, spreading from there to the Peruvian Andes, northern South America, Central America, and Mexico. These findings present the first evidence of guava domestication in the Americas, contributing to a deeper understanding of its evolutionary history.

A high-density linkage map construction in guava (Psidium guajava L.) using genotyping by sequencing and identification of QTLs for leaf, peel, and pulp color in an intervarietal mapping population

Psidium guajava L. is an important fruit crop in the tropical and subtropical regions of the world. The advanced breeding methods are not employed for important commercial traits like peel and pulp color, seed hardiness, fruit size, etc., due to the scarcity of genome-wide molecular markers and high-density linkage maps. In this study, we employed single-nucleotide polymorphism (SNP) markers and identified quantitative trait loci (QTL) regions that are associated with color traits of leaf, peel, and pulp in the guava intervarietal mapping population. The mapping population was developed from the contrasting genotypes of fruit and leaf color. Variations in color among the segregating hybrids were recorded both visually and using a Color reader. A high-density linkage map of guava was constructed using the SNP markers from genotyping by sequencing (GBS) of 150 hybrid individuals of the cross 'Arka Poorna' (green) x 'Purple Local' (purple). The integrated linkage map consisted of 1426 SNPs mapped on 11 linkage groups (LG), spanning a total distance of around 730 cM with an average of 129.6 markers per LG. Through QTL analysis for color traits, a minor QTL region was identified for visually scored leaf color and peel color on LG1, whereas a major QTL was detected for pulp color in LG4. The Hunter color values (L* and, a*) also had major QTLs with overlapping marker intervals for leaf and peel colors, establishing the association of SNP markers to the trait. The QTLs harbored genes and transcription factors involved in lycopene and anthocyanin pigment biosynthesis. This is the first report of a high-density linkage map based on SNP markers in guava and QTL mapping for color characters in leaf, fruit peel and pulp. The genotyping information generated in this study can aid in genetic engineering and marker-assisted breeding in guava.

Assessing Genetic Variation among Strychnos spinosa Lam. Morphotypes Using Simple Sequence Repeat Markers

Strychnos spinosa Lam., commonly known as green monkey orange, is a highly valued indigenous fruit tree in South Africa with potential for domestication and commercialization. However, no study has reported on the molecular diversity of Strychnos spinosa morphotypes. Therefore, this study aimed to determine genetic variation among 32 Strychnos spinosa morphotypes using simple sequence repeat (SSR) markers. Fourteen amplified SSR markers produced 159 alleles, with a mean of 5.68 per locus. The polymorphic information content (PIC) values ranged from 0.22 (Ssp_1) to 0.84 (Ssp_6). Morphotypes were clustered in a biplot based on their genetic distances. The dendrogram chiefly discriminated morphotypes according to variation of pericarp texture. The population structure had the highest delta value K = 3, thus the 32 morphotypes were divided into three subpopulations based on the Bayesian approach. The affinities produced by the population structure agreed with the genetic distance of closely related morphotypes. This study is the first to report on SSR marker development and their successful use for genetic diversity and population structure studies of Strychnos spinosa. It provides insights into the molecular characterisation of Strychnos spinosa. This can lead to breeding programs and crop improvement programs, particularly in varietal developmental programs, which can contribute to alleviating food security challenges.

Construction of a genetic linkage map and QTL mapping of fruit quality traits in guava (Psidium guajava L.)

Guava (Psidium guajava L.) is an important fruit crop of the Indian sub-continent, with potential for improvements in quality and yield. The goal of the present study was to construct a genetic linkage map in an intraspecific cross between the elite cultivar 'Allahabad Safeda' and the Purple Guava landrace to identify the genomic regions responsible for important fruit quality traits, viz., total soluble solids, titratable acidity, vitamin C, and sugars. This population was phenotyped in field trials (as a winter crop) for three consecutive years, and showed moderate-to-high values of heterogeneity coefficients along with higher heritability (60.0%-97.0%) and genetic-advance-over-mean values (13.23%-31.17%), suggesting minimal environmental influence on the expression of fruit-quality traits and indicating that these traits can be improved by phenotypic selection methods. Significant correlations and strong associations were also detected among fruit physico-chemical traits in segregating progeny. The constructed linkage map consisted of 195 markers distributed across 11 chromosomes, spanning a length of 1,604.47 cM (average inter-loci distance of 8.80 markers) and with 88.00% coverage of the guava genome. Fifty-eight quantitative trait loci (QTLs) were detected in three environments with best linear unbiased prediction (BLUP) values using the composite interval mapping algorithm of the BIP (biparental populations) module. The QTLs were distributed on seven different chromosomes, explaining 10.95%-17.77% of phenotypic variance, with the highest LOD score being 5.96 for qTSS.AS.pau-6.2. Thirteen QTLs detected across multiple environments with BLUPs indicate stability and utility in a future breeding program for guava. Furthermore, seven QTL clusters with stable or common individual QTLs affecting two or more different traits were located on six linkage groups (LGs), explaining the correlation among fruit-quality traits. Thus, the multiple environmental evaluations conducted here have increased our understanding of the molecular basis of phenotypic variation, providing the basis for future high-resolution fine-mapping and paving the way for marker-assisted breeding of fruit-quality traits.

Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit “apple of the tropics”, and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

Construction of Genetic Linkage Map and Mapping QTL Specific to Leaf Anthocyanin Colouration in Mapping Population 'Allahabad Safeda' × 'Purple Guava (Local)' of Guava (Psidium guajava L.)

In the present investigation, F1 hybrids were developed in guava (Psidium guajava L.) by crossing high leaf-anthocyanin reflective-index (ARI1) content cultivars purple guava (local) 'PG', 'CISH G-1' and low leaf-ARI1 content cultivar Seedless 'SL' with Allahabad Safeda 'AS'. On the basis of phenotypic observations, high ARI1 content was observed in the cross 'AS' × 'PG' (0.214). Further, an SSR-markers-based genetic linkage map was developed from a mapping population of 238 F1 individuals derived from cross 'AS' × 'PG'. The linkage map comprised 11 linkage groups (LGs), spanning 1601.7 cM with an average marker interval distance of 29.61 cM between adjacent markers. Five anthocyanin-content related gene-specific markers from apple were tested for parental polymorphism in the genotypes 'AS' and 'PG'. Subsequently, a marker, viz., 'MdMYB10F1', revealed a strong association with leaf anthocyanin content in the guava mapping population. QTL (qARI-6-1) on LG6 explains much of the variation (PVE = 11.51% with LOD = 4.67) in levels of leaf anthocyanin colouration. This is the first report of amplification/utilization of apple anthocyanin-related genes in guava. The genotypic data generated from the genetic map can be further exploited in future for the enrichment of linkage maps and for identification of complex quantitative trait loci (QTLs) governing economically important fruit quality traits in guava.

A Pipeline for the Development of Microsatellite Markers using Next Generation Sequencing Data

Background: Also known as Simple Sequence Repetitions (SSRs), microsatellites are profoundly informative molecular markers and powerful tools in genetics and ecology studies on plants. Objective: This research presents a workflow for developing microsatellite markers using genome skimming. Methods: The pipeline was proposed in several stages that must be performed sequentially: obtaining DNA sequences, identifying microsatellite regions, designing primers, and selecting candidate microsatellite regions to develop the markers. Our pipeline efficiency was analyzed using Illumina sequencing data from the non-model tree species Pterodon emarginatus Vog. Results: The pipeline revealed 4,382 microsatellite regions and drew 7,411 pairs of primers for P. emarginatus. However, a much larger number of microsatellite regions with the potential to develop markers were discovered from our pipeline. We selected 50 microsatellite regions with high potential for developing markers and organized 29 microsatellite regions in sets for multiplex PCR. Conclusion: The proposed pipeline is a powerful tool for fast and efficient development of microsatellite markers on a large scale in several species, especially nonmodel plant species.

Genomic Approaches for Improvement of Tropical Fruits: Fruit Quality, Shelf Life and Nutrient Content

The breeding of tropical fruit trees for improving fruit traits is complicated, due to the long juvenile phase, generation cycle, parthenocarpy, polyploidy, polyembryony, heterozygosity and biotic and abiotic factors, as well as a lack of good genomic resources. Many molecular techniques have recently evolved to assist and hasten conventional breeding efforts. Molecular markers linked to fruit development and fruit quality traits such as fruit shape, size, texture, aroma, peel and pulp colour were identified in tropical fruit crops, facilitating Marker-assisted breeding (MAB). An increase in the availability of genome sequences of tropical fruits further aided in the discovery of SNP variants/Indels, QTLs and genes that can ascertain the genetic determinants of fruit characters. Through multi-omics approaches such as genomics, transcriptomics, metabolomics and proteomics, the identification and quantification of transcripts, including non-coding RNAs, involved in sugar metabolism, fruit development and ripening, shelf life, and the biotic and abiotic stress that impacts fruit quality were made possible. Utilizing genomic assisted breeding methods such as genome wide association (GWAS), genomic selection (GS) and genetic modifications using CRISPR/Cas9 and transgenics has paved the way to studying gene function and developing cultivars with desirable fruit traits by overcoming long breeding cycles. Such comprehensive multi-omics approaches related to fruit characters in tropical fruits and their applications in breeding strategies and crop improvement are reviewed, discussed and presented here.

Development of Genome-Wide Functional Markers Using Draft Genome Assembly of Guava (Psidium guajava L.) cv. Allahabad Safeda to Expedite Molecular Breeding

Guava (Psidium guajava L.), a rich source of nutrients, is an important tropical and subtropical fruit of the Myrtaceae family and exhibits magnificent diversity. Genetic diversity analysis is the first step toward the identification of parents for hybridization, genetic mapping, and molecular breeding in any crop species. A diversity analysis based on whole-genome functional markers increases the chances of identifying genetic associations with agronomically important traits. Therefore, here, we sequenced the genome of guava cv. Allahabad Safeda on an Illumina platform and generated a draft assembly of ~304 MB. The assembly of the Allahabad Safeda genome constituted >37.95% repeat sequences, gene prediction with RNA-seq data as evidence identified 14,115 genes, and BLAST n/r, Interproscan, PfamScan, BLAST2GO, and KEGG annotated 13,957 genes. A comparative protein transcript analysis of tree species revealed the close relatedness of guava with Eucalyptus. Comparative transcriptomics-based SSR/InDel/SNP-PCR ready genome-wide markers in greenish-yellow skinned and white fleshed-Allahabad Safeda to four contrasting cultivars viz apple-color-skinned and white-fleshed-Lalima, greenish-yellow-skinned and pink-fleshed-Punjab Pink, purple-black-skinned and purple-fleshed-Purple Local and widely used rootstock-Lucknow-49 were developed. The molecular markers developed here revealed a high level of individual heterozygosity within genotypes in 22 phenotypically diverse guava cultivars. Principal coordinate, STRUCTURE clustering, and neighbor-joining-based genetic diversity analysis identified distinct clusters associated with fruit skin and flesh color. The genome sequencing of guava, functional annotation, comparative transcriptomics-based genome-wide markers, and genetic diversity analysis will expand the knowledge of genomes of climacteric fruits, facilitating trait-based molecular breeding and diversifying the nutritional basket.