Elucidating the biology of transcription factor-DNA interaction for accurate identification of cis-regulatory elements

Targeting the RUNX3-miR-186-3p-DAT-IGF1R axis as a therapeutic strategy in a Parkinson's disease model

With the increasing age of the population worldwide, the incidence rate of Parkinson's disease (PD) is increasing annually. Currently, the treatment strategy for PD only improves clinical symptoms. No effective treatment strategy can slow down the progression of the disease. In the present study, whole transcriptome sequencing was used to obtain the mRNA and miRNA expression profiles in a PD mouse model, which revealed the pathogenesis of PD. The transcription factor RUNX3 upregulated the miR-186-3p expression in the PD model. Furthermore, the high miR-186-3p expression in PD can be targeted to inhibit the DAT expression, resulting in a decrease in the dopamine content of dopaminergic neurons. Moreover, miR-186-3p can be targeted to inhibit the IGF1R expression and prevent the activation of the IGF1R-P-PI3K-P-AKT pathway, thus increasing the apoptosis of dopaminergic neurons by regulating the cytochrome c-Bax-cleaved caspase-3 pathway. Our research showed that the RUNX3-miR-186-3p-DAT-IGF1R axis plays a key role in the pathogenesis of PD, and miR-186-3p is a potential target for the treatment of PD.

Molecular mechanism analysis of nontuberculous mycobacteria infection in patients with cystic fibrosis

Aim: This study aims to explore the molecular mechanisms of cystic fibrosis (CF) complicated with nontuberculous mycobacteria (NTM) infection. Materials & methods: Expression profiles of CF with NTM-infected patients were downloaded from GEO database. Intersection analysis yielded 78 genes associated with CF with NTM infection. The protein-protein interaction (PPI) network and the functions of hub genes were investigated. Results: Five hub genes (PIK3R1, IL1A, CXCR4, ACTN1, PFN1) were identified, which were primarily enriched in actin-related biological processes and pathways. Transcription factors RELA, JUN, NFKB1 and FOS that regulated hub genes modulated IL1A expression, while 21 other transcription factors regulated CXCR4 expression. Conclusion: In summary, this study may provide new insights into the mechanisms of CF with NTM infection.

Identification of Senescence-Related Biomarkers and Regulatory Networks in Intracerebral Hemorrhage

OBJECTIVES

Intracerebral hemorrhage (ICH) is a severe neurological disorder with substantial societal implications. Cellular senescence plays a critical role in ICH pathogenesis. This study aims to identify senescence-related biomarkers in ICH for diagnostic and therapeutic purposes.

METHODS

Raw data from GSE24265 in Gene Expression Omnibus was downloaded. Senescence-related genes were acquired from CellAge. Differential gene analysis was done between patients with ICH and controls. The intersection of ICH differentially expressed genes and senescence-related genes for senescence-related ICH genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed. Protein-protein interaction network was constructed through the Search Tool for the Retrieval of Interacting Genes. Single sample gene set enrichment analysis was done for immune cell infiltration and function evaluation in control and ICH groups. miRWalk2.0 database was used for microRNA predictions targeting ICH biomarkers. Transcriptional regulatory relationships unraveled by sentence-based text mining database was employed to predict transcription factors regulating identified biomarkers.

RESULTS

Thirteen senescence-related ICH genes were identified. They were primarily enriched in the positive regulation of angiogenesis and the Advanced Glycation End Product -Receptor for AGE signaling pathway in diabetic complications. Validation in the GSE149317 data set and receiver operating characteristic analysis highlighted Caveolin 1, C-X-C Motif Chemokine Ligand 1, ETS proto-oncogene 1, transcription factor, and Serpin Family E Member 1 as potential ICH biomarkers. Single sample gene set enrichment analysis revealed increased Type 2 T helper cell 2_cells, Treg cells, and immune functions like Antigen-presenting cells_co_stimulation in patients with ICH. Fourteen microRNA, including has-miR-6728-3p, were predicted to regulate these biomarkers. transcription factors such as PPARG, RARA, HMGA1, and NFKB1 were identified as potential regulators of the ICH biomarkers.

CONCLUSION

Caveolin 1, C-X-C Motif Chemokine Ligand 1, ETS proto-oncogene 1, transcription factor, and Serpin Family E Member 1 may serve as valuable biomarkers in ICH. Targeting these genes could contribute to ICH prevention and treatment.

Transcription factor binding site divergence across maize inbred lines drives transcriptional and phenotypic variation

Regulatory elements are important constituents of plant genomes that have shaped ancient and modern crops. Their identification, function, and diversity in crop genomes however are poorly characterized, thus limiting our ability to harness their power for further agricultural advances using induced or natural variation. Here, we use DNA affinity purification-sequencing (DAP-seq) to map transcription factor (TF) binding events for 200 maize TFs belonging to 30 distinct families and heterodimer pairs in two distinct inbred lines historically used for maize hybrid plant production, providing empirical binding site annotation for 5.3% of the maize genome. TF binding site comparison in B73 and Mo17 inbreds reveals widespread differences, driven largely by structural variation, that correlate with gene expression changes. TF binding site presence-absence variation helps clarify complex QTL such as vgt1, an important determinant of maize flowering time, and DICE, a distal enhancer involved in herbivore resistance. Modification of TF binding regions via CRISPR-Cas9 mediated editing alters target gene expression and phenotype. Our functional catalog of maize TF binding events enables collective and comparative TF binding analysis, and highlights its value for agricultural improvement.

Genomic Exploration: Unraveling the Intricacies of Indica Rice Oryza sativa L. Germin-Like Protein Gene 12-3 (OsGLP12-3) Promoter via Cloning, Sequencing, and In Silico Analysis

Germin and Germin-like proteins (GLPs) are a class of plant proteins that are part of the Cupins superfamily, found in several plant organs including roots, seeds, leaves, and nectar glands. They play a crucial role in plant defense against pathogens and environmental stresses. Herein, this study focused on the promoter analysis of OsGLP12-3 in rice cultivar Swat-1 to elucidate its regulation and functions. The region (1863bp) of the OsGLP12-3 promoter from Swat-1 genomic DNA was amplified, purified, quantified, and cloned using Topo cloning technology, followed by sequencing. Further in silico comparative analysis was conducted between the OsGLP12-3 promoters from Nipponbare and Swat-1 using the Plant CARE database, identifying 24 cis-acting regulatory elements with diverse functions. These elements exhibited distinct distribution patterns in the 2 rice varieties. The OsGLP12-3 promoter revealed an abundance of regulatory elements associated with biotic and abiotic stress responses. Computational tools were employed to analyze the regulatory features of this region. In silico expression analysis of OsGLP12-3, considering various developmental stages, stress conditions, hormones, and expression timing, was performed using the TENOR tool. Pairwise alignment indicated 86% sequence similarity between Nipponbare and Swat-1. Phylogenetic analysis was conducted to explore the evolutionary relationship between the OsGLP12-3 and other plant GLPs. Additionally, 2 unique regulatory elements were modeled and docked, GARE and MBS to understand their hydrogen bonding interactions in gene regulation. The study highlights the importance of OsGLP12-3 in plant defense against biotic and abiotic stresses, supported by its expression patterns in response to various stressors and the presence of specific regulatory elements within its promoter region.

Comparative analysis of PLATZ transcription factors in six poplar species and analysis of the role of PtrPLATZ14 in leaf development

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are a class of plant-specific transcription factor that play a crucial role in plant growth, development, and stress response. However, the evolutionary relationship of the PLATZ gene family across the Populus genus and the biological functions of the PLATZ protein require further investigation. In this study, we identified 133 PLATZ genes from six Populus species belonging to four Populus sections. Synteny analysis of the PLATZ gene family indicated that whole genome duplication events contributed to the expansion of the PLATZ family. Among the nine paralogous pairs, the protein structure of PtrPLATZ14/18 pair exhibited significant differences with others. Through gene expression patterns and co-expression networks, we discovered divergent expression patterns and sub-networks, and found that the members of pair PtrPLATZ14/18 might play different roles in the regulation of macromolecule biosynthesis and modification. Furthermore, we found that PtrPLATZ14 regulates poplar leaf development by affecting cell size control genes PtrGRF/GIF and PtrTCP. In conclusion, our study provides a theoretical foundation for exploring the evolution relationships and functions of the PLATZ gene family within Populus species and provides insights into the function and potential mechanism of PtrPLATZ14 in leaf morphology that were diverse across the Populus genus.

Genome-wide analysis of plant specific YABBY transcription factor gene family in carrot (Dacus carota) and its comparison with Arabidopsis

YABBY gene family is a plant-specific transcription factor with DNA binding domain involved in various functions i.e. regulation of style, length of flowers, and polarity development of lateral organs in flowering plants. Computational methods were utilized to identify members of the YABBY gene family, with Carrot (Daucus carota) 's genome as a foundational reference. The structure of genes, location of the chromosomes, protein motifs and phylogenetic investigation, syntony and transcriptomic analysis, and miRNA targets were analyzed to unmask the hidden structural and functional characteristics YABBY gene family in Carrots. In the following research, it has been concluded that 11 specific YABBY genes irregularly dispersed on all 9 chromosomes and proteins assembled into five subgroups i.e. AtINO, AtCRC, AtYAB5, AtAFO, and AtYAB2, which were created on the well-known classification of Arabidopsis. The wide ranges of YABBY genes in carrots were dispersed due to segmental duplication, which was detected as prevalent when equated to tandem duplication. Transcriptomic analysis showed that one of the DcYABBY genes was highly expressed during anthocyanin pigmentation in carrot taproots. The cis-regulatory elements (CREs) analysis unveiled elements that particularly respond to light, cell cycle regulation, drought induce ability, ABA hormone, seed, and meristem expression. Furthermore, a relative study among Carrot and Arabidopsis genes of the YABBY family indicated 5 sub-families sharing common characteristics. The comprehensive evaluation of YABBY genes in the genome provides a direction for the cloning and understanding of their functional properties in carrots. Our investigations revealed genome-wide distribution and role of YABBY genes in the carrots with best-fit comparison to Arabidopsis thaliana.

Comprehensive in silico characterization of NAC transcription factor family of Pinellia ternata and functional analysis of PtNAC66 under high-temperature tolerance in transgenic Arabidopsis thaliana

Pinellia ternata, a valuable Chinese herb, suffers yield reduction due to "sprout tumble" under high temperatures. However, the mechanisms underlying its high-temperature stress remain poorly understood. NAM, ATAF1/2, and CUC2 (NAC) transcription factors regulate plant tissue growth and abiotic stress. Hence, there has been no comprehensive research conducted on NAC transcription factors in P. ternata. We identified 98 PtNAC genes unevenly distributed across 13 chromosomes, grouped into 15 families via phylogenetic analysis. Gene expression analysis revealed diverse expression patterns of PtNAC genes in different tissue types. Further studies revealed that PtNAC5/7/17/35/43/47/57/66/86 genes were highly expressed in various tissues of P. ternata and induced by heat stress, among which PtNAC66 was up-regulated at the highest folds induced by heat temperature. PtNAC66 is a nuclear protein that can selectively bind to the cis-responsive region NACRS but lacks the ability to activate transcription in yeast. For further research, PtNAC66 was cloned and transgenic Arabidopsis was obtained. PtNAC66 overexpression increased high-temperature tolerance compared to wild-type plants. Transcriptome profiling demonstrated that overexpression of PtNAC66 led to significant modification of genes responsible for regulating binding, catalytic activity, transcription regulator activity and transporter activity response genes. Additionally, PtNAC66 was found to bind the promoters of CYP707A3, MYB102 and NAC055, respectively, and inhibited their expression, affecting the high-temperature stress response in Arabidopsis. Our research established the foundation for functional studies of PtNAC genes in response to high-temperature forcing by characterizing the P. ternata NAC gene family and examining the biological role of PtNAC66 in plant high-temperature tolerance.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

Extreme restructuring of cis -regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis -regulatory sequences. It remains unclear how such drastic cis -regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis- regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ∼125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 ( CLV3 ) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis -regulatory sequences. Our results demonstrate remarkable malleability in cis -regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis -regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis -regulation to effectively engineer trait variation from conserved productivity genes in crops.

Author summary

We investigated the evolution of cis -regulatory elements (CREs) and their interactions in the regulation of a plant stem cell regulator gene, CLAVATA3 (CLV3) , in Arabidopsis and tomato. Despite diverging ∼125 million years ago, the function and expression of CLV3 is conserved in these species; however, cis -regulatory sequences upstream and downstream have drastically diverged, preventing identification of conserved non-coding sequences between them. We used CRISPR-Cas9 to engineer dozens of mutations within the cis -regulatory regions of Arabidopsis and tomato CLV3. In tomato, our results show that tomato CLV3 function primarily relies on interactions among CREs in the 5' non-coding region, unlike Arabidopsis CLV3 , which depends on a more balanced distribution of functional CREs between the 5' and 3' regions. Therefore, despite a high degree of functional conservation, our study demonstrates divergent regulatory strategies between two distantly related CLV3 orthologs, with substantial alterations in regulatory sequences, their spatial arrangement, and their relative effects on CLV3 regulation. These results suggest that regulatory regions are not only extremely robust to mutagenesis, but also that the sequences underlying this robustness can be lineage-specific for conserved genes, due to the complex and often redundant interactions among CREs that ensure proper gene function amidst large-scale sequence turnover.

Sea-ATI unravels novel vocabularies of plant active cistrome

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants.

Whole-genome identification and multiple abiotic stresses expression pattern profiling analysis of PLATZ transcription factor family members in Pecan (Carya illinoensis)

Plant AT-rich sequence and zinc-binding (PLATZ), as a plant-specific transcription factor, have been identified and studied in a variety of plants. However, there are no reports about PLATZ proteins in Carya illinoensis (pecan). Here, 24 C. illinoensis CiPLATZs have been identified and divided into 4 groups. Gene structure, motif composition, conserved domain and cis-acting elements analysis indicated that the PLATZ gene family was highly conserved. Transcriptome data combined with qRT-PCR analysis revealed that CiPLATZ6, CiPLATZ12, CiPLATZ13, CiPLATZ14 and CiPLATZ23 were highly expressed in multiple tissues of C. illinoensis and strongly responded to drought, salt and heat stress. Among them, CiPLATZ6, CiPLATZ12 and CiPLATZ23 were all located in the nucleus and had no transcriptional autoactivation ability in yeast cells, and acted as transcriptional suppressors in plants. In addition, the CiPLATZ23-overexpressing transgenic Arabidopsis thaliana showed enhanced tolerance to drought. Measurements of physiological indicators and analysis of stress-related genes expression levels in transgenic A. thaliana were used to support this conclusion. The results of this study are helpful to understand the structural feature and function of CiPLATZs, and provide candidate genes for molecular breeding of drought tolerance of C. illinoensis.

Multi-Omics Pipeline and Omics-Integration Approach to Decipher Plant's Abiotic Stress Tolerance Responses

The present day's ongoing global warming and climate change adversely affect plants through imposing environmental (abiotic) stresses and disease pressure. The major abiotic factors such as drought, heat, cold, salinity, etc., hamper a plant's innate growth and development, resulting in reduced yield and quality, with the possibility of undesired traits. In the 21st century, the advent of high-throughput sequencing tools, state-of-the-art biotechnological techniques and bioinformatic analyzing pipelines led to the easy characterization of plant traits for abiotic stress response and tolerance mechanisms by applying the 'omics' toolbox. Panomics pipeline including genomics, transcriptomics, proteomics, metabolomics, epigenomics, proteogenomics, interactomics, ionomics, phenomics, etc., have become very handy nowadays. This is important to produce climate-smart future crops with a proper understanding of the molecular mechanisms of abiotic stress responses by the plant's genes, transcripts, proteins, epigenome, cellular metabolic circuits and resultant phenotype. Instead of mono-omics, two or more (hence 'multi-omics') integrated-omics approaches can decipher the plant's abiotic stress tolerance response very well. Multi-omics-characterized plants can be used as potent genetic resources to incorporate into the future breeding program. For the practical utility of crop improvement, multi-omics approaches for particular abiotic stress tolerance can be combined with genome-assisted breeding (GAB) by being pyramided with improved crop yield, food quality and associated agronomic traits and can open a new era of omics-assisted breeding. Thus, multi-omics pipelines together are able to decipher molecular processes, biomarkers, targets for genetic engineering, regulatory networks and precision agriculture solutions for a crop's variable abiotic stress tolerance to ensure food security under changing environmental circumstances.

Comprehensive analysis of HSF genes from celery (Apium graveolens L.) and functional characterization of AgHSFa6-1 in response to heat stress

High temperature stress is regarded as one of the significant abiotic stresses affecting the composition and distribution of natural habitats and the productivity of agriculturally significant plants worldwide. The HSF family is one of the most important transcription factors (TFs) families in plants and capable of responding rapidly to heat and other abiotic stresses. In this study, 29 AgHSFs were identified in celery and classified into three classes (A, B, and C) and 14 subgroups. The gene structures of AgHSFs in same subgroups were conserved, whereas in different classes were varied. AgHSF proteins were predicted to be involved in multiple biological processes by interacting with other proteins. Expression analysis revealed that AgHSF genes play a significant role in response to heat stress. Subsequently, AgHSFa6-1, which was significantly induced by high temperature, was selected for functional validation. AgHSFa6-1 was identified as a nuclear protein, and can upregulate the expression of certain downstream genes (HSP98.7, HSP70-1, BOB1, CPN60B, ADH2, APX1, GOLS1) in response to high-temperature treatment. Overexpression of AgHSFa6-1 in yeast and Arabidopsis displayed higher thermotolerance, both morphologically and physiologically. In response to heat stress, the transgenic plants produced considerably more proline, solute protein, antioxidant enzymes, and less MDA than wild-type (WT) plants. Overall, this study revealed that AgHSF family members perform a key role in response to high temperature, and AgHSFa6-1 acts as a positive regulator by augmenting the ROS-scavenging system to maintain membrane integrity, reducing stomatal apertures to control water loss, and upregulating the expression level of heat-stress sensitive genes to improve celery thermotolerance.

Genome-wide identification of PLATZ genes related to cadmium tolerance in Populus trichocarpa and characterization of the role of PtPLATZ3 in phytoremediation of cadmium

Plant AT-rich sequence and zinc-binding (PLATZ) proteins are a class of plant-specific zinc finger transcription factors that perform critical functions in plant development and resistance. However, the function of PLATZs in heavy metal tolerance has not yet been investigated. Moreover, only a few PLATZ proteins have been functionally characterized in tree species. In this study, we identified 18 PtPLATZ genes in Populus trichocarpa, an important woody model plant, and classified them into five groups. PtPLATZ genes attributed to the same clade usually possess similar exon-intron structures containing two or three introns, as well as a similar motif composition. Furthermore, chromosomal location analysis indicated an uneven distribution of PtPLATZ genes on 13 of the 19 Populus chromosomes. Promoter cis-acting element prediction and gene expression analysis showed that PtPLATZ genes were highly responsive to heavy metal stress. Heterologous yeast expression revealed that PtPLATZ1, PtPLATZ2, PtPLATZ3, PtPLATZ4, PtPLATZ8 and PtPLATZ9 are significantly involved in Cd tolerance. In addition, transgenic expression of PtPLATZ3 significantly enhanced Cd tolerance and accumulation, slowed the decline in chlorophyll content, maintained membrane integrity in Populus, and increased the expression of genes related to Cd tolerance and accumulation. In conclusion, our results suggest the potential of PtPLATZ3 to improve Cd tolerance and accumulation in Populus, which is of great significance for phytoremediation.

Cis-Regulation by NACs: A Promising Frontier in Wheat Crop Improvement

Crop traits are controlled by multiple genes; however, the complex spatio-temporal transcriptional behavior of genes cannot be fully understood without comprehending the role of transcription factors (TFs) and the underlying mechanisms of the binding interactions of their cis-regulatory elements. NAC belongs to one of the largest families of plant-specific TFs and has been associated with the regulation of many traits. This review provides insight into the cis-regulation of genes by wheat NACs (TaNACs) for the improvement in yield-related traits, including phytohormonal homeostasis, leaf senescence, seed traits improvement, root modulation, and biotic and abiotic stresses in wheat and other cereals. We also discussed the current potential, knowledge gaps, and prospects of TaNACs.

A Novel Strategy for Identifying NSCLC MicroRNA Biomarkers and Their Mechanism Analysis Based on a Brand-New CeRNA-Hub-FFL Network

Finding reliable miRNA markers and revealing their potential mechanisms will play an important role in the diagnosis and treatment of NSCLC. Most existing computational methods for identifying miRNA biomarkers only consider the expression variation of miRNAs or rely heavily on training sets. These deficiencies lead to high false-positive rates. The independent regulatory model is an important complement to traditional models of co-regulation and is more impervious to the dataset. In addition, previous studies of miRNA mechanisms in the development of non-small cell lung cancer (NSCLC) have mostly focused on the post-transcriptional level and did not distinguish between NSCLC subtypes. For the above problems, we improved mainly in two areas: miRNA identification based on both the NOG network and biological functions of miRNA target genes; and the construction of a 4-node directed competitive regulatory network to illustrate the mechanisms. NSCLC was classified as lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) in this work. One miRNA biomarker of LUAD (miR-708-5p) and four of LUSC (miR-183-5p, miR-140-5p, miR-766-5p, and miR-766-3p) were obtained. They were validated using literature and external datasets. The ceRNA-hub-FFL involving transcription factors (TFs), microRNAs (miRNAs), mRNAs, and long non-coding RNAs (lncRNAs) was constructed. There were multiple interactions among these components within the net at the transcriptional, post-transcriptional, and protein levels. New regulations were revealed by the network. Meanwhile, the network revealed the reasons for the previous conflicting conclusions on the roles of CD44, ACTB, and ITGB1 in NSCLC, and demonstrated the necessity of typing studies on NSCLC. The novel miRNA markers screening method and the 4-node directed competitive ceRNA-hub-FFL network constructed in this work can provide new ideas for screening tumor markers and understanding tumor development mechanisms in depth.