protein-DNA Complexes Research Articles

EMSA-SELEX-seq method for analysis of binding site sequences in DNA-protein complexes

The BOB1 protein (OBF1, OCA-B) is a transcriptional coactivator of two POU domain proteins — OCT1, expressed in all cells, and lymphoid-specific OCT2. The interaction of BOB1 with OCT1/2 plays an important role in the regulation of immune responses in both physiological and pathological contexts. BOB1 is known to form a ternary complex with OCT1/2 bound to DNA in monomeric and certain dimeric configurations, changing the sequence specificity of the binding. To analyze DNA sequences from these complexes, in this work we proposed the EMSA-SELEX-seq method, based on the separation of OCT/BOB1 complexes of various compositions in a non-denaturing polyacrylamide gel (EMSA) followed by the isolation and amplification of the oligonucleotides that they contain (SELEX). Based on several rounds of the enrichment followed by the NGS sequencing and bioinformatics analysis, the DNA sequences were determined and the relevance of this approach was confirmed. Thus, the proposed EMSA-SELEX-seq method allows the analysis of DNA sequences in DNA-protein complexes with varying dimensions of its protein components.

Comprehensive in silico analysis of single nucleotide polymorphism and molecular dynamics simulation of human GATA6 protein in ventricular septal defect

Congenital heart disease (CHD) represents nearly one-third of congenital birth defects annually, with ventricular septal defect (VSD) being the most common type. The aim of this study was to explore the role of specific GATA binding protein 6 gene (GATA6) mutations as a potential etiological factor in the development of VSD through an in silico approach. Data were collected from the human gene databases: DisGeNET and GeneCards, with protein-protein interaction networks constructed via STRING and Cytoscape. Gene ontology and pathway enrichment analyses were conducted using DAVID, with data analysis in R with significance set at FDR p<0.05. Target single nucleotide polymorphisms (SNPs) of GATA6 were obtained from NCBI dbSNP, and non-synonymous single nucleotide polymorphism (nsSNP) effects were predicted using SIFT, PolyPhen-2, I-Mutant 2.0, Fathmm, MutPred 2.0, SNP&GO, and PON-P2. Conserved regions of GATA6 were analyzed using ConSurf, with functional classification, variant conservation, and stability changes evaluated in Google Colab. Multiple sequence alignment was performed using ClustalW. Mutation modeling and molecular dynamics analysis, using GROMACS, revealed that among 87 intersecting genes, 16 proteins were interconnected with GATA6, showing a centrality value of 0.4378. Gene ontology analysis highlighted atrioventricular canal development, protein-DNA complexes, and transcription factor regulation as key processes for cardiac development, especially in the ventricular septum. NsSNP and molecular dynamics analyses identified rs387906818 and rs387906820 as having the highest pathogenic potential for VSD due to amino acid structural changes.

Upstream regulation of microRNA-9 through a complex cellular machinery during neurogenesis

SummaryWhile microRNAs (miRs) like miR-9 are crucial for neurogenesis and neuronal differentiation, their regulatory mechanisms are not well understood. miR-9 is highly expressed in the brain and plays a significant role in neurogenesis. Using the Collaborative Cross resource, we identified significant quantitative trait loci (QTL) through genetic analyses. We then characterized over 130 candidate genes within these QTL regions using RNA interference, qPCR, and neuronal differentiation assays, narrowing them down to 13 promising candidates. Among these, Panx2, Polr1c, and Mgea5 were found to colocalize in the neurogenic niches of the SVZ and DG regions, as shown by immunofluorescence. Further ChIP-seq and Co-IP analyses revealed their interaction and binding to the miR-9 locus, forming a DNA-protein regulatory complex we termed ’miRSome-9.’ A 3C/ChIP-loop assay confirmed the chromatin organization of miRSome-9 at the miR-9 locus, shedding light on the upstream mechanisms regulating miR-9 expression during neurogenesis.

Expression and Purification of pAG-MNase for CUT&RUN.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a chromatin profiling strategy that releases specific DNA-protein complexes from in situ chromatin to the supernatant. CUT&RUN employs the hybrid protein pAG-MNase, composed of IgG binding proteins A and G, followed by micrococcal nuclease. Here, I describe the expression and purification of recombinant pAG-MNase as well as the assessment of its enzymatic activity. Purified pAG-MNase is hyperactive and enables large-scale analyses.

DNAproDB: an updated database for the automated and interactive analysis of protein-DNA complexes.

DNAproDB (https://dnaprodb.usc.edu/) is a database, visualization tool, and processing pipeline for analyzing structural features of protein-DNA interactions. Here, we present a substantially updated version of the database through additional structural annotations, search, and user interface functionalities. The update expands the number of pre-analyzed protein-DNA structures, which are automatically updated weekly. The analysis pipeline identifies water-mediated hydrogen bonds that are incorporated into the visualizations of protein-DNA complexes. Tertiary structure-aware nucleotide layouts are now available. New file formats and external database annotations are supported. The website has been redesigned, and interacting with graphs and data is more intuitive. We also present a statistical analysis on the updated collection of structures revealing salient patterns in protein-DNA interactions.

Molecular dynamics and small-angle x-ray scattering: a comparison computational and experimental approaches to studying the structure of biological complexes

The results of studying DNA-protein complexes using two independent structural methods – molecular dynamics (MD) and small-angle X-ray scattering (SAXS) – are compared. MD is a computational method that allows visualization of macromolecule behavior in real environmental conditions based on the laws of physics but suffers from numerous simplifications. SAXS is an X-ray method that allows the reconstruction of the three-dimensional structure of an object in solution based on the one-dimensional profile of small-angle scattering, which presents the problem of ambiguity in solving inverse problems. The use of structural characteristics of complexes obtained by the SAXS method for validating 3D structural models obtained in MD experiments has significantly reduced the ambivalence of theoretical predictions and demonstrated the effectiveness of combining MD and SAXS methods for solving structural biology problems.

WRKY Transcription Factors (TFs) as Key Regulators of Plant Resilience to Environmental Stresses: Current Perspective

Plants encounter various stresses in their natural environments and can effectively respond to only one stress at a time. Through a complex gene network, transcription factors (TFs) such as WRKY TFs regulate a diverse array of stress responses. The clarification of the structural characteristics of WRKY proteins, along with recent advancements in molecular dynamics simulations, has shed light on the formation, stability, and interactions of DNA–protein complexes. This has provided a novel viewpoint regarding the control of WRKY TFs. The investigation of superfamilies, encompassing their historical development, diversity, and evolutionary patterns, has become feasible due to the transcriptome approach’s capacity to provide extensive and comprehensive transcripts. The significance of WRKY TFs lies in their pivotal role within several signaling cascades and regulatory networks that influence plant defense responses. The present review summarizes the functional aspects of the high-volume sequence data of WRKY TFs from different species studied to date. Moreover, a comparative analysis approach was utilized to determine the functions of the identified WRKY TFs in response to both abiotic and biotic stresses, as revealed through numerous studies on different plant species. The results of this review will be pivotal in understanding evolutionary events and the significance of WRKY TFs in the context of climate change, incorporating new scientific evidence to propose an innovative viewpoint.

Charging of DNA Complexes in Positive-Mode Native Electrospray Ionization Mass Spectrometry.

Native mass spectrometry (nMS) provides insights into the structures and dynamics of biomacromolecules in their native-like states by preserving noncovalent interactions through "soft" electrospray ionization (ESI). For native proteins, the number of charges that are acquired scales with the surface area and mass. Here, we explore the effect of highly negatively charged DNA on the ESI charge of protein complexes and find a reduction of the mass-to-charge ratio as well as a greater variation. The charge state distributions of pure DNA assemblies show a lower mass-to-charge ratio than proteins due to their greater density in the gas phase, whereas the charge of protein-DNA complexes can additionally be influenced by the distribution of the ESI charges, ion pairing events, and collapse of the DNA components. Our findings suggest that structural features of protein-DNA complexes can result in lower charge states than expected for proteins.

The Impact of Disease Severity on the Serum Levels of Significant Neutrophil Extracellular Trap (NET) Proteins in Patients with Psoriasis.

Psoriasis is an inflammatory skin disease with various symptoms of differing severities and with the reported prominent involvement of neutrophil extracellular traps (NETs). The excitation of neutrophils, e.g., by interleukin 8 (IL-8) or lipopolysaccharide (LPS), leads to the citrullination of histones and the release of protein-DNA complexes into the extracellular space, where they are digested by DNases. Our aim was to explore data on the levels of protein-complexed DNAs neutrophil elastase-DNA (NE-DNA) and myeloperoxidase-DNA (MPO-DNA), citrullinated histones (citH2, citH3, citH4), and NET-degrading enzyme DNase I in the serum of psoriatic patients with varying severities of clinical symptoms assessed with the Psoriasis Area Severity Index (PASI), Body Surface Area (BSA), and Dermatology Life Quality Index (DLQI) scores. The levels of factors were detected in 52 patients with psoriasis and 22 healthy volunteers by the enzyme-linked immunosorbent assay (ELISA). The results showed the elevated levels of NE-DNA, MPO-DNA, citH3, and DNase I in the patients with psoriasis compared to healthy volunteers (p < 0.05). Additionally, changes were noticed in the levels of NE-DNA, citH3, and DNase I, depending on the severity of symptoms (p < 0.05). In mild psoriasis (PASI < 10, BSA < 10, DLQI < 10), the suppressing activity of the enzyme caused the impaired ability to remove the physiological level of NETs, whereas in moderate to severe psoriasis (PASI ≥ 10, BSA ≥ 10, DLQI ≥ 10), the enhanced activity of DNase I failed to remove NETs due to the observed overexpression. It may, thus, be concluded that the mechanism of action of NETs, which play an undeniable role in psoriatic diseases, seem to follow two different paths depending on the severity of disease, which may be crucial in selecting potential anti-NET treatment methods.

Identification of sanguinarine as c-MYC transcription inhibitor through enhancing the G-quadruplex-NM23-H2 interactions

c-MYC is a proto-oncogene ubiquitously overexpressed in various cancers. The formation of G-quadruplex (G4) structures within the c-MYC promoter region can regulate its transcription by interfering with protein binding. Consequently, small molecules targeting c-MYC G4 have emerged as promising anticancer agents. Herein, we report that sanguinarine (SG) and its analogs exhibit a high affinity for c-MYC G4 and potently modulate G4-protein interactions within a natural product library. Notably, SG uniquely enhances NM23-H2 binding to c-MYC G4, both in vitro and in cellular contexts, leading to c-MYC transcriptional repression and subsequent inhibition of cancer cell growth in an NM23-H2-dependent manner. Mechanistic studies and molecular modeling suggest that SG binds to the c-MYC G4/NM23-H2 interface, acting as an orthosteric stabilizer of the DNA-protein complex and preventing c-MYC transcription. Our findings identify SG as a potent c-MYC transcription inhibitor and provide a novel strategy for developing G4-targeting anticancer therapeutics through modulation of G4-protein interactions.

Nucleosomal DNA unwinding pathway through canonical and non-canonical histone disassembly

The nucleosome including H2A.B, a mammalian-specific H2A variant, plays pivotal roles in spermatogenesis, embryogenesis, and oncogenesis, indicating unique involvement in transcriptional regulation distinct from canonical H2A nucleosomes. Despite its significance, the exact regulatory mechanism remains elusive. This study utilized solid-state nanopores to investigate DNA unwinding dynamics, applying local force between DNA and histones. Comparative analysis of canonical H2A and H2A.B nucleosomes demonstrated that the H2A.B variant required a lower voltage for complete DNA unwinding. Furthermore, synchronization analysis and molecular dynamics simulations indicate that the H2A.B variant rapidly unwinds DNA, causing the H2A-H2B dimer to dissociate from DNA immediately upon disassembly of the histone octamer. In contrast, canonical H2A nucleosomes unwind DNA at a slower rate, suggesting that the H2A-H2B dimer undergoes a state of stacking at the pore. These findings suggest that nucleosomal DNA in the H2A.B nucleosomes undergoes a DNA unwinding process involving histone octamer disassembly distinct from that of canonical H2A nucleosomes, enabling smoother unwinding. The integrated approach of MD simulations and nanopore measurements is expected to evolve into a versatile tool for studying molecular interactions, not only within nucleosomes but also through the forced dissociation of DNA-protein complexes.

Chemical crosslinking extends and complements UV crosslinking in analysis of RNA/DNA nucleic acid-protein interaction sites by mass spectrometry.

UV (ultra-violet) crosslinking with mass spectrometry (XL-MS) has been established for identifying RNA-and DNA-binding proteins along with their domains and amino acids involved. Here, we explore chemical XL-MS for RNA-protein, DNA-protein, and nucleotide-protein complexes in vitro and in vivo . We introduce a specialized nucleotide-protein-crosslink search engine, NuXL, for robust and fast identification of such crosslinks at amino acid resolution. Chemical XL-MS complements UV XL-MS by generating different crosslink species, increasing crosslinked protein yields in vivo almost four-fold and thus it expands the structural information accessible via XL-MS. Our workflow facilitates integrative structural modelling of nucleic acid-protein complexes and adds spatial information to the described RNA-binding properties of enzymes, for which crosslinking sites are often observed close to their cofactor-binding domains. In vivo UV and chemical XL-MS data from E. coli cells analysed by NuXL establish a comprehensive nucleic acid-protein crosslink inventory with crosslink sites at amino acid level for more than 1500 proteins. Our new workflow combined with the dedicated NuXL search engine identified RNA crosslinks that cover most RNA-binding proteins, with DNA and RNA crosslinks detected in transcriptional repressors and activators.

Geometric deep learning of protein–DNA binding specificity

Predicting protein–DNA binding specificity is a challenging yet essential task for understanding gene regulation. Protein–DNA complexes usually exhibit binding to a selected DNA target site, whereas a protein binds, with varying degrees of binding specificity, to a wide range of DNA sequences. This information is not directly accessible in a single structure. Here, to access this information, we present Deep Predictor of Binding Specificity (DeepPBS), a geometric deep-learning model designed to predict binding specificity from protein–DNA structure. DeepPBS can be applied to experimental or predicted structures. Interpretable protein heavy atom importance scores for interface residues can be extracted. When aggregated at the protein residue level, these scores are validated through mutagenesis experiments. Applied to designed proteins targeting specific DNA sequences, DeepPBS was demonstrated to predict experimentally measured binding specificity. DeepPBS offers a foundation for machine-aided studies that advance our understanding of molecular interactions and guide experimental designs and synthetic biology.

The Dynamic Plasticity of P. aeruginosa CueR Copper Transcription Factor upon Cofactor and DNA Binding.

Bacteria use specialized proteins, like transcription factors, to rapidly control metal ion balance. CueR is a Gram-negative bacterial copper regulator. The structure of E. coli CueR complexed with Cu(I) and DNA was published, since then many studies have shed light on its function. However, P. aeruginosa CueR, which shows high sequence similarity to E. coli CueR, has been less studied. Here, we applied room-temperature electron paramagnetic resonance (EPR) measurements to explore changes in dynamics of P. aeruginosa CueR in dependency of copper concentrations and interaction with two different DNA promoter regions. We showed that P. aeruginosa CueR is less dynamic than the E. coli CueR protein and exhibits much higher sensitivity to DNA binding as compared to its E. coli CueR homolog. Moreover, a difference in dynamical behavior was observed when P. aeruginosa CueR binds to the copZ2 DNA promoter sequence compared to the mexPQ-opmE promoter sequence. Such dynamical differences may affect the expression levels of CopZ2 and MexPQ-OpmE proteins in P. aeruginosa. Overall, such comparative measurements of protein-DNA complexes derived from different bacterial systems reveal insights about how structural and dynamical differences between two highly homologous proteins lead to quite different DNA sequence-recognition and mechanistic properties.

Exploring protein-mediated compaction of DNA by coarse-grained simulations and unsupervised learning

Protein-DNA interactions and protein-mediated DNA compaction play key roles in a range of biological processes. The length scales typically involved in DNA bending, bridging, looping, and compaction (≥1 kbp) are challenging to address experimentally or by all-atom molecular dynamics simulations, making coarse-grained simulations a natural approach. Here, we present a simple and generic coarse-grained model for DNA-protein and protein-protein interactions and investigate the role of the latter in the protein-induced compaction of DNA. Our approach models the DNA as a discrete worm-like chain. The proteins are treated in the grand canonical ensemble, and the protein-DNA binding strength is taken from experimental measurements. Protein-DNA interactions are modeled as an isotropic binding potential with an imposed binding valency without specific assumptions about the binding geometry. To systematically and quantitatively classify DNA-protein complexes, we present an unsupervised machine learning pipeline that receives a large set of structural order parameters as input, reduces the dimensionality via principal-component analysis, and groups the results using a Gaussian mixture model. We apply our method to recent data on the compaction of viral genome-length DNA by HIV integrase and find that protein-protein interactions are critical to the formation of looped intermediate structures seen experimentally. Our methodology is broadly applicable to DNA-binding proteins and protein-induced DNA compaction and provides a systematic and semi-quantitative approach for analyzing their mesoscale complexes.

Computational insights into intrinsically disordered regions in protein-nucleic acid complexes

We study transitions in intrinsically disordered regions (IDRs) upon complex formation, utilizing X-ray-solved structural dataset of protein-DNA and protein-RNA complexes, along with their available unbound protein forms. The identified IDRs are categorized into three classes: Disordered-to-Ordered (D-O), Disordered-to-Partial Ordered (D-PO) and Disordered-to-Disordered (D-D) after comparing them in unbound and complex forms. In the D-O class, IDRs form secondary structures like coils, helices, and strands upon binding to nucleic acids. Though a majority of these IDRs are present at the surface of the complexes, a significant number of IDRs are also observed at the interfaces and are involved in polar interactions. The hydrogen bonds made by the interface IDRs (B_IDRs) with phosphates and bases of nucleic acids are comparatively more than those formed with sugars. B_IDRs form more H-bonds with the ribose in protein-RNA than with the deoxyribose in protein-DNA. Among the B_IDRs, Arg and Lys prefer to interact with the major and minor grooves of DNA and RNA, respectively. Ser, however, prefers the minor groove in both the nucleic acids. Interestingly, we report 61 and 48 IDRs in 31 protein-DNA and 22 protein-RNA complexes, respectively, suggesting nucleic acid binding to proteins may also result in ordered-to-disordered transitions.

Comparison of the Amyloid Plaque Proteome in Down Syndrome, Early-Onset Alzheimer's Disease and Late-Onset Alzheimer's Disease.

Down syndrome (DS) is strongly associated with Alzheimer's disease (AD), attributable to APP overexpression. DS exhibits Amyloid-β (Aβ) and Tau pathology similar to early-onset AD (EOAD) and late-onset AD (LOAD). The study aimed to evaluate the Aβ plaque proteome of DS, EOAD and LOAD. Using unbiased localized proteomics, we analyzed amyloid plaques and adjacent plaque-devoid tissue ('non-plaque') from post-mortem paraffin-embedded tissues in four cohorts (n = 20/group): DS (59.8 ± 4.99 y/o), EOAD (63 ± 4.07 y/o), LOAD (82.1 ± 6.37 y/o) and controls (66.4 ± 13.04). We assessed functional associations using Gene Ontology (GO) enrichment and protein interaction networks. We identified differentially abundant Aβ plaque proteins vs. non-plaques (FDR < 5%, fold-change > 1.5) in DS (n = 132), EOAD (n = 192) and in LOAD (n = 128); there were 43 plaque-associated proteins shared between all groups. Positive correlations (p < 0.0001) were observed between plaque-associated proteins in DS and EOAD (R2 = 0.77), DS and LOAD (R2 = 0.73), and EOAD vs. LOAD (R2 = 0.67). Top Biological process (BP) GO terms (p < 0.0001) included lysosomal transport for DS, immune system regulation for EOAD, and lysosome organization for LOAD. Protein networks revealed a plaque enriched signature across all cohorts involving APP metabolism, immune response, and lysosomal functions. In DS, EOAD and LOAD non-plaque vs. control tissue, we identified 263, 269, and 301 differentially abundant proteins, including 65 altered non-plaque proteins across all cohorts. Differentially abundant non-plaque proteins in DS showed a significant (p < 0.0001) but weaker positive correlation with EOAD (R2 = 0.59) and LOAD (R2 = 0.33) compared to the stronger correlation between EOAD and LOAD (R2 = 0.79). The top BP GO term for all groups was chromatin remodeling (DS p = 0.0013, EOAD p = 5.79×10- 9, and LOAD p = 1.69×10- 10). Additional GO terms for DS included extracellular matrix (p = 0.0068), while EOAD and LOAD were associated with protein-DNA complexes and gene expression regulation (p < 0.0001). We found strong similarities among the Aβ plaque proteomes in individuals with DS, EOAD and LOAD, and a robust association between the plaque proteomes and lysosomal and immune-related pathways. Further, non-plaque proteomes highlighted altered pathways related to chromatin structure and extracellular matrix (ECM), the latter particularly associated with DS. We identified novel Aβ plaque proteins, which may serve as biomarkers or therapeutic targets.

Characterization of the DNA accessibility of chloroplast genomes in grasses

Although the chloroplast genome (cpDNA) of higher plants is known to exist as a large protein-DNA complex called ‘plastid nucleoid’, researches on its DNA state and regulatory elements are limited. In this study, we performed the assay for transposase-accessible chromatin sequencing (ATAC-seq) on five common tissues across five grasses, and found that the accessibility of different regions in cpDNA varied widely, with the transcribed regions being highly accessible and accessibility patterns around gene start and end sites varying depending on the level of gene expression. Further analysis identified a total of 3970 putative protein binding footprints on cpDNAs of five grasses. These footprints were enriched in intergenic regions and co-localized with known functional elements. Footprints and their flanking accessibility varied dynamically among tissues. Cross-species analysis showed that footprints in coding regions tended to overlap non-degenerate sites and contain a high proportion of highly conserved sites, indicating that they are subject to evolutionary constraints. Taken together, our results suggest that the accessibility of cpDNA has biological implications and provide new insights into the transcriptional regulation of chloroplasts.

HNRNPA2B1 stabilizes NFATC3 levels to potentiate its combined actions with FOSL1 to mediate vasculogenic mimicry in GBM cells

BackgroundVasculogenic mimicry (VM) is an enigmatic physiological feature that influences blood supply within glioblastoma (GBM) tumors for their sustained growth. Previous studies identify NFATC3, FOSL1 and HNRNPA2B1 as significant mediators of VEGFR2, a key player in vasculogenesis, and their molecular relationships may be crucial for VM in GBM.AimsThe aim of this study was to understand how NFATC3, FOSL1 and HNRNPA2B1 collectively influence VM in GBM.MethodsWe have investigated the underlying gene regulatory mechanisms for VM in GBM cell lines U251 and U373 in vitro and in vivo. In vitro cell-based assays were performed to explore the role of NFATC3, FOSL1 and HNRNPA2B1 in GBM cell proliferation, VM and migration, in the context of RNA interference (RNAi)-mediated knockdown alongside corresponding controls. Western blotting and qRT-PCR assays were used to examine VEGFR2 expression levels. CO-IP was employed to detect protein–protein interactions, ChIP was used to detect DNA–protein complexes, and RIP was used to detect RNA–protein complexes. Histochemical staining was used to detect VM tube formation in vivo.ResultsFocusing on NFATC3, FOSL1 and HNRNPA2B1, we found each was significantly upregulated in GBM and positively correlated with VM-like cellular behaviors in U251 and U373 cell lines. Knockdown of NFATC3, FOSL1 or HNRNPA2B1 each resulted in decreased levels of VEGFR2, a key growth factor gene that drives VM, as well as the inhibition of proliferation, cell migration and extracorporeal VM activity. Chromatin immunoprecipitation (ChIP) studies and luciferase reporter gene assays revealed that NFATC3 binds to the promoter region of VEGFR2 to enhance VEGFR2 gene expression. Notably, FOSL1 interacts with NFATC3 as a co-factor to potentiate the DNA-binding capacity of NFATC3, resulting in enhanced VM-like cellular behaviors. Also, level of NFATC3 protein in cells was enhanced through HNRNPA2B1 binding of NFATC3 mRNA. Furthermore, RNAi-mediated silencing of NFATC3, FOSL1 and HNRNPA2B1 in GBM cells reduced their capacity for tumor formation and VM-like behaviors in vivo.ConclusionTaken together, our findings identify NFATC3 as an important mediator of GBM tumor growth through its molecular and epistatic interactions with HNRNPA2B1 and FOSL1 to influence VEGFR2 expression and VM-like cellular behaviors.Graphical 1. NFATC3 binds to the promoter region of VEGFR2 to enhance VEGFR2 gene expression which leads to an increase in VM of GBM.2. FOSL1 interacts with NFATC3 to further facilitate VEGFR2 gene expression and VM.3. HNRNPA2B1 enhances NFATC3 mRNA stability to increase VEGFR2 expression and VM.

Cla4A, a Novel Regulator of Gene Expression Networks Required for Asexual and Insect-Pathogenic Lifecycles of Beauveria bassiana.

Cla4, an orthologous p21-activated kinase crucial for non-entomopathogenic fungal lifestyles, has two paralogs (Cla4A/B) functionally unknown in hypocrealean entomopathogens. Here, we report a regulatory role of Cla4A in gene expression networks of Beauveria bassiana required for asexual and entomopathogenic lifecycles while Cla4B is functionally redundant. The deletion of cla4A resulted in severe growth defects, reduced stress tolerance, delayed conidiation, altered conidiation mode, impaired conidial quality, and abolished pathogenicity through cuticular penetration, contrasting with no phenotype affected by cla4B deletion. In ∆cla4A, 5288 dysregulated genes were associated with phenotypic defects, which were restored by targeted gene complementation. Among those, 3699 genes were downregulated, including more than 1300 abolished at the transcriptomic level. Hundreds of those downregulated genes were involved in the regulation of transcription, translation, and post-translational modifications and the organization and function of the nuclear chromosome, chromatin, and protein-DNA complex. DNA-binding elements in promoter regions of 130 dysregulated genes were predicted to be targeted by Cla4A domains. Samples of purified Cla4A extract were proven to bind promoter DNAs of 12 predicted genes involved in multiple stress-responsive pathways. Therefore, Cla4A acts as a novel regulator of genomic expression and stability and mediates gene expression networks required for insect-pathogenic fungal adaptations to the host and environment.