The water lily genome and the early evolution of flowering plants.

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1-3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. colorata genome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. colorata might support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. colorata genome sheds light on the early evolution of angiosperms.

EVLncRNAs 3.0: an updated comprehensive database for manually curated functional long non-coding RNAs validated by low-throughput experiments.

Long noncoding RNAs (lncRNAs) have emerged as crucial regulators across diverse biological processes and diseases. While high-throughput sequencing has enabled lncRNA discovery, functional characterization remains limited. The EVLncRNAs database is the first and exclusive repository for all experimentally validated functional lncRNAs from various species. After previous releases in 2018 and 2021, this update marks a major expansion through exhaustive manual curation of nearly 25 000 publications from 15 May 2020, to 15 May 2023. It incorporates substantial growth across all categories: a 154% increase in functional lncRNAs, 160% in associated diseases, 186% in lncRNA-disease associations, 235% in interactions, 138% in structures, 234% in circular RNAs, 235% in resistant lncRNAs and 4724% in exosomal lncRNAs. More importantly, it incorporated additional information include functional classifications, detailed interaction pathways, homologous lncRNAs, lncRNA locations, COVID-19, phase-separation and organoid-related lncRNAs. The web interface was substantially improved for browsing, visualization, and searching. ChatGPT was tested for information extraction and functional overview with its limitation noted. EVLncRNAs 3.0 represents the most extensive curated resource of experimentally validated functional lncRNAs and will serve as an indispensable platform for unravelling emerging lncRNA functions. The updated database is freely available at https://www.sdklab-biophysics-dzu.net/EVLncRNAs3/.

Haplotype-resolved genomes of octoploid species in Phyllanthaceae family reveal a critical role for polyploidization and hybridization in speciation.

The Phyllanthaceae family comprises a diverse range of plants with medicinal, edible, and ornamental value, extensively cultivated worldwide. Polyploid species commonly occur in Phyllanthaceae. Due to the rather complex genomes and evolutionary histories, their speciation process has been still lacking in research. In this study, we generated chromosome-scale haplotype-resolved genomes of two octoploid species (Phyllanthus emblica and Sauropus spatulifolius) in Phyllanthaceae family. Combined with our previously reported one tetraploid (Sauropus androgynus) and one diploid species (Phyllanthus cochinchinensis) from the same family, we explored their speciation history. The three polyploid species were all identified as allopolyploids with subgenome A/B. Each of their two distinct subgenome groups from various species was uncovered to independently share a common diploid ancestor (Ancestor-AA and Ancestor-BB). Via different evolutionary routes, comprising various scenarios of bifurcating divergence, allopolyploidization (hybrid polyploidization), and autopolyploidization, they finally evolved to the current tetraploid S. androgynus, and octoploid S. spatulifolius and P. emblica, respectively. We further discuss the variations in copy number of alleles and the potential impacts within the two octoploids. In addition, we also investigated the fluctuation of metabolites with medical values and identified the key factor in its biosynthesis process in octoploids species. Our study reconstructed the evolutionary history of these Phyllanthaceae species, highlighting the critical roles of polyploidization and hybridization in their speciation processes. The high-quality genomes of the two octoploid species provide valuable genomic resources for further research of evolution and functional genomics.

EVlncRNA-Dpred: improved prediction of experimentally validated lncRNAs by deep learning.

Long non-coding RNAs (lncRNAs) played essential roles in nearly every biological process and disease. Many algorithms were developed to distinguish lncRNAs from mRNAs in transcriptomic data and facilitated discoveries of more than 600 000 of lncRNAs. However, only a tiny fraction (<1%) of lncRNA transcripts (~4000) were further validated by low-throughput experiments (EVlncRNAs). Given the cost and labor-intensive nature of experimental validations, it is necessary to develop computational tools to prioritize those potentially functional lncRNAs because many lncRNAs from high-throughput sequencing (HTlncRNAs) could be resulted from transcriptional noises. Here, we employed deep learning algorithms to separate EVlncRNAs from HTlncRNAs and mRNAs. For overcoming the challenge of small datasets, we employed a three-layer deep-learning neural network (DNN) with a K-mer feature as the input and a small convolutional neural network (CNN) with one-hot encoding as the input. Three separate models were trained for human (h), mouse (m) and plant (p), respectively. The final concatenated models (EVlncRNA-Dpred (h), EVlncRNA-Dpred (m) and EVlncRNA-Dpred (p)) provided substantial improvement over a previous model based on support-vector-machines (EVlncRNA-pred). For example, EVlncRNA-Dpred (h) achieved 0.896 for the area under receiver-operating characteristic curve, compared with 0.582 given by sequence-based EVlncRNA-pred model. The models developed here should be useful for screening lncRNA transcripts for experimental validations. EVlncRNA-Dpred is available as a web server at https://www.sdklab-biophysics-dzu.net/EVlncRNA-Dpred/index.html, and the data and source code can be freely available along with the web server.

DNA methylation remodeled amino acids biosynthesis regulates flower senescence in carnation (Dianthus caryophyllus).

Dynamic DNA methylation regulatory networks are involved in many biological processes. However, how DNA methylation patterns change during flower senescence and their relevance with gene expression and related molecular mechanism remain largely unknown. Here, we used whole genome bisulfite sequencing to reveal a significant increase of DNA methylation in the promoter region of genes during natural and ethylene-induced flower senescence in carnation (Dianthus caryophyllus L.), which was correlated with decreased expression of DNA demethylase gene DcROS1. Silencing of DcROS1 accelerated while overexpression of DcROS1 delayed carnation flower senescence. Moreover, among the hypermethylated differentially expressed genes during flower senescence, we identified two amino acid biosynthesis genes, DcCARA and DcDHAD, with increased DNA methylation and reduced expression in DcROS1 silenced petals, and decreased DNA methylation and increased expression in DcROS1 overexpression petals, accompanied by decreased or increased amino acids content. Silencing of DcCARA and DcDHAD accelerates carnation flower senescence. We further showed that adding corresponding amino acids could largely rescue the senescence phenotype of DcROS1, DcCARA and DcDHAD silenced plants. Our study not only demonstrates an essential role of DcROS1-mediated remodeling of DNA methylation in flower senescence but also unravels a novel epigenetic regulatory mechanism underlying DNA methylation and amino acid biosynthesis during flower senescence.

First Report of Coptis Wilt Disease Caused by Plectosphaerella cucumerina in China.

Coptis (Coptis chinensis) belongs to the Ranunculaceae family, the rhizomes used in traditional Chinese medicine. Since 2021, an uncommon stem and leaf wilt disease, with an average disease incidence of 70%~90%, has been observed in Guangdong and Guangxi provinces. The early wilt symptoms were observed on older leaves and stems, and the whole seedling wilted and died. The rhizome of the diseased seedlings changed in color, became necrotic, and rotted. Symptomatic roots and stems were surface-sterilized with 70% ethanol for 30 s, followed by 0.2% NaClO for 2-3 min, rinsed in sterile water three times, and then placed on potato dextrose agar (PDA) at 25℃for 14 days. Fungal growth was observed, and six isolates with similar morphology were obtained. The 14-day-old colonies on PDA were buff with few aerial hyphae and slimy surfaces. Aerial hyphae were sparse with simple or branched conidiophores. Conidia were hyaline, smooth, ovoid, septate or aseptate, and 5.77 to 9.53 × 2.15 to 3.32 µm (n = 50). Three of the six isolates were subjected to further analysis. The genomic DNA of three isolates (CCF1-1, CCF1-2, CCF1-3) was extracted using Axygen MAG-FRAG- I-50 (Axygen Bio-Tek) for molecular identification. Partial sequences of the internal transcribed spacer of rDNA (ITS) and large subunit rDNA (LSU) were amplified using the primers ITS1/4 and LR5F/LROR, respectively (Vilgalys and Hester 1990). Their sequences were aligned by MEGA X (Kumar et al., 2018), and the sequences of each region showed 100% sequence similarity among our isolates. A BLAST search of ITS and LSU sequences (accession nos. ON377369, ON428244) showed that both regions had the highest nucleotide similarities (99.43 to 99.89%) to the Plectosphaerella cucumerina strains. Based on morphological and molecular analyses, the isolates were identified as P. cucumerina (Palm et al. 1995). The pathogenicity of our isolates CCF1-1, CCF1-2, CCF1-3 was tested on ten 2-month-old healthy seedlings of coptis, respectively. For the seedlings, 30 ml of fungal conidial suspension (1×106 conidia/ml) or sterile water, as control, were poured into their root area. Conidia suspension were prepared from 14-day-old cultures on PDA by eluting with sterilized water. The seedlings were incubated at 25°C and 75% relative humidity under a 12-h/12-h light/dark cycle. The test was repeated three times. After 20 days, only seedlings inoculated with P. cucumerina exhibited symptoms similar to those diseased seedlings in the field. The control seedling had no symptoms. The morphologically similar fungus was re-isolated from the tested seedlings, thus fulfilling Koch's postulates. Based on molecular, morphological, and pathogenic properties, P. cucumerina is the causal fungal pathogen of coptis wilt disease. Previously, P. cucumerina has been related to wilt disease in strawberry and Chinese cabbage (Yang et al. 2023; Gao et al. 2022), but to our knowledge, this is the first report of P. cucumerina causing wilt disease on coptis in China. Coptis wilt disease tends to occur in a warm and rainy environment, and strengthening the detection and quarantine of seedlings is the key to preventing the occurrence and spread of the disease.

First Report of Anthracnose on Wurfbainia villosa var. villosa Caused by Colletotrichum gloeosporioides in China.

Wurfbainia villosa var. villosa is a traditional Chinese herbal medicine under the family Zingiberaceae, and its ripe fruits (called Fructus Amomi) are widely used clinically for the treatment of gastrointestinal disorders (Yang et al. 2023; Chen et al. 2023). In September 2023, plants of W. villosa var. villosa exhibited anthracnose-like symptoms on leaf with a disease incidence of 35% (n = 100 investigated plants) in an approximately 90 m2 field in Guangning, China (N23°42'51.70″, E112°26'35.75″). Light yellowish-green spots (~2 mm diameter) initially appeared on the infected leaves, gradually formed sub-circular or irregular spots, then fused and expanded, resulting in wilting of the leaves. To identify the causal agent, 10 symptomatic leaves were collected and transferred to the laboratory. The symptomatic leaf samples were surface sterilized in 0.5% NaClO for 2 min, and in 70% ethanol for 30 s, then washed three times with sterile water and air-dried on sterile filter paper. The leaf tissues were placed on potato dextrose agar (PDA) medium containing 100 µg mL-1 of ampicillin (Sigma-Aldrich, St. Louis, MO) and incubated for 7 days at 28°C in darkness. Nine isolates with similar colony morphology were isolated from the 10 plated leaves. Three representative isolates (GNAF03, GNAF06, GNAF09 with approximately 3.5 cm in diameter after 3 days of incubation) appeared gray to dark brown with dense aerial hyphae at the front and gray to black colonies on the reverse of the plates. Conidia were cylindrical and measured 21.2 to 29.3 µm long × 7.1 to 9.6 µm wide (n = 50). Appressoria were formed by the tips of germ tubes or hyphae and were brown, ellipsoid, thick-walled, and smooth-margined, measuring 10.2 to 12.3 µm long × 6.4 to 8.2 µm wide (n = 50). Morphologically, the fungal isolates resembled Colletotrichum sp. (Weir et al. 2012). For molecular analysis, genomic DNA was extracted from fresh mycelia of the three isolates, and the primers ACT-512F/ACT-783R, CL1/CL2A, GDF/GDR, and ITS1/ITS4 were used to amplify partial regions of rDNA-ITS, actin (ACT), calmodulin (CAL), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) regions, respectively (Weir et al. 2012). The resulting sequences with more than 99% nucleotide identity to C. gloeosporioides were submitted to GenBank (accession numbers PP552725, PP552726, and OR827444 for ACT; PP552727, PP552728, and OR827443 for CAL; PP552729, PP552730, and OR827445 for GAPDH; PP549996, PP549999, and OR841394 for ITS). A phylogenetic tree was generated by the maximum likelihood method using the concatenated sequences of ACT, CAL, GADPH, and ITS by Polysuite software (Damm et al. 2020). Based on morphological and molecular analysis, the three isolates were characterized as C. gloeosporioides. The pathogenicity of the GNAF09 isolate was assessed on W. villosa var. villosa seedling leaves inoculated by spraying with 40 µL of conidial suspension at 106 conidia mL-1 or wounded with a sterile toothpick then inoculated with mycelial agar plugs (5 mm diameter). Control leaves were inoculated with 40 µL of sterile distilled water or agar plugs without mycelia. The inoculated plants were placed in a humid chamber at 28°C with 80% humidity and a 12 h light-dark photoperiod. Symptoms similar to those seen on naturally infected leaves were observed on all inoculated leaves after 7 days inoculation. Re-isolation was performed from 80% of the inoculated leaves and isolates were confirmed as C. gloeosporioides morphologically, confirming Koch's postulates, and by sequencing the ACT, CAL, GADPH, and ITS regions. The control groups remained asymptomatic. In previous studies, C. gloeosporioides has also caused anthracnose on Chinese medicinal plants, including Baishao (Radix paeoniae alba) (Zhang et al. 2017) and Rubia cordifolia L. (Tang et al. 2020). To our knowledge, this is the first report of C. gloeosporioides causing anthracnose on W. villosa var. villosa in China. The results of our report serve as valuable references for further research on this disease.

Effects of a mixture of glycerol monolaurate and cinnamaldehyde supplementation on laying performance, egg quality, and antioxidant status in laying hens.

BACKGROUND: This study aimed to determine the effects of a mixture of glycerol monolaurate and cinnamaldehyde (GCM) supplementation on the laying performance, egg quality, antioxidant capacity, and serum parameters of laying hens. A total of 1120 14-week-old Jingfen-1 strain laying hens with similar performance were randomly allocated to four dietary treatments: control, and GCM groups supplemented with 250, 500, or 1000 mg kg-1 for 12 weeks. RESULTS: Compared with the control group, GCM-supplemented groups significantly reduced (P < 0.05) the rate of unqualified eggs of laying hens aged 17-24 weeks. Supplementation of GCM significantly increased (P < 0.05) yolk color and serum glutathione peroxidase (GSH-Px) activity but decreased (P < 0.05) the hydrogen peroxide (H2 O2 ) content in the serum of laying hens at the age of 20 weeks. Furthermore, groups supplemented with GCM showed a significant increase (P < 0.05) in Haugh unit, yolk color, activities of total superoxide dismutase and GSH-Px, and the glucose content in serum, and a decrease (P < 0.05) in the content of urea nitrogen and H2 O2 and malondialdehyde in serum of laying hens at the age of 24 weeks. 500 mg kg-1 GCM supplementation significantly increased (P < 0.05) the number of large white follicles and 1000 mg kg-1 GCM supplementation decreased the number of large yellow follicles in 28-week-old laying hens. CONCLUSION: These results indicated that GCM supplementation has positive effects on reducing egg loss and improving egg quality in the early laying period of laying hens. © 2023 Society of Chemical Industry.

[Genome-wide identification and expression analysis of TCP gene family of Andrographis paniculata under abiotic stress].

Andrographis paniculata is an important medicinal plant in the Lingnan region of China, which has the functions of clearing heat, removing toxins, and resisting bacteria and inflammation. The TCP gene family is a class of transcription factors that regulate plant growth, development, and stress response. In order to analysis the role of the TCP gene family under abiotic stress in A. paniculata, this study identified the TCP gene family of A. paniculata at the genome-wide level and analyzed its expression pattern in response to abiotic stress. The results showed that the A. paniculata TCP gene family had 23 members, with length of amino acid ranging from 136 to 508, the relative molecular mass between 14 854.71 and 55 944.90 kDa, and the isoelectric point between 5.67 and 10.39. All members were located in the nucleus and unevenly distributed on 13 chromosomes. Phylogenetic analysis classified them into three subfamilies: PCF, CIN and CYC/TB1. Gene structure and conserved motif analysis showed that most members of the TCP gene family contained motif 1, motif 2, motif 3 in the same order and 1-3 CDS. The analysis of promoter cis-acting elements showed that the transcriptional expression of the TCP gene family in A. paniculata might be induced by light, hormones, and adversity stress. In light of the expression pattern analysis and qRT-PCR verification, the expression of ApTCP4, ApTCP5, ApTCP6, and ApTCP11 involved in response by various abiotic stresses such as drought, high temperature, and MeJA. This study lays the foundation for in-depth exploration of the functions of A. paniculata TCP genes in response to abiotic stress.

Genomic-wide identification and expression analysis of R2R3-MYB transcription factors related to flavonol biosynthesis in Morinda officinalis.

BACKGROUND: The R2R3-MYB transcription factors are a crucial and extensive gene family in plants, which participate in diverse processes, including development, metabolism, defense, differentiation, and stress response. In the Lingnan region of China, Morinda officinalis is extensively grown and is renowned for its use as both a medicinal herb and food source. However, there are relatively few reports on the R2R3-MYB transcription factor family in M.officinalis. RESULTS: In this study, we identified 97 R2R3-MYB genes in the genome of Morinda officinalis and classified them into 32 subgroups based on phylogenetic comparison with Arabidopsis thaliana. The lack of recent whole-genome duplication events in M.officinalis may be the reason for the relatively few members of the R2R3-MYB family. We also further analyzed the physical and chemical characteristics, conserved motifs, gene structure, and chromosomal location. Gene duplication events found 21 fragment duplication pairs and five tandem duplication event R2R3-MYB genes in M.officinalis may also affect gene family expansion. Based on phylogenetic analysis, cis-element analysis, co-expression analysis and RT-qPCR, we concluded that MoMYB33 might modulate flavonol levels by regulating the expression of 4-coumarate-CoA ligase Mo4CL2, chalcone isomerase MoCHI3, and flavonol synthase MoFLS4/11/12. MoMYB33 and AtMYB111 showed the highest similarity of 79% and may be involved in flavonol synthase networks by the STRING database. Moreover, we also identified MoMYB genes that respond to methyl Jasmonate (MeJA) and abscisic acid (ABA) stress by RT-qPCR. CONCLUSIONS: This study offers a thorough comprehension of R2R3-MYB in M.officinalis, which lays the foundation for the regulation of flavonol synthesis and the response of MoMYB genes to phytohormones in M.officinalis.

Genome-wide identification and expression pattern analysis of R2R3-MYB transcription factor gene family involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii.

BACKGROUND: R2R3-MYB transcription factors regulate secondary metabolism, stress responses and development in various plants. Puerarin is a bioactive ingredient and most abundant secondary metabolite isolated from Pueraria lobata. The biosynthesis of puerarin proceeds via the phenylpropanoid pathway and isoflavonoids pathway, in which 9 key enzymes are involved. The expression of these structural genes is under control of specific PtR2R3-MYB genes in different plant tissues. However, how PtR2R3-MYB genes regulates structural genes in puerarin biosynthesis remains elusive. This study mined the PtR2R3-MYB genes involved in puerarin biosynthesis and response to hormone in Pueraria lobata var. thomsonii. RESULTS: A total of 209 PtR2R3-MYB proteins were identified, in which classified into 34 subgroups based on the phylogenetic topology and the classification of the R2R3-MYB superfamily in Arabidopsis thaliana. Furtherly physical and chemical characteristics, gene structure, and conserved motif analysis were also used to further analyze PtR2R3-MYBs. Combining puerarin content and RNA-seq data, speculated on the regulated puerarin biosynthesis of PtR2R3-MYB genes and structural genes, thus 21 PtR2R3-MYB genes and 25 structural genes were selected for validation gene expression and further explore its response to MeJA and GSH treatment by using qRT-PCR analysis technique. Correlation analysis and cis-acting element analysis revealed that 6 PtR2R3-MYB genes (PtMYB039, PtMYB057, PtMYB080, PtMYB109, PtMYB115 and PtMYB138) and 7 structural genes (PtHID2, PtHID9, PtIFS3, PtUGT069, PtUGT188, PtUGT286 and PtUGT297) were directly or indirectly regulation of puerarin biosynthesis in ZG11. It is worth noting that after MeJA and GSH treatment for 12-24 h, the expression changes of most candidate genes were consistent with the correlation of puerarin biosynthesis, which also shows that MeJA and GSH have the potential to mediate puerarin biosynthesis by regulating gene expression in ZG11. CONCLUSIONS: Overall, this study provides a comprehensive understanding of the PtR2R3-MYB and will paves the way to reveal the transcriptional regulation of puerarin biosynthesis and response to phytohormone of PtR2R3-MYB genes in Pueraria lobata var. thomsonii.

A novel genome sequence of Jasminum sambac helps uncover the molecular mechanism underlying the accumulation of jasmonates.

Jasminum sambac is a well-known plant for its attractive and exceptional fragrance, the flowers of which are used to produce scented tea. Jasmonate (JA), an important plant hormone was first identified in Jasminum species. Jasmine plants contain abundant JA naturally, of which the molecular mechanisms of synthesis and accumulation are not clearly understood. Here, we report a telomere-to-telomere consensus assembly of a double-petal J. sambac genome along with two haplotype-resolved genomes. We found that gain-and-loss, positive selection, and allelic specific expression of aromatic volatile-related genes contributed to the stronger flower fragrance in double-petal J. sambac compared with single- and multi-petal jasmines. Through comprehensive comparative genomic, transcriptomic, and metabolomic analyses of double-petal J. sambac, we revealed the genetic basis of the production of aromatic volatiles and salicylic acid (SA), and the accumulation of JA under non-stress conditions. We identified several key genes associated with JA biosynthesis, and their non-stress related activities lead to extraordinarily high concentrations of JA in tissues. High JA synthesis coupled with low degradation in J. sambac results in accumulation of high JA under typical environmental conditions, similar to the accumulation mechanism of SA. This study offers important insights into the biology of J. sambac, and provides valuable genomic resources for further utilization of natural products.

EVLncRNAs 2.0: an updated database of manually curated functional long non-coding RNAs validated by low-throughput experiments.

Long non-coding RNAs (lncRNAs) play important functional roles in many diverse biological processes. However, not all expressed lncRNAs are functional. Thus, it is necessary to manually collect all experimentally validated functional lncRNAs (EVlncRNA) with their sequences, structures, and functions annotated in a central database. The first release of such a database (EVLncRNAs) was made using the literature prior to 1 May 2016. Since then (till 15 May 2020), 19 245 articles related to lncRNAs have been published. In EVLncRNAs 2.0, these articles were manually examined for a major expansion of the data collected. Specifically, the number of annotated EVlncRNAs, associated diseases, lncRNA-disease associations, and interaction records were increased by 260%, 320%, 484% and 537%, respectively. Moreover, the database has added several new categories: 8 lncRNA structures, 33 exosomal lncRNAs, 188 circular RNAs, and 1079 drug-resistant, chemoresistant, and stress-resistant lncRNAs. All records have checked against known retraction and fake articles. This release also comes with a highly interactive visual interaction network that facilitates users to track the underlying relations among lncRNAs, miRNAs, proteins, genes and other functional elements. Furthermore, it provides links to four new bioinformatics tools with improved data browsing and searching functionality. EVLncRNAs 2.0 is freely available at https://www.sdklab-biophysics-dzu.net/EVLncRNAs2/.

Chromosome-level genome assembly and population genetic analysis of a critically endangered rhododendron provide insights into its conservation.

Rhododendrons are woody plants, famous throughout the world as having high horticultural value. However, many wild species are currently threatened with extinction. Here, we report for the first time a high-quality, chromosome-level genome of Rhododendron griersonianum, which has contributed to approximately 10% of all horticultural rhododendron varieties but which in its wild form has been evaluated as critically endangered. The final genome assembly, which has a contig N50 size of approximately 34 M and a total length of 677 M, is the highest-quality genome sequenced within the genus to date, in part due to its low heterozygosity (0.18%). Identified repeats constitute approximately 57% of the genome, and 38 280 protein-coding genes were predicted with high support. We further resequenced 31 individuals of R. griersonianum as well as 30 individuals of its widespread relative R. delavayi, and performed additional conservation genomic analysis. The results showed that R. griersonianum had lower genetic diversity (θ = 2.58e-3; π = 1.94e-3) when compared not only to R. delavayi (θ = 11.61e-3, π = 12.97e-3), but also to most other woody plants. Furthermore, three severe genetic bottlenecks were detected using both the Stairway plot and fastsimcoal2 analysis, which are thought to have occurred in the late Middle Pleistocene and the Last Glacial Maximum (LGM) period. After these bottlenecks, R. griersonianum recovered and maintained a constant effective population size (>25 000) until now. Intriguingly, R. griersonianum has accumulated significantly more deleterious mutations in the homozygous state than R. delavayi, and several deleterious mutations (e.g., in genes involved in the response to heat stress) are likely to have harmed the adaptation of this plant to its surroundings. This high-quality, chromosome-level genome and the population genomic analysis of the critically endangered R. griersonianum will provide an invaluable resource as well as insights for future study in this species to facilitate conservation and in the genus Rhododendron in general.

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

MAIN CONCLUSION: 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.

The complete chloroplast genome sequence of the medicinal plant Abrus pulchellus subsp. cantoniensis: genome structure, comparative and phylogenetic relationship analysis.

Abrus pulchellus subsp. cantoniensis, an endemic medicinal plant in southern China, is clinically used to treat jaundice hepatitis, cholecystitis, stomachache and breast carbuncle. Here, we assembled and analyzed the first complete chloroplast (cp) genome of A. pulchellus subsp. cantoniensis. The A. pulchellus subsp. cantoniensis cp genome size is 156,497 bp with 36.5% GC content. The cp genome encodes 130 genes, including 77 protein-coding genes, 30 tRNA genes and four rRNA genes, of which 19 genes are duplicated in the inverted repeats (IR) regions. A total of 30 codons exhibited codon usage bias with A/U-ending. Moreover, 53 putative RNA editing sites were predicted in 20 genes, all of which were cytidine to thymine transitions. Repeat sequence analysis identified 45 repeat structures and 125 simple-sequence repeats (SSRs) in A. pulchellus subsp. cantoniensis cp genome. In addition, 19 mononucleotides (located in atpB, trnV-UAC, ycf3, atpF, rps16, rps18, clpP, rpl16, trnG-UCC and ndhA) and three compound SSRs (located in ndhA, atpB and rpl16) showed species specificity between A. pulchellus subsp. cantoniensis and Abrus precatorius, which might be informative sources for developing molecular markers for species identification. Furthermore, phylogenetic analysis inferred that A. pulchellus subsp. cantoniensis was closely related to A. precatorius, and the genus Abrus formed a subclade with Canavalia in the Millettioid/Phaseoloid clade. These data provide a valuable resource to facilitate the evolutionary relationship and species identification of this species.

MicroRNAs in Medicinal Plants.

Medicinal plant microRNAs (miRNAs) are an endogenous class of small RNA central to the posttranscriptional regulation of gene expression. Biosynthetic research has shown that the mature miRNAs in medicinal plants can be produced from either the standard messenger RNA splicing mechanism or the pre-ribosomal RNA splicing process. The medicinal plant miRNA function is separated into two levels: (1) the cross-kingdom level, which is the regulation of disease-related genes in animal cells by oral intake, and (2) the intra-kingdom level, which is the participation of metabolism, development, and stress adaptation in homologous or heterologous plants. Increasing research continues to enrich the biosynthesis and function of medicinal plant miRNAs. In this review, peer-reviewed papers on medicinal plant miRNAs published on the Web of Science were discussed, covering a total of 78 species. The feasibility of the emerging role of medicinal plant miRNAs in regulating animal gene function was critically evaluated. Staged progress in intra-kingdom miRNA research has only been found in a few medicinal plants, which may be mainly inhibited by their long growth cycle, high demand for growth environment, immature genetic transformation, and difficult RNA extraction. The present review clarifies the research significance, opportunities, and challenges of medicinal plant miRNAs in drug development and agricultural production. The discussion of the latest results furthers the understanding of medicinal plant miRNAs and helps the rational design of the corresponding miRNA/target genes functional modules.

Genome-Wide Identification of SnRK1 Catalytic α Subunit and FLZ Proteins in Glycyrrhiza inflata Bat. Highlights Their Potential Roles in Licorice Growth and Abiotic Stress Responses.

Sucrose non-fermenting-1-related protein kinase-1 (SnRK1) and its scaffolding proteins, FCS-like zinc finger proteins (FLZs), are well conserved in land plants and involved in various processes of plant growth and stress responses. Glycyrrhiza inflata Bat. is a widely used licorice species with strong abiotic stress resistance, in which terpenoids and flavonoids are the major bioactive components. Here, we identified 2 SnRK1 catalytic α subunit encoding genes (GiSnRK1α1 and GiSnRK1α2) and 21 FLZ genes in G. inflata. Polygenetic analysis showed that the 21 GiFLZs could be divided into three groups. A total of 10 representative GiFLZ proteins interact with GiSnRK1α1, and they display overlapped subcellular localization (mainly in the nucleus and the cytoplasm) when transiently expressed in Nicotiana benthamiana leaf cells. Coinciding with the existence of various phytohormone-responsive and stress-responsive cis-regulatory elements in the GiSnRK1α and GiFLZ gene promoters, GiFLZs are actively responsive to methyl jasmonic acid (MeJA) and abscisic acid (ABA) treatments, and several GiFLZs and GiSnRK1α1 are regulated by drought and saline-alkaline stresses. Interestingly, GiSnRK1α and 20 of 21 GiFLZs (except for GiFLZ2) show higher expression in the roots than in the leaves. These data provide comprehensive information on the SnRK1 catalytic α subunit and the FLZ proteins in licorice for future functional characterization.

[Identification and expression analysis of R2R3-MYB gene family in Andrographis paniculata].

The plant growth, development, and secondary metabolism are regulated by R2 R3-MYB transcription factors. This study identified the R2 R3-MYB genes in the genome of Andrographis paniculata and analyzed the chromosomal localization, gene structure, and conserved domains, phylogenetic relationship, and promoter cis-acting elements of these R2 R3-MYB genes. Moreover, the gene expression profiles of R2 R3-MYB genes under abiotic stress and hormone treatments were generated by RNA-seq and validated by qRT-PCR. The results showed that A. paniculata contained 73 R2 R3-MYB genes on 21 chromosomes. These members belonged to 34 subfamilies, 19 of which could be classified into the known subfamilies in Arabidopsis thaliana. The 73 R2 R3-MYB members included 36 acidic proteins and 37 basic proteins, with the lengths of 148-887 aa. The domains, motifs, and gene structures of R2 R3-MYBs in A. paniculata were conserved. The promoter regions of these genes contains a variety of cis-acting elements related to the responses to environmental factors and plant hormones including light, ABA, MeJA, and drought. Based on the similarity of functions of R2 R3-MYBs in the same subfamily and the transcription profiles, ApMYB13/21/35/67/73(S22) may regulate drought stress through ABA pathway; ApMYB20(S11) and ApMYB55(S2) may play a role in the response of A. paniculata to high temperature and UV-C stress; ApMYB5(S7) and ApMYB33(S20) may affect the accumulation of andrographolide by regulating the expression of key enzymes in the MEP pathway. This study provides theoretical reference for further research on the functions of R2 R3-MYB genes in A. paniculata and breeding of A. paniculata varieties with high andrographolide content.

The Rhododendron Plant Genome Database (RPGD): a comprehensive online omics database for Rhododendron.

BACKGROUND: The genus Rhododendron L. has been widely cultivated for hundreds of years around the world. Members of this genus are known for great ornamental and medicinal value. Owing to advances in sequencing technology, genomes and transcriptomes of members of the Rhododendron genus have been sequenced and published by various laboratories. With increasing amounts of omics data available, a centralized platform is necessary for effective storage, analysis, and integration of these large-scale datasets to ensure consistency, independence, and maintainability. RESULTS: Here, we report our development of the Rhododendron Plant Genome Database (RPGD; http://bioinfor.kib.ac.cn/RPGD/ ), which represents the first comprehensive database of Rhododendron genomics information. It includes large amounts of omics data, including genome sequence assemblies for R. delavayi, R. williamsianum, and R. simsii, gene expression profiles derived from public RNA-Seq data, functional annotations, gene families, transcription factor identification, gene homology, simple sequence repeats, and chloroplast genome. Additionally, many useful tools, including BLAST, JBrowse, Orthologous Groups, Genome Synteny Browser, Flanking Sequence Finder, Expression Heatmap, and Batch Download were integrated into the platform. CONCLUSIONS: RPGD is designed to be a comprehensive and helpful platform for all Rhododendron researchers. Believe that RPGD will be an indispensable hub for Rhododendron studies.