Harnessing the power of microbes: Enhancing soybean growth in an acidic soil through AMF inoculation rather than P-fertilization.

The low phosphorus (P) availability of acidic soils severely limits leguminous plant growth and productivity. Improving the soil P nutritional status can be achieved by increasing the P-content through P-fertilization or stimulating the mineralization of organic P via arbuscular mycorrhizal fungi (AMF) application; however, their corresponding impacts on plant and soil microbiome still remain to be explored. Here, we examined the effects of AMF-inoculation and P-fertilization on the growth of soybean with different P-efficiencies, as well as the composition of rhizo-microbiome in an acidic soil. The growth of recipient soybean NY-1001, which has a lower P-efficiency, was not significantly enhanced by AMF-inoculation or P-fertilization. However, the plant biomass of higher P-efficiency transgenic soybean PT6 was significantly increased by 46.74%-65.22% through AMF-inoculation. Although there was no discernible difference in plant biomass between PT6 and NY-1001 in the absence of AMF-inoculation and P-fertilization, PT6 had approximately 1.9-2.5 times the plant biomass of NY-1001 after AMF-inoculation. Therefore, the growth advantage of higher P-efficiency soybean was achieved through the assistance of AMF rather than P-fertilization in available P-deficient acidic soil. Most nitrogen (N)-fixing bacteria and some functional genes related to N-fixation were abundant in endospheric layer, as were the P-solubilizing Pseudomonas plecoglossicida, and annotated P-metabolism genes. These N-fixing and P-solubilizing bacteria were positive correlated with each other. Lastly, the two most abundant phytopathogenic fungi species accumulated in endospheric layer, they exhibited positive correlations with N-fixing bacteria, but displayed negative interactions with the majority of the other dominant non-pathogenic genera with potential antagonistic activity.

PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data.

Complete mitochondrial genomes (mitogenomes) of plants are valuable resources for nucleocytoplasmic interactions, plant evolution, and plant cytoplasmic male sterile line breeding. However, the complete assembly of plant mitogenomes is challenging due to frequent recombination events and horizontal gene transfers. Previous studies have adopted Illumina, PacBio, and Nanopore sequencing data to assemble plant mitogenomes, but the poor assembly completeness, low sequencing accuracy, and high cost limit the sampling capacity. Here, we present an efficient assembly toolkit (PMAT) for de novo assembly of plant mitogenomes using low-coverage HiFi sequencing data. PMAT has been applied to the de novo assembly of 13 broadly representative plant mitogenomes, outperforming existing organelle genome assemblers in terms of assembly accuracy and completeness. By evaluating the assembly of plant mitogenomes from different sequencing data, it was confirmed that PMAT only requires 1× HiFi sequencing data to obtain a complete plant mitogenome. The source code for PMAT is available at https://github.com/bichangwei/PMAT. The developed PMAT toolkit will indeed accelerate the understanding of evolutionary variation and breeding application of plant mitogenomes.

The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees.

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

A modified natural small molecule inhibits triple-negative breast cancer growth by interacting with Tubb3.

BACKGROUND: Triple-negative breast cancer (TNBC) is a malignant tumor without specific therapeutic targets and a poor prognosis. Chemotherapy is currently the first-line therapeutic option for TNBC. However, due to the heterogeneity of TNBC, not all of TNBC patients are responsive to chemotherapeutic agents. Therefore, the demand for new targeted agents is critical. ß-tubulin isotype III (Tubb3) is a prognostic factor associated with cancer progression, including breast cancer, and targeting Tubb3 may lead to improve TNBC disease control. Shikonin, the active compound in the roots of Lithospermun erythrorhizon suppresses the growth of various types of tumors, and its efficacy can be improved by altering its chemical structure. PURPOSE: In this work, the anti-TNBC effect of a shikonin derivative (PMMB276) was investigated, and its mechanism was also investigated. STUDY DESIGN/METHODS: This study combines flow cytometry, immunofluorescence staining, immunoblotting, immunoprecipitation, siRNA silencing, and the iTRAQ proteomics assay to analyze the inhibition potential of PMMB276 on TNBC. In vivo study was performed, Balb/c female murine models with or without the small molecule treatments. RESULTS: Herein, we screened 300 in-house synthesized analogs of shikonin against TNBC and identified a novel small molecule, PMMB276; it suppressed cell proliferation, induced apoptosis, and arrested the cell cycle at the G2/M phase, suggesting that it could have a tumor suppressive role in TNBC. Tubb3 was identified as the target of PMMB276 using proteomic and biological activity analyses. Meanwhile, PMMB276 regulated microtubule dynamics in vitro by inducing microtubule depolymerization and it could act as a tubulin stabilizer by a different process than that of paclitaxel. Moreover, suppressing or inhibiting Tubb3 with PMMB276 reduced the growth of breast cancer in an experimental mouse model, indicating that Tubb3 plays a significant role in TNBC progression. CONCLUSION: The findings support the therapeutic potential of PMMB276, a Tubb3 inhibitor, as a treatment for TNBC. Our findings might serve as a foundation for the utilization of shikonin and its derivatives in the development of anti-TNBC.

Genome-Wide Comparison and Functional Characterization of HMGR Gene Family Associated with Shikonin Biosynthesis in Lithospermum erythrorhizon.

3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), as the rate-limiting enzyme in the mevalonate pathway, is essential for the biosynthesis of shikonin in Lithospermum erythrorhizon. However, in the absence of sufficient data, the principles of a genome-wide in-depth evolutionary exploration of HMGR family members in plants, as well as key members related to shikonin biosynthesis, remain unidentified. In this study, 124 HMGRs were identified and characterized from 36 representative plants, including L. erythrorhizon. Vascular plants were found to have more HMGR family genes than nonvascular plants. The phylogenetic tree revealed that during lineage and species diversification, the HMGRs evolved independently and intronless LerHMGRs emerged from multi-intron HMGR in land plants. Among them, Pinus tabuliformis and L. erythrorhizon had the most HMGR gene duplications, with 11 LerHMGRs most likely expanded through WGD/segmental and tandem duplications. In seedling roots and M9 cultured cells/hairy roots, where shikonin biosynthesis occurs, LerHMGR1 and LerHMGR2 were expressed significantly more than other genes. The enzymatic activities of LerHMGR1 and LerHMGR2 further supported their roles in catalyzing the conversion of HMG-CoA to mevalonate. Our findings provide insight into the molecular evolutionary properties and function of the HMGR family in plants and a basis for the genetic improvement of efficiently produced secondary metabolites in L. erythrorhizon.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Design, synthesis and biological evaluation of novel modified dual-target shikonin derivatives for colorectal cancer treatment.

Warburg effect provides energy and material essential for tumor proliferation, the reverse of Warburg effect provides insights into the development of a novel anti-cancer strategy. Pyruvate kinase 2 (PKM2) and pyruvate dehydrogenase kinase 1 (PDK1) are two key enzymes in tumor glucose metabolism pathway that not only contribute to the Warburg effect through accelerating aerobic glycolysis, but also serve as druggable target for colorectal cancer (CRC). Considering that targeting PKM2 or PDK1 alone does not seem to be sufficient to remodel abnormal glucose metabolism and achieve significant antitumor activity, a series of novel benzenesulfonyl shikonin derivatives were designed to regulate PKM2 and PDK1 simultaneously. By means of molecular docking and antiproliferative screen, we found that compound Z10 could act as the combination of PKM2 activator and PDK1 inhibitor, thereby significantly inhibited glycolysis that reshaping tumor metabolism. Moreover, Z10 could inhibit proliferation, migration and induce apoptosis in CRC cell HCT-8. Finally, the in vivo anti-tumor activity of Z10 was evaluated in a colorectal cancer cell xenograft model in nude mice and the results demonstrated that Z10 induced tumor cell apoptosis and inhibited tumor cell proliferation with lower toxicity than shikonin. Our findings indicated that it is feasible to alter tumor energy metabolism through multi-target synergies, and the dual-target benzenesulfonyl shikonin derivative Z10 could be a potential anti-CRC agent.

The Metasequoia genome and evolutionary relationships among redwoods.

Redwood trees (Sequoioideae), including Metasequoia glyptostroboides (dawn redwood), Sequoiadendron giganteum (giant sequoia), and Sequoia sempervirens (coast redwood), are threatened and widely recognized iconic tree species. Genomic resources for redwood trees could provide clues to their evolutionary relationships. Here, we report the 8-Gb reference genome of M. glyptostroboides and a comparative analysis with two related species. More than 62% of the M. glyptostroboides genome is composed of repetitive sequences. Clade-specific bursts of long terminal repeat retrotransposons may have contributed to genomic differentiation in the three species. The chromosomal synteny between M. glyptostroboides and S. giganteum is extremely high, whereas there has been significant chromosome reorganization in S. sempervirens. Phylogenetic analysis of marker genes indicates that S. sempervirens is an autopolyploid, and more than 48% of the gene trees are incongruent with the species tree. Results of multiple analyses suggest that incomplete lineage sorting (ILS) rather than hybridization explains the inconsistent phylogeny, indicating that genetic variation among redwoods may be due to random retention of polymorphisms in ancestral populations. Functional analysis of ortholog groups indicates that gene families of ion channels, tannin biosynthesis enzymes, and transcription factors for meristem maintenance have expanded in S. giganteum and S. sempervirens, which is consistent with their extreme height. As a wetland-tolerant species, M. glyptostroboides shows a transcriptional response to flooding stress that is conserved with that of analyzed angiosperm species. Our study offers insights into redwood evolution and adaptation and provides genomic resources to aid in their conservation and management.

The sex determination gene of Populus deltoides, PdFERR, interacts with ERF96 to promote the development of female flower organs.

The female-specifically expressed response regulator (PdFERR) gene in Populus deltoides, a sex determination gene (an orthologous gene of ARR17 in Populus tremula), was found to promote femaleness in heterologous expression lines of Arabidopsis. None of the genes in the Arabidopsis genome seem to be orthologous to PdFERR. Although originating from two evolutionarily distant plants, the dioecious poplar FERR might promote femaleness in the hermaphroditic Arabidopsis through an evolutionary consistent regulatory pathway. However, there is no molecular evidence to support this viewpoint. In this study, to identify the shared downstream orthologous gene of PdFERR, we used yeast two-hybrid assay to screen potential interactors of PdFERR in Arabidopsis. We identified the ethylene response factor 96 (AtERF96) and confirmed the interaction via in vivo and in vitro assays. The ERF96 orthologous gene in P. deltoides was also experimentally confirmed to interact with PdFERR. PdFERR could then promote femaleness in poplar or Arabidopsis through interactions with ERF96, which provide a new perspective for understanding the PdFERR gene regulating sex differentiation.

Genome-Wide Analysis of the Expansin Gene Family in Populus and Characterization of Expression Changes in Response to Phytohormone (Abscisic Acid) and Abiotic (Low-Temperature) Stresses.

Expansins are a group of cell wall enzyme proteins that help to loosen cell walls by breaking hydrogen bonds between cellulose microfibrils and hemicellulose. Expansins are essential plant proteins that are involved in several key processes, including seed germination, the growth of pollen tubes and root hairs, fruit ripening and abscission processes. Currently, there is a lack of knowledge concerning the role of expansins in woody plants. In this study, we analyzed expansin genes using Populus genome as the study target. Thirty-six members of the expansin gene family were identified in Populus that were divided into four subfamilies (EXPA, EXPB, EXLA and EXLB). We analyzed the molecular structure, chromosome localization, evolutionary relationships and tissue specificity of these genes and investigated expression changes in responses to phytohormone and abiotic stresses of the expansin genes of Populus tremula L. (PtEXs). Molecular structure analysis revealed that each PtEX protein had several conserved motifs and all of the PtEXs genes had multiple exons. Chromosome structure analysis showed that the expansin gene family is distributed on 14 chromosomes. The PtEXs gene family expansion patterns showed segmental duplication. Transcriptome data of Populus revealed that 36 PtEXs genes were differently expressed in different tissues. Cis-element analysis showed that the PtEXs were closely associated with plant development and responses to phytohormone and abiotic stress. Quantitative real-time PCR showed that abscisic acid (ABA) and low-temperature treatment affected the expression of some PtEXs genes, suggesting that these genes are involved in responses to phytohormone and abiotic stress. This study provides a further understanding of the expansin gene family in Populus and forms a basis for future functional research studies.

Function and molecular mechanism of a poplar placenta limited MIXTA gene in regulating differentiation of plant epidermal cells.

The placenta in fruits of most plants either desiccate and shrink as the fruits mature or develop further to form the fleshy tissues. In poplars, placental epidermal cells protrude collectively to produce catkin fibers. In this study, three carpel limited MIXTA genes, PdeMIXTA02, PdeMIXTA03, PdeMIXTA04, were find to specifically expressed in carpel immediately after pollination. Heterologous expression of the three genes in Arabidopsis demonstrated that PdeMIXTA04 significantly promoted trichomes density and could restore trichomes in the trichomeless mutant. By contrast, such functions were not observed with PdeMIXTA02, PdeMIXTA03. In situ hybridization revealed that PdeMIXTA04 was explicitly expressed in poplar placental epidermal cells. We also confirmed trichome-specific expression of the PdeMIXTA04 promoter. Multiple experimental proofs have confirmed the interaction between PdeMIXTA04, PdeMYC and PdeWD40, indicating PdeMIXTA04 functioned through the MYB-bHLH-WD40 ternary complex. Our work provided distinctive understanding of the molecular mechanism triggering differentiation of poplar catkins.

The proposed role of MSL-lncRNAs in causing sex lability of female poplars.

Labile sex expression is frequently observed in dioecious plants, but the underlying genetic mechanism remains largely unknown. Sex plasticity is also observed in many Populus species. Here we carried out a systematic study on a maleness-promoting gene, MSL, detected in the Populus deltoides genome. Our results showed that both strands of MSL contained multiple cis-activating elements, which generated long non-coding RNAs (lncRNAs) promoting maleness. Although female P. deltoides did not have the male-specific MSL gene, a large number of partial sequences with high sequence similarity to this gene were detected in the female poplar genome. Based on sequence alignment, the MSL sequence could be divided into three partial sequences, and heterologous expression of these partial sequences in Arabidopsis confirmed that they could promote maleness. Since activation of the MSL sequences can only result in female sex lability, we propose that MSL-lncRNAs might play a role in causing sex lability of female poplars.

Mycorrhizae Enhance Soybean Plant Growth and Aluminum Stress Tolerance by Shaping the Microbiome Assembly in an Acidic Soil.

Strongly acidic soils are characterized by high aluminum (Al) toxicity and low phosphorus (P) availability, which suppress legume plant growth and nodule development. Arbuscular mycorrhizal fungi (AMF) stimulate rhizobia and enhance plant P uptake. However, it is unclear how this symbiotic soybean-AMF-rhizobial trio promotes soybean growth in acidic soils. We examined the effects of AMF and rhizobium addition on the growth of two soybean genotypes, namely, Al-tolerant and Al-sensitive soybeans as well as their associated bacterial and fungal communities in an acidic soil. With and without rhizobial addition, AMF significantly increased the fresh shoot and root biomass of Al-tolerant soybean by 47%/87% and 37%/24%, respectively. This increase in plant biomass corresponded to the enrichment of four plant growth-promoting rhizobacteria (PGPR) in the rhizospheric soil, namely, Chitinophagaceae bacterium 4GSH07, Paraburkholderia soli, Sinomonas atrocyanea, and Aquincola tertiaricarbonis. For Al-sensitive soybean, AMF addition increased the fresh shoot and root biomass by 112%/64% and 30%/217%, respectively, with/without rhizobial addition. Interestingly, this significant increase coincided with a decrease in the pathogenic fungus Nigrospora oryzae as well as an increase in S. atrocyanea, A. tertiaricarbonis, and Talaromyces verruculosus (a P-solubilizing fungus) in the rhizospheric soil. Lastly, the compartment niche along the soil-plant continuum shaped microbiome assembly, with pathogenic/saprotrophic microbes accumulating in the rhizospheric soil and PGPR related to nitrogen fixation or stress resistance (e.g., Rhizobium leguminosarum and Sphingomonas azotifigens) accumulating in the endospheric layer. IMPORTANCE Taken together, this study examined the effects of arbuscular mycorrhizal fungi (AMF) and rhizobial combinations on the growth of Al-tolerant and Al-sensitive soybeans as well as their associated microbial communities in acidic soils and concluded that AMF enhances soybean growth and Al stress tolerance by recruiting PGPR and altering the root-associated microbiome assembly in a host-dependent manner. In the future, these findings will help us better understand the impacts of AMF on rhizosphere microbiome assembly and will contribute to the development of soybean breeding techniques for the comprehensive use of PGPR in sustainable agriculture.

PKM2/PDK1 dual-targeted shikonin derivatives restore the sensitivity of EGFR-mutated NSCLC cells to gefitinib by remodeling glucose metabolism.

Pyruvate kinase 2 (PKM2) and pyruvate dehydrogenase kinase 1 (PDK1) are two key enzymes in tumor glucose metabolism pathway that not only promote tumor growth and proliferation through accelerating aerobic glycolysis, but also contribute to drug resistance of non-small cell lung cancer (NSCLC). Considering that targeting PKM2 or PDK1 alone seems insufficient to remodel abnormal glucose metabolism to achieve significant antitumor activity, we proposed a "two-step approach" that regulates PKM2 and PDK1 synchronously. Firstly, we found that the combination of ML265 (PKM2 activator) and AZD7545 (PDK1 inhibitor) could synergistically inhibit proliferation and induce apoptosis in H1299 cells. Base on this, we designed a series of novel shikonin (SK) thioether derivatives as PKM2/PDK1 dual-target agents, among which the most potent compound E5 featuring a 2-methyl substitution on the benzene ring exerted significantly increased inhibitory activity toward EGFR mutant NSCLC cell H1975 (IC50 = 1.51 µmol/L), which was 3 and 17-fold more active than the lead compound SK (IC50 = 4.56 µmol/L) and the positive control gefitinib (IC50 = 25.56 µmol/L), respectively. Additionally, E5 also showed good anti-tumor activity in xenografted mouse models, with significantly lower toxicity side effects than SK. Moreover, E5 also inhibited the entry of PKM2 into nucleus to regulate the transcriptional activation of oncogenes, thus restoring the sensitivity of H1975 cell to gefitinib. Collectively, these data demonstrate that E5, a dual inhibitor of PKM2/PDK1, may be a promising adjunct to gefitinib in the treatment of EGFR-TKIs resistant NSCLC, deserving further investigation.

Discrepancies in rhizobacterial assembly caused by glyphosate application and herbicide-tolerant soybean Co-expressing GAT and EPSPS.

There are concerns that the innovation of genetically modified herbicide-tolerant (GMHT) plants, as well as the application of herbicide to such GMHT plants, could have an impact on ecological interactions and unintentionally harm non-targeted organisms. Consequently, we intend to use full-length 16 S rDNA amplicon sequencing to examine changes in the bacterial community in the rhizosphere of GMHT soybean (Z106) harboring 5-enolpyruvylshikimate-3-phosphate synthase and Glyphosate N-acetyltransferase genes and GMHT soybean treated with glyphosate (Z106G). Glyphosate application significantly impacted bacterial alpha diversity (species richness, and Shannon diversity). Permutational multivariate analysis of variance of beta diversity demonstrated that soil compartments and growth stages had a substantial impact on soybean rhizobacterial communities (soil compartments, growth stages, P = 0.001). Community composition revealed that Z106G soils were abundant in Taibaiella and Arthrobacter pascens at maturity, while Chryseobacterium joostei and Stenotrophomonas maltophilia predominated in Z106 soils during flowering. Nitrogen-fixing and phosphate-solubilizing microbes were found in higher proportions in the rhizosphere than in bulk soil, with Sinorhizobium being more abundant in Z106 and Bacillus and Stenotrophomonas being more prevalent in Z106G rhizosphere soils. Collectively, our findings suggest glyphosate application and glyphosate-tolerant soybean as potential regulators of soybean rhizobacterial composition.

Genome-Wide Comparative Analysis of the Fasciclin-like Arabinogalactan Proteins (FLAs) in Salicacea and Identification of Secondary Tissue Development-Related Genes.

Fasciclin-like arabinogalactan proteins (FLAs) are a subclass of arabinogalactan proteins (AGPs) containing both AGP-like glycated domains and fasciclin (FAS) domains, which are involved in plant growth and development and synthesis of the cell wall. However, these proteins have not been identified or analyzed in willow, Salix, the sister genus of Populus. In this study, we performed a whole genome study of the FLA gene family of Salix suchowensis and compared it with the FLA gene family of Populus deltoides. The results showed the presence of 40 and 46 FLA genes in P. deltoides and S. suchowensis, distributed on 17 and 16 chromosomes, respectively. Four pairs of tandem repeat genes were found in willow, while poplar had no tandem repeat genes. Twelve and thirteen pairs of duplicated gene fragments were identified in poplar and willow, respectively. The multispecies phylogenetic tree showed that the FLA gene family could be divided into four groups (I-IV), with Group 1 showing significant expansion in woody plants. A gene expression analysis showed that PdeFLA19/27 in Group I of poplar was highly expressed, specifically during the secondary growth period of the stem and the rapid elongation of seed hairs. In the Group I genes of S. suchowensis, SsuFLA25/26/28 was also highly expressed during the secondary growth period, whereas increased expression of SsuFLA35 was associated with seed hair tissue. These results provide important clues about the differences in the FLA gene family during the evolution of herbs and woody plants, and suggest that the FLA gene family may play an essential role in regulating the secondary growth of woody plants. It also provides a reference for further studies on the regulation of secondary growth and seed hair development by FLA genes in poplar and willow.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis.

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

The host niches of soybean rather than genetic modification or glyphosate application drive the assembly of root-associated microbial communities.

Plant roots significantly influence soil microbial diversity, and soil microorganisms play significant roles in both natural and agricultural ecosystems. Although the genetically modified (GM) crops with enhanced insect and herbicide resistance are thought to have unmatched yield and stress resistance advantages, thorough and in-depth case studies still need to be carried out in a real-world setting due to the potential effects of GM plants on soil microbial communities. In this study, three treatments were used: a recipient soybean variety Jack, a triple transgenic soybean line JD321, and the glyphosate-treated JD321 (JD321G). Three sampling stages (flowering, seed filling and maturing), as well as three host niches of soybean rhizosphere [intact roots (RT), rhizospheric soil (RS) and surrounding soil (SS)] were established. In comparison to Jack, the rhizospheric soil of JD321G had higher urease activity and lower nitrite reductase at the flowering stage. Different treatments and different sampling stages existed no significant effects on the compositions of microbial communities at different taxonomic levels. However, at the genus level, the relative abundance of three plant growth-promoting fungal genera (i.e. Mortierella, Chaetomium and Pseudombrophila) increased while endophytic bacteria Chryseobacterium and pathogenic bacteria Streptomyces decreased from the inside to the outside of the roots (i.e. RT → RS → SS). Moreover, two bacterial genera, Bradyrhizobium and Ensifer were more abundant in RT than in RS and SS, as well as three species, Agrobacterium radiobacter, Ensifer fredii and Ensifer meliloti, which are closely related to nitrogen-fixation. Furthermore, five clusters of orthologous groups (COGs) associated to nitrogen-fixation genes were higher in RT than in RS, whereas only one COG annotated as dinitrogenase iron-molybdenum cofactor biosynthesis protein was lower. Overall, the results imply that the rhizosphere host niches throughout the soil-plant continuum largely control the composition and function of the root-associated microbiome of triple transgenic soybean.

Genome-Wide Identification of LeBAHDs in Lithospermum erythrorhizon and In Vivo Transgenic Studies Confirm the Critical Roles of LeBAHD1/LeSAT1 in the Conversion of Shikonin to Acetylshikonin.

The BAHD acyltransferase family is a unique class of plant proteins that acylates plant metabolites and participates in plant secondary metabolic processes. However, the BAHD members in Lithospermum erythrorhizon remain unknown and uncharacterized. Although the heterologously expressed L. erythrorhizon BAHD family member LeSAT1 in Escherichia coli has been shown to catalyze the conversion of shikonin to acetylshikonin in vitro, its in vivo role remains unknown. In this study, the characterization, evolution, expression patterns, and gene function of LeBAHDs in L. erythrorhizon were explored by bioinformatics and transgenic analysis. We totally identified 73 LeBAHDs in the reference genome of L. erythrorhizon. All LeBAHDs were phylogenetically classified into five clades likely to perform different functions, and were mainly expanded by dispersed and WGD/segmental duplication. The in vivo functional investigation of the key member LeBAHD1/LeSAT1 revealed that overexpression of LeBAHD1 in hairy roots significantly increased the content of acetylshikonin as well as the conversion rate of shikonin to acetylshikonin, whereas the CRISPR/Cas9-based knockout of LeBAHD1 in hairy roots displayed the opposite trend. Our results not only confirm the in vivo function of LeBAHD1/LeSAT1 in the biosynthesis of acetylshikonin, but also provide new insights for the biosynthetic pathway of shikonin and its derivatives.

Novel shikonin derivatives suppress cell proliferation, migration and induce apoptosis in human triple-negative breast cancer cells via regulating PDK1/PDHC axis.

AIMS: PDK1 is one of the key enzymes in the glucose metabolism pathway, which is abnormally high expressed in breast cancer tissues and can promote tumor proliferation and metastasis. PDK1 and the PDHC/PDK axis are important targets for regulating glucose metabolism and anti-tumor activity. In this study, we evaluated the anti-tumor activities of a series of semi-synthesized shikonin (SK) derivatives against human breast cancer cells. MAIN METHODS: The anti-proliferation activity of SK derivatives against human breast cancer cell lines was tested by CCK-8 and EdU assay. Flow cytometry was utilized to evaluate cell apoptosis, reactive oxygen species and cell cycle distribution. Cell migration ability was determined by wound healing and trans-well assay. PDK1 targeting effect was confirmed by western bolting, molecular docking, bio-layer interferometry and PDK1 enzyme activity assay. Nude-mouse transplanted tumor model was used to evaluate their anti-tumor effect in vivo. KEY FINDINGS: Findings revealed that SK derivatives had good anti-proliferation ability against MDA-MB-231 cell. They induced cell apoptosis by regulating the mitochondrial apoptosis and death receptor pathway. They also inhibited cell migration by suppressing EMT progression. Molecular docking, PDK1 affinity and enzyme activity demonstrated their PDK1 targeting. In vivo antitumor experiment showed that E2 could significantly inhibit tumor growth with lower side-effect on mice than SK. SIGNIFICANCE: In conclusion, the novel SK derivatives E2 and E5 inhibited tumor glycolysis by targeting PDK1 and ultimately induced apoptosis. Our data demonstrated that E2 would be a good lead compound for the treatment of human TNBC as a novel PDK1 inhibitor.