Integration of genomics, transcriptomics and metabolomics identifies candidate loci underlying fruit weight in loquat

Chloroplast genomes of Eriobotrya elliptica and an unknown wild loquat "YN-1"

The chloroplast genomes of wild loquat can help to determine their place in the history of evolution. Here, we sequenced and assembled two novel wild loquat's chloroplast genomes, one is Eriobotrya elliptica, and the other is an unidentified wild loquat, which we named "YN-1". Their sizes are 159,471 bp and 159,399 bp, respectively. We also assembled a cultivated loquat named 'JFZ', its chloroplast genome size is 159,156 bp. A comparative study was conducted with six distinct species of loquats, including five wild loquats and one cultivated loquat. The results showed that both E. elliptica and "YN-1" have 127 genes, one gene more than E. fragrans, which is psbK. Regions trnF-GAA-ndhJ, petG-trnP-UGG, and rpl32-trnL-UAG were found to exhibit high variability. It was discovered that there was a positive selection on rpl22 and rps12. RNA editing analysis found several chilling stress-specific RNA editing sites, especially in rpl2 gene. Phylogenetic analysis results showed that "YN-1" is closely related to E. elliptica, E. obovata and E. henryi.

Transcriptomics reveal useful resources for examining fruit development and variation in fruit size in Coccinia grandis

Introduction

The Cucurbitaceae family comprises many agronomically important members, that bear nutritious fruits and vegetables of great economic importance. Coccinia grandis, commonly known as Ivy gourd, belongs to this family and is widely consumed as a vegetable. Members of this family are known to display an impressive range of variation in fruit morphology. Although there have been studies on flower development in Ivy gourd, fruit development remains unexplored in this crop.

Methods

In this study, comparative transcriptomics of two Ivy gourd cultivars namely "Arka Neelachal Kunkhi" (larger fruit size) and "Arka Neelachal Sabuja" (smaller fruit size) differing in their average fruit size was performed. A de novo transcriptome assembly for Ivy gourd was developed by collecting fruits at different stages of development (5, 10, 15, and 20 days after anthesis i.e. DAA) from these two varieties. The transcriptome was analyzed to identify differentially expressed genes, transcription factors, and molecular markers.

Results

The transcriptome of Ivy gourd consisted of 155205 unigenes having an average contig size of 1472bp. Unigenes were annotated on publicly available databases to categorize them into different biological functions. Out of these, 7635 unigenes were classified into 38 transcription factor (TF) families, of which Trihelix TFs were most abundant. A total of 11,165 unigenes were found to be differentially expressed in both the varieties and the in silico expression results were validated through real-time PCR. Also, 98768 simple sequence repeats (SSRs) were identified in the transcriptome of Ivy gourd.

Discussion

This study has identified a number of genes, including transcription factors, that could play a crucial role in the determination of fruit shape and size in Ivy gourd. The presence of polymorphic SSRs indicated a possibility for marker-assisted selection for crop breeding in Ivy gourd. The information obtained can help select candidate genes that may be implicated in regulating fruit development and size in other fruit crops.

Multi-omics analysis provides new insights into the changes of important nutrients and fructose metabolism in loquat bud sport mutant

Bud sport is a common and stable somatic variation in perennial fruit trees, and often leads to significant modification of fruit traits and affects the breeding value. To investigate the impact of bud sport on the main metabolites in the fruit of white-fleshed loquat, we conducted a multi-omics analysis of loquat fruits at different developmental stages of a white-fleshed bud sport mutant of Dongting loquat (TBW) and its wild type (TBY). The findings from the detection of main fruit quality indices and metabolites suggested that bud sport resulted in a reduction in the accumulation of carotenoids, fructose, titratable acid and terpenoids at the mature stage of TBW, while leading to the accumulation of flavonoids, phenolic acids, amino acids and lipids. The comparably low content of titratable acid further enhances the balanced and pleasent taste profile of TBW. Expression patterns of differentially expressed genes involved in fructose metabolism exhibited a significant increase in the expression level of S6PDH (EVM0006243, EVM0044405) prior to fruit maturation. The comparison of protein sequences and promoter region of S6PDH between TBY and TBW revealed no structural variations that would impact gene function or expression, indicating that transcription factors may be responsible for the rapid up-regulation of S6PDH before maturation. Furthermore, correlation analysis helped to construct a comprehensive regulatory network of fructose metabolism in loquat, including 23 transcription factors, six structural genes, and nine saccharides. Based on the regulatory network and existing studies, it could be inferred that transcription factors such as ERF, NAC, MYB, GRAS, and bZIP may promote fructose accumulation in loquat flesh by positively regulating S6PDH. These findings improve our understanding of the nutritional value and breeding potential of white-fleshed loquat bud sport mutant, as well as serve as a foundation for exploring the genes and transcription factors that regulate fructose metabolism in loquat.

Construction of a High-Density Paulownia Genetic Map and QTL Mapping of Important Phenotypic Traits Based on Genome Assembly and Whole-Genome Resequencing

Quantitative trait locus (QTL) mapping based on a genetic map is a very effective method of marker-assisted selection in breeding, and whole-genome resequencing is one of the useful methods to obtain high-density genetic maps. In this study, the hybrid assembly of Illumina, PacBio, and chromatin interaction mapping data was used to construct high-quality chromosomal genome sequences of Paulownia fortunei, with a size of 476.82 Mb, a heterozygosity of 0.52%, and a contig and scaffold N50s of 7.81 Mb and 21.81 Mb, respectively. Twenty scaffolds with a total length of 437.72 Mb were assembled into 20 pseudochromosomes. Repeat sequences with a total length of 243.96 Mb accounted for 51.16% of the entire genome. In all, 26,903 protein-coding gene loci were identified, and 26,008 (96.67%) genes had conserved functional motifs. Further comparative genomics analysis preliminarily showed that the split of P. fortunei with Tectona grandis likely occurred 38.8 (33.3-45.1) million years ago. Whole-genome resequencing was used to construct a merged genetic map of 20 linkage groups, with 2993 bin markers (3,312,780 SNPs), a total length of 1675.14 cm, and an average marker interval of 0.56 cm. In total, 73 QTLs for important phenotypic traits were identified (19 major QTLs with phenotypic variation explained ≥ 10%), including 10 for the diameter at breast height, 7 for the main trunk height, and 56 for branch-related traits. These results not only enrich P. fortunei genomic data but also form a solid foundation for fine QTL mapping and key marker/gene mining of Paulownia, which is of great significance for the directed genetic improvement of these species.

Comparative Anatomical and Transcriptomics Reveal the Larger Cell Size as a Major Contributor to Larger Fruit Size in Apricot

Fruit size is one of the essential quality traits and influences the economic value of apricots. To explore the underlying mechanisms of the formation of differences in fruit size in apricots, we performed a comparative analysis of anatomical and transcriptomics dynamics during fruit growth and development in two apricot cultivars with contrasting fruit sizes (large-fruit Prunus armeniaca 'Sungold' and small-fruit P. sibirica 'F43'). Our analysis identified that the difference in fruit size was mainly caused by the difference in cell size between the two apricot cultivars. Compared with 'F43', the transcriptional programs exhibited significant differences in 'Sungold', mainly in the cell expansion period. After analysis, key differentially expressed genes (DEGs) most likely to influence cell size were screened out, including genes involved in auxin signal transduction and cell wall loosening mechanisms. Furthermore, weighted gene co-expression network analysis (WGCNA) revealed that PRE6/bHLH was identified as a hub gene, which interacted with 1 TIR1, 3 AUX/IAAs, 4 SAURs, 3 EXPs, and 1 CEL. Hence, a total of 13 key candidate genes were identified as positive regulators of fruit size in apricots. The results provide new insights into the molecular basis of fruit size control and lay a foundation for future breeding and cultivation of larger fruits in apricot.

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Integration of multi-omics technologies for crop improvement: Status and prospects

With the rapid development of next-generation sequencing (NGS), multi-omics techniques have been emerging as effective approaches for crop improvement. Here, we focus mainly on addressing the current status and future perspectives toward omics-related technologies and bioinformatic resources with potential applications in crop breeding. Using a large amount of omics-level data from the functional genome, transcriptome, proteome, epigenome, metabolome, and microbiome, clarifying the interaction between gene and phenotype formation will become possible. The integration of multi-omics datasets with pan-omics platforms and systems biology could predict the complex traits of crops and elucidate the regulatory networks for genetic improvement. Different scales of trait predictions and decision-making models will facilitate crop breeding more intelligent. Potential challenges that integrate the multi-omics data with studies of gene function and their network to efficiently select desirable agronomic traits are discussed by proposing some cutting-edge breeding strategies for crop improvement. Multi-omics-integrated approaches together with other artificial intelligence techniques will contribute to broadening and deepening our knowledge of crop precision breeding, resulting in speeding up the breeding process.

Combination of Genomics, Transcriptomics Identifies Candidate Loci Related to Cold Tolerance in Dongxiang Wild Rice

Rice, a cold-sensitive crop, is a staple food for more than 50% of the world’s population. Low temperature severely compromises the growth of rice and challenges China’s food safety. Dongxiang wild rice (DXWR) is the most northerly common wild rice in China and has strong cold tolerance, but the genetic basis of its cold tolerance is still unclear. Here, we report quantitative trait loci (QTLs) analysis for seedling cold tolerance (SCT) using a high-density single nucleotide polymorphism linkage map in the backcross recombinant inbred lines that were derived from a cross of DXWR, and an indica cultivar, GZX49. A total of 10 putative QTLs were identified for SCT under 4 °C cold treatment, each explaining 2.0–6.8% of the phenotypic variation in this population. Furthermore, transcriptome sequencing of DXWR seedlings before and after cold treatment was performed, and 898 and 3413 differentially expressed genes (DEGs) relative to 0 h in cold-tolerant for 4 h and 12 h were identified, respectively. Gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) analysis were performed on these DEGs. Using transcriptome data and genetic linkage analysis, combined with qRT-PCR, sequence comparison, and bioinformatics, LOC_Os08g04840 was putatively identified as a candidate gene for the major effect locus qSCT8. These findings provided insights into the genetic basis of SCT for the improvement of cold stress potential in rice breeding programs.

Roles of Auxin in the Growth, Development, and Stress Tolerance of Horticultural Plants

Auxin, a plant hormone, regulates virtually every aspect of plant growth and development. Many current studies on auxin focus on the model plant Arabidopsis thaliana, or on field crops, such as rice and wheat. There are relatively few studies on what role auxin plays in various physiological processes of a range of horticultural plants. In this paper, recent studies on the role of auxin in horticultural plant growth, development, and stress response are reviewed to provide novel insights for horticultural researchers and cultivators to improve the quality and application of horticultural crops.

High-density genetic linkage mapping reveals low stability of QTLs across environments for economic traits in Eucalyptus

Introduction

Eucalyptus urophylla, E. tereticornis and their hybrids are the most important commercial forest tree species in South China where they are grown for pulpwood and solid wood production. Construction of a fine-scale genetic linkage map and detecting quantitative trait loci (QTL) for economically important traits linked to these end-uses will facilitate identification of the main candidate genes and elucidate the regulatory mechanisms.

Method

A high-density consensus map (a total of 2754 SNPs with 1359.18 cM) was constructed using genotyping by sequencing (GBS) on clonal progenies of E. urophylla × tereticornis hybrids. QTL mapping of growth and wood property traits were conducted in three common garden experiments, resulting in a total of 108 QTLs. A total of 1052 candidate genes were screened by the efficient combination of QTL mapping and transcriptome analysis.

Results

Only ten QTLs were found to be stable across two environments, and only one (qSG10Stable mapped on chromosome 10, and associated with lignin syringyl-to-guaiacyl ratio) was stable across all three environments. Compared to other QTLs, qSG10Stable explained a very high level of phenotypic variation (18.4-23.6%), perhaps suggesting that QTLs with strong effects may be more stably inherited across multiple environments. Screened candidate genes were associated with some transcription factor families, such as TALE, which play an important role in the secondary growth of plant cell walls and the regulation of wood formation.

Discussion

While QTLs such as qSG10Stable, found to be stable across three sites, appear to be comparatively uncommon, their identification is likely to be a key to practical QTL-based breeding. Further research involving clonally-replicated populations, deployed across multiple target planting sites, will be required to further elucidate QTL-by-environment interactions.