Unveiling orphan receptor-like kinases in plants: novel client discovery using high-confidence library predictions in the Kinase-Client (KiC) assay.

Plants are remarkable in their ability to adapt to changing environments, with receptor-like kinases (RLKs) playing a pivotal role in perceiving and transmitting environmental cues into cellular responses. Despite extensive research on RLKs from the plant kingdom, the function and activity of many kinases, i.e., their substrates or "clients", remain uncharted. To validate a novel client prediction workflow and learn more about an important RLK, this study focuses on P2K1 (DORN1), which acts as a receptor for extracellular ATP (eATP), playing a crucial role in plant stress resistance and immunity. We designed a Kinase-Client (KiC) assay library of 225 synthetic peptides, incorporating previously identified P2K phosphorylated peptides and novel predictions from a deep-learning phosphorylation site prediction model (MUsite) and a trained hidden Markov model (HMM) based tool, HMMER. Screening the library against purified P2K1 cytosolic domain (CD), we identified 46 putative substrates, including 34 novel clients, 27 of which may be novel peptides, not previously identified experimentally. Gene Ontology (GO) analysis among phosphopeptide candidates revealed proteins associated with important biological processes in metabolism, structure development, and response to stress, as well as molecular functions of kinase activity, catalytic activity, and transferase activity. We offer selection criteria for efficient further in vivo experiments to confirm these discoveries. This approach not only expands our knowledge of P2K1's substrates and functions but also highlights effective prediction algorithms for identifying additional potential substrates. Overall, the results support use of the KiC assay as a valuable tool in unraveling the complexities of plant phosphorylation and provide a foundation for predicting the phosphorylation landscape of plant species based on peptide library results.

The Auto-Regulation of ATL2 E3 Ubiquitin Ligase Plays an Important Role in the Immune Response against Alternaria brassicicola in Arabidopsis thaliana.

The ubiquitin/26S proteasome system is a crucial regulatory mechanism that governs various cellular processes in plants, including signal transduction, transcriptional regulation, and responses to biotic and abiotic stressors. Our study shows that the RING-H2-type E3 ubiquitin ligase, Arabidopsis Tóxicos en Levadura 2 (ATL2), is involved in response to fungal pathogen infection. Under normal growth conditions, the expression of the ATL2 gene is low, but it is rapidly and significantly induced by exogenous chitin. Additionally, ATL2 protein stability is markedly increased via chitin treatment, and its degradation is prolonged when 26S proteasomal function is inhibited. We found that an atl2 null mutant exhibited higher susceptibility to Alternaria brassicicola, while plants overexpressing ATL2 displayed increased resistance. We also observed that the hyphae of A. brassicicola were strongly stained with trypan blue staining, and the expression of A. brassicicola Cutinase A (AbCutA) was dramatically increased in atl2. In contrast, the hyphae were weakly stained, and AbCutA expression was significantly reduced in ATL2-overexpressing plants. Using bioinformatics, live-cell confocal imaging, and cell fractionation analysis, we revealed that ATL2 is localized to the plasma membrane. Further, it is demonstrated that the ATL2 protein possesses E3 ubiquitin ligase activity and found that cysteine 138 residue is critical for its function. Moreover, ATL2 is necessary to successfully defend against the A. brassicicola fungal pathogen. Altogether, our data suggest that ATL2 is a plasma membrane-integrated protein with RING-H2-type E3 ubiquitin ligase activity and is essential for the defense response against fungal pathogens in Arabidopsis.

Involvement of Peptidoglycan Receptor Proteins in Mediating the Growth-Promoting Effects of Bacillus pumilus TUAT1 in Arabidopsis thaliana.

Bacillus pumilus TUAT1 acts as plant growth-promoting rhizobacteria (PGPR) for various plants like rice and Arabidopsis. Under stress conditions, B. pumilus TUAT1 forms spores with a thick peptidoglycan (PGN) cell wall. Previous research showed that spores were significantly more effective than vegetative cells in enhancing plant growth. In Arabidopsis, the lysin-motif proteins LYM1, LYM3, and CERK1 are required for recognizing bacterial peptidoglycans (PGNs) to mediate immunity. Here, we examined the involvement of PGN receptor proteins in the PGP effects of B. pumilus TUAT1 using Arabidopsis mutants defective in PGNs receptors. Root growth of wild-type, cerk1-1, lym1-1 and lym1-2 mutant plants was significantly increased by TUAT1 inoculation, but this was not the case for lym3-1 and lym3-2 mutant plants. RNA-seq analysis revealed that the expression of a number of defense-related genes was upregulated in lym3 mutant plants. These results suggested that B. pumilus TUAT1 may act to reduce the defense response, which is dependent on a functional LYM3. The expression of the defense-responsive gene, WRKY29, was significantly induced by the elicitor flg-22, both in wild-type and lym3 mutant plants, while this induction was significantly reduced by treatment with B. pumilus TUAT1 and PGNs in wild-type, but not in lym3 mutant plants. These findings suggest that the PGNs of B. pumilus TUAT1 may be recognized by the LYM3 receptor protein, suppressing the defense response, which results in plant growth promotion in a trade-off between defense and growth.

High-Throughput f-LAESI-IMS-MS for Mapping Biological Nitrogen Fixation One Cell at a Time.

For the characterization of the metabolic heterogeneity of cell populations, high-throughput single-cell analysis platforms are needed. In this study, we utilized mass spectrometry (MS) enhanced with ion mobility separation (IMS) and coupled with an automated sampling platform, fiber-based laser ablation electrospray ionization (f-LAESI), for in situ high-throughput single-cell metabolomics in soybean (Glycine max) root nodules. By fully automating the in situ sampling platform, an overall sampling rate of 804 cells/h was achieved for high numbers (>500) of tissue-embedded plant cells. This is an improvement by a factor of 13 compared to the previous f-LAESI-MS configuration. By introducing IMS, the molecular coverage improved, and structural isomers were separated on a millisecond time scale. The enhanced f-LAESI-IMS-MS platform produced 259 sample-related peaks/cell, almost twice as much as the 131 sample-related peaks/cell produced by f-LAESI-MS without IMS. Using the upgraded system, two types of metabolic heterogeneity characterization methods became possible. For unimodal metabolite abundance distributions, the metabolic noise reported on the metabolite level variations within the cell population. For bimodal distributions, the presence of metabolically distinct subpopulations was established. Discovering these latent cellular phenotypes could be linked to the presence of different cell states, e.g., proliferating bacteria in partially occupied plant cells and quiescent bacteroids in fully occupied cells in biological nitrogen fixation, or spatial heterogeneity due to altered local environments.

Enhancing powdery mildew resistance in soybean by targeted mutation of MLO genes using the CRISPR/Cas9 system.

BACKGROUND: Powdery mildew is a major disease that causes great losses in soybean yield and seed quality. Disease-resistant varieties, which are generated by reducing the impact of susceptibility genes through mutation in host plants, would be an effective approach to protect crops from this disease. The Mildew Locus O (MLO) genes are well-known susceptibility genes for powdery mildew in plant. In this study, we utilized the CRISPR/Cas9 system to induce targeted mutations in the soybean GmMLO genes to improve powdery mildew resistance. RESULTS: A dual-sgRNA CRISPR/Cas9 construct was designed and successfully transferred into the Vietnamese soybean cultivar DT26 through Agrobacterium tumefaciens-mediated transformation. Various mutant forms of the GmMLO genes including biallelic, chimeric and homozygous were found at the T0 generation. The inheritance and segregation of CRISPR/Cas9-induced mutations were confirmed and validated at the T1 and T2 generations. Out of six GmMLO genes in the soybean genome, we obtained the Gmmlo02/Gmmlo19/Gmmlo23 triple and Gmmlo02/Gmmlo19/Gmmlo20/Gmmlo23 quadruple knockout mutants at the T2 generation. When challenged with Erysiphe diffusa, a fungus that causes soybean powdery mildew, all mutant plants showed enhanced resistance to the pathogen, especially the quadruple mutant. The powdery mildew severity in the mutant soybeans was reduced by up to 36.4% compared to wild-type plants. In addition, no pleiotropic effect on soybean growth and development under net-house conditions was observed in the CRISPR/Cas9 mutants. CONCLUSIONS: Our results indicate the involvement of GmMLO02, GmMLO19, GmMLO20 and GmMLO23 genes in powdery mildew susceptibility in soybean. Further research should be conducted to investigate the roles of individual tested genes and the involvement of other GmMLO genes in this disease infection mechanism. Importantly, utilizing the CRISPR/Cas9 system successfully created the Gmmlo transgene-free homozygous mutant lines with enhanced resistance to powdery mildew, which could be potential materials for soybean breeding programs.

Phosphorylation-mediated regulation of integrin-linked kinase 5 by purinoreceptor P2K2.

Extracellular ATP (eATP) in plants plays a crucial role as a ligand for purinoreceptors, mediating purinergic signaling and regulating diverse biological functions, including responses to abiotic and biotic stresses. DORN1/P2K1 (LecRK I.9) was the first identified plant purinoreceptor. P2K2 (LecRK I.5) was subsequently identified as an additional plant purinoreceptor and shown to directly interact with P2K1. Recently, we reported that P2K1 interacts with Integrin-linked kinase 5 (ILK5), a Raf-like MAPKKK protein, and phosphorylates ILK5 to regulate purinergic signaling in relation to plant innate immunity. Here, we report that P2K2 also interacts with the ILK5 protein in planta. Furthermore, we demonstrate that P2K2 phosphorylates ILK5 in the presence of [γ-32P] ATP, similar to P2K1. However, unlike P2K1, P2K2 exhibits strong phosphorylation even when the Serine 192 residue of ILK5 is mutated to Alanine (ILK5S192A), suggesting the possibility of phosphorylation of other residues to fully regulate ILK5 protein function.

Acetylene Reduction Assay: A Measure of Nitrogenase Activity in Plants and Bacteria.

Nitrogen is one of the most abundant elements in the biosphere, but its gaseous form is not biologically available to many organisms, including plants and animals. Diazotrophic microorganisms can convert atmospheric nitrogen into ammonia, a form that can be absorbed by plants in a process called biological nitrogen fixation (BNF). BNF is catalyzed by the enzyme nitrogenase, which not only reduces N2 to NH3 , but also reduces other substrates such as acetylene. The acetylene reduction assay (ARA) can be used to measure nitrogenase activity in diazotrophic organisms, either in symbiotic associations or in their free-living state. The technique uses gas chromatography to measure the reduction of acetylene to ethylene by nitrogenase in a simple, quick, and inexpensive manner. Here, we demonstrate how to: prepare nodulated soybean plants and culture free-living Azospirillum brasilense for the ARA, use the gas chromatograph to detect the ethylene formed, and calculate the nitrogenase activity based on the peaks generated by the chromatograph. The methods shown here using example organisms can be easily adapted to other nodulating plants and diazotrophic bacteria. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Acetylene reduction assay in root nodules Basic Protocol 2: Acetylene reduction assay using diazotrophic bacteria Basic Protocol 3: Calculation of nitrogenase activity Support Protocol 1: Production of acetylene from calcium carbide Support Protocol 2: Calibration of the gas chromatograph Support Protocol 3: Total protein quantification.

The Raf-like MAPKKK INTEGRIN-LINKED KINASE 5 regulates purinergic receptor-mediated innate immunity in Arabidopsis.

Mitogen-activated protein (MAP) kinase signaling cascades play important roles in eukaryotic defense against various pathogens. Activation of the extracellular ATP (eATP) receptor P2K1 triggers MAP kinase 3 and 6 (MPK3/6) phosphorylation, which leads to an elevated plant defense response. However, the mechanism by which P2K1 activates the MAPK cascade is unclear. In this study, we show that in Arabidopsis thaliana, P2K1 phosphorylates the Raf-like MAP kinase kinase kinase (MAPKKK) INTEGRIN-LINKED KINASE 5 (ILK5) on serine 192 in the presence of eATP. The interaction between P2K1 and ILK5 was confirmed both in vitro and in planta and their interaction was enhanced by ATP treatment. Similar to P2K1 expression, ILK5 expression levels were highly induced by treatment with ATP, flg22, Pseudomonas syringae pv. tomato DC3000, and various abiotic stresses. ILK5 interacts with and phosphorylates the MAP kinase MKK5. Moreover, phosphorylation of MPK3/6 was significantly reduced upon ATP treatment in ilk5 mutant plants, relative to wild-type (WT). The ilk5 mutant plants showed higher susceptibility to P. syringae pathogen infection relative to WT plants. Plants expressing only the mutant ILK5S192A protein, with decreased kinase activity, did not activate the MAPK cascade upon ATP addition. These results suggest that eATP activation of P2K1 results in transphosphorylation of the Raf-like MAPKKK ILK5, which subsequently triggers the MAPK cascade, culminating in activation of MPK3/6 associated with an elevated innate immune response.

Super-Resolution Imaging of Plant Receptor-Like Kinases Uncovers Their Colocalization and Coordination with Nanometer Resolution.

Plant cell signaling often relies on the cellular organization of receptor-like kinases (RLKs) within membrane nanodomains to enhance signaling specificity and efficiency. Thus, nanometer-scale quantitative analysis of spatial organizations of RLKs could provide new understanding of mechanisms underlying plant responses to environmental stress. Here, we used stochastic optical reconstruction fluorescence microscopy (STORM) to quantify the colocalization of the flagellin-sensitive-2 (FLS2) receptor and the nanodomain marker, remorin, within Arabidopsis thaliana root hair cells. We found that recovery of FLS2 and remorin in the plasma membrane, following ligand-induced internalization by bacterial-flagellin-peptide (flg22), reached ~85% of their original membrane density after ~90 min. The pairs colocalized at the membrane at greater frequencies, compared with simulated randomly distributed pairs, except for directly after recovery, suggesting initial uncoordinated recovery followed by remorin and FLS2 pairing in the membrane. The purinergic receptor, P2K1, colocalized with remorin at similar frequencies as FLS2, while FLS2 and P2K1 colocalization occurred at significantly lower frequencies, suggesting that these RLKs mostly occupy distinct nanodomains. The chitin elicitor receptor, CERK1, colocalized with FLS2 and remorin at much lower frequencies, suggesting little coordination between CERK1 and FLS2. These findings emphasize STORM's capacity to observe distinct nanodomains and degrees of coordination between plant cell receptors, and their respective immune pathways.

Mapping Genetic Variation in Arabidopsis in Response to Plant Growth-Promoting Bacterium Azoarcus olearius DQS-4T.

Plant growth-promoting bacteria (PGPB) can enhance plant health by facilitating nutrient uptake, nitrogen fixation, protection from pathogens, stress tolerance and/or boosting plant productivity. The genetic determinants that drive the plant-bacteria association remain understudied. To identify genetic loci highly correlated with traits responsive to PGPB, we performed a genome-wide association study (GWAS) using an Arabidopsis thaliana population treated with Azoarcus olearius DQS-4T. Phenotypically, the 305 Arabidopsis accessions tested responded differently to bacterial treatment by improving, inhibiting, or not affecting root system or shoot traits. GWA mapping analysis identified several predicted loci associated with primary root length or root fresh weight. Two statistical analyses were performed to narrow down potential gene candidates followed by haplotype block analysis, resulting in the identification of 11 loci associated with the responsiveness of Arabidopsis root fresh weight to bacterial inoculation. Our results showed considerable variation in the ability of plants to respond to inoculation by A. olearius DQS-4T while revealing considerable complexity regarding statistically associated loci with the growth traits measured. This investigation is a promising starting point for sustainable breeding strategies for future cropping practices that may employ beneficial microbes and/or modifications of the root microbiome.

Computational prediction and in vitro analysis of the potential ligand binding site within the extracellular ATP receptor, P2K2.

The plant extracellular ATP (eATP) receptor, P2K2, binds eATP with strong ligand affinity through its extracellular lectin domain. Ligand binding activates the intracellular kinase domain of P2K2 resulting in a variety of intracellular responses and, ultimately, increased plant immunity to invading fungal and bacterial pathogens. Here, using a computational prediction approach, we developed a tertiary structure model of the P2K2 extracellular lectin domain. In silico target docking of ATP to the P2K2-binding site predicted interaction with several residues through hydrophobic interactions and hydrogen bonding. Our confirmation of the modeling was obtained by showing that H99, R144, and S256 are key residues essential for in vitro binding of ATP by P2K2.

Critical role for uricase and xanthine dehydrogenase in soybean nitrogen fixation and nodule development.

De novo purine biosynthesis is required for the incorporation of fixed nitrogen in ureide exporting nodules, as formed on soybean [Glycine max (L.) Merr.] roots. However, in many cases, the enzymes involved in this pathway have been deduced strictly from genome annotations with little direct genetic evidence, such as mutant studies, to confirm their biochemical function or importance to nodule development. While efforts to develop large mutant collections of soybean are underway, research on this plant is still hampered by the inability to obtain mutations in any specific gene of interest. Using a forward genetic approach, as well as CRISPR/Cas9 gene editing via Agrobacterium rhizogenes-mediated hairy root transformation, we identified and characterized the role of GmUOX (Uricase) and GmXDH (Xanthine Dehydrogenase) in nitrogen fixation and nodule development in soybean. The gmuox knockout soybean mutants displayed nitrogen deficiency chlorosis and early nodule senescence, as exemplified by the reduced nitrogenase (acetylene reduction) activity in nodules, the internal greenish-white internal appearance of nodules, and diminished leghemoglobin production. In addition, gmuox1 nodules showed collapsed infected cells with degraded cytoplasm, aggregated bacteroids with no discernable symbiosome membranes, and increased formation of poly-ß-hydroxybutyrate granules. Similarly, knockout gmxdh mutant nodules, generated with the CRISPR/Cas9 system, also exhibited early nodule senescence. These genetic studies confirm the critical role of the de novo purine metabolisms pathway not only in the incorporation of fixed nitrogen but also in the successful development of a functional, nitrogen-fixing nodule. Furthermore, these studies demonstrate the great utility of the CRISPR/Cas9 system for studying root-associated gene traits when coupled with hairy root transformation.

Near Infrared Emitting Semiconductor Polymer Dots for Bioimaging and Sensing.

Semiconducting polymer dots (Pdots) are rapidly becoming one of the most studied nanoparticles in fluorescence bioimaging and sensing. Their small size, high brightness, and resistance to photobleaching make them one of the most attractive fluorophores for fluorescence imaging and sensing applications. This paper highlights our recent advances in fluorescence bioimaging and sensing with nanoscale luminescent Pdots, specifically the use of organic dyes as dopant molecules to modify the optical properties of Pdots to enable deep red and near infrared fluorescence bioimaging applications and to impart sensitivity of dye doped Pdots towards selected analytes. Building on our earlier work, we report the formation of secondary antibody-conjugated Pdots and provide Cryo-TEM evidence for their formation. We demonstrate the selective targeting of the antibody-conjugated Pdots to FLAG-tagged FLS2 membrane receptors in genetically engineered plant leaf cells. We also report the formation of a new class of luminescent Pdots with emission wavelengths of around 1000 nm. Finally, we demonstrate the formation and utility of oxygen sensing Pdots in aqueous media.

Spatial Mapping of Plant N-Glycosylation Cellular Heterogeneity Inside Soybean Root Nodules Provided Insights Into Legume-Rhizobia Symbiosis.

Although ubiquitously present, information on the function of complex N-glycan posttranslational modification in plants is very limited and is often neglected. In this work, we adopted an enzyme-assisted matrix-assisted laser desorption/ionization mass spectrometry imaging strategy to visualize the distribution and identity of N-glycans in soybean root nodules at a cellular resolution. We additionally performed proteomics analysis to probe the potential correlation to proteome changes during symbiotic rhizobia-legume interactions. Our ion images reveal that intense N-glycosylation occurs in the sclerenchyma layer, and inside the infected cells within the infection zone, while morphological structures such as the cortex, uninfected cells, and cells that form the attachment with the root are fewer N-glycosylated. Notably, we observed different N-glycan profiles between soybean root nodules infected with wild-type rhizobia and those infected with mutant rhizobia incapable of efficiently fixing atmospheric nitrogen. The majority of complex N-glycan structures, particularly those with characteristic Lewis-a epitopes, are more abundant in the mutant nodules. Our proteomic results revealed that these glycans likely originated from proteins that maintain the redox balance crucial for proper nitrogen fixation, but also from enzymes involved in N-glycan and phenylpropanoid biosynthesis. These findings indicate the possible involvement of Lewis-a glycans in these critical pathways during legume-rhizobia symbiosis.

High-Resolution Translatome Analysis Reveals Cortical Cell Programs During Early Soybean Nodulation.

Nodule organogenesis in legumes is regulated temporally and spatially through gene networks. Genome-wide transcriptome, proteomic, and metabolomic analyses have been used previously to define the functional role of various plant genes in the nodulation process. However, while significant progress has been made, most of these studies have suffered from tissue dilution since only a few cells/root regions respond to rhizobial infection, with much of the root non-responsive. To partially overcome this issue, we adopted translating ribosome affinity purification (TRAP) to specifically monitor the response of the root cortex to rhizobial inoculation using a cortex-specific promoter. While previous studies have largely focused on the plant response within the root epidermis (e.g., root hairs) or within developing nodules, much less is known about the early responses within the root cortex, such as in relation to the development of the nodule primordium or growth of the infection thread. We focused on identifying genes specifically regulated during early nodule organogenesis using roots inoculated with Bradyrhizobium japonicum. A number of novel nodulation gene candidates were discovered, as well as soybean orthologs of nodulation genes previously reported in other legumes. The differential cortex expression of several genes was confirmed using a promoter-GUS analysis, and RNAi was used to investigate gene function. Notably, a number of differentially regulated genes involved in phytohormone signaling, including auxin, cytokinin, and gibberellic acid (GA), were also discovered, providing deep insight into phytohormone signaling during early nodule development.

Hydroporphyrin-Doped Near-Infrared-Emitting Polymer Dots for Cellular Fluorescence Imaging.

Near-infrared (NIR) fluorescent semiconductor polymer dots (Pdots) have shown great potential for fluorescence imaging due to their exceptional chemical and photophysical properties. This paper describes the synthesis of NIR-emitting Pdots with great control and tunability of emission peak wavelength. The Pdots were prepared by doping poly[(9,9-dioctylfluorenyl-2,7-diyl)-alt-co-(1,4-benzo-(2,1',3)-thiadiazole)] (PFBT), a semiconducting polymer commonly used as a host polymer in luminescent Pdots, with a series of chlorins and bacteriochlorins with varying functional groups. Chlorins and bacteriochlorins are ideal dopants due to their high hydrophobicity, which precludes their use as molecular probes in aqueous biological media but on the other hand prevents their leakage when doped into Pdots. Additionally, chlorins and bacteriochlorins have narrow deep red to NIR-emission bands and the wide array of synthetic modifications available for modifying their molecular structure enables tuning their emission predictably and systematically. Transmission electron microscopy (TEM) and dynamic light scattering (DLS) measurements show the chlorin- and bacteriochlorin-doped Pdots to be nearly spherical with an average diameter of 46 ± 12 nm. Efficient energy transfer between PFBT and the doped chlorins or bacteriochlorins decreases the PFBT donor emission to near baseline level and increases the emission of the doped dyes that serve as acceptors. The chlorin- and bacteriochlorin-doped Pdots show narrow emission bands ranging from 640 to 820 nm depending on the doped dye. The paper demonstrates the utility of the systematic chlorin and bacteriochlorin synthesis approach by preparing Pdots of varying emission peak wavelength, utilizing them to visualize multiple targets using wide-field fluorescence microscopy, binding them to secondary antibodies, and determining the binding of secondary antibody-conjugated Pdots to primary antibody-labeled receptors in plant cells. Additionally, the chlorin- and bacteriochlorin-doped Pdots show a blinking behavior that could enable their use in super-resolution imaging methods like STORM.

A Multi-Level Iterative Bi-Clustering Method for Discovering miRNA Co-regulation Network of Abiotic Stress Tolerance in Soybeans.

Although growing evidence shows that microRNA (miRNA) regulates plant growth and development, miRNA regulatory networks in plants are not well understood. Current experimental studies cannot characterize miRNA regulatory networks on a large scale. This information gap provides an excellent opportunity to employ computational methods for global analysis and generate valuable models and hypotheses. To address this opportunity, we collected miRNA-target interactions (MTIs) and used MTIs from Arabidopsis thaliana and Medicago truncatula to predict homologous MTIs in soybeans, resulting in 80,235 soybean MTIs in total. A multi-level iterative bi-clustering method was developed to identify 483 soybean miRNA-target regulatory modules (MTRMs). Furthermore, we collected soybean miRNA expression data and corresponding gene expression data in response to abiotic stresses. By clustering these data, 37 MTRMs related to abiotic stresses were identified, including stress-specific MTRMs and shared MTRMs. These MTRMs have gene ontology (GO) enrichment in resistance response, iron transport, positive growth regulation, etc. Our study predicts soybean MTRMs and miRNA-GO networks under different stresses, and provides miRNA targeting hypotheses for experimental analyses. The method can be applied to other biological processes and other plants to elucidate miRNA co-regulation mechanisms.

Large-Scale Integrative Analysis of Soybean Transcriptome Using an Unsupervised Autoencoder Model.

Plant tissues are distinguished by their gene expression patterns, which can help identify tissue-specific highly expressed genes and their differential functional modules. For this purpose, large-scale soybean transcriptome samples were collected and processed starting from raw sequencing reads in a uniform analysis pipeline. To address the gene expression heterogeneity in different tissues, we utilized an adversarial deconfounding autoencoder (AD-AE) model to map gene expressions into a latent space and adapted a standard unsupervised autoencoder (AE) model to help effectively extract meaningful biological signals from the noisy data. As a result, four groups of 1,743, 914, 2,107, and 1,451 genes were found highly expressed specifically in leaf, root, seed and nodule tissues, respectively. To obtain key transcription factors (TFs), hub genes and their functional modules in each tissue, we constructed tissue-specific gene regulatory networks (GRNs), and differential correlation networks by using corrected and compressed gene expression data. We validated our results from the literature and gene enrichment analysis, which confirmed many identified tissue-specific genes. Our study represents the largest gene expression analysis in soybean tissues to date. It provides valuable targets for tissue-specific research and helps uncover broader biological patterns. Code is publicly available with open source at https://github.com/LingtaoSu/SoyMeta.