Construction of transcript regulation mechanism prediction models based on binding motif environment of transcription factor AoXlnR in Aspergillus oryzae.

DNA-binding transcription factors (TFs) play a central role in transcriptional regulation mechanisms, mainly through their specific binding to target sites on the genome and regulation of the expression of downstream genes. Therefore, a comprehensive analysis of the function of these TFs will lead to the understanding of various biological mechanisms. However, the functions of TFs in vivo are diverse and complicated, and the identified binding sites on the genome are not necessarily involved in the regulation of downstream gene expression. In this study, we investigated whether DNA structural information around the binding site of TFs can be used to predict the involvement of the binding site in the regulation of the expression of genes located downstream of the binding site. Specifically, we calculated the structural parameters based on the DNA shape around the DNA binding motif located upstream of the gene whose expression is directly regulated by one TF AoXlnR from Aspergillus oryzae, and showed that the presence or absence of expression regulation can be predicted from the sequence information with high accuracy ([Formula: see text]-1.0) by machine learning incorporating these parameters.

Deciphering the genomes of motility-deficient mutants of Vibrio alginolyticus 138-2.

The motility of Vibrio species plays a pivotal role in their survival and adaptation to diverse environments and is intricately associated with pathogenicity in both humans and aquatic animals. Numerous mutant strains of Vibrio alginolyticus have been generated using UV or EMS mutagenesis to probe flagellar motility using molecular genetic approaches. Identifying these mutations promises to yield valuable insights into motility at the protein structural physiology level. In this study, we determined the complete genomic structure of 4 reference specimens of laboratory V. alginolyticus strains: a precursor strain, V. alginolyticus 138-2, two strains showing defects in the lateral flagellum (VIO5 and YM4), and one strain showing defects in the polar flagellum (YM19). Subsequently, we meticulously ascertained the specific mutation sites within the 18 motility-deficient strains related to the polar flagellum (they fall into three categories: flagellar-deficient, multi-flagellar, and chemotaxis-deficient strains) by whole genome sequencing and mapping to the complete genome of parental strains VIO5 or YM4. The mutant strains had an average of 20.6 (±12.7) mutations, most of which were randomly distributed throughout the genome. However, at least two or more different mutations in six flagellar-related genes were detected in 18 mutants specifically selected as chemotaxis-deficient mutants. Genomic analysis using a large number of mutant strains is a very effective tool to comprehensively identify genes associated with specific phenotypes using forward genetics.

Identification of plb1 mutation that extends longevity via activating Sty1 MAPK in Schizosaccharomyces pombe.

To understand the lifespan of higher organisms, including humans, it is important to understand lifespan at the cellular level as a prerequisite. So, fission yeast is a good model organism for the study of lifespan. To identify the novel factors involved in longevity, we are conducting a large-scale screening of long-lived mutant strains that extend chronological lifespan (cell survival in the stationary phase) using fission yeast. One of the newly acquired long-lived mutant strains (No.98 mutant) was selected for analysis and found that the long-lived phenotype was due to a missense mutation (92Phe → Ile) in the plb1+ gene. plb1+ gene in fission yeast is a nonessential gene encoding a homolog of phospholipase B, but its functions under normal growth conditions, as well as phospholipase B activity, remain unresolved. Our analysis of the No.98 mutant revealed that the plb1 mutation reduces the integrity of the cellular membrane and cell wall and activates Sty1 via phosphorylation.

The GGDEF protein Dgc2 suppresses both motility and biofilm formation in the filamentous cyanobacterium Leptolyngbya boryana.

Colony pattern formations of bacteria with motility manifest complicated morphological self-organization phenomena. Leptolyngbya boryana is a filamentous cyanobacterium, which has been used as a genetic model organism for studying metabolism including photosynthesis and nitrogen fixation. A widely used type strain [wild type (WT) in this article] of this species has not been reported to show any motile activity. However, we isolated a spontaneous mutant strain that shows active motility (gliding activity) to give rise to complicated colony patterns, including comet-like wandering clusters and disk-like rotating vortices on solid media. Whole-genome resequencing identified multiple mutations in the genome of the mutant strain. We confirmed that inactivation of the candidate gene dgc2 (LBDG_02920) in the WT background was sufficient to give rise to motility and morphologically complex colony patterns. This gene encodes a protein containing the GGDEF motif which is conserved at the catalytic domain of diguanylate cyclase (DGC). Although DGC has been reported to be involved in biofilm formation, the dgc2 mutant significantly facilitated biofilm formation, suggesting a role for the dgc2 gene in suppressing both gliding motility and biofilm formation. Thus, Leptolyngbya is expected to be an excellent genetic model for studying dynamic colony pattern formation and to provide novel insights into the role of DGC family genes in biofilm formation. IMPORTANCE Self-propelled bacteria often exhibit complex collective behaviors, such as formation of dense-moving clusters, which are exemplified by wandering comet-like and rotating disk-like colonies; however, the molecular details of how these structures are formed are scant. We found that a strain of the filamentous cyanobacterium Leptolyngbya deficient in the GGDEF protein gene dgc2 elicits motility and complex and dynamic colony pattern formation, including comet-like and disk-like clusters. Although c-di-GMP has been reported to activate biofilm formation in some bacterial species, disruption of dgc2 unexpectedly enhanced it, suggesting a novel role for this GGDEF protein for inhibiting both colony pattern formation and biofilm formation.

Two Trk/Ktr/HKT-type potassium transporters, TrkG and TrkH, perform distinct functions in Escherichia coli K-12.

Escherichia coli K-12 possesses two versions of Trk/Ktr/HKT-type potassium ion (K+) transporters, TrkG and TrkH. The current paradigm is that TrkG and TrkH have largely identical characteristics, and little information is available regarding their functional differences. Here, we show using cation uptake experiments with K+ transporter knockout mutants that TrkG and TrkH have distinct ion transport activities and physiological roles. K+-transport by TrkG required Na+, whereas TrkH-mediated K+ uptake was not affected by Na+. An aspartic acid located five residues away from a critical glycine in the third pore-forming region might be involved in regulation of Na+-dependent activation of TrkG. In addition, we found that TrkG but not TrkH had Na+ uptake activity. Our analysis of K+ transport mutants revealed that TrkH supported cell growth more than TrkG; however, TrkG was able to complement loss of TrkH-mediated K+ uptake in E. coli. Furthermore, we determined that transcription of trkG in E. coli was downregulated but not completely silenced by the xenogeneic silencing factor H-NS (histone-like nucleoid structuring protein or heat-stable nucleoid-structuring protein). Taken together, the transport function of TrkG is clearly distinct from that of TrkH, and TrkG seems to have been accepted by E. coli during evolution as a K+ uptake system that coexists with TrkH.

Existence of two substates in the O intermediate of the bacteriorhodopsin photocycle.

The proton pumping cycle of bacteriorhodopsin (bR) is initiated when the retinal chromophore with the 13-trans configuration is photo-isomerized into the 13-cis configuration. To understand the recovery processes of the initial retinal configuration that occur in the late stage of the photocycle, we have performed a comprehensive analysis of absorption kinetics data collected at various pH levels and at different salt concentrations. The result of analysis revealed the following features of the late stages of the trans photocycle. i) Two substates occur in the O intermediate. ii) The visible absorption band of the first substate (O1) appears at a much shorter wavelength than that of the late substate (O2). iii) O1 is in rapid equilibrium with the preceding state (N), but O1 becomes less stable than N when an ionizable residue (X1) with a pKa value of 6.5 (in 2 M KCl) is deprotonated. iv) At a low pH and at a low salt concentration, the decay time constant of O2 is longer than those of the preceding states, but the relationship between these time constants is altered when the medium pH or the salt concentration is increased. On the basis of the present observations and previous studies on the structure of the chromophore in O, we suspect that the retinal chromophore in O1 takes on a distorted 13-cis configuration and the O1-to-O2 transition is accompanied by cis-to-trans isomerization about C13C14 bond.

Two substates in the O intermediate of the light-driven proton pump archaerhodopsin-2.

The proton pumping cycle of archaerhodopsin-2 (aR2) was investigated over a wide pH range and at different salt concentrations. We have found that two substates, which are spectroscopically and kinetically distinguishable, occur in the O intermediate. The first O-intermediate (O1) absorbs maximumly at ~580 nm, whereas the late O-intermediate (O2) absorbs maximumly at 605 nm. At neutral pH, O1 is in rapid equilibrium with the N intermediate. When the medium pH is increased, O1 becomes less stable than N and, in proportion to the amount of O1 in the dynamic equilibrium between N and O1, the formation rate of O2 decreases. By contrast, the decay rate of O2 increases ~100 folds when the pH of a low-salt membrane suspension is increased from 5.5 to 7.5 or when the salt concentration is increased to 2 M KCl. Together with our recent study on two substates in the O intermediate of bacteriorhodopsin (bR), the present study suggests that the thermally activated re-isomerization of the retinylidene chromophore into the initial all-trans configuration takes place in the O1-to-O2 transition; that is, O1 contains a distorted 13-cis chromophore. It is also found that the pKa value of the key ionizable residue (Asp101aR2, Asp96bR) in the proton uptake channel is elevated in the O1 state of aR2 as compared to the O1 state of bR. This implies that the structural property of O1 in the aR2 photocycle can be investigated over a wide pH range.

Mutations in the stator protein PomA affect switching of rotational direction in bacterial flagellar motor.

The flagellar motor rotates bi-directionally in counter-clockwise (CCW) and clockwise (CW) directions. The motor consists of a stator and a rotor. Recent structural studies have revealed that the stator is composed of a pentameric ring of A subunits and a dimer axis of B subunits. Highly conserved charged and neighboring residues of the A subunit interacts with the rotor, generating torque through a gear-like mechanism. The rotational direction is controlled by chemotaxis signaling transmitted to the rotor, with less evidence for the stator being involved. In this study, we report novel mutations that affect the switching of the rotational direction at the putative interaction site of the stator to generate rotational force. Our results highlight an aspect of flagellar motor function that appropriate switching of the interaction states between the stator and rotor is critical for controlling the rotational direction.

Identification of ksg1 mutation showing long-lived phenotype in fission yeast.

Fission yeast is a good model organism for the study of lifespan. To elucidate the mechanism, we screened for long-lived mutants. We found a nonsense mutation in the ksg1+ gene, which encodes an ortholog of mammalian PDK1 (phosphoinositide-dependent protein kinase). The mutation was in the PH domain of Ksg1 and caused defect in membrane localization and protein stability. Analysis of the ksg1 mutant revealed that the reduced amounts and/or activity of the Ksg1 protein are responsible for the increased lifespan. Ksg1 is essential for growth and known to phosphorylate multiple substrates, but the substrate responsible for the long-lived phenotype of ksg1 mutation is not yet known. Genetic analysis showed that deletion of pck2 suppressed the long-lived phenotype of ksg1 mutant, suggesting that Pck2 might be involved in the lifespan extension caused by ksg1 mutation.

Identification of sur2 mutation affecting the lifespan of fission yeast.

Yeast is a suitable model system to analyze the mechanism of lifespan. In this study, to identify novel factors involved in chronological lifespan, we isolated a mutant with a long chronological lifespan and found a missense mutation in the sur2+ gene, which encodes a homolog of Saccharomyces cerevisiae sphingolipid C4-hydroxylase in fission yeast. Characterization of the mutant revealed that loss of sur2 function resulted in an extended chronological lifespan. The effect of caloric restriction, a well-known signal for extending lifespan, is thought to be dependent on the sur2+ gene.

A slight bending of an α-helix in FliM creates a counterclockwise-locked structure of the flagellar motor in Vibrio.

Many bacteria swim by rotating flagella. The chemotaxis system controls the direction of flagellar rotation. Vibrio alginolyticus, which has a single polar flagellum, swims smoothly by rotating the flagellar motor counterclockwise (CCW) in response to attractants. In response to repellents, the motor frequently switches its rotational direction between CCW and clockwise (CW). We isolated a mutant strain that swims with a CW-locked rotation of the flagellum, which pulls rather than pushes the cell. This CW phenotype arises from a R49P substitution in FliM, which is the component in the C-ring of the motor that binds the chemotaxis signalling protein, phosphorylated CheY. However, this phenotype is independent of CheY, indicating that the mutation produces a CW conformation of the C-ring in the absence of CheY. The crystal structure of FliM with the R49P substitution showed a conformational change in the N-terminal α-helix of the middle domain of FliM (FliMM). This helix should mediates FliM-FliM interaction. The structural models of wild type and mutant C-ring showed that the relatively small conformational change in FliMM induces a drastic rearrangement of the conformation of the FliMM domain that generates a CW conformation of the C-ring.

Loss of cell wall integrity genes cpxA and mrcB causes flocculation in Escherichia coli.

Flocculation has been recognized for hundreds of years as an important phenomenon in brewing and wastewater treatment. However, the underlying molecular mechanisms remain elusive. The lack of a distinct phenotype to differentiate between slow-growing mutants and floc-forming mutants prevents the isolation of floc-related gene by conventional mutant screening. To overcome this, we performed a two-step Escherichia coli mutant screen. The initial screen of E. coli for mutants conferring floc production during high salt treatment yielded a mutant containing point mutations in 61 genes. The following screen of the corresponding single-gene mutants identified two genes, mrcB, encoding a peptidoglycan-synthesizing enzyme and cpxA, encoding a histidine kinase of a two-component signal transduction system that contributed to salt tolerance and flocculation prevention. Both single mutants formed flocs during high salt shock, these flocs contained cytosolic proteins. ΔcpxA exhibited decreased growth with increasing floc production and addition of magnesium to ΔcpxA suppressed floc production effectively. In contrast, the growth of ΔmrcB was inconsistent under high salt conditions. In both strains, flocculation was accompanied by the release of membrane vesicles containing inner and outer membrane proteins. Of 25 histidine kinase mutants tested, ΔcpxA produced the highest amount of proteins in floc. Expression of cpxP was up-regulated by high salt in ΔcpxA, suggesting that high salinity and activation of CpxR might promote floc formation. The finding that ΔmrcB or ΔcpxA conferred floc production indicates that cell envelope stress triggered by unfavorable environmental conditions cause the initiation of flocculation in E. coli.

Haloarcula mannanilytica sp. nov., a galactomannan-degrading haloarchaeon isolated from commercial salt.

A mannan-degrading halophilic archaeal strain, MD130-1T, was isolated from a commercial salt sample. Cells were motile, rod-shaped, and stained Gram-negative. Colonies were pink pigmented. Strain MD130-1T was able to grow at 1.5-4.6 M NaCl (optimum, 3.6 M) at pH 6.0-8.0 (optimum, pH 7.0) and at 25-50 °C (optimum, 40 °C). The DNA G+C content was 62.1 mol% (genome). The orthologous 16S rRNA gene sequence showed the highest similarity (99.4 %) to those of Haloarcula japonica JCM 7785T and Haloarcula hispanica JCM 8911T. The values of genome relatedness between strain MD130-1T and Haloarcula species were 84.33-85.96 % in ANIb and 30.4-32.9 % using GGDC formula 2. The polar lipids of strain MD130-1T were phosphatidylglycerol, phosphatidylglycerol phosphate methyl ester and triglycosyl diether-2. Based on the results of phenotypic and phylogenetic analyses, the strain represents a new species of the genus Haloarcula, for which the name Haloarcula mannanilytica sp. nov. is proposed. The type strain is MD130-1T (=JCM 33835T=KCTC 4287T) isolated from commercial salt made in Ishikawa prefecture, Japan.

Energetics and kinetics of substrate analog-coupled staphylococcal nuclease folding revealed by a statistical mechanical approach.

Protein conformational changes associated with ligand binding, especially those involving intrinsically disordered proteins, are mediated by tightly coupled intra- and intermolecular events. Such reactions are often discussed in terms of two limiting kinetic mechanisms, conformational selection (CS), where folding precedes binding, and induced fit (IF), where binding precedes folding. It has been shown that coupled folding/binding reactions can proceed along both CS and IF pathways with the flux ratio depending on conditions such as ligand concentration. However, the structural and energetic basis of such complex reactions remains poorly understood. Therefore, we used experimental, theoretical, and computational approaches to explore structural and energetic aspects of the coupled-folding/binding reaction of staphylococcal nuclease in the presence of the substrate analog adenosine-3',5'-diphosphate. Optically monitored equilibrium and kinetic data, combined with a statistical mechanical model, gave deeper insight into the relative importance of specific and Coulombic protein-ligand interactions in governing the reaction mechanism. We also investigated structural aspects of the reaction at the residue level using NMR and all-atom replica-permutation molecular dynamics simulations. Both approaches yielded clear evidence for accumulation of a transient protein-ligand encounter complex early in the reaction under IF-dominant conditions. Quantitative analysis of the equilibrium/kinetic folding revealed that the ligand-dependent CS-to-IF shift resulted from stabilization of the compact transition state primarily by weakly ligand-dependent Coulombic interactions with smaller contributions from specific binding energies. At a more macroscopic level, the CS-to-IF shift was represented as a displacement of the reaction "route" on the free energy surface, which was consistent with a flux analysis.

The role of ROC75 as a daytime component of the circadian oscillator in Chlamydomonas reinhardtii.

The circadian clocks in chlorophyte algae have been studied in two model organisms, Chlamydomonas reinhardtii and Ostreococcus tauri. These studies revealed that the chlorophyte clocks include some genes that are homologous to those of the angiosperm circadian clock. However, the genetic network architectures of the chlorophyte clocks are largely unknown, especially in C. reinhardtii. In this study, using C. reinhardtii as a model, we characterized RHYTHM OF CHLOROPLAST (ROC) 75, a clock gene encoding a putative GARP DNA-binding transcription factor similar to the clock proteins LUX ARRHYTHMO (LUX, also called PHYTOCLOCK 1 [PCL1]) and BROTHER OF LUX ARRHYTHMO (BOA, also called NOX) of the angiosperm Arabidopsis thaliana. We observed that ROC75 is a day/subjective day-phase-expressed nuclear-localized protein that associates with some night-phased clock genes and represses their expression. This repression may be essential for the gating of reaccumulation of the other clock-related GARP protein, ROC15, after its light-dependent degradation. The restoration of ROC75 function in an arrhythmic roc75 mutant under constant darkness leads to the resumption of circadian oscillation from the subjective dawn, suggesting that the ROC75 restoration acts as a morning cue for the C. reinhardtii clock. Our study reveals a part of the genetic network of C. reinhardtii clock that could be considerably different from that of A. thaliana.

Presynaptic MAST kinase controls opposing postsynaptic responses to convey stimulus valence in Caenorhabditis elegans.

Presynaptic plasticity is known to modulate the strength of synaptic transmission. However, it remains unknown whether regulation in presynaptic neurons can evoke excitatory and inhibitory postsynaptic responses. We report here that the Caenorhabditis elegans homologs of MAST kinase, Stomatin, and Diacylglycerol kinase act in a thermosensory neuron to elicit in its postsynaptic neuron an excitatory or inhibitory response that correlates with the valence of thermal stimuli. By monitoring neural activity of the valence-coding interneuron in freely behaving animals, we show that the alteration between excitatory and inhibitory responses of the interneuron is mediated by controlling the balance of two opposing signals released from the presynaptic neuron. These alternative transmissions further generate opposing behavioral outputs necessary for the navigation on thermal gradients. Our findings suggest that valence-encoding interneuronal activity is determined by a presynaptic mechanism whereby MAST kinase, Stomatin, and Diacylglycerol kinase influence presynaptic outputs.

gas1 mutation extends chronological lifespan via Pmk1 and Sty1 MAPKs in Schizosaccharomyces pombe.

In the longevity research by using yeasts, chronological lifespan is defined as the survival time after entry into stationary phase. Previously, screening for long lived mutants of Schizosaccharomyces pombe was performed to identify the novel factors involved in longevity. From this screening, one long lived mutant called as No.36 was obtained. In this study, we identified the mutation caused in gas1+, which encodes glucanosyltransferase (gas1-287 mutation) is responsible for the longevity of No.36 mutant. Through the analysis of this mutant, we found that cell wall perturbing agent micafungin also extends chronological lifespan in fission yeast. This lifespan extension depended on both Pmk1 and Sty1 MAP kinases, and longevity caused by the gas1-287 mutation also depended on these kinases. In summary, we propose that the gas1-287 mutation causes longevity as the similar mechanism as cell wall stress depending on Pmk1 and Sty1 MAPK pathways.

The mechanosensitive channel YbdG from Escherichia coli has a role in adaptation to osmotic up-shock.

Mechanosensitive channels play an important role in the adaptation of cells to hypo-osmotic shock. Among members of this channel family in Escherichia coli, the exact function and physiological role of the mechanosensitive channel homolog YbdG remain unclear. Characterization of YbdG's physiological role has been hampered by its lack of measurable transport activity. Using a nitrosoguanidine mutagenesis-aided screen in combination with next-generation sequencing, here we isolated a mutant with a point mutation in ybdG This mutation (resulting in a I167T change) conferred sensitivity to high osmotic stress, and the mutant cells differed from WT cells in morphology during hyperosmotic stress at alkaline pH. Interestingly, unlike the cells containing the I167T variant, a null-ybdG mutant did not exhibit this sensitivity and phenotype. Although I167T was located near the putative ion-conducting pore in a transmembrane region of YbdG, no change in ion channel activities of YbdG-I167T was detected. Of note, introduction of the WT C-terminal cytosolic region of YbdG into the I167T variant complemented the osmo-sensitive phenotype. Co-precipitation of proteins interacting with the C-terminal YbdG region led to the isolation of HldD and FbaA, whose overexpression in cells containing the YbdG-I167T variant partially rescued the osmo-sensitive phenotype. This study indicates that YbdG functions as a component of a mechanosensing system that transmits signals triggered by external osmotic changes to intracellular factors. The cellular role of YbdG uncovered here goes beyond its predicted function as an ion or solute transport protein.

Cholecystokinin induces crowing in chickens.

Animals that communicate using sound are found throughout the animal kingdom. Interestingly, in contrast to human vocal learning, most animals can produce species-specific patterns of vocalization without learning them from their parents. This phenomenon is called innate vocalization. The underlying molecular basis of both vocal learning in humans and innate vocalization in animals remains unknown. The crowing of a rooster is also innately controlled, and the upstream center is thought to be localized in the nucleus intercollicularis (ICo) of the midbrain. Here, we show that the cholecystokinin B receptor (CCKBR) is a regulatory gene involved in inducing crowing in roosters. Crowing is known to be a testosterone (T)-dependent behavior, and it follows that roosters crow but not hens. Similarly, T-administration induces chicks to crow. By using RNA-sequencing to compare gene expression in the ICo between the two comparison groups that either crow or do not crow, we found that CCKBR expression was upregulated in T-containing groups. The expression of CCKBR and its ligand, cholecystokinin (CCK), a neurotransmitter, was observed in the ICo. We also showed that crowing was induced by intracerebroventricular administration of an agonist specific for CCKBR. Our findings therefore suggest that the CCK system induces innate vocalization in roosters.

Comprehensive investigation of the gene expression system regulated by an Aspergillus oryzae transcription factor XlnR using integrated mining of gSELEX-Seq and microarray data.

BACKGROUND: Transcription factors (TFs) specifically bind to DNA sequences and control the expression of target genes. AoXlnR is a key TF involved in the expression of xylanolytic and cellulolytic enzymes in the filamentous fungi, Aspergillus oryzae. Genomic SELEX-Seq (gSELEX-Seq) can reveal the in vitro binding sites of a TF in a genome. To date, the gene expression network controlled by AoXlnR in A. oryzae is not fully explored. In this study, the data from gSELEX-Seq analysis and data mining were applied toward a comprehensive investigation of the AoXlnR-regulated transcriptional network in A. oryzae. RESULTS: Around 2000 promoters were selected as AoXlnR-binding DNAs using gSELEX-Seq, consequently identifying the genes downstream of them. On the other hand, 72 differentially expressed genes (DEGs) related to AoXlnR had been determined by microarray analysis. The intersecting set of genes, that were found using the gSELEX-Seq and the microarray analysis, had 51 genes. Further, the canonical AoXlnR-binding motifs, 5'-GGCT(A/G) A-3', were successfully identified in gSELEX-Seq. The motif numbers in each promoter of the DEGs and differential expression levels were correlated by in silico analysis. The analysis showed that the presence of both 5'-GGCTAA-3' and 5'-GGCTGA-3' motif has significantly high correlation with the differential expression levels of the genes. CONCLUSIONS: Genes regulated directly by AoXlnR were identified by integrated mining of data obtained from gSELEX-Seq and microarray. The data mining of the promoters of differentially expressed genes revealed the close relation between the presence of the AoXlnR-binding motifs and the expression levels of the downstream genes. The knowledge obtained in this study can contribute greatly to the elucidation of AoXlnR-mediated cellulose and xylan metabolic network in A. oryzae. The pipeline, which is based on integrated mining of data consisting of both in vitro characterization of the DNA-binding sites and TF phenotype, can be a robust platform for comprehensive analysis of the gene expression network via the TFs.