Genomic Resources to Guide Improvement of the Shea Tree

Genome-wide association study of fat content and fatty acid composition of shea tree (Vitellaria paradoxa C.F. Gaertn subsp. paradoxa)

BACKGROUND

Fat content (FC) and fatty acids (FA) are the most important traits in shea tree breeding, controlled by several genes with relatively small effects. Therefore, determining the genes involved in the biosynthesis of such traits is crucial for improving oil quantity and quality and for the domestication process of the species. To identify the quantitative trait nucleotides (QTNs) controlling FC and FA, we conducted a multi-locus genome-wide association study (GWAS) using six multi-locus GWAS methods for FC and FA in 122 superior shea trees (SSTs). SSTs were genotyped using DArTseq, resulting in 7,559 non-redundant single nucleotide polymorphism markers.

RESULTS

Fat content varied from 36 to 58% with a mean of 50%. Fatty acid composition was 51.26 ± 4.21, 38.76 ± 4.67, 6.45 ± 0.76 and 3.53 ± 0.52% for oleic, stearic, linoleic and palmitic acids, respectively. A very high negative correlation coefficient (-0.98) was found between stearic and oleic acids. A total of 47 significant QTNs associated with fat-related traits were detected by the GWAS methods. Among these QTNs, 25 were identified as common QTNs based on their detection by multiple GWAS methods. Using the superior allele information of the 4 common QTNs associated with fat content in 17 high-fat and 21 low-fat SSTs, we found a higher percentage of superior alleles in SSTs with high FC (47.1%) than in SSTs with low FC (14.3%). Pathway analysis of the common QTNs identified 24 potential candidate genes likely involved in the biosynthesis of FC and FA composition in shea tree seeds.

CONCLUSIONS

These findings will contribute to the discovery of the polygenic networks controlling FC in shea tree, improve our understanding of the genetic basis and regulation of FC, and be useful for molecular breeding of high-fat shea tree cultivars.

Divergent evolutionary paces among eudicot plants revealed by simultaneously duplicated genes produced billions of years ago

Polyploidization often occurs more than once along an evolutionary lineage to form extant plants. Major core eudicot plants share a whole-genome triplication (ceWGT), through which thousands of simultaneously duplicated genes are retained in extant genomes, providing a valuable starting line to check the difference in their evolutionary paces. Here, by characterizing the synonymous nucleotide substitutions (Ks) between these duplicates from 28 representative plants from 21 families, we checked the various evolutionary rates among plants among plants subjected to different rounds of extra polyploidization events. We found up to 68.04% difference in evolutionary rates among the selected plants. A statistical correlation analysis (correlation coefficient =0.57, at significant level = 0.01) indicated that plants affected by extra polyploidies have evolved faster than plants without such extra polyploidies showing that (additional) polyploidization has resulted in elevated genetic diversity. Comparing the plants affected by additional polyploidization and plants without it, the duplicated genes produced by the ceWGT and retained in extant genomes have gathered 4.75% more nucleotide substitutions in the former plants. By identifying the fast- and slowly evolving genes, we showed that genes evolving at divergent rates were often related to different evolutionary paths. By performing correction to evolutionary rates using a genome-scale approach, we revised the estimated timing of key evolutionary events. The present effort exploited the simultaneously duplicated genes produced by the shared polyploidization and help deepen the understanding of the role of polyploidization, especially its long-term effect in plant evolution and biological innovation.

Genetic diversity and population structure of superior shea trees (Vitellaria paradoxa subsp. paradoxa) using SNP markers for the establishment of a core collection in Côte d'Ivoire

BACKGROUND

The shea tree is a well-known carbon sink in Africa that requires a sustainable conservation of its gene pool. However, the genetic structure of its population is not well studied, especially in Côte d'Ivoire. In this study, 333 superior shea tree genotypes conserved in situ in Côte d'Ivoire were collected and genotyped with the aim of investigating its genetic diversity and population structure to facilitate suitable conservation and support future breeding efforts to adapt to climate change effects.

RESULTS

A total of 7,559 filtered high-quality single nucleotide polymorphisms (SNPs) were identified using the genotyping by sequencing technology. The gene diversity (HE) ranged between 0.1 to 0.5 with an average of 0.26, while the polymorphism information content (PIC) value ranged between 0.1 to 0.5 with an average of 0.24, indicating a moderate genetic diversity among the studied genotypes. The population structure model classified the 333 genotypes into three genetic groups (GP1, GP2, and GP3). GP1 contained shea trees that mainly originated from the Poro, Tchologo, and Hambol districts, while GP2 and GP3 contained shea trees collected from the Bagoué district. Analysis of molecular variance (AMOVA) identified 55% variance within populations and 45% variance within individuals, indicating a very low genetic differentiation (or very high gene exchange) between these three groups (FST = 0.004, gene flow Nm = 59.02). Morphologically, GP1 displayed spreading tree growth habit, oval nut shape, higher mean nut weight (10.62 g), wide leaf (limb width = 4.63 cm), and small trunk size (trunk circumference = 133.4 cm). Meanwhile, GP2 and GP3 showed similar morphological characteristics: erect and spreading tree growth habit, ovoid nut shape, lower mean nut weight (GP2: 8.89 g; GP3: 8.36 g), thin leaf (limb width = 4.45 cm), and large trunk size ( GP2: 160.5 cm, GP3: 149.1 cm). A core set of 100 superior shea trees, representing 30% of the original population size and including individuals from all four study districts, was proposed using the "maximum length sub-tree function" in DARwin v. 6.0.21.

CONCLUSION

These findings provide new knowledge of the genetic diversity and population structure of Ivorian shea tree genetic resources for the design of effective collection and conservation strategies for the efficient use of inbreeding.

The Primula edelbergii S-locus is an example of a jumping supergene

Research on supergenes, non-recombining genomic regions housing tightly linked genes that control complex phenotypes, has recently gained prominence in genomics. Heterostyly, a floral heteromorphism promoting outcrossing in several angiosperm families, is controlled by the S-locus supergene. The S-locus has been studied primarily in closely related Primula species and, more recently, in other groups that independently evolved heterostyly. However, it remains unknown whether genetic architecture and composition of the S-locus are maintained among species that share a common origin of heterostyly and subsequently diverged across larger time scales. To address this research gap, we present a chromosome-scale genome assembly of Primula edelbergii, a species that shares the same origin of heterostyly with Primula veris (whose S-locus has been characterized) but diverged from it 18 million years ago. Comparative genomic analyses between these two species allowed us to show, for the first time, that the S-locus can 'jump' (i.e. translocate) between chromosomes maintaining its function in controlling heterostyly. Additionally, we found that four S-locus genes were conserved but reshuffled within the supergene, seemingly without affecting their expression, thus we could not detect changes explaining the lack of self-incompatibility in P. edelbergii. Furthermore, we confirmed that the S-locus is not undergoing genetic degeneration. Finally, we investigated P. edelbergii evolutionary history within Ericales in terms of whole genome duplications and transposable element accumulation. In summary, our work provides a valuable resource for comparative analyses aimed at investigating the genetics of heterostyly and the pivotal role of supergenes in shaping the evolution of complex phenotypes.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

An innovative optimized protocol for high-quality genomic DNA extraction from recalcitrant Shea tree (Vitellaria paradoxa, C.F. Gaertn) plant and its suitability for downstream applications

BACKGROUND

It is not always easy to find a universal protocol for the extraction of genomic DNA (gDNA) from plants. Extraction of gDNA from plants such as shea with a lot of polysaccharides in their leaves is done in two steps: a first step to remove the polysaccharides and a second step for the extraction of the gDNA. In this work, we designed a protocol for extracting high-quality gDNA from shea tree and demonstrate its suitability for downstream molecular applications.

METHODS

Fifty milligrams of leaf and root tissues were used to test the efficiency of our protocol. The quantity of gDNA was measured with the NanoDrop spectrometer and the quality was checked on agarose gel. Its suitability for use in downstream applications was tested with restriction enzymes, SSRs and RAPD polymerase chain reactions and Sanger sequencing.

RESULTS

The average yield of gDNA was 5.17; 3.96; 2.71 and 2.41 µg for dry leaves, dry roots, fresh leaves and fresh roots respectively per 100 mg of tissue. Variance analysis of the yield showed significant difference between all tissue types. Leaf gDNA quality was better compared to root gDNA at the absorbance ratio A260/280 and A260/230. The minimum amplifiable concentration of leaf gDNA was 1 pg/µl while root gDNA remained amplifiable at 10 pg/µl. Genomic DNA obtained was also suitable for sequencing.

CONCLUSION

This protocol provides an efficient, convenient and cost effective DNA extraction method suitable for use in various vitellaria paradoxa genomic studies.

Genome-wide diversity analysis suggests divergence among Upper Guinea and the Dahomey Gap populations of the Sisrè berry (Syn: miracle fruit) plant (Synsepalum dulcificum [Schumach. & Thonn.] Daniell) in West Africa

Although Synsepalum dulcificum is viewed as one of the most economically promising orphan tree crops worldwide, its genetic improvement and sustainable conservation are hindered by a lack of understanding of its evolutionary history and current population structure. Here, we report for the first time the application of genome-wide single nucleotide polymorphism genotyping to a diverse panel of S. dulcificum accessions to depict the genetic diversity and population structure of the species in the Dahomey Gap (DG) and Upper Guinea (UG) regions to infer its evolutionary history. Our findings suggest low overall genetic diversity but strong population divergence within the species. Neighbor-joining analysis detected two genetic groups in the UG and DG regions, while STRUCTURE distinguished three genetic groups, corresponding to the UG, Western DG, and Central DG regions. Application of Monmonier's algorithm revealed the existence of a barrier disrupting connectivity between the UG and DG groups. The Western DG group consistently exhibited the highest levels of nucleotide and haplotype diversities, while that of the Central DG exhibited the lowest. Analyses of Tajima's D, Fu's Fs, and Achaz Y* statistics suggest that while both UG and Central DG groups likely experienced recent expansions, the Western DG group is at equilibrium. These findings suggest a geographical structuring of genetic variation which supports the conclusion of differential evolutionary histories among West African groups of S. dulcificum. These results provide foundational insights to guide informed breeding population development and design sustainable conservation strategies for this species.

Complex genome assembly based on long-read sequencing

High-quality genome chromosome-scale sequences provide an important basis for genomics downstream analysis, especially the construction of haplotype-resolved and complete genomes, which plays a key role in genome annotation, mutation detection, evolutionary analysis, gene function research, comparative genomics and other aspects. However, genome-wide short-read sequencing is difficult to produce a complete genome in the face of a complex genome with high duplication and multiple heterozygosity. The emergence of long-read sequencing technology has greatly improved the integrity of complex genome assembly. We review a variety of computational methods for complex genome assembly and describe in detail the theories, innovations and shortcomings of collapsed, semi-collapsed and uncollapsed assemblers based on long reads. Among the three methods, uncollapsed assembly is the most correct and complete way to represent genomes. In addition, genome assembly is closely related to haplotype reconstruction, that is uncollapsed assembly realizes haplotype reconstruction, and haplotype reconstruction promotes uncollapsed assembly. We hope that gapless, telomere-to-telomere and accurate assembly of complex genomes can be truly routinely achieved using only a simple process or a single tool in the future.