The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis

Genome-Wide Identification and Expression Profile Analysis of the Phenylalanine Ammonia-Lyase Gene Family in Hevea brasiliensis

The majority of the world's natural rubber comes from the rubber tree (Hevea brasiliensis). As a key enzyme for synthesizing phenylpropanoid compounds, phenylalanine ammonia-lyase (PAL) has a critical role in plant satisfactory growth and environmental adaptation. To clarify the characteristics of rubber tree PAL family genes, a genome-wide characterization of rubber tree PALs was conducted in this study. Eight PAL genes (HbPAL1-HbPAL8), which spread over chromosomes 3, 7, 8, 10, 12, 13, 14, 16, and 18, were found to be present in the genome of H. brasiliensis. Phylogenetic analysis classified HbPALs into groups I and II, and the group I HbPALs (HbPAL1-HbPAL6) displayed similar conserved motif compositions and gene architectures. Tissue expression patterns of HbPALs quantified by quantitative real-time PCR (qPCR) proved that distinct HbPALs exhibited varying tissue expression patterns. The HbPAL promoters contained a plethora of cis-acting elements that responded to hormones and stress, and the qPCR analysis demonstrated that abiotic stressors like cold, drought, salt, and H2O2-induced oxidative stress, as well as hormones like salicylic acid, abscisic acid, ethylene, and methyl jasmonate, controlled the expression of HbPALs. The majority of HbPALs were also regulated by powdery mildew, anthracnose, and Corynespora leaf fall disease infection. In addition, HbPAL1, HbPAL4, and HbPAL7 were significantly up-regulated in the bark of tapping panel dryness rubber trees relative to that of healthy trees. Our results provide a thorough comprehension of the characteristics of HbPAL genes and set the groundwork for further investigation of the biological functions of HbPALs in rubber trees.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.

Heterozygosity analysis of spontaneous 2n female gametes and centromere mapping of the diploid Hevea brasiliensis based on full-sib triploid populations

KEY MESSAGE

Unreduced megagametophytes via second-division restitution were confirmed through heterozygosity analysis, and four candidate physical centromeres of rubber were located for the first time. The evaluation of maternal heterozygosity restitution (MHR) is vital in identifying the mechanism of 2n gametogenesis and assessing the utilization value of 2n gametes. In this study, three full-sib triploid populations were employed to evaluate the MHR of 2n female gametes of rubber tree clone GT1 and to confirm their genetic derivation. The 2n female gametes of GT1 were derived from second-division restitution (SDR) and transmitted more than half of the parental heterozygosity. In addition, low recombination frequency markers were developed, and four candidate physical centromeres of rubber tree were located for the first time. The confirmation that 2n female gametes of rubber tree clone GT1 are derived from SDR provides insights into the molecular mechanisms of 2n gametogenesis. In addition, the identified centromere location will aid in the development of centromeric markers for the rapid identification of the 2n gametogenesis mechanism.

Genome-wide identification of rubber tree pathogenesis-related 10 (PR-10) proteins with biological relevance to plant defense

Pathogenesis-related 10 (PR-10) is a group of small intracellular proteins that is one of 17 subclasses of pathogenesis-related proteins in plants. The PR-10 proteins have been studied extensively and are well-recognized for their contribution to host defense against phytopathogens in several plant species. Interestingly, the accumulation of PR-10 proteins in the rubber tree, one of the most economically important crops worldwide, after being infected by pathogenic organisms has only recently been reported. In this study, the homologous proteins of the PR-10 family were systemically identified from the recently available rubber tree genomes in the NCBI database. The sequence compositions, structural characteristics, protein physical properties, and phylogenetic relationships of identified PR-10 proteins in rubber trees support their classification into subgroups, which mainly consist of Pru ar 1-like major allergens and major latex-like (MLP) proteins. The rubber tree PR10-encoding genes were majorly clustered on chromosome 15. The potential roles of rubber tree PR-10 proteins are discussed based on previous reports. The homologous proteins in the PR-10 family were identified in the recent genomes of rubber trees and were shown to be crucial in host responses to biotic challenges. The genome-wide identification conducted here will accelerate the future study of rubber tree PR-10 proteins. A better understanding of these defense-related proteins may contribute to alternative ways of developing rubber tree clones with desirable traits in the future.

TCOD: an integrated resource for tropical crops

Tropical crops are vital for tropical agriculture, with resource scarcity, functional diversity and extensive market demand, providing considerable economic benefits for the world's tropical agriculture-producing countries. The rapid development of sequencing technology has promoted a milestone in tropical crop research, resulting in the generation of massive amount of data, which urgently needs an effective platform for data integration and sharing. However, the existing databases cannot fully satisfy researchers' requirements due to the relatively limited integration level and untimely update. Here, we present the Tropical Crop Omics Database (TCOD, https://ngdc.cncb.ac.cn/tcod), a comprehensive multi-omics data platform for tropical crops. TCOD integrates diverse omics data from 15 species, encompassing 34 chromosome-level de novo assemblies, 1 255 004 genes with functional annotations, 282 436 992 unique variants from 2048 WGS samples, 88 transcriptomic profiles from 1997 RNA-Seq samples and 13 381 germplasm items. Additionally, TCOD not only employs genes as a bridge to interconnect multi-omics data, enabling cross-species comparisons based on homology relationships, but also offers user-friendly online tools for efficient data mining and visualization. In short, TCOD integrates multi-species, multi-omics data and online tools, which will facilitate the research on genomic selective breeding and trait biology of tropical crops.

Comprehensive Analysis of the Catalase (CAT) Gene Family and Expression Patterns in Rubber Tree (Hevea brasiliensis) under Various Abiotic Stresses and Multiple Hormone Treatments

Catalase (CAT) is one of the key enzymes involved in antioxidant defense systems and mainly scavenges H2O2 and plays a vital role in plant growth, development, and various adverse stresses. To date, a systematic study of the CAT gene family in rubber tree has not been reported. In this study, five HbCAT gene family members were identified from the rubber tree genome, and these were mainly clustered into two subfamilies. Gene structure and motif analysis showed that exon-intron and motif patterns were conserved across different plant species. Sequence analysis revealed that HbCAT proteins contain one active catalytic site, one heme-ligand signature sequence, three conserved amino acid residues (His, Tyr, and Asn), and one peroxisome-targeting signal 1 (PTS1) sequence. Fragment duplication is a selection pressure for the evolution of the HbCAT family based on Ka/Ks values. Analysis of cis-acting elements in the promoters indicated that HbCAT gene expression might be regulated by abscisic acid (ABA), salicylic acid (SA), and MYB transcription factors; furthermore, these genes might be involved in plant growth, development, and abiotic stress responses. A tissue-specific expression analysis showed that HbCATs gradually increased with leaf development and were highly expressed in mature leaves. Gene expression profiling exhibited the differential expression of the HbCATs under cold, heat, drought, and NaCl stresses. Our results provide comprehensive information about the HbCAT gene family, laying the foundation for further research on its function in rubber tree.

Identification and characterization of two closely related virga-like viruses latently infecting rubber trees (Hevea brasiliensis)

A novel virga-like virus, provisionally named Rubber tree latent virus 2 (RTLV2), was identified from rubber tree (Hevea brasiliensis). It is a close relative of the previously reported Rubber tree latent virus 1 (RTLV1). The complete genomes of RTLV1 and RTLV2 were sequenced and comparatively analyzed in terms of genome organization, putative gene products and phylogenetic relationship. Both RTLV1 and RTLV2 have positive-sense single-stranded RNA genomes that encode seven open reading frames (ORFs), forming a similar genomic layout. In phylogenetic analyses based on replicase and coat protein amino acid sequences, RTLV1 and RTLV2 were clustered with unclassified virga-like viruses. They are distinct from currently recognized plant virus families. RTLV1 and RTLV2 can be distinguished from members of Virgaviridae by the presence of a putative coat protein duplex and a poly(A) tail at the 3'-terminus. The authenticity of RTLV1 and RTLV2 as infectious viruses was confirmed through field investigations and transmissibility assays. In conclusion, RTLV1 and RTLV2 represent a novel plant virus group that does not readily fit into current virus families.

Advances in Genome Sequencing and Natural Rubber Biosynthesis in Rubber-Producing Plants

Natural rubber (cis-1,4-polyisoprene, NR) is an important raw material utilized widely in the manufacturing of medical, agricultural, and industrial products. Rubber tree (Hevea brasiliensis) and several alternative rubber-producing plants (Taraxacum kok-saghyz, Lactuca sativa, and Parthenium argentatum) have the capability to produce high-quality NR. With the progress of genome sequencing, similar rubber biosynthesis pathways have been discovered among different rubber-producing plant species. NR is synthesized and stored in rubber particles, which are specialized organelles comprising a hydrophobic NR core surrounded by a lipid monolayer and membrane-bound proteins. The rubber transferase complex is considered to be the pivotal enzyme involved in catalyzing NR biosynthesis. However, the exact compositions of the RT complex in rubber-producing plants remain elusive and poorly understood. Here, we review the progress of genome sequencing, natural rubber biosynthesis, and the components of the RT complex in rubber-producing plants. We emphasize that identifying the detailed components of the RT complex holds great significance for exploring the mechanism of NR biosynthesis and accelerating molecular breeding in rubber-producing plants.

Telomere-to-telomere assembly of cassava genome reveals the evolution of cassava and divergence of allelic expression

Cassava is a crucial crop that makes a significant contribution to ensuring human food security. However, high-quality telomere-to-telomere cassava genomes have not been available up to now, which has restricted the progress of haploid molecular breeding for cassava. In this study, we constructed two nearly complete haploid resolved genomes and an integrated, telomere-to-telomere gap-free reference genome of an excellent cassava variety, 'Xinxuan 048', thereby providing a new high-quality genomic resource. Furthermore, the evolutionary history of several species within the Euphorbiaceae family was revealed. Through comparative analysis of haploid genomes, it was found that two haploid genomes had extensive differences in linear structure, transcriptome features, and epigenetic characteristics. Genes located within the highly divergent regions and differentially expressed alleles are enriched in the functions of auxin response and the starch synthesis pathway. The high heterozygosity of cassava 'Xinxuan 048' leads to rapid trait segregation in the first selfed generation. This study provides a theoretical basis and genomic resource for molecular breeding of cassava haploids.

Plant exosome nanovesicles (PENs): green delivery platforms

Natural plants have been attracting increasing attention in biomedical research due to their numerous benefits. Plant exosome-derived vesicles, some of the plant's components, are small nanoscale vesicles secreted by plant cells. These vesicles are rich in bioactive substances and play significant roles in intercellular communication, information transfer, and maintaining homeostasis in organisms. They also hold promise for treating diseases, and their vesicular structures make them suitable carriers for drug delivery, with large-scale production feasible. Therefore, this paper aims to provide an overview of nanovesicles from different plant sources and their extraction methods. We also outline the biological activities of nanovesicles, including their anti-inflammatory, anti-viral, and anti-tumor properties, and systematically introduce their applications in drug delivery. These applications include transdermal delivery, targeted drug delivery, gene delivery, and their potential use in the modern food industry. This review provides new ideas and methods for future research on plant exosomes, including their empowerment by artificial intelligence and gene editing, as well as their potential application in the biomedicine, food, and agriculture industries.

Phylogenetic Analysis of the PR-4 Gene Family in Euphorbiaceae and Its Expression Profiles in Tung Tree (Vernicia fordii)

Pathogenesis-related protein-4 (PR-4) is generally believed to be involved in physiological processes. However, a comprehensive investigation of this protein in tung tree (Vernicia fordii) has yet to be conducted. In this study, we identified 30 PR-4 genes in the genomes of Euphorbiaceae species and investigated their domain organization, evolution, promoter cis-elements, expression profiles, and expression profiles in the tung tree. Sequence and structural analyses indicated that VF16136 and VF16135 in the tung tree could be classified as belonging to Class II and I, respectively. Phylogenetic and Ka/Ks analyses revealed that Hevea brasiliensis exhibited a significantly expanded number of PR-4 genes. Additionally, the analysis of promoter cis-elements suggested that two VfPR-4 genes may play a role in the response to hormones and biotic and abiotic stress of tung trees. Furthermore, the expression patterns of VfPR-4 genes and their responses to 6-BA, salicylic acid, and silver nitrate in inflorescence buds of tung trees were evaluated using qRT-PCR. Notably, the expression of two VfPR-4 genes was found to be particularly high in leaves and early stages of tung seeds. These results suggest that VF16136 and VF16135 may have significant roles in the development of leaves and seeds in tung trees. Furthermore, these genes were found to be responsive to 6-BA, salicylic acid, and silver nitrate in the development of inflorescence buds. This research provides valuable insights for future investigation into the functions of PR-4 genes in tung trees.

Multiomics studies with co-transformation reveal microRNAs via miRNA-TF-mRNA network participating in wood formation in Hevea brasiliensis

Introduction

MicroRNAs (miRNAs) are small endogenous non-coding RNAs that play an important role in wood formation in plants. However, the significance of the link between miRNAs and their target transcripts in wood formation remains unclear in rubber tree (Hevea brasiliensis).

Methods

In this study, we induced the formation of reaction wood by artificially bending rubber trees for 300 days and performed small RNA sequencing and transcriptome deep sequencing (RNA-seq) to describe the complement of miRNAs and their targets contributing to this process.

Results and discussion

We identified 5, 11, and 2 differentially abundant miRNAs in normal wood (NW) compared to tension wood (TW), in NW relative to opposite wood (OW), and between TW and OW, respectively. We also identified 12 novel miRNAs and 39 potential miRNA-mRNA pairs with different accumulation patterns in NW, TW, and OW. We noticed that many miRNAs targeted transcription factor genes, which were enriched in KEGG pathways associated with phenylpropanoid biosynthesis, phenylalanine metabolism, and pyruvate metabolism. Thus, miRNA-TF-mRNA network involved in wood formation via tension wood model were constructed. We validated the differential accumulation of miRNAs and their targets by RT-qPCR analysis and overexpressed miRNA in Nicotiana benthamiana with its potential target gene. These results will provide a reference for a deep exploration of growth and development in rubber tree.

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker-trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.

Chromosomal-Scale Genome Assemblies of Two Coastal Plant Species, Scaevola taccada and S. hainanensis-Insight into Adaptation Outside of the Common Range

While most of the species in Goodeniaceae family, excluding the Scaevola genus, are endemic to Australasia, S. taccada and S. hainanensis have expanded their distribution range to the tropical coastlines of the Atlantic and Indian Oceans. S. taccada appears to be highly adapted to coastal sandy lands and cliffs, and it has become invasive in places. S. hainanensis is found mainly in salt marshes near mangrove forests, and is at risk of extinction. These two species provide a good system to investigate adaptive evolution outside the common distribution range of this taxonomic group. Here, we report their chromosomal-scale genome assemblies with the objective of probing their genomic mechanisms related to divergent adaptation after leaving Australasia. The scaffolds were assembled into eight chromosome-scale pseudomolecules, which covered 90.12% and 89.46% of the whole genome assembly for S. taccada and S. hainanensis, respectively. Interestingly, unlike many mangroves, neither species has undergone whole-genome duplication. We show that private genes, specifically copy-number expanded genes are essential for stress response, photosynthesis, and carbon fixation. The gene families that are expanded in S. hainanensis and contracted in S. taccada might have facilitated adaptation to high salinity in S. hainanensis. Moreover, the genes under positive selection in S. hainanensis have contributed to its response to stress and its tolerance of flooding and anoxic environments. In contrast, compared with S. hainanensis, the more drastic copy number expansion of FAR1 genes in S. taccada might have facilitated its adaptation to the stronger light radiation present in sandy coastal lands. In conclusion, our study of the chromosomal-scale genomes of S. taccada and S. hainanensis provides novel insights into their genomic evolution after leaving Australasia.

Functional characterization of a bark-specific monoterpene synthase potentially involved in wounding- and methyl jasmonate-induced linalool emission in rubber (Hevea brasiliensis)

Rubber (Hevea brasiliensis) is a latex-producing plant that often encounters mechanical wounding, as well as pathogen and pest attacks through wound sites during and after tapping. Terpenoids play an important role in the ecological interactions of many plant species, and their diversity is mainly generated by enzymes known as terpene synthases (TPS). In this study, one cDNA sequence encoding a putative terpene synthase, HbTPS20, was obtained from the bark tissues of H. brasiliensis. The encoded protein contains 610 amino acids with a putative N-terminal plastid transit peptide of approximately 70 residues. It belongs to the TPS-b subfamily. Further phylogenetic analysis showed that HbTPS20 formed a separate branch that diverged from the progenitor of all other potentially functional terpene synthases of the rubber TPS-b subfamily. The truncated HbTPS20 without the signal peptide coding sequence was successfully expressed in E. coli and in vitro enzymatic assays with geranyl diphosphate (GPP) or neryl diphosphate (NPP) as a substrate defined HbTPS20 as an active linalool synthase (HbLIS) with the ability to produce linalool as the principal product. RT-qPCR analysis showed that the highest transcript levels of HbTPS20 were found in barks, and this gene was expressed at 2.26- and 250-fold greater levels in the bark tissues of wounded and MeJA-treated plants, respectively, than in those of the control plants. This indicates that this gene may be involved in the induced stress responses of rubber.

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Analysis of Genetic Diversity and Resistance to Foliar Pathogens Based on Genotyping-by-Sequencing of a Para Rubber Diversity Panel and Progeny of an Interspecific Cross

Para rubber trees (Hevea brasiliensis) are the largest major source of natural rubber in the world. Its major pathogens are Phytophthora spp., Corynespora cassiicola, and Colletotrichum spp. A rubber diversity panel of 116 clones using over 12,000 single nucleotide polymorphisms (SNPs) from DArTSeq genotyping revealed clear phylogenetic differences in clones that originated from different geographical regions of the world. An integrated linkage map constructed with an F₁ progeny of 86 from an interspecific cross between H. brasiliensis and H. benthamiana using 23,978 markers [10,323 SNPs and 13,655 SilicoDArTs] spanned 3947.83 cM with 0.83 cM average marker-interval. The genome scaffolds that were anchored to the linkage map, covering 1.44 Gb of H. brasiliensis reference genome, revealed a high level of collinearity between the genetic map and reference genome. Association analysis identified 12 SNPs significantly associated with the resistance against Phytophthora, Corynespora, and Colletotrichum in six linkage groups: 2, 6, 12, 14, 17, and 18. Kompetitive Allele-Specific PCR marker assays were developed for those 12 SNPs, screened with 178 individuals, and detected clear separation between two genotypes. Within the proximity to those SNPs, 41 potentially key genes that have previously been reported to associate with plant disease resistance were predicted with high confidence.

Improving Plant Photosynthesis through Light-Harvesting Upconversion Nanoparticles

Nanotechnology is considered as an emerging effective means to augment plant photosynthesis. However, there is still a lot of work to be done in this field. Here, we applied the upconversion nanoparticles (UCNPs) on lettuce leaves and found that the UCNPs were able to transport into the lettuce body and colocalize with the chloroplasts. It was proved that UCNPs could harvest the near-infrared light of sunlight and increase the electron transfer rate in the photosynthesis process, thus increasing the photosynthesis rate. The gene expression analysis showed that more than 90% of gene expression in photosynthesis was upregulated. After spraying the UCNP solution on the leaves of lettuce and placing the lettuce under sunlight for 1 week, the wet/dry weight of the leaves increased by 53.33% and 45.71%, respectively. This nanoengineering of light-harvesting UCNPs may have great potential for applications in agriculture.

The first high-quality chromosome-level genome assembly of Phyllanthaceae (Phyllanthus cochinchinensis) provides insights into flavonoid biosynthesis

We report the genome assembly of P. cochinchinensis, as the first high-quality chromosome-level genome of Phyllanthaceae which is rich in medicinal plants. Phyllanthus cochinchinensis, a member of the Phyllanthaceae, is one of the famous medicinal plants in South China. Here, we report a de novo chromosome-level genome assembly for P. cochinchinensis using a combination of Nanopore and Illumina sequencing technologies. In total, the assembled genome consists of 284.88 Mb genomic sequences with a contig N50 of 10.32 Mb, representing ~ 95.49% of the estimated genome size. By applying Hi-C data, 13 pseudochromosomes of P. cochinchinensis were constructed, covering ~ 99.87% of the assembled sequences. The genome is annotated with 59.12% repetitive sequences and 20,836 protein-coding genes. Whole-genome duplication of P. cochinchinensis is likely shared with Ricinus communis as well as Vitis vinifera. Homologous genes within the flavonoid pathway for P. cochinchinensis were identified and copy numbers and expression level of related genes revealed potential critical genes involved in flavonoid biosynthesis. This study provides the first whole-genome sequence for the Phyllanthaceae, confirms the evolutionary status of Phyllanthus from the genomic level, and provides foundations for accelerating functional genomic research of species from Phyllanthus.

A divide-and-conquer approach for genomic prediction in rubber tree using machine learning

Rubber tree (Hevea brasiliensis) is the main feedstock for commercial rubber; however, its long vegetative cycle has hindered the development of more productive varieties via breeding programs. With the availability of H. brasiliensis genomic data, several linkage maps with associated quantitative trait loci have been constructed and suggested as a tool for marker-assisted selection. Nonetheless, novel genomic strategies are still needed, and genomic selection (GS) may facilitate rubber tree breeding programs aimed at reducing the required cycles for performance assessment. Even though such a methodology has already been shown to be a promising tool for rubber tree breeding, increased model predictive capabilities and practical application are still needed. Here, we developed a novel machine learning-based approach for predicting rubber tree stem circumference based on molecular markers. Through a divide-and-conquer strategy, we propose a neural network prediction system with two stages: (1) subpopulation prediction and (2) phenotype estimation. This approach yielded higher accuracies than traditional statistical models in a single-environment scenario. By delivering large accuracy improvements, our methodology represents a powerful tool for use in Hevea GS strategies. Therefore, the incorporation of machine learning techniques into rubber tree GS represents an opportunity to build more robust models and optimize Hevea breeding programs.

Divergence of flowering-related genes to control flowering in five Euphorbiaceae genomes

Reproductive growth and vegetative growth are a pair of main contradictions in the process of plant growth. Flowering, as part of reproductive growth, is a key switch in the life cycle of higher plants, which affects the yield and economic benefits of plants to a certain extent. The Euphorbiaceae species, including castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), cassava (Manihot esculenta), and rubber tree (Hevea brasiliensis), have important economic values because they are raw materials for the production of biodiesel, rubber, etc. The flowering mechanisms are still excluded in the Euphorbiaceae species. The flowering-related genes of Arabidopsis thaliana (Arabidopsis) were used as a reference to determine the orthologs of these genes in Euphorbiaceae genomes. The result showed that 146, 144, 114, 114, and 149 of 207 A. thaliana genes were respectively matched to R. communis, V. fordii, J. curcas, H. brasiliensis, and M. esculenta. These identified genes were clustered into seven pathways including gibberellins, floral meristem identity (FMI), vernalization, photoperiod, floral pathway integrators (FPIs), and autonomous pathways. Then, some key numbers of flowering-related genes are widely conserved in the Euphorbiaceae genomes including but not limited to FPI genes LFY, SOC1, FT, and FMI genes AG, CAL, and FUL. However, some genes, including FRI, FLC, and GO, were missing in several or all five Euphorbiaceae species. In this study, we proposed the putative mechanisms of flowering-related genes to control flowering and provided new candidate flowering genes for using marker-assisted breeding to improve variety quality.

Insights to the superoxide dismutase genes and its roles in Hevea brasiliensis under abiotic stress

The superoxide dismutase (SOD) protein significantly influences the development and growth of plants and their reaction to abiotic stresses. However, little is known about the characteristics of rubber tree SOD genes and their expression changes under abiotic stresses. The present study recognized 11 SOD genes in the rubber tree genome, including 7 Cu/ZnSODs, 2 MnSODs, and 2 FeSODs. Except for HbFSD1, SODs were scattered on five chromosomes. The phylogenetic analysis of SOD proteins in rubber trees and a few other plants demonstrated that the SOD proteins contained three major subgroups. Moreover, the genes belonging to the same clade contained similar gene structures, which confirmed their classification further. The extension of the SOD gene family in the rubber tree was mainly induced by the segmental duplication events. The cis-acting components analysis showed that HbSODs were utilized in many biological procedures. The transcriptomics data indicated that the phosphorylation of the C-terminal domain of RNA polymerase II might control the cold response genes through the CBF pathway and activate the SOD system to respond to cold stress. The qRT-PCR results showed that the expression of HbCSD1 was significantly downregulated under drought and salt stresses, which might dominate the adaption capability to different stresses. Additionally, salt promoted the expression levels of HbMSD1 and HbMSD2, exhibiting their indispensable role in the salinity reaction. The study results will provide a theoretical basis for deep research on HbSODs in rubber trees.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03328-7.

Characterization of the complete chloroplast genome of Hevea pauciflora (Euphorbiaceae), an important wild relative of the rubber tree

Hevea pauciflora belongs to the Euphorbiaceae family, an important wild relative of the rubber tree. This study sequenced, assembled, and annotated the complete chloroplast genome of H. pauciflora. The complete chloroplast genome is 161,123 bp with a canonical quadripartite structure containing a large single-copy (LSC) region (89,109 bp), a small single-copy (SSC) region (18,376 bp), and two inverted repeat regions (IRa and IRb) (26,819 bp, each). A total of 134 genes were annotated, including 86 protein-coding genes, four pseudogenes, 36 tRNA genes, and eight rRNA genes. The 134 genes include four major groups: ‘self-replication’, ‘photosynthesis’, ‘unknown function’, and ‘others’. A phylogenetic analysis clustered H. pauciflora, H. brasiliensis, H. camargoana, and H. benthamiana into one clade, consistent with traditional taxonomy. This study provides useful data for further studies of Hevea genus and the phylogenetic relationships of Euphorbiaceae species.

Evolution and function of ubiquitin-specific proteases (UBPs): Insight into seed development roles in tung tree (Vernicia fordii)

The tung oil produced by the tung tree (Vernicia fordii) provides resources for the manufacture of biodiesel. Ubiquitin-specific proteases (UBPs) are the largest group of deubiquitinases and play key roles in regulating development and stress responses. Here, 21 UBPs were identified in V. fordii, roughly one-half the number found in Manihot esculenta and Hevea brasiliensis. Most UBP duplications are produced from whole-genome duplication (WGD), and significant differences in gene retention existed among Euphorbiaceae. The great majority of UBP-containing blocks in V. fordii, V. montana, Ricinus communis, and Jatropha curcas exhibited extensive conservation with the duplicated regions of M. esculenta and H. brasiliensis. These blocks formed 14 orthologous groups, indicating they shared WGD with UBPs in M. esculenta and H. brasiliensis, but most of these UBPs copies were lost. The UBP orthologs contained significant functional divergence which explained the susceptibility of V. fordii to Fusarium wilt and the resistance of V. montana to Fusarium wilt. The expression patterns and experiments suggested that Vf03G1417 could affect the seed-related traits and positively regulate the seed oil accumulation. This study provided important insights into the evolution of UBPs in Euphorbiaceae and identified important candidate VfUBPs for marker-assisted breeding in V. fordii.

The chromosome-scale genome sequence of Triadica sebifera provides insight into fatty acids and anthocyanin biosynthesis

The Chinese tallow tree (Triadica sebifera) can produce oil with high content of unsaturated fatty acids in seeds and shows attractive leaf color in autumn and winter. Here, the 739 Mb chromosome-scale genome sequence of the Chinese tallow tree was assembled and it reveals the Chinese tallow tree is a tetraploid. Numerous genes related to nutrition assimilation, energy utilization, biosynthesis of secondary metabolites and resistance significantly expanded or are specific to the Chinese tallow tree. These genes would enable the Chinese tallow tree to obtain high adaptability. More genes in fatty acids biosynthesis in its genome, especially for unsaturated fatty acids biosynthesis, and higher expression of these genes in seeds would be attributed to its high content of unsaturated fatty acids. Cyanidin 3-O-glucoside was identified as the major component of anthocyanin in red leaves. All structural genes in anthocyanin biosynthesis show significantly higher expression in red leaves than in green leaves. Transcription factors, seven MYB and one bHLH, were predicted to regulate these anthocyanin biosynthesis genes. Collectively, we provided insight into the polyploidization, high adaptability and biosynthesis of the high content of unsaturated fatty acids in seeds and anthocyanin in leaves for the Chinese tallow tree.

A chromosome-level assembly of the economically-important Chinese tallow tree is presented alongside functional analyses into biosynthetic products in the leaves and seeds.

Comparative genomics reveals a dynamic genome evolution in the ectomycorrhizal milk-cap (Lactarius) mushrooms

Ectomycorrhizal fungi play a key role in forests by establishing mutualistic symbioses with woody plants. Genome analyses have identified conserved symbiosis-related traits among ectomycorrhizal fungal species, but the molecular mechanisms underlying host specificity remain poorly known. We sequenced and compared the genomes of seven species of milk-cap fungi (Lactarius, Russulales) with contrasting host specificity. We also compared these genomes with those of symbiotic and saprotrophic Russulales species, aiming to identify genes involved in their ecology and host specificity. The size of Lactarius genomes is significantly larger than other Russulales species, owing to a massive accumulation of transposable elements and duplication of dispensable genes. As expected, their repertoire of genes coding for plant cell wall-degrading enzymes is restricted, but they retained a substantial set of genes involved in microbial cell wall degradation. Notably, Lactarius species showed a striking expansion of genes encoding proteases, such as secreted ectomycorrhiza-induced sedolisins. A high copy number of genes coding for small secreted LysM proteins and Lactarius-specific lectins were detected, which may be linked to host specificity. This study revealed a large diversity in the genome landscapes and gene repertoires within Russulaceae. The known host specificity of Lactarius symbionts may be related to mycorrhiza-induced species-specific genes, including secreted sedolisins.

The Chromosome-Scale Reference Genome of Macadamia tetraphylla Provides Insights Into Fatty Acid Biosynthesis

Macadamia is an evergreen tree belonging to the Proteaceae family. The two commercial macadamia species, Macadamia integrifolia and M. tetraphylla, are highly prized for their edible kernels. The M. integrifolia genome was recently sequenced, but the genome of M. tetraphylla has to date not been published, which limits the study of biological research and breeding in this species. This study reports a high-quality genome sequence of M. tetraphylla based on the Oxford Nanopore Technologies technology and high-throughput chromosome conformation capture techniques (Hi-C). An assembly of 750.87 Mb with 51.11 Mb N50 length was generated, close to the 740 and 758 Mb size estimates by flow cytometry and k-mer analysis, respectively. Genome annotation indicated that 61.42% of the genome is composed of repetitive sequences and 34.95% is composed of long terminal repeat retrotransposons. Up to 31,571 protein-coding genes were predicted, of which 92.59% were functionally annotated. The average gene length was 6,055 bp. Comparative genome analysis revealed that the gene families associated with defense response, lipid transport, steroid biosynthesis, triglyceride lipase activity, and fatty acid metabolism are expanded in the M. tetraphylla genome. The distribution of fourfold synonymous third-codon transversion showed a recent whole-genome duplication event in M. tetraphylla. Genomic and transcriptomic analysis identified 187 genes encoding 33 crucial oil biosynthesis enzymes, depicting a comprehensive map of macadamia lipid biosynthesis. Besides, the 55 identified WRKY genes exhibited preferential expression in root as compared to that in other tissues. The genome sequence of M. tetraphylla provides novel insights for breeding novel varieties and genetic improvement of agronomic traits.

Genome-Wide Identification, Evolution and Expression Analysis of the Glutathione S-Transferase Supergene Family in Euphorbiaceae

Euphorbiaceae, a family of plants mainly grown in the tropics and subtropics, is also widely distributed all over the world and is well known for being rich in rubber, oil, medicinal materials, starch, wood and other economically important plant products. Glutathione S-transferases (GSTs) constitute a family of proteins encoded by a large supergene family and are widely expressed in animals, bacteria, fungi and plants, but with few reports of them in Euphorbiaceae plants. These proteins participate in and regulate the detoxification and oxidative stress response of heterogeneous organisms, resistance to stress, growth and development, signal transduction and other related processes. In this study, we identified and analyzed the whole genomes of four species of Euphorbiaceae, namely Ricinus communis, Jatropha curcas, Hevea brasiliensis, and Manihot esculenta, which have high economic and practical value. A total of 244 GST genes were identified. Based on their sequence characteristics and conserved domain types, the GST supergene family in Euphorbiaceae was classified into 10 subfamilies. The GST supergene families of Euphorbiaceae and Arabidopsis have been found to be highly conserved in evolution, and tandem repeats and translocations in these genes have made the greatest contributions to gene amplification here and have experienced strong purification selection. An evolutionary analysis showed that Euphorbiaceae GST genes have also evolved into new subtribes (GSTO, EF1BG, MAPEG), which may play a specific role in Euphorbiaceae. An analysis of expression patterns of the GST supergene family in Euphorbiaceae revealed the functions of these GSTs in different tissues, including resistance to stress and participation in herbicide detoxification. In addition, an interaction analysis was performed to determine the GST gene regulatory mechanism. The results of this study have laid a foundation for further analysis of the functions of the GST supergene family in Euphorbiaceae, especially in stress and herbicide detoxification. The results have also provided new ideas for the study of the regulatory mechanism of the GST supergene family, and have provided a reference for follow-up genetics and breeding work.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

Large structural variations in the haplotype-resolved African cassava genome

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

HbSnRK2.6 Functions in ABA-Regulated Cold Stress Response by Promoting HbICE2 Transcriptional Activity in Hevea brasiliensis

Low temperature remarkably limits rubber tree (Hevea brasiliensis Muell. Arg.) growth, latex production, and geographical distribution, but the underlying mechanisms of Hevea brasiliensis cold stress response remain elusive. Here, we identified HbSnRK2.6 as a key component in ABA signaling functions in phytohormone abscisic acid (ABA)-regulated cold stress response in Hevea brasiliensis. Exogenous application of ABA enhances Hevea brasiliensis cold tolerance. Cold-regulated (COR) genes in the CBF pathway are upregulated by ABA. Transcript levels of all five HbSnRK2.6 members are significantly induced by cold, while HbSnRK2.6A, HbSnRK2.6B, and HbSnRK2.6C can be further activated by ABA under cold conditions. Additionally, HbSnRK2.6s are localized in the cytoplasm and nucleus, and can physically interact with HbICE2, a crucial positive regulator in the cold signaling pathway. Overexpression of HbSnRK2.6A or HbSnRK2.6B in Arabidopsis extensively enhances plant responses to ABA and expression of COR genes, leading to increased cold stress tolerance. Furthermore, HbSnRK2.6A and HbSnRK2.6B can promote transcriptional activity of HbICE2, thus, increasing the expression of HbCBF1. Taken together, we demonstrate that HbSnRK2.6s are involved in ABA-regulated cold stress response in Hevea brasiliensis by regulating transcriptional activity of HbICE2.

RNA Sequencing Reveals Phenylpropanoid Biosynthesis Genes and Transcription Factors for Hevea brasiliensis Reaction Wood Formation

Given the importance of wood in many industrial applications, much research has focused on wood formation, especially lignin biosynthesis. However, the mechanisms governing the regulation of lignin biosynthesis in the rubber tree (Hevea brasiliensis) remain to be elucidated. Here, we gained insight into the mechanisms of rubber tree lignin biosynthesis using reaction wood (wood with abnormal tissue structure induced by gravity or artificial mechanical treatment) as an experimental model. We performed transcriptome analysis of rubber tree mature xylem from tension wood (TW), opposite wood (OW), and normal wood (NW) using RNA sequencing (RNA-seq). A total of 214, 1,280, and 32 differentially expressed genes (DEGs) were identified in TW vs. NW, OW vs. NW, and TW vs. OW, respectively. GO and KEGG enrichment analysis of DEGs from different comparison groups showed that zeatin biosynthesis, plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interaction pathways may play important roles in reaction wood formation. Sixteen transcripts involved in phenylpropanoid biosynthesis and 129 transcripts encoding transcription factors (TFs) were used to construct a TF-gene regulatory network for rubber tree lignin biosynthesis. Among them, MYB, C2H2, and NAC TFs could regulate all the DEGs involved in phenylpropanoid biosynthesis. Overall, this study identified candidate genes and TFs likely involved in phenylpropanoid biosynthesis and provides novel insights into the mechanisms regulating rubber tree lignin biosynthesis.

Genome-Wide Identification and Analysis of Small Nucleolar RNAs and Their Roles in Regulating Latex Regeneration in the Rubber Tree (Hevea brasiliensis)

Small nucleolar RNAs (snoRNAs) are a class of conserved nuclear RNAs that play important roles in the modification of ribosomal RNAs (rRNAs) in plants. In rubber trees, rRNAs are run off with latex flow during tapping and need to be regenerated for maintaining the functions of the laticifer cells. SnoRNAs are expected to play essential roles in the regeneration of rRNAs. However, snoRNAs in the rubber tree have not been sufficiently characterized thus far. In this study, we performed nuclear RNA sequencing (RNA-seq) to identify snoRNAs globally and investigate their roles in latex regeneration. We identified a total of 3,626 snoRNAs by computational prediction with nuclear RNA-seq data. Among these snoRNAs, 50 were highly expressed in latex; furthermore, the results of reverse transcription polymerase chain reaction (RT-PCR) showed the abundant expression of 31 of these snoRNAs in latex. The correlation between snoRNA expression and adjusted total solid content (TSC/C) identified 13 positively yield-correlated snoRNAs. To improve the understanding of latex regeneration in rubber trees, we developed a novel insulated tapping system (ITS), which only measures the latex regenerated in specific laticifers. Using this system, a laticifer-abundant snoRNA, HbsnoR28, was found to be highly correlated with latex regeneration. To the best of our knowledge, this is the first report to globally identify snoRNAs that might be involved in latex regeneration regulation and provide new clues for unraveling the mechanisms underlying the regulation of latex regeneration.

High-quality genome assembly of an important biodiesel plant, Euphorbia lathyris L

Caper spurge, Euphorbia lathyris L., is an important energy crop and medicinal crop. Here, we generated a high-quality, chromosome-level genome assembly of caper spurge using Oxford Nanopore sequencing, Illumina sequencing, and Hi-C technology. The final genome assembly was ∼988.9 Mb in size, 99.8% of which could be grouped into 10 pseudochromosomes, with contig and scaffold N50 values of 32.6 and 95.7 Mb, respectively. A total of 651.4 Mb repetitive sequences and 36,342 protein-coding genes were predicted in the genome assembly. Comparative genomic analysis showed that caper spurge and castor bean clustered together. We found that no independent whole-genome duplication event had occurred in caper spurge after its split from castor bean, and recent substantial amplification of LTR retrotransposons (LTR-RTs) has contributed significantly to its genome expansion. Furthermore, based on gene homology searching, we identified a number of candidate genes involved in the biosynthesis of fatty acids and triacylglycerols. The reference genome presented here will be highly useful for the further study of the genetics, genomics, and breeding of this high-value crop, as well as for evolutionary studies of spurge family and angiosperms.

Gapless indica rice genome reveals synergistic contributions of active transposable elements and segmental duplications to rice genome evolution

The ultimate goal of genome assembly is a high-accuracy gapless genome. Here, we report a new assembly pipeline that is used to produce a gapless genome for the indica rice cultivar Minghui 63. The resulting 397.71-Mb final assembly is composed of 12 contigs with a contig N50 size of 31.93 Mb. Each chromosome is represented by a single contig and the genomic sequences of all chromosomes are gapless. Quality evaluation of this gapless genome assembly showed that gene regions in our assembly have the highest completeness compared with the other 15 reported high-quality rice genomes. Further comparison with the japonica rice genome revealed that the gapless indica genome assembly contains more transposable elements (TEs) and segmental duplications (SDs), the latter of which produce many duplicated genes that can affect agronomic traits through dose effect or sub-/neo-functionalization. The insertion of TEs can also affect the expression of duplicated genes, which may drive the evolution of these genes. Furthermore, we found the expansion of nucleotide-binding site with leucine-rich repeat disease-resistance genes and cis-zeatin-O-glucosyltransferase growth-related genes in SDs in the gapless indica genome assembly, suggesting that SDs contribute to the adaptive evolution of rice disease resistance and developmental processes. Collectively, our findings suggest that active TEs and SDs synergistically contribute to rice genome evolution.

Deciphering the roles of leucine-rich repeat receptor-like protein kinases (LRR-RLKs) in response to Fusarium wilt in the Vernicia fordii (Tung tree)

Leucine-rich repeat receptor-like protein kinases (LRR-RLKs) are vital for plant growth and development, signal transduction, immunity, and play diverse roles in plant defense responses. However, the LRR-RLK genes have not been systematically studied in Vernicia fordii (tung tree), especially its response to Fusarium wilt. Here, we carried out an integrative analysis of LRR-RLKs among five Euphorbiaceae species: Hevea brasiliensis (rubber tree), Manihot esculenta (cassava), Jatropha curcas (physic nut), Ricinus communis (castor bean), and V. fordii, which contained 223, 311, 186, 138, and 167 LRR-RLKs, respectively. Maximum-likelihood tree was estimated using LRR-RLKs of Arabidopsis thaliana as a template, and they allowed us to divide Euphorbiaceae LRR-RLKs into 22 groups. There are 126 segmental and 30 tandem duplications in these Euphorbiaceae genomes by synteny analysis. The tissue-specific expression patterns revealed that V. fordii LRR-RLKs (VfLRR-RLKs) were differentially expressed in various tissues, and some of them exhibited specific expression in meristems tissues, which suggested their potential functions during organ formation and cell fate specification. Two VfLRR-RLK pairs (Vf01G2125 and Vf03G1740, Vf06G2687 and Vf10G1659), which generated by tandem duplication events, were associated with possible resistance to Fusarium wilt infection. The qRT-PCR confirmed these four VfLRR-RLKs contained opposite expression profiles during pathogen infection in V. fordii and V. montana. Taken together, our data systematically analyzed the LRR-RLK family in Euphorbiaceae genomes for the first time. We highlight the putative roles of VfLRR-RLKs in response to Fusarium wilt infection, and VfLRR-RLKs may be further applied in marker-assisted breeding to control Fusarium wilt in V. fordii.

Tobacco rattle virus-induced gene silencing in Hevea brasiliensis

Virus-induced gene silencing (VIGS) is a powerful gene-silencing tool that has been intensively applied in plants. To data, the application of VIGS in rubber tree has not yet been reported. In this study, we described the efficient gene silencing in rubber tree by VIGS. The gene encoding Hevea brasiliensis phytoene desaturase (HbPDS) was identified in rubber tree genome. Small interfering RNAs from HbPDS and the silencing gene fragment were predicted and a length of 399 bp was selected to be tested. We showed that the tobacco rattle virus (TRV)-VIGS could induce effective HbPDS silencing in rubber tree. This study was the first to report VIGS in rubber tree. The present TRV-VIGS method could be used to perform reverse genetic approaches to identify unknown gene functions and might be further applied to produce gene silenced rubber tree plants, to advance functional gene of rubber tree.

Optimizing imputation of marker data from genotyping-by-sequencing (GBS) for genomic selection in non-model species: Rubber tree (Hevea brasiliensis) as a case study

Genotyping-by-sequencing (GBS) provides the marker density required for genomic predictions (GP). However, GBS gives a high proportion of missing SNP data which, for species without a chromosome-level genome assembly, must be imputed without knowing the SNP physical positions. Here, we compared GP accuracy with seven map-independent and two map-dependent imputation approaches, and when using all SNPs against the subset of genetically mapped SNPs. We used two rubber tree (Hevea brasiliensis) datasets with three traits. The results showed that the best imputation approaches were LinkImputeR, Beagle and FImpute. Using the genetically mapped SNPs increased GP accuracy by 4.3%. Using LinkImputeR on all the markers allowed avoiding genetic mapping, with a slight decrease in GP accuracy. LinkImputeR gave the highest level of correctly imputed genotypes and its performances were further improved by its ability to define a subset of SNPs imputed optimally. These results will contribute to the efficient implementation of genomic selection with GBS. For Hevea, GBS is promising for rubber yield improvement, with GP accuracies reaching 0.52.

Comparison of Morphological Characteristics and Determination of Different Patterns for Rubber Particles in Dandelion and Different Rubber Grass Varieties

Russian dandelion Taraxacum kok-saghyz (TKS) is one promising alternative crop for natural rubber production. However, it is easily confused with other dandelions. In this study, we performed a systematical comparison of the morphological characteristics for different TKS varieties and common dandelion Taraxacum officinale (TO). Our results demonstrated that several obvious differences in morphology can be found between TKS and TO. TO leaf is a pinnate shape, its margin is heavily jagged and its base is cuneate, but TKS leaf is more cuneate and its leaf margin is nearly smooth and round. There are obvious differences for the outer bracts of TO and TKS flower buds. TKS bracts are oblanceolate, apex obtuse, margin smooth and sinuate, and its outer layer of flower buds and faceplate involucre sepal is buckled inward to form a certain angle. TKS is self-incompatible, and its seeds are spindle-shaped achene and show upright plumpness. A large amount of laticifer cells and rubber particles can be detected from many TKS tissues, and dry roots of TKS contain high contents of natural rubber. Laticifer cells and rubber particles can only be examined in the vein, stem, and roots of TKS. Our statical results also revealed that the numbers of laticifer cells and rubber particles have a positive relationship with the rubber content in TKS roots. These morphological features can help us to easily distinguish TKS from common dandelion and approximately estimate the rubber content in the roots of different TKS varieties for TKS breeding in future.

High-quality de novo assembly of the Eucommia ulmoides haploid genome provides new insights into evolution and rubber biosynthesis

We report the acquisition of a high-quality haploid chromosome-scale genome assembly for the first time in a tree species, Eucommia ulmoides, which is known for its rubber biosynthesis and medicinal applications. The assembly was obtained by applying PacBio and Hi-C technologies to a haploid that we specifically generated. Compared to the initial genome release, this one has significantly improved assembly quality. The scaffold N50 (53.15 MB) increased 28-fold, and the repetitive sequence content (520 Mb) increased by 158.24 Mb, whereas the number of gaps decreased from 104,772 to 128. A total of 92.87% of the 26,001 predicted protein-coding genes identified with multiple strategies were anchored to the 17 chromosomes. A new whole-genome duplication event was superimposed on the earlier γ paleohexaploidization event, and the expansion of long terminal repeats contributed greatly to the evolution of the genome. The more primitive rubber biosynthesis of this species, as opposed to that in Hevea brasiliensis, relies on the methylerythritol-phosphate pathway rather than the mevalonate pathway to synthesize isoprenyl diphosphate, as the MEP pathway operates predominantly in trans-polyisoprene-containing leaves and central peels. Chlorogenic acid biosynthesis pathway enzymes were preferentially expressed in leaves rather than in bark. This assembly with higher sequence contiguity can foster not only studies on genome structure and evolution, gene mapping, epigenetic analysis and functional genomics but also efforts to improve E. ulmoides for industrial and medical uses through genetic engineering.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Complete chloroplast genome of Hevea benthamiana, a SALB-resistant relative wild species of rubber tree

Hevea benthamiana is a SALB-resistant wild species of H. brasiliensis, the only source of mass production of high quality natural rubber. This study sequenced and analyzed the chloroplast genome of H. benthamiana. The chloroplast genome of H. benthamiana contains 161,124 bp and consists of 51,495 bp of A (31.96%), 52,022 bp of T (32.29%), 28,915 bp of G (17.95%), and 28,692 bp of C (17.81%). The ring-shaped genome includes four regions: a large single-copy region (LSC, 89,110 bp), a small single copy (SSC, 18,376 bp) region, and two inverted repeat regions (IRs, 26,819 bp). A total of 134 genes were annotated, of which 86 encode proteins; four are pseudogenes; 36 are tRNA genes, and eight are rRNA genes. Phylogenetic analyses showed that H. benthamiana is very closely related to H. Brasiliensis, this result indicates that H. benthamiana is highly valuable for the breeding of SALB-resistant varieties of rubber trees.