Epigenetic map and genetic map basis of complex traits in cassava population

Recent Advances in Studies of Genomic DNA Methylation and Its Involvement in Regulating Drought Stress Response in Crops

As global arid conditions worsen and groundwater resources diminish, drought stress has emerged as a critical impediment to plant growth and development globally, notably causing declines in crop yields and even the extinction of certain cultivated species. Numerous studies on drought resistance have demonstrated that DNA methylation dynamically interacts with plant responses to drought stress by modulating gene expression and developmental processes. However, the precise mechanisms underlying these interactions remain elusive. This article consolidates the latest research on the role of DNA methylation in plant responses to drought stress across various species, focusing on methods of methylation detection, mechanisms of methylation pattern alteration (including DNA de novo methylation, DNA maintenance methylation, and DNA demethylation), and overall responses to drought conditions. While many studies have observed significant shifts in genome-wide or gene promoter methylation levels in drought-stressed plants, the identification of specific genes and pathways involved remains limited. This review aims to furnish a reference for detailed research into plant responses to drought stress through epigenetic approaches, striving to identify drought resistance genes regulated by DNA methylation, specific signaling pathways, and their molecular mechanisms of action.

Current status and impending progress for cassava structural genomics

We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Identification and validation of QTLs for tuber quality related traits in greater yam (Dioscorea alata L.)

Two Dioscorea alata populations were generated by hand pollination between contrasted diploid genitors. Population A (74F × Kabusa) was composed of 121 progenies while population B (74F × 14M) involved 193 progenies. These two populations were assessed over two consecutive years regarding important tuber quality traits. Analysis of variance showed that the genotype had the greatest influence on the phenotypic scores. Also for some traits, effect of the year_replicate was strong. The heritabilities of most traits were high. Based on these data and a reference high-density genetic map of greater yam, a total of 34 quantitative trait loci (QTLs) were detected on 8 of the 20 yam chromosomes. They corresponded to five of each of the following traits: tuber size, shape regularity, tubercular roots, skin texture, tuber flesh oxidation, six for oxidation ratio and three for flesh colour. The fraction of total phenotypic variance attributable to a single QTL ranged from 11.1 to 43.5%. We detected significant correlations between traits and QTL colocalizations that were consistent with these correlations. A majority of QTLs (62%) were found on linkage group LG16, indicating that this chromosome could play a major role in genetic control of the investigated traits. In addition, an inversion involving this chromosome was detected in the Kabusa male. Nine QTLs were validated on a diversity panel, including three for tuber size, three for shape regularity, two for skin texture and one for tubercular roots. The approximate physical localization of validated QTLs allowed the identification of various candidates genes. The validated QTLs should be useful for breeding programs using marker-assisted selection to improve yam tuber quality.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

Cassava breeding and agronomy in Asia: 50 years of history and future directions

In Asia, cassava (Manihot esculenta) is cultivated by more than 8 million farmers, driving the rural economy of many countries. The International Center for Tropical Agriculture (CIAT), in partnership with national agricultural research institutes (NARIs), instigated breeding and agronomic research in Asia, 1983. The breeding program has successfully released high-yielding cultivars resulting in an average yield increase from 13.0 t ha-1 in 1996 to 21.3 t ha-1 in 2016, with significant economic benefits. Following the success in increasing yields, cassava breeding has turned its focus to higher-value traits, such as waxy cassava, to reach new market niches. More recently, building resistance to invasive pests and diseases has become a top priority due to the emergent threat of cassava mosaic disease (CMD). The agronomic research involves driving profitability with advanced technologies focusing on better agronomic management practices thereby maintaining sustainable production systems. Remote sensing technologies are being tested for trait discovery and large-scale field evaluation of cassava. In summary, cassava breeding in Asia is driven by a combination of food and market demand with technological innovations to increase the productivity. Further, exploration in the potential of data-driven agriculture is needed to empower researchers and producers for sustainable advancement.

Construction of an ultrahigh-density genetic linkage map for Jatropha curcas L. and identification of QTL for fruit yield

BACKGROUND

As an important biofuel plant, the demand for higher yield Jatropha curcas L. is rapidly increasing. However, genetic analysis of Jatropha and molecular breeding for higher yield have been hampered by the limited number of molecular markers available.

RESULTS

An ultrahigh-density linkage map for a Jatropha mapping population of 153 individuals was constructed and covered 1380.58 cM of the Jatropha genome, with average marker density of 0.403 cM. The genetic linkage map consisted of 3422 SNP and indel markers, which clustered into 11 linkage groups. With this map, 13 repeatable QTLs (reQTLs) for fruit yield traits were identified. Ten reQTLs, qNF-1, qNF-2a, qNF-2b, qNF-2c, qNF-3, qNF-4, qNF-6, qNF-7a, qNF-7b and qNF-8, that control the number of fruits (NF) mapped to LGs 1, 2, 3, 4, 6, 7 and 8, whereas three reQTLs, qTWF-1, qTWF-2 and qTWF-3, that control the total weight of fruits (TWF) mapped to LGs 1, 2 and 3, respectively. It is interesting that there are two candidate critical genes, which may regulate Jatropha fruit yield. We also identified three pleiotropic reQTL pairs associated with both the NF and TWF traits.

CONCLUSION

This study is the first to report an ultrahigh-density Jatropha genetic linkage map construction, and the markers used in this study showed great potential for QTL mapping. Thirteen fruit-yield reQTLs and two important candidate genes were identified based on this linkage map. This genetic linkage map will be a useful tool for the localization of other economically important QTLs and candidate genes for Jatropha.

Genome-Wide Association Studies of 11 Agronomic Traits in Cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is a major tuberous crop produced worldwide. In this study, we sequenced 158 diverse cassava varieties and identified 349,827 single-nucleotide polymorphisms (SNPs) and indels. In each chromosome, the number of SNPs and the physical length of the respective chromosome were in agreement. Population structure analysis indicated that this panel can be divided into three subgroups. Genetic diversity analysis indicated that the average nucleotide diversity of the panel was 1.21 × 10^-4 for all sampled landraces. This average nucleotide diversity was 1.97 × 10^-4, 1.01 × 10^-4, and 1.89 × 10^-4 for subgroups 1, 2, and 3, respectively. Genome-wide linkage disequilibrium (LD) analysis demonstrated that the average LD was about ∼8 kb. We evaluated 158 cassava varieties under 11 different environments. Finally, we identified 36 loci that were related to 11 agronomic traits by genome-wide association analyses. Four loci were associated with two traits, and 62 candidate genes were identified in the peak SNP sites. We found that 40 of these genes showed different expression profiles in different tissues. Of the candidate genes related to storage roots, Manes.13G023300, Manes.16G000800, Manes.02G154700, Manes.02G192500, and Manes.09G099100 had higher expression levels in storage roots than in leaf and stem; on the other hand, of the candidate genes related to leaves, Manes.05G164500, Manes.05G164600, Manes.04G057300, Manes.01G202000, and Manes.03G186500 had higher expression levels in leaves than in storage roots and stem. This study provides basis for research on genetics and the genetic improvement of cassava.

Identification, Expression, and Functional Analysis of the Fructokinase Gene Family in Cassava

Fructokinase (FRK) proteins play important roles in catalyzing fructose phosphorylation and participate in the carbohydrate metabolism of storage organs in plants. To investigate the roles of FRKs in cassava tuber root development, seven FRK genes (MeFRK1-7) were identified, and MeFRK1-6 were isolated. Phylogenetic analysis revealed that the MeFRK family genes can be divided into α (MeFRK1, 2, 6, 7) and β (MeFRK3, 4, 5) groups. All the MeFRK proteins have typical conserved regions and substrate binding residues similar to those of the FRKs. The overall predicted three-dimensional structures of MeFRK1-6 were similar, folding into a catalytic domain and a β-sheet ''lid" region, forming a substrate binding cleft, which contains many residues involved in the binding to fructose. The gene and the predicted three-dimensional structures of MeFRK3 and MeFRK4 were the most similar. MeFRK1-6 displayed different expression patterns across different tissues, including leaves, stems, tuber roots, flowers, and fruits. In tuber roots, the expressions of MeFRK3 and MeFRK4 were much higher compared to those of the other genes. Notably, the expression of MeFRK3 and MeFRK4 as well as the enzymatic activity of FRK were higher at the initial and early expanding tuber stages and were lower at the later expanding and mature tuber stages. The FRK activity of MeFRK3 and MeFRK4 was identified by the functional complementation of triple mutant yeast cells that were unable to phosphorylate either glucose or fructose. The gene expression and enzymatic activity of MeFRK3 and MeFRK4 suggest that they might be the main enzymes in fructose phosphorylation for regulating the formation of tuber roots and starch accumulation at the tuber root initial and expanding stages.