Transcriptome-based investigation of cirrus development and identifying microsatellite markers in rattan (Daemonorops jenkinsiana)

Analyzing lignin biosynthesis pathways in rattan using improved co-expression networks of NACs and MYBs

BACKGROUND

The rattan is a valuable plant resource with multiple applications in tropical forests. Calamus simplicifolius and Daemonorops jenkinsiana are the two most representative rattan species, supplying over 95% of the raw materials for the rattan industry. Hence, the wood properties of both rattans have always attracted researchers' attention.

RESULTS

We re-annotated the genomes, obtained 81 RNA-Seq datasets, and developed an improved pipeline to increase the reliability of co-expression networks of both rattans. Based on the data and pipeline, co-expression relationships were detected in 11 NACs, 49 MYBs, and 86 lignin biosynthesis genes in C. simplicifolius and four NACs, 59 MYBs, and 76 lignin biosynthesis genes in D. jenkinsiana, respectively. Among these co-expression pairs, several genes had a close relationship to the development of wood properties. Additionally, we detected the enzyme gene on the lignin biosynthesis pathway was regulated by either NAC or MYB, while LACCASES was regulated by both NAC and MYB. For D. jenkinsiana, the lignin biosynthesis regulatory network was characterized by positive regulation, and MYB possible negatively regulate non-expressed lignin biosynthesis genes in stem tissues. For C. simplicifolius, NAC may positively regulate highly expressed genes and negatively regulate non-expressed lignin biosynthesis genes in stem tissues. Furthermore, we established core regulatory networks of NAC and MYB for both rattans.

CONCLUSIONS

This work improved the accuracy of rattan gene annotation by integrating an efficient co-expression network analysis pipeline, enhancing gene coverage and accuracy of the constructed network, and facilitating an understanding of co-expression relationships among NAC, MYB, and lignin biosynthesis genes in rattan and other plants.

Draft genome of Korthalsia laciniosa (Griff.) Mart., a climbing rattan elucidates its phylogenetic position

Korthalsia laciniosa (Griff.) Mart. is a climbing rattan used as a source of durable and flexible cane. In the present study, the draft genome of K. laciniosa was sequenced, de novo assembled and annotated. Genome-wide identification of MADS-Box transcription factors revealed loss of Mβ, and Mγ genes belonging to Type I subclass in the rattan lineage. Mining of the genome revealed presence of 13 families of lignin biosynthetic pathway genes and expression profiling of nine major genes documented relatively lower level of expression in cirrus when compared to leaflet and petiole. The chloroplast genome was re-constructed and analysis revealed the phylogenetic relatedness of this genus to Eugeissona, in contrast with its present taxonomic position. The genomic resource generated in the present study will accelerate population structure analysis, genetic resource conservation, phylogenomics and facilitate understanding the unique developmental processes like gender expression at molecular level.

Coexpression Analysis Reveals Dynamic Modules Regulating the Growth and Development of Cirri in the Rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Rattan is regarded as one of the major non-timber forest products, second only to wood and bamboo, worldwide. Although the published genomes of Calamus simplicifolius and Daemonorops jenkinsiana have facilitated genome-wide gene functional analyses, coexpression networks (CENs) provide more comprehensive and complete annotations of gene function at the transcriptome level. Thus, we analyzed the CENs of the two rattans, C. simplicifolius and D. jenkinsiana, by integrating the genome sequences and analyzing in-house transcriptome data from different development stages of their cirri using a well-developed strategy. A total of 3,504 and 3,027 functional modules were identified in C. simplicifolius and D. jenkinsiana, respectively, based on a combination of CENs, gene family classification, and function enrichment tools. These modules covered the major developmental processes, including photosynthesis, lignin biosynthesis, flavonoid biosynthesis, and phenylpropanoid biosynthesis. Reference annotations were refined using CENs and functional modules. Moreover, we obtained novel insights into the regulation of cirrus growth and development in rattans. Furthermore, Rattan-NET (http://rattan.bamboogdb.org/), an online database with analysis tools for gene set enrichment analysis, module enrichment, network comparison analysis, and cis-element analysis, was constructed for the easy analysis of gene function and regulation modules involved in the growth and development of cirri in rattans.

The chromosome-level genome assemblies of two rattans (Calamus simplicifolius and Daemonorops jenkinsiana)

Background

Calamus simplicifolius and Daemonorops jenkinsiana are two representative rattans, the most significant material sources for the rattan industry. However, the lack of reference genome sequences is a major obstacle for basic and applied biology on rattan.

Findings

We produced two chromosome-level genome assemblies of C. simplicifolius and D. jenkinsiana using Illumina, Pacific Biosciences, and Hi-C sequencing data. A total of ∼730 Gb and ∼682 Gb of raw data covered the predicted genome lengths (∼1.98 Gb of C. simplicifolius and ∼1.61 Gb of D. jenkinsiana) to ∼372 × and ∼426 × read depths, respectively. The two de novo genome assemblies, ∼1.94 Gb and ∼1.58 Gb, were generated with scaffold N50s of ∼160 Mb and ∼119 Mb in C. simplicifolius and D. jenkinsiana, respectively. The C. simplicifolius and D. jenkinsiana genomes were predicted to harbor  51,235 and  53,342 intact protein-coding gene models, respectively. Benchmarking Universal Single-Copy Orthologs evaluation demonstrated that genome completeness reached 96.4% and 91.3% in the C. simplicifolius and D. jenkinsiana genomes, respectively. Genome evolution showed that four Arecaceae plants clustered together, and the divergence time between the two rattans was ∼19.3 million years ago. Additionally, we identified 193 and 172 genes involved in the lignin biosynthesis pathway in the C. simplicifolius and D. jenkinsiana genomes, respectively.

Conclusions

We present the first de novo assemblies of two rattan genomes (C. simplicifolius and D. jenkinsiana). These data will not only provide a fundamental resource for functional genomics, particularly in promoting germplasm utilization for breeding, but also serve as reference genomes for comparative studies between and among different species.