Best known as licorice, Glycyrrhiza Linn., a genus of herbaceous perennial legume, has been used as a traditional herbal medicine in Asia and a flavoring agent for tobacco and food industry in Europe and America. Abiotic stresses and hormonal treatments can significantly impact the development and metabolism of secondary metabolites in Glycyrrhiza. To better understand the biosynthesis of the trace-amount bioactive compounds, we first screened for the suitable reference genes for quantitative real-time reverse transcription PCR (qRT-PCR) analysis in Glycyrrhiza. The expression profiles of 14 candidate reference genes, including Actin1 (ACT), Clathrin complex AP1 (CAC), Cyclophilin (CYP), Heat-shock protein 40 (DNAJ), Dehydration responsive element binding gene (DREB), Translation elongation factor1 (EF1), Ras related protein (RAN), Translation initiation factor (TIF1), β-Tubulin (TUB), Ubiquitin-conjugating enzyme E2 (UBC2), ATP binding-box transpoter 2 (ABCC2), COP9 signal compex subunit 3 (COPS3), Citrate synthase (CS), and R3H domain protein 2 (R3HDM2) from two congeneric species, Glycyrrhiza uralensis F. and Glycyrrhiza inflata B., were examined under abiotic stresses (osmotic and salinity) and hormonal treatments (Abscisic acid (ABA) and methyl jasmonic acid (MeJA)) using a panel of software, including geNorm, NormFinder, BestKeeper, and Delta CT. The overall stability, however, was provided by RefFinder, a comprehensive ranking system integrating inputs from all four algorithms. In G. uralensis, the most stable reference genes under osmotic stress, salt stress, ABA treatment, and MeJA treatment were TIF1, DNAJ, CS, and ABCC2 for leaves and DNAJ, DREB, CAC, and CAC for roots, respectively. In comparison, the top ranked genes were TUB, CAC, UBC2, and RAN for leaves and TIF1, ABCC2, CAC, and UBC2 for roots, respectively, under stress and hormonal treatments in G. inflata. ACT and TIF1, on the other hand, were the least stable genes under the most experimental conditions in the two congeneric species. Finally, our survey of the reference genes in legume shows that EF, ACT, UBC2, and TUB were the top choices for the abiotic stresses while EF, UBC2, CAC, and ABCC2 were recommended for the hormonal treatments in Leguminosae. Our combined results provide reliable normalizers for accurate gene quantifications in Glycyrrhiza species, which will allow us to exploit its medicinal potential in general and antiviral activities in particular.

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.

Hygrophila polysperma is a type of amphibious plant that originates from Acanthaceae. Here, we report its first complete chloroplast (cp) genome. The complete cp genome is 146,675 bp in length with 38.3% of GC content. There are 130 genes including 86 protein coding genes, 36 tRNA genes, and 8 rRNA genes in this genome. Simple short sequence (SSR) analysis found 30 SSRs, 24 of which are located in a large single-copy region. Nucleotide diversity identified six most divergent sequences (trns-GCU, psaA-pafI, psaI-pafII, ycf2, rpl32, and ycf1) among 3 close-related species, H. polysperma, H. ringens, and Asteracantha longifolia. A phylogenetic tree among H. polysperma and another 30 related species was constructed based on the common coding sequence of the cp genome and showed that H. polysperma is most closely related to H. ringens (both belong to subtribe Hygrophilinae) and, together, they form a clade that is sister to A. longifolia. This study provides a basis for systemic and evolution studies as well as the development of molecular markers for species identification and genetic breeding.