Genome-wide identification and expression analysis of the WRKY transcription factor family in flax (Linum usitatissimum L.).

BACKGROUND: Members of the WRKY protein family, one of the largest transcription factor families in plants, are involved in plant growth and development, signal transduction, senescence, and stress resistance. However, little information is available about WRKY transcription factors in flax (Linum usitatissimum L.). RESULTS: In this study, comprehensive genome-wide characterization of the flax WRKY gene family was conducted that led to prediction of 102 LuWRKY genes. Based on bioinformatics-based predictions of structural and phylogenetic features of encoded LuWRKY proteins, 95 LuWRKYs were classified into three main groups (Group I, II, and III); Group II LuWRKYs were further assigned to five subgroups (IIa-e), while seven unique LuWRKYs (LuWRKYs 96-102) could not be assigned to any group. Most LuWRKY proteins within a given subgroup shared similar motif compositions, while a high degree of motif composition variability was apparent between subgroups. Using RNA-seq data, expression patterns of the 102 predicted LuWRKY genes were also investigated. Expression profiling data demonstrated that most genes associated with cellulose, hemicellulose, or lignin content were predominantly expressed in stems, roots, and less in leaves. However, most genes associated with stress responses were predominantly expressed in leaves and exhibited distinctly higher expression levels in developmental stages 1 and 8 than during other stages. CONCLUSIONS: Ultimately, the present study provides a comprehensive analysis of predicted flax WRKY family genes to guide future investigations to reveal functions of LuWRKY proteins during plant growth, development, and stress responses.

First Report of Leaf Spot on Industrial Hemp (Cannabis sativa L.) Caused by Alternaria alternata in China.

Industrial hemp is an economically important plant with traditional uses for textiles, paper, building materials, food and medicine (Li 1974; Russo et al. 2008; Zlas et al. 1993). In August 2020, an estimated 80% of the industrial hemp plants with leaf spots were observed in greenhouse in Minzhu town, Harbin City, Heilongjiang Province, China (45.8554°N, 126.8167°E), resulting in yield losses of 20%. Leaf symptoms began as small spots on the upper surface of leaves and gradually developed into brown spots with light yellow halos. These irregular spots expanded gradually and eventually covered the entire leaf; the center of the spots was easily perforated. To identify the pathogen, 20 diseased leaves were collected, and small sections of (3 to 5 mm) were taken from the margins of lesions of infected leaves. The pieces were sterilized with 75% alcohol for 30 s, a 0.1% mercuric chloride solution for 1 min, and then rinsed three times with sterile water. Samples were then cultured on potato dextrose agar at 28℃ in darkness for 4 days. A single-spore culture was obtained by monosporic isolation. Conidiophores were simple or branched, straight or flexuous, brown, and measured 22 to 61 µm long × 4 to 5 µm wide (n = 50). Conidia were solitary or in chains, brown or dark brown, obclavate, obpyriform or ellipsoid. Conidia ranged from 23 to 55 µm long × 10 to 15 µm wide (n = 50) with one to eight transverse and several longitudinal septa. For molecular identification (Jayawardena et al. 2019), genomic DNA of pathogenic isolate (MZ1287) was extracted by a cetyltrimethylammonium bromide protocol. Four gene regions including the rDNA internal transcribed spacer (ITS), glyceraldehyde-3-phosplate dehydrogenase (GAPDH), translation elongation factor 1-alpha (TEF1) and RNA polymerase II beta subunit (RPB2) were amplified with primers ITS1/ITS4, GDF1/GDR1, EF1-728F/EF1-986R and RPB2-5F/RPB2-7cR, respectively (White et al. 1990). Resulting sequences were deposited in GenBank with accession numbers of MW272539.1, MW303956.1, MW415414.1 and MW415413.1, respectively. A BLASTn analysis showed 100% homology with A. alternata (GenBank accession nos. MN615420.1, MH926018.1, MN615423.1 and KP124770.1), respectively. A neighbor-joining phylogenetic tree was constructed by combining all sequenced loci in MEGA7. The isolate MZ1287 clustered in the A. alternata clade with 100% bootstrap support. Thus, based on morphological (Simmons 2007) and molecular characteristics, the pathogen was identified as A. alternata. To test pathogenicity, leaves of ten healthy, 2-month-old potted industrial hemp plants were sprayed using a conidial suspension (1×106 spores/ml). Control plants were sprayed with sterile water. All plants were incubated in a greenhouse at 25℃ for a 16 h light and 8 h dark period at 90% relative humidity. The experiment was repeated three times. After two weeks, leaf spots of industrial hemp developed on the inoculated leaves while the control plants remained asymptomatic. The A. alternata pathogen was re-isolated from the diseased leaves on inoculated plants, fulfilling Koch's postulates. Based on morphology, sequencing, and pathogenicity test, the pathogen was identified as A. alternata. To our knowledge, this is the first report of A. alternata causing leaf spot disease of industrial hemp (Cannabis sativa L.) in China and is worthy of our attention for the harm it may cause to industrial hemp production.

Exogenous Application of Phytohormones Promotes Growth and Regulates Expression of Wood Formation-Related Genes in Populus simonii × P. nigra.

Although phytohormones are known to be important signal molecules involved in wood formation, their roles are still largely unclear. Here, Populus simonii × P. nigra seedlings were treated with different concentrations of exogenous phytohormones, indole-3-acetic acid (IAA), gibberellin (GA3), and brassinosteroid (BR), and the effects of phytohormones on growth were investigated. Next, 27 genes with known roles in wood formation were selected for qPCR analysis to determine tissue-specificity and timing of responses to phytohormone treatments. Compared to the control, most IAA, GA3, and BR concentrations significantly increased seedling height. Meanwhile, IAA induced significant seedling stem diameter and cellulose content increases that peaked at 3 and 30 mg·L-1, respectively. Significant increase in cellulose content was also observed in seedlings treated with 100 mg·L-1 GA3. Neither stem diameter nor cellulose content of seedlings were affected by BR treatment significantly, although slight effects were observed. Anatomical measurements demonstrated improved xylem, but not phloem, development in IAA- and BR-treated seedlings. Most gene expression patterns induced by IAA, GA3, and BR differed among tissues. Many IAA response genes were also regulated by GA3, while BR-induced transcription was weaker and slower in Populus than for IAA and GA3. These results reveal the roles played by phytohormones in plant growth and lay the foundation for exploring molecular regulatory mechanisms of wood formation in Populus.

Genomic variations and association study of agronomic traits in flax.

BACKGROUND: Flax (Linum usitatissimum. L) is an ancient oilseed and natural fiber crop. It could be divided into three categories by use, namely oil flax, fiber flax and oil-fiber dual purpose (OF). Cultivated flax is widely used in the food and textile industry. It is of great significance to elucidate the genetic characteristics of flax collections for accelerating the process of breeding improvement in this dual purpose crop. With the development of next-generation sequencing, we can use new methods, such as SLAF-seq (specific-locus amplified fragment sequencing), to decode unknown genomes of species. In this study, a high-through sequencing of flax collections using SLAF-seq was conducted. The evolutionary tendency was defined and candidate genes associated with agronomic traits of flax species were identified by Genome-Wide Association Studying (GWAS). RESULTS: A flax collection consisting of 224 varieties were sequenced by SLAF-seq. In total, 346,639 SLAF tags were developed from all accessions, with an average sequencing depth of 7.19 for each accession. A total of 584,987 SNPs (single nucleotide polymorphism) with an MAF > 0.05 were identified from these SLAFs. The population structure division and phylogenetic analysis indicated a strong divergence among three kinds of flax groups. The genome-wide variation uncovered that oil flax had the highest genetic diversity and was considered to be the ancestor of fiber flax and oil-fiber flax. Sixteen associated peak SNPs for six traits were obtained by GWAS of oil-related traits using EMMAX (efficient mixed-model association eXpedited). Candidate genes and their related pathway were evaluated. A new GWAS was developed for fiber properties using the GLM (General linear model) model and a number of loci were identified. CONCLUSIONS: To our knowledge, this is the first study on discovery multiple loci for important agronomic traits of flax species using GWAS strategy. These results will provide the highest possibility of incorporating both high fiber and good oil traits in a single variety.

Advancement of Research Progress on Synthesis Mechanism of Cannabidiol (CBD).

Cannabis sativa L. is a multipurpose crop with high value for food, textiles, and other industries. Its secondary metabolites, including cannabidiol (CBD), have potential for broad application in medicine. With the CBD market expanding, traditional production may not be sufficient. Here we review the potential for the production of CBD using biotechnology. We describe the chemical and biological synthesis of cannabinoids, the associated enzymes, and the application of metabolic engineering, synthetic biology, and heterologous expression to increasing production of CBD.

Identification of differentially expressed genes in flax (Linum usitatissimum L.) under saline-alkaline stress by digital gene expression.

The salinization and alkalization of soil are widespread environmental problems, and alkaline salt stress is more destructive than neutral salt stress. Therefore, understanding the mechanism of plant tolerance to saline-alkaline stress has become a major challenge. However, little attention has been paid to the mechanism of plant alkaline salt tolerance. In this study, gene expression profiling of flax was analyzed under alkaline-salt stress (AS2), neutral salt stress (NSS) and alkaline stress (AS) by digital gene expression. Three-week-old flax seedlings were placed in 25 mM Na2CO3 (pH11.6) (AS2), 50mM NaCl (NSS) and NaOH (pH11.6) (AS) for 18 h. There were 7736, 1566 and 454 differentially expressed genes in AS2, NSS and AS compared to CK, respectively. The GO category gene enrichment analysis revealed that photosynthesis was particularly affected in AS2, carbohydrate metabolism was particularly affected in NSS, and the response to biotic stimulus was particularly affected in AS. We also analyzed the expression pattern of five categories of genes including transcription factors, signaling transduction proteins, phytohormones, reactive oxygen species proteins and transporters under these three stresses. Some key regulatory gene families involved in abiotic stress, such as WRKY, MAPKKK, ABA, PrxR and ion channels, were differentially expressed. Compared with NSS and AS, AS2 triggered more differentially expressed genes and special pathways, indicating that the mechanism of AS2 was more complex than NSS and AS. To the best of our knowledge, this was the first transcriptome analysis of flax in response to saline-alkaline stress. These data indicate that common and diverse features of saline-alkaline stress provide novel insights into the molecular mechanisms of plant saline-alkaline tolerance and offer a number of candidate genes as potential markers of tolerance to saline-alkaline stress.