Genes and regulatory networks involved in persistence of Mycobacterium tuberculosis

Regulatory role of Mycobacterium tuberculosis MtrA on dormancy/resuscitation revealed by a novel target gene-mining strategy

Introduction

The unique dormancy of Mycobacterium tuberculosis plays a significant role in the major clinical treatment challenge of tuberculosis, such as its long treatment cycle, antibiotic resistance, immune escape, and high latent infection rate.

Methods

To determine the function of MtrA, the only essential response regulator, one strategy was developed to establish its regulatory network according to high-quality genome-wide binding sites.

Results and discussion

The complex modulation mechanisms were implied by the strong bias distribution of MtrA binding sites in the noncoding regions, and 32.7% of the binding sites were located inside the target genes. The functions of 288 potential MtrA target genes predicted according to 294 confirmed binding sites were highly diverse, and DNA replication and damage repair, lipid metabolism, cell wall component biosynthesis, cell wall assembly, and cell division were the predominant pathways. Among the 53 pathways shared between dormancy/resuscitation and persistence, which accounted for 81.5% and 93.0% of the total number of pathways, respectively, MtrA regulatory genes were identified not only in 73.6% of their mutual pathways, but also in 75.4% of the pathways related to dormancy/resuscitation and persistence respectively. These results suggested the pivotal roles of MtrA in regulating dormancy/resuscitation and the apparent relationship between dormancy/resuscitation and persistence. Furthermore, the finding that 32.6% of the MtrA regulons were essential in vivo and/or in vitro for M. tuberculosis provided new insight into its indispensability. The findings mentioned above indicated that MtrA is a novel promising therapeutic target for tuberculosis treatment since the crucial function of MtrA may be a point of weakness for M. tuberculosis.

Revisiting tuberculosis screening: An insight to complementary diagnosis and prospective molecular approaches for the recognition of the dormant TB infection in human and cattle hosts

Tuberculosis (TB) is defined as a chronic infection in both human and cattle hosts and many subclinical cases remain undetected. After the pathogen is inhaled by a host, phagocyted bacilli can persist inside macrophages surviving intracellularly. Hosts develop granulomatous lesions in the lungs or lymph nodes, limiting infection. However, bacilli become persister cells. Immunological diagnosis of TB is performed basically by routine tuberculin skin test (TST), and in some cases, by ancillary interferon-gamma release assay (IGRA). The concept of human latent TB infection (LTBI) by M. tuberculosis is recognized in cohorts without symptoms by routine clinical diagnostic tests, and nowadays IGRA tests are used to confirm LTBI with either active or latent specific antigens of M. tuberculosis. On the other hand, dormant infection in cattle by M. bovis has not been described by TST or IGRA testing as complications occur by cross-reactive immune responses to homolog antigens of environmental mycobacteria or a false-negative test by anergic states of a wained bovine immunity, evidencing the need for deciphering more specific biomarkers by new-generation platforms of analysis for detection of M. bovis dormant infection. The study and description of bovine latent TB infection (boLTBI) would permit the recognition of hidden animal infection with an increase in the sensitivity of routine tests for an accurate estimation of infected dairy cattle. Evidence of immunological and experimental analysis of LTBI should be taken into account to improve the study and the description of the still neglected boLTBI.

TREASURE: Text Mining Algorithm Based on Affinity Analysis and Set Intersection to Find the Action of Tuberculosis Drugs against Other Pathogens

Tuberculosis (TB) is one of the top causes of death in the world. Though TB is known as the world’s most infectious killer, it can be treated with a combination of TB drugs. Some of these drugs can be active against other infective agents, in addition to TB. We propose a framework called TREASURE (Text mining algoRithm basEd on Affinity analysis and Set intersection to find the action of tUberculosis dRugs against other pathogEns), which particularly focuses on the extraction of various drug–pathogen relationships in eight different TB drugs, namely pyrazinamide, moxifloxacin, ethambutol, isoniazid, rifampicin, linezolid, streptomycin and amikacin. More than 1500 research papers from PubMed are collected for each drug. The data collected for this purpose are first preprocessed, and various relation records are generated for each drug using affinity analysis. These records are then filtered based on the maximum co-occurrence value and set intersection property to obtain the required inferences. The inferences produced by this framework can help the medical researchers in finding cures for other bacterial diseases. Additionally, the analysis presented in this model can be utilized by the medical experts in their disease and drug experiments.

Identification and in silico functional prediction of lineage-specific SNPs distributed in DosR-related proteins and resuscitation-promoting factor proteins of Mycobacterium tuberculosis

One-third of the world population is infected by Mycobacterium tuberculosis, which may persist in the latent or dormant state. Bacteria can shift to dormancy when encountering harsh conditions such as low oxygen, nutrient starvation, high acidity and host immune defenses. Genes related to the dormancy survival regulator (DosR) regulon are responsible for the inhibition of aerobic respiration and replication, which is required to enter dormancy. Conversely, resuscitation-promoting factor (rpf) proteins participate in reactivation from dormancy and the development of active tuberculosis (TB). Many DosR regulon and rpf proteins are immunodominant T cell antigens that are highly expressed in latent TB infection. They could serve as TB vaccine candidates and be used for diagnostic development. We explored the genetic polymorphisms of 50 DosR-related genes and 5 rpf genes among 1,170 previously sequenced clinical M. tuberculosis genomes. Forty-three lineage- or sublineage-specific nonsynonymous single nucleotide polymorphisms (nsSNPs) were identified. Ten nsSNPs were specific to all Mtb isolates belonging to lineage 1 (L1). Two common sublineages, the Beijing family (L2.2) and EAI2 (L1.2.1), differed at as many as 26 lineage- or sublineage-specific SNPs. DosR regulon genes related to membrane proteins and the rpf family possessed mean dN/dS ratios greater than one, suggesting that they are under positive selection. Although the T cell epitope regions of DosR-related and rpf antigens were quite conserved, we found that the epitopes in L1 had higher rates of genetic polymorphisms than the other lineages. Some mutations in immunogenic epitopes of the antigens were specific to particular M. tuberculosis lineages. Therefore, the genetic diversity of the DosR regulon and rpf proteins might impact the adaptation of M. tuberculosis to the dormant state and the immunogenicity of latency antigens, which warrants further investigation.

Truncation of C-Terminal Intrinsically Disordered Region of Mycobacterial Rv1915 Facilitates Production of “Difficult-to-Purify” Recombinant Drug Target

Availability of purified drug target is a prerequisite for its structural and functional characterization. However, aggregation of recombinant protein as inclusion bodies (IBs) is a common problem during the large scale production of overexpressed protein in heterologous host. Such proteins can be recovered from IB pool using some mild solubilizing agents such as low concentration of denaturants or detergents, alcohols and osmolytes. This study reports optimization of solubilization buffer for recovery of soluble and biologically active recombinant mycobacterial Rv1915/ICL2a from IBs. Even though the target protein could be solubilized successfully with mild agents (sarcosine and βME) without using denaturants, it failed to bind on Ni-NTA resin. The usual factors such as loss of His6-tag due to proteolysis, masking of the tag due to its location or protein aggregation were investigated, but the actual explanation, provided through bioinformatics analysis, turned out to be presence of intrinsically disordered protein regions (IDPRs) at the C-terminus. These regions due to their inability to fold into ordered structure may lead to non-specific protein aggregation and hence reduced binding to Ni-NTA affinity matrix. With this rationale, 90 residues from the C-terminal of Rv1915/ICL2 were truncated, the variant successfully purified and characterized for ICL and MICL activities, supporting the disordered nature of Rv1915/ICL2a C-terminal. When a region that has definite structure associated in some mycobaterial strains such as CDC 1551 and disorder in others for instance Mycobacterium tuberculosis H37Rv, it stands to reason that larger interface in the later may have implication in binding to other cellular partner.

Impact of the Polymorphism of the PACRG and CD80 Genes on the Development of the Different Stages of Tuberculosis Infection

BACKGROUND

Tuberculosis (TB) is one of the most significant health-care problems worldwide. The host's genetics play an important role in the development of TB in humans. The disease progresses through several stages, each of which can be under the control of different genes. The precise genes influencing the different stages of the disease are not yet identified. The aim of the current study was to determine the associations between primary and secondary TB and the polymorphisms of novel candidate genes for TB susceptibility, namely CD79A, HCST, CXCR4, CD4, CD80, CP, PACRG, and CD69.

METHODS

A total of 357 patients with TB (130 cases with primary TB and 227 cases with secondary TB) from the Siberian region of Russia as well as 445 healthy controls were studied. The study was performed at the Research Institute of Medical Genetics, Tomsk NRMC, Tomsk, Russia, between July 2015 and November 2016. Genotyping was carried out using MALDI-TOF mass spectrometry and PCR-RFLP. The associations between the single-nucleotide polymorphisms and TB were assessed using logistic regression adjusting for covariates (age and gender). Multiple testing was addressed via the experiment-wise permutation approach. The statistical significance threshold was a P value less than 0.05 for the permutation P values. The analyses were done in R 3.2 statistical software.

RESULTS

An association was established between the rs1880661 variant of the CD80 gene and secondary TB and the rs10945890 variant of the PACRG gene and both primary and secondary TB. However, the same allele of PACRG appeared to be both a risk factor for reactivation (secondary TB) and a protector against primary infection.

CONCLUSION

The results suggested that the CD80 and PACRG genes were associated with susceptibility to different forms of TB infection in the Russian population.

Ribosomal maturation factor (RimP) is essential for survival of nontuberculous mycobacteria Mycobacterium fortuitum under in vitro acidic stress conditions

Mycobacterium fortuitum is an important human pathogenic NTM, which resists stress conditions inside macrophages by exploitation of specific genes. TnphoA-based transposon mutagenesis was employed to identify membrane genes responsible for survival of M. fortuitum under such stress conditions. A library of about 450 mutants was constructed after electroporation of vector pRT291 into wild-type M. fortuitum. On the basis of blue color development and alkaline phosphatase assay, 20 mutants were shortlisted to screen for growth and survival under acidic stress at pH 6.5, 5.5, 4.5, and 3.5. Mutant MT727 showed reduced growth and survival under acidic stress. The acid susceptible mutant MT727 was subjected to other in vitro stress conditions prevalent inside macrophages including oxidative, nutrient starvation and nitrosative stress. However, the mutant showed no appreciable difference in growth behavior under oxidative, nutrient starvation and nitrosative stress conditions in comparison to the wild type. Genomic and bioinformatics analysis of MT727 led to identification of putative ribosomal maturation factor RimP of M. fortuitum to be affected by mutagenesis, showing closest homology to M. abscessus RimP. In silico functional interaction of RimP protein using STRING database showed its interaction with proteins of ribosomal assembly and maturation. Results indicate role of rimP gene in survival of M. fortuitum under acidic stress conditions which may be further explored for use as a potential drug target against M. fortuitum and other mycobacterial infections.

Site-specific His/Asp phosphoproteomic analysis of prokaryotes reveals putative targets for drug resistance

BACKGROUND

Phosphorylation of amino acid residues on proteins is an important and common post-translational modification in both eukaryotes and prokaryotes. Most research work has been focused on phosphorylation of serine, threonine or tyrosine residues, whereas phosphorylation of other amino acids are significantly less clear due to the controversy on their stability under standard bioanalytical conditions.

RESULTS

Here we applied a shotgun strategy to analyze the histidine and aspartate phosphorylations in different microbes. Our results collectively indicate that histidine and aspartate phosphorylations frequently occur also in proteins that are not part of the two-component systems. Noticeably, a number of the modified proteins are pathogenesis-related or essential for survival in host. These include the zinc ion periplasmic transporter ZnuA in Acinetobacter baumannii SK17, the multidrug and toxic compound extrusion (MATE) channel YeeO in Klebsiella pneumoniae NTUH-K2044, branched amino acid transporter AzlC in Vibrio vulnificus and the RNA-modifying pseudouridine synthase in Helicobacter pylori.

CONCLUSIONS

In summary, histidine and aspartate phosphorylation is likely to be ubiquitous and to take place in proteins of various functions. This work also sheds light into how these functionally important proteins and potential drug targets might be regulated at a post-translational level.

Significant under expression of the DosR regulon in M. tuberculosis complex lineage 6 in sputum

Mycobacterium africanum lineage (L) 6 is an important pathogen in West Africa, causing up to 40% of pulmonary tuberculosis (TB). The biology underlying the clinical differences between M. africanum and M. tuberculosis sensu stricto remains poorly understood. We performed ex vivo expression of 2179 genes of the most geographically dispersed cause of human TB, M. tuberculosis L4 and the geographically restricted, M. africanum L6 directly from sputa of 11 HIV-negative TB patients from The Gambia who had not started treatment. The DosR regulon was the most significantly decreased category in L6 relative to L4. Further, we identified nonsynonymous mutations in major DosR regulon genes of 44 L6 genomes of TB patients from The Gambia and Ghana. Using Lebek's test, we assessed differences in oxygen requirements for growth. L4 grew only at the aerobic surface while L6 grew throughout the medium. In the host, the DosR regulon is critical for M. tuberculosis in adaptation to oxygen limitation. However, M. africanum L6 appears to have adapted to growth under hypoxic conditions or to different biological niches. The observed under expression of DosR in L6 fits with the genomic changes in DosR genes, microaerobic growth and the association with extrapulmonary disease.

Combined Systems Approaches Reveal a Multistage Mode of Action of a Marine Antimicrobial Peptide against Pathogenic Escherichia coli and Its Protective Effect against Bacterial Peritonitis and Endotoxemia

A marine arenicin-3 derivative, N4, displayed potent antibacterial activity against Gram-negative bacteria, but its antibacterial mode of action remains elusive. The mechanism of action of N4 against pathogenic Escherichia coli was first researched by combined cytological and transcriptomic techniques in this study. The N4 peptide permeabilized the outer membrane within 1 min, disrupted the plasma membrane after 0.5 h, and localized in the cytoplasm within 5 min. Gel retardation and circular dichroism (CD) spectrum analyses demonstrated that N4 bound specifically to DNA and disrupted the DNA conformation from the B type to the C type. N4 inhibited 21.1% of the DNA and 20.6% of the RNA synthesis within 15 min. Several hallmarks of apoptosis-like cell death were exhibited by N4-induced E. coli, such as cell cycle arrest in the replication (R) and division(D) phases, reactive oxygen species production, depolarization of the plasma membrane potential, and chromatin condensation within 0.5 h. Deformed cell morphology, disappearance of the plasma membrane, leakage of the contents, and ghost cell formation were demonstrated by transmission electron microscopy, and nearly 100% of the bacteria were killed by N4. A total of 428 to 663 differentially expressed genes are involved in the response to N4, which are associated mainly with membrane biogenesis (53.9% to 56.7%) and DNA binding (13.3% to 14.9%). N4-protected mice that were lethally challenged with lipopolysaccharide (LPS) exhibited reduced levels of interleukin-6 (IL-6), IL-1β, and tumor necrosis factor alpha (TNF-α) in serum and protected the lungs from LPS-induced injury. These data facilitate an enhanced understanding of the mechanisms of marine antimicrobial peptides (AMPs) against Gram-negative bacteria and provide guidelines in developing and applying novel multitarget AMPs in the field of unlimited marine resources as therapeutics.

Identification of Lysine Acetylation in Mycobacterium abscessus Using LC-MS/MS after Immunoprecipitation

Mycobacterium abscessus (MAB), which manifests in the pulmonary system, is one of the neglected causes of nontuberculous mycobacteria (NTM) infection. Treatment against MAB is difficult, characterized by its intrinsic antibiotic drug resistance. Lysine acetylation can alter the physiochemical property of proteins in living organisms. This study aimed to determine if this protein post-translational modification (PTM) exists in a clinical isolate M. abscessus GZ002. We used the antiacetyl-lysine immunoprecipitation to enrich the low-abundant PTM proteins, followed by the LC-MS/MS analysis. The lysine acetylome of M. abscessus GZ002 was determined. There were 459 lysine acetylation sites found in 289 acetylated proteins. Lysine acetylation occurred in 5.87% of the M. abscessus GZ002 proteome, and at least 25% of them were growth essential. Aerobic respiration and carbohydrate metabolic pathways of M. abscessus GZ002 were enriched with lysine acetylation. Through bioinformatics analysis, we identified four major acetyl motif logos (K(ac)Y, K(ac)F, K(ac)H, and DK(ac)). Further comparison of the reported M. tuberculosis (MTB) acetylomes and that of MAB GZ002 revealed several common features between these two species. The lysine residues of several antibiotic-resistance, virulence, and persistence-related proteins were acetylated in both MAB GZ002 and MTB. There were 51 identical acetylation sites in 37 proteins found in common between MAB GZ002 and MTB. Overall, we demonstrate a profile of lysine acetylation in MAB GZ002 proteome that shares similarities with MTB. Interventions that target at these conserved sections may be valuable as anti-NTM or anti-TB therapies.

Comparative analyses of the proteins from Mycobacterium tuberculosis and human genomes: Identification of potential tuberculosis drug targets

Tuberculosis, one of the major infectious diseases affecting human beings is caused by the bacillus Mycobacterium tuberculosis. Increased resistance to known drugs commonly used for the treatment of tuberculosis has created an urgent need to identify new targets for validation and to develop drugs. In this study, we have used various bioinformatics tools, to compare the protein sequences from twenty-three M. tuberculosis genome strains along with the known human protein sequences, in order to identify the 'conserved' M. tuberculosis proteins absent in human. Further, based on the analysis of protein interaction networks, we selected one-hundred and forty proteins that were predicted as potential M. tuberculosis drug targets and prioritized according to the ranking of 'clusters' of interacting proteins. Comparison of the predicted 140 TB targets with literature indicated that 46 of them were previously reported, thereby increasing the confidence in our predictions of the remaining 94 targets too. The analyses of the structures and functions corresponding to the predicted potential TB drug targets indicated a diverse range of proteins that included ten 'druggable' targets with some of the known drugs.

Towards understanding the biological function of the unusual chaperonin Cpn60.1 (GroEL1) of Mycobacterium tuberculosis

The 60 kDa heat shock proteins, also known as Cpn60s (GroELs) are components of the essential protein folding machinery of the cell, but are also dominant antigens in many infectious diseases. Although generally essential for cellular survival, in some organisms such as Mycobacterium tuberculosis, one or more paralogous Cpn60s are known to be dispensable. In M. tuberculosis, Cpn60.2 (GroEL2) is essential for cell survival, but the biological role of the non-essential Cpn60.1 (GroEL1) is still elusive. To understand the relevance of Cpn60.1 (GroEL1) in M. tuberculosis physiology, detailed transcriptomic analyses for the wild type H37Rv and cpn60.1 knockout (groEL1-KO) were performed under in vitro stress conditions: stationary phase, cold shock, low aeration, mild cold shock and low pH. Additionally, the survival of the groEL1-KO was assessed in macrophages at multiplicity of infection (MOI) of 1:1 and 1:5. We observed that survival under low aeration was significantly compromised in the groEL1-KO. Further, the gene expression analyses under low aeration showed change in expression of several key virulence factors like two component system PhoP/R and MprA/B, sigma factors SigM and C and adversely affected known hypoxia response regulators Rv0081, Rv0023 and DosR. Our work is therefore suggestive of an important role of Cpn60.1 (GroEL1) for survival under low aeration by affecting the expression of genes known for hypoxia response.

Transcription factor regulatory network for early lung immune response to tuberculosis in mice

Numerous transcription factors (TFs) have been suggested to have a role in Mycobacterium tuberculosis infection; however, the TFs involved in the early immune response of lung cells remains to be fully elucidated. The present study aimed to identify TFs which may have a role in the early immune response to tuberculosis and the gene regulatory networks in which they are involved. Gene expression data obtained from microarray analysis of the early lung immune response to tuberculosis (Gene Expression Omnibus; accession no. GSE23014) was integrated with data for TF binding sites and protein-protein interactions in order to construct a TF regulatory network. The role of TFs in protein complexes, active modules, topology of the network and regulation of immune processes were investigated. The results demonstrated that the constructed gene regulatory network harbored 1,270 differentially expressed (DE) genes with 4,070 regulatory and protein-protein interactions. In addition, it was revealed that 17 DE TFs were involved in the positive regulation of numerous immunological and biological processes, including T cell activation, T cell proliferation and tuberculosis-associated gene expression, in the constructed regulatory network. Signal transducer and activator of transcription 4, interferon regulatory factor 8, spleen focus-forming virus proviral integration 1, enhancer of zeste homolog 2 and kruppel-like factor 4 were predicted to be the primary TFs regulating the DE genes during early lung infection by M. tuberculosis, as determined through various analyses of the gene regulatory network. In conclusion, the present study identified novel TFs involved in the early response to M. tuberculosis infection, which may enhance current understanding of the molecular mechanism underlying tuberculosis infection and introduced potential targets for novel tuberculosis therapies.

Latent tuberculosis infection: myths, models, and molecular mechanisms

The aim of this review is to present the current state of knowledge on human latent tuberculosis infection (LTBI) based on clinical studies and observations, as well as experimental in vitro and animal models. Several key terms are defined, including "latency," "persistence," "dormancy," and "antibiotic tolerance." Dogmas prevalent in the field are critically examined based on available clinical and experimental data, including the long-held beliefs that infection is either latent or active, that LTBI represents a small population of nonreplicating, "dormant" bacilli, and that caseous granulomas are the haven for LTBI. The role of host factors, such as CD4(+) and CD8(+) T cells, T regulatory cells, tumor necrosis factor alpha (TNF-α), and gamma interferon (IFN-γ), in controlling TB infection is discussed. We also highlight microbial regulatory and metabolic pathways implicated in bacillary growth restriction and antibiotic tolerance under various physiologically relevant conditions. Finally, we pose several clinically important questions, which remain unanswered and will serve to stimulate future research on LTBI.

Integrated gene co-expression network analysis in the growth phase of Mycobacterium tuberculosis reveals new potential drug targets

We have carried out weighted gene co-expression network analysis of Mycobacterium tuberculosis to gain insights into gene expression architecture during log phase growth. The differentially expressed genes between at least one pair of 11 different M. tuberculosis strains as source of biological variability were used for co-expression network analysis. This data included genes with highest coefficient of variation in expression. Five distinct modules were identified using topological overlap based clustering. All the modules together showed significant enrichment in biological processes: fatty acid biosynthesis, cell membrane, intracellular membrane bound organelle, DNA replication, Quinone biosynthesis, cell shape and peptidoglycan biosynthesis, ribosome and structural constituents of ribosome and transposition. We then extracted the co-expressed connections which were supported either by transcriptional regulatory network or STRING database or high edge weight of topological overlap. The genes trpC, nadC, pitA, Rv3404c, atpA, pknA, Rv0996, purB, Rv2106 and Rv0796 emerged as top hub genes. After overlaying this network on the iNJ661 metabolic network, the reactions catalyzed by 15 highly connected metabolic genes were knocked down in silico and evaluated by Flux Balance Analysis. The results showed that in 12 out of 15 cases, in 11 more than 50% of reactions catalyzed by genes connected through co-expressed connections also had altered fluxes. The modules 'Turquoise', 'Blue' and 'Red' also showed enrichment in essential genes. We could map 152 of the previously known or proposed drug targets in these modules and identified 15 new potential drug targets based on their high degree of co-expressed connections and strong correlation with module eigengenes.

Latent Tuberculosis: Models, Computational Efforts and the Pathogen's Regulatory Mechanisms during Dormancy

Latent tuberculosis is a clinical syndrome that occurs after an individual has been exposed to the Mycobacterium tuberculosis (Mtb) Bacillus, the infection has been established and an immune response has been generated to control the pathogen and force it into a quiescent state. Mtb can exit this quiescent state where it is unresponsive to treatment and elusive to the immune response, and enter a rapid replicating state, hence causing infection reactivation. It remains a gray area to understand how the pathogen causes a persistent infection and it is unclear whether the organism will be in a slow replicating state or a dormant non-replicating state. The ability of the pathogen to adapt to changing host immune response mechanisms, in which it is exposed to hypoxia, low pH, nitric oxide (NO), nutrient starvation, and several other anti-microbial effectors, is associated with a high metabolic plasticity that enables it to metabolize under these different conditions. Adaptive gene regulatory mechanisms are thought to coordinate how the pathogen changes their metabolic pathways through mechanisms that sense changes in oxygen tension and other stress factors, hence stimulating the pathogen to make necessary adjustments to ensure survival. Here, we review studies that give insights into latency/dormancy regulatory mechanisms that enable infection persistence and pathogen adaptation to different stress conditions. We highlight what mathematical and computational models can do and what they should do to enhance our current understanding of TB latency.

A transcriptional dynamic network during Arabidopsis thaliana pollen development

Background

To understand transcriptional regulatory networks (TRNs), especially the coordinated dynamic regulation between transcription factors (TFs) and their corresponding target genes during development, computational approaches would represent significant advances in the genome-wide expression analysis. The major challenges for the experiments include monitoring the time-specific TFs' activities and identifying the dynamic regulatory relationships between TFs and their target genes, both of which are currently not yet available at the large scale. However, various methods have been proposed to computationally estimate those activities and regulations. During the past decade, significant progresses have been made towards understanding pollen development at each development stage under the molecular level, yet the regulatory mechanisms that control the dynamic pollen development processes remain largely unknown. Here, we adopt Networks Component Analysis (NCA) to identify TF activities over time couse, and infer their regulatory relationships based on the coexpression of TFs and their target genes during pollen development.

Results

We carried out meta-analysis by integrating several sets of gene expression data related to Arabidopsis thaliana pollen development (stages range from UNM, BCP, TCP, HP to 0.5 hr pollen tube and 4 hr pollen tube). We constructed a regulatory network, including 19 TFs, 101 target genes and 319 regulatory interactions. The computationally estimated TF activities were well correlated to their coordinated genes' expressions during the development process. We clustered the expression of their target genes in the context of regulatory influences, and inferred new regulatory relationships between those TFs and their target genes, such as transcription factor WRKY34, which was identified that specifically expressed in pollen, and regulated several new target genes. Our finding facilitates the interpretation of the expression patterns with more biological relevancy, since the clusters corresponding to the activity of specific TF or the combination of TFs suggest the coordinated regulation of TFs to their target genes.

Conclusions

Through integrating different resources, we constructed a dynamic regulatory network of Arabidopsis thaliana during pollen development with gene coexpression and NCA. The network illustrated the relationships between the TFs' activities and their target genes' expression, as well as the interactions between TFs, which provide new insight into the molecular mechanisms that control the pollen development.