Novel tuberculosis susceptibility candidate genes revealed by the reconstruction and analysis of associative networks

An Accurate and Efficient Approach to Knowledge Extraction from Scientific Publications Using Structured Ontology Models, Graph Neural Networks, and Large Language Models

The rapid growth of biomedical literature makes it challenging for researchers to stay current. Integrating knowledge from various sources is crucial for studying complex biological systems. Traditional text-mining methods often have limited accuracy because they don't capture semantic and contextual nuances. Deep-learning models can be computationally expensive and typically have low interpretability, though efforts in explainable AI aim to mitigate this. Furthermore, transformer-based models have a tendency to produce false or made-up information-a problem known as hallucination-which is especially prevalent in large language models (LLMs). This study proposes a hybrid approach combining text-mining techniques with graph neural networks (GNNs) and fine-tuned large language models (LLMs) to extend biomedical knowledge graphs and interpret predicted edges based on published literature. An LLM is used to validate predictions and provide explanations. Evaluated on a corpus of experimentally confirmed protein interactions, the approach achieved a Matthews correlation coefficient (MCC) of 0.772. Applied to insomnia, the approach identified 25 interactions between 32 human proteins absent in known knowledge bases, including regulatory interactions between MAOA and 5-HT2C, binding between ADAM22 and 14-3-3 proteins, which is implicated in neurological diseases, and a circadian regulatory loop involving RORB and NR1D1. The hybrid GNN-LLM method analyzes biomedical literature efficiency to uncover potential molecular interactions for complex disorders. It can accelerate therapeutic target discovery by focusing expert verification on the most relevant automatically extracted information.

Microbiome Alteration in Lung Tissues of Tuberculosis Patients Revealed by Metagenomic Next-Generation Sequencing and Immune-Related Transcriptional Profile Identified by Transcriptome Sequencing

This study explored alterations in the respiratory microbiome and transcriptome after Mycobacterium tuberculosis infection in tuberculosis (TB) patients. Metagenomic next-generation sequencing (mNGS) was adopted to reveal the microbiome in lung tissues from 110 TB and 25 nontuberculous (NonTB) patients. Transcriptome sequencing was performed in TB tissues (n = 3), tissues adjacent to TB (ParaTB, n = 3), and NonTB tissues (n = 3) to analyze differentially expressed genes (DEGs) and functional pathways. The microbial β diversity (p = 0.01325) in TB patients differed from that in the NonTB group, with 17 microbial species distinctively distributed. Eighty-three co-up-regulated DEGs were identified in the TB versus NonTB and the TB versus ParaTB comparison groups, and six were associated with immune response to Mtb. These DEGs were significantly enriched in the signaling pathways such as immune response, NF-κB, and B cell receptor. Data in the lung tissue microbiome and transcriptome in TB patients offer a sufficient understanding of the pathogenesis of TB.

Reconstruction of the regulatory hypermethylation network controlling hepatocellular carcinoma development during hepatitis C viral infection

Hepatocellular carcinoma (HCC) has been associated with hepatitis C viral (HCV) infection as a potential risk factor. Nonetheless, the precise genetic regulatory mechanisms triggered by the virus, leading to virus-induced hepatocarcinogenesis, remain unclear. We hypothesized that HCV proteins might modulate the activity of aberrantly methylated HCC genes through regulatory pathways. Virus-host regulatory pathways, interactions between proteins, gene expression, transport, and stability regulation, were reconstructed using the ANDSystem. Gene expression regulation was statistically significant. Gene network analysis identified four out of 70 HCC marker genes whose expression regulation by viral proteins may be associated with HCC: DNA-binding protein inhibitor ID - 1 (ID1), flap endonuclease 1 (FEN1), cyclin-dependent kinase inhibitor 2A (CDKN2A), and telomerase reverse transcriptase (TERT). It suggested the following viral protein effects in HCV/human protein heterocomplexes: HCV NS3(p70) protein activates human STAT3 and NOTC1; NS2-3(p23), NS5B(p68), NS1(E2), and core(p21) activate SETD2; NS5A inhibits SMYD3; and NS3 inhibits CCN2. Interestingly, NS3 and E1(gp32) activate c-Jun when it positively regulates CDKN2A and inhibit it when it represses TERT. The discovered regulatory mechanisms might be key areas of focus for creating medications and preventative therapies to decrease the likelihood of HCC development during HCV infection.

Human alveolar macrophage response to Mycobacterium tuberculosis: immune characteristics underlying large inter-individual variability

Background

Mycobacterium tuberculosis (M.tb), the causative bacterium of tuberculosis (TB), establishes residence and grows in human alveolar macrophages (AMs). Inter-individual variation in M.tb-human AM interactions can indicate TB risk and the efficacy of therapies and vaccines; however, we currently lack an understanding of the gene and protein expression programs that dictate this variation in the lungs.

Results

Herein, we systematically analyze interactions of a virulent M.tb strain H37Rv with freshly isolated human AMs from 28 healthy adult donors, measuring host RNA expression and secreted candidate proteins associated with TB pathogenesis over 72h. A large set of genes possessing highly variable inter-individual expression levels are differentially expressed in response to M.tb infection. Eigengene modules link M.tb growth rate with host transcriptional and protein profiles at 24 and 72h. Systems analysis of differential RNA and protein expression identifies a robust network with IL1B, STAT1, and IDO1 as hub genes associated with M.tb growth. RNA time profiles document stimulation towards an M1-type macrophage gene expression followed by emergence of an M2-type profile. Finally, we replicate these results in a cohort from a TB-endemic region, finding a substantial portion of significant differentially expressed genes overlapping between studies.

Conclusions

We observe large inter-individual differences in bacterial uptake and growth, with tenfold variation in M.tb load by 72h.The fine-scale resolution of this work enables the identification of genes and gene networks associated with early M.tb growth dynamics in defined donor clusters, an important step in developing potential biological indicators of individual susceptibility to M.tb infection and response to therapies.

Apoptosis Genes as a Key to Identification of Inverse Comorbidity of Huntington's Disease and Cancer

Cancer and neurodegenerative disorders present overwhelming challenges for healthcare worldwide. Epidemiological studies showed a decrease in cancer rates in patients with neurodegenerative disorders, including the Huntington disease (HD). Apoptosis is one of the most important processes for both cancer and neurodegeneration. We suggest that genes closely connected with apoptosis and associated with HD may affect carcinogenesis. We applied reconstruction and analysis of gene networks associated with HD and apoptosis and identified potentially important genes for inverse comorbidity of cancer and HD. The top 10 high-priority candidate genes included APOE, PSEN1, INS, IL6, SQSTM1, SP1, HTT, LEP, HSPA4, and BDNF. Functional analysis of these genes was carried out using gene ontology and KEGG pathways. By exploring genome-wide association study results, we identified genes associated with neurodegenerative and oncological disorders, as well as their endophenotypes and risk factors. We used publicly available datasets of HD and breast and prostate cancers to analyze the expression of the identified genes. Functional modules of these genes were characterized according to disease-specific tissues. This integrative approach revealed that these genes predominantly exert similar functions in different tissues. Apoptosis along with lipid metabolism dysregulation and cell homeostasis maintenance in the response to environmental stimulus and drugs are likely key processes in inverse comorbidity of cancer in patients with HD. Overall, the identified genes represent the promising targets for studying molecular relations of cancer and HD.

Primary and Secondary micro-RNA Modulation the Extrinsic Pathway of Apoptosis in Hepatocellular Carcinoma

Abstract-One of the most common malignant liver diseases is hepatocellular carcinoma, which has a high recurrence rate and a low five-year survival rate. It is very heterogeneous both in structure and between patients, which complicates the diagnosis, prognosis and response to treatment. In this regard, an individualized, patient-centered approach becomes important, in which the use of mimetics and hsa-miRNA inhibitors involved in the pathogenesis of the disease may be determinative. From this point of view hsa-miRNAs are of interest, their aberrant expression is associated with poor prognosis for patients and is associated with tumor progression due to dysregulation of programmed cell death (apoptosis). However, the effect of hsa-miRNA on tumor development depends not only on its direct effect on expression of genes, the primary targets, but also on secondary targets mediated by regulatory pathways. While the former are actively studied, the role of secondary targets of these hsa-miRNAs in modulating apoptosis is still unclear. The present work summarizes data on hsa-miRNAs whose primary targets are key genes of the extrinsic pathway of apoptosis. Their aberrant expression is associated with early disease relapse and poor patient outcome. For these hsa-miRNAs, using the software package ANDSystem, we reconstructed the regulation of the expression of secondary targets and analyzed their impact on the activity of the extrinsic pathway of apoptosis. The potential effect of hsa-miRNAs mediated by action on secondary targets is shown to negatively correlate with the number of primary targets. It is also shown that hsa-miR-373, hsa-miR-106b and hsa-miR-96 have the highest priority as markers of hepatocellular carcinoma, whose action on secondary targets enhances their anti-apoptotic effect.

Computer analysis of regulation of hepatocarcinoma marker genes hypermethylated by HCV proteins

Hepatitis C virus (HCV) is a risk factor that leads to hepatocellular carcinoma (HCC) development. Epigenetic changes are known to play an important role in the molecular genetic mechanisms of virus-induced oncogenesis. Aberrant DNA methylation is a mediator of epigenetic changes that are closely associated with the HCC pathogenesis and considered a biomarker for its early diagnosis. The ANDSystem software package was used to reconstruct and evaluate the statistical significance of the pathways HCV could potentially use to regulate 32 hypermethylated genes in HCC, including both oncosuppressor and protumorigenic ones identified by genome-wide analysis of DNA methylation. The reconstructed pathways included those affecting protein-protein interactions (PPI), gene expression, protein activity, stability, and transport regulations, the expression regulation pathways being statistically significant. It has been shown that 8 out of 10 HCV proteins were involved in these pathways, the HCV NS3 protein being implicated in the largest number of regulatory pathways. NS3 was associated with the regulation of 5 tumor-suppressor genes, which may be the evidence of its central role in HCC pathogenesis. Analysis of the reconstructed pathways has demonstrated that following the transcription factor inhibition caused by binding to viral proteins, the expression of a number of oncosuppressors (WT1, MGMT, SOCS1, P53) was suppressed, while the expression of others (RASF1, RUNX3, WIF1, DAPK1) was activated. Thus, the performed gene-network reconstruction has shown that HCV proteins can influence not only the methylation status of oncosuppressor genes, but also their transcriptional regulation. The results obtained can be used in the search for pharmacological targets to develop new drugs against HCV-induced HCC.

Plasma metabolomics and gene regulatory networks analysis reveal the role of nonstructural SARS-CoV-2 viral proteins in metabolic dysregulation in COVID-19 patients

Metabolomic analysis of blood plasma samples from COVID-19 patients is a promising approach allowing for the evaluation of disease progression. We performed the metabolomic analysis of plasma samples of 30 COVID-19 patients and the 19 controls using the high-performance liquid chromatography (HPLC) coupled with tandem mass spectrometric detection (LC-MS/MS). In our analysis, we identified 103 metabolites enriched in KEGG metabolic pathways such as amino acid metabolism and the biosynthesis of aminoacyl-tRNAs, which differed significantly between the COVID-19 patients and the controls. Using ANDSystem software, we performed the reconstruction of gene networks describing the potential genetic regulation of metabolic pathways perturbed in COVID-19 patients by SARS-CoV-2 proteins. The nonstructural proteins of SARS-CoV-2 (orf8 and nsp5) and structural protein E were involved in the greater number of regulatory pathways. The reconstructed gene networks suggest the hypotheses on the molecular mechanisms of virus-host interactions in COVID-19 pathology and provide a basis for the further experimental and computer studies of the regulation of metabolic pathways by SARS-CoV-2 proteins. Our metabolomic analysis suggests the need for nonstructural protein-based vaccines and the control strategy to reduce the disease progression of COVID-19.

Locally Adjust Networks Based on Connectivity and Semantic Similarities for Disease Module Detection

For studying the pathogenesis of complex diseases, it is important to identify the disease modules in the system level. Since the protein-protein interaction (PPI) networks contain a number of incomplete and incorrect interactome, most existing methods often lead to many disease proteins isolating from disease modules. In this paper, we propose an effective disease module identification method IDMCSS, where the used human PPI networks are obtained by adding some potential missing interactions from existing PPI networks, as well as removing some potential incorrect interactions. In IDMCSS, a network adjustment strategy is developed to add or remove links around disease proteins based on both topological and semantic information. Next, neighboring proteins of disease proteins are prioritized according to a suggested similarity between each of them and disease proteins, and the protein with the largest similarity with disease proteins is added into a candidate disease protein set one by one. The stopping criterion is set to the boundary of the disease proteins. Finally, the connected subnetwork having the largest number of disease proteins is selected as a disease module. Experimental results on asthma demonstrate the effectiveness of the method in comparison to existing algorithms for disease module identification. It is also shown that the proposed IDMCSS can obtain the disease modules having crucial biological processes of asthma and 12 targets for drug intervention can be predicted.

Functional Expression, Purification and Identification of Interaction Partners of PACRG

PACRG (Parkin co-regulated gene) shares a bi-directional promoter with the Parkinson’s disease-associated gene Parkin, but the physiological roles of PACRG have not yet been fully elucidated. Recombinant expression methods are indispensable for protein structural and functional studies. In this study, the coding region of PACRG was cloned to a conventional vector pQE80L, as well as two cold-shock vectors pCold II and pCold-GST, respectively. The constructs were transformed into Escherichia coli (DE3), and the target proteins were overexpressed. The results showed that the cold-shock vectors are more suitable for PACRG expression. The soluble recombinant proteins were purified with Ni²⁺ chelating column, glutathione S-transferase (GST) affinity chromatography and gel filtration. His₆ pull down assay and LC-MS/MS were carried out for identification of PACRG-binding proteins in HEK293T cell lysates, and a total number of 74 proteins were identified as potential interaction partners of PACRG. GO (Gene ontology) enrichment analysis (FunRich) of the 74 proteins revealed multiple molecular functions and biological processes. The highest proportion of the 74 proteins functioned as transcription regulator and transcription factor activity, suggesting that PACRG may play important roles in regulation of gene transcription.

ANDDigest: a new web-based module of ANDSystem for the search of knowledge in the scientific literature

BACKGROUND

The rapid growth of scientific literature has rendered the task of finding relevant information one of the critical problems in almost any research. Search engines, like Google Scholar, Web of Knowledge, PubMed, Scopus, and others, are highly effective in document search; however, they do not allow knowledge extraction. In contrast to the search engines, text-mining systems provide extraction of knowledge with representations in the form of semantic networks. Of particular interest are tools performing a full cycle of knowledge management and engineering, including automated retrieval, integration, and representation of knowledge in the form of semantic networks, their visualization, and analysis. STRING, Pathway Studio, MetaCore, and others are well-known examples of such products. Previously, we developed the Associative Network Discovery System (ANDSystem), which also implements such a cycle. However, the drawback of these systems is dependence on the employed ontologies describing the subject area, which limits their functionality in searching information based on user-specified queries.

RESULTS

The ANDDigest system is a new web-based module of the ANDSystem tool, permitting searching within PubMed by using dictionaries from the ANDSystem tool and sets of user-defined keywords. ANDDigest allows performing the search based on complex queries simultaneously, taking into account many types of objects from the ANDSystem's ontology. The system has a user-friendly interface, providing sorting, visualization, and filtering of the found information, including mapping of mentioned objects in text, linking to external databases, sorting of data by publication date, citations number, journal H-indices, etc. The system provides data on trends for identified entities based on dynamics of interest according to the frequency of their mentions in PubMed by years.

CONCLUSIONS

The main feature of ANDDigest is its functionality, serving as a specialized search for information about multiple associative relationships of objects from the ANDSystem's ontology vocabularies, taking into account user-specified keywords. The tool can be applied to the interpretation of experimental genetics data, the search for associations between molecular genetics objects, and the preparation of scientific and analytical reviews. It is presently available at https://anddigest.sysbio.ru/ .

Clinical relevance of the lnc-HNF1B-3:1 genetic polymorphisms in Western Chinese tuberculosis patients

BACKGROUND

Tuberculosis remains a global public health problem. Genetic polymorphisms may affect the susceptibility, clinical characteristics, and adverse drug reactions of patients with TB. The present study aimed to examine the association of single nucleotide polymorphisms of lncRNA-HNF1B-3:1 with the clinical manifestation of TB in a Western Chinese population.

METHOD

A total of 526 tuberculosis patients and 561 healthy subjects were recruited in Western China. The correlation between lnc-HNF1B-3:1 polymorphism and tuberculosis susceptibility was investigated. Moreover, the influence on adverse drug reactions following treatment was explored. A total of 7 SNPs within the lnc-HNF1B-3:1 locus was genotyped by the improved multiplex ligation detection reaction method.

RESULTS

No significant associations were noted between TB susceptibility and the presence of all 7 SNPs of the lnc-HNF1B-3:1 as determined by single-locus analysis (All P > .05). The AA genotype of rs12939622 (in the dominant model) and the AA genotype of rs4262994 (in the recessive model) caused increased susceptibility of the subjects to fever (P < .001 and P = .008, respectively). The Rs2542670 G allele was associated with increased risk of thrombocytopenia, leukopenia, and chronic kidney damage following drug administration (P = .007, .029, .003, respectively).

CONCLUSION

The present study reported for the first time that the rs12939622, rs4262994 and rs2542670 genotypes in lnc-HNF1B-3:1 locus may influence the clinical manifestations of tuberculosis.

A Survey of Gene Prioritization Tools for Mendelian and Complex Human Diseases

Modern high-throughput experiments provide us with numerous potential associations between genes and diseases. Experimental validation of all the discovered associations, let alone all the possible interactions between them, is time-consuming and expensive. To facilitate the discovery of causative genes, various approaches for prioritization of genes according to their relevance for a given disease have been developed. In this article, we explain the gene prioritization problem and provide an overview of computational tools for gene prioritization. Among about a hundred of published gene prioritization tools, we select and briefly describe 14 most up-to-date and user-friendly. Also, we discuss the advantages and disadvantages of existing tools, challenges of their validation, and the directions for future research.

Genome-Wide Association Study in Mexican Holstein Cattle Reveals Novel Quantitative Trait Loci Regions and Confirms Mapped Loci for Resistance to Bovine Tuberculosis

Bovine tuberculosis (bTB) is a disease of cattle that represents a risk to public health and causes severe economic losses to the livestock industry. Recently, genetic studies, like genome-wide association studies (GWAS) have greatly improved the investigation of complex diseases identifying thousands of disease-associated genomic variants. Here, we present evidence of genetic variants associated with resistance to TB in Mexican dairy cattle using a case-control approach with a selective DNA pooling experimental design. A total of 154 QTLRs (quantitative trait loci regions) at 10% PFP (proportion of false positives), 42 at 5% PFP and 5 at 1% PFP have been identified, which harbored 172 annotated genes. On BTA13, five new QTLRs were identified in the MACROD2 and KIF16B genes, supporting their involvement in resistance to bTB. Six QTLRs harbor seven annotated genes that have been previously reported as involved in immune response against Mycobacterium spp: BTA (Bos taurus autosome) 1 (CD80), BTA3 (CTSS), BTA 3 (FCGR1A), BTA 23 (HFE), BTA 25 (IL21R), and BTA 29 (ANO9 and SIGIRR). We identified novel QTLRs harboring genes involved in Mycobacterium spp. immune response. This is a first screening for resistance to TB infection on Mexican dairy cattle based on a dense SNP (Single Nucleotide Polymorphism) chip.

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.

Sequence variant analysis of RNA sequences in severe equine asthma

Background

Severe equine asthma is a chronic inflammatory disease of the lung in horses similar to low-Th2 late-onset asthma in humans. This study aimed to determine the utility of RNA-Seq to call gene sequence variants, and to identify sequence variants of potential relevance to the pathogenesis of asthma.

Methods

RNA-Seq data were generated from endobronchial biopsies collected from six asthmatic and seven non-asthmatic horses before and after challenge (26 samples total). Sequences were aligned to the equine genome with Spliced Transcripts Alignment to Reference software. Read preparation for sequence variant calling was performed with Picard tools and Genome Analysis Toolkit (GATK). Sequence variants were called and filtered using GATK and Ensembl Variant Effect Predictor (VEP) tools, and two RNA-Seq predicted sequence variants were investigated with both PCR and Sanger sequencing. Supplementary analysis of novel sequence variant selection with VEP was based on a score of <0.01 predicted with Sorting Intolerant from Tolerant software, missense nature, location within the protein coding sequence and presence in all asthmatic individuals. For select variants, effect on protein function was assessed with Polymorphism Phenotyping 2 and screening for non-acceptable polymorphism 2 software. Sequences were aligned and 3D protein structures predicted with Geneious software. Difference in allele frequency between the groups was assessed using a Pearson’s Chi-squared test with Yates’ continuity correction, and difference in genotype frequency was calculated using the Fisher’s exact test for count data.

Results

RNA-Seq variant calling and filtering correctly identified substitution variants in PACRG and RTTN. Sanger sequencing confirmed that the PACRG substitution was appropriately identified in all 26 samples while the RTTN substitution was identified correctly in 24 of 26 samples. These variants of uncertain significance had substitutions that were predicted to result in loss of function and to be non-neutral. Amino acid substitutions projected no change of hydrophobicity and isoelectric point in PACRG, and a change in both for RTTN. For PACRG, no difference in allele frequency between the two groups was detected but a higher proportion of asthmatic horses had the altered RTTN allele compared to non-asthmatic animals.

Discussion

RNA-Seq was sensitive and specific for calling gene sequence variants in this disease model. Even moderate coverage (<10–20 counts per million) yielded correct identification in 92% of samples, suggesting RNA-Seq may be suitable to detect sequence variants in low coverage samples. The impact of amino acid alterations in PACRG and RTTN proteins, and possible association of the sequence variants with asthma, is of uncertain significance, but their role in ciliary function may be of future interest.

Genome-wide association study identifies loci associated with milk leukocyte phenotypes following experimental challenge with Streptococcus uberis

Mastitis is a detrimental disease in the dairy industry that decreases milk quality and costs upwards of $2 billion annually. Often, mastitis results from bacteria entering the gland through the teat opening. Streptococcus uberis is responsible for a high percentage of subclinical and clinical mastitis. Following an intramammary experimental challenge with S. uberis on Holstein cows (n = 40), milk samples were collected and somatic cell counts (SCC) were determined by the Dairy Herd Improvement Association Laboratory. Traditional genome-wide association studies (GWAS) have utilized test day SCC or SCC lactation averages to identify loci of interest. Our approach utilizes SCC collected following a S. uberis experimental challenge to generate three novel phenotypes: (1) area under the curve (AUC) of SCC for 0-7 days and (2) 0-28 days post-challenge; and (3) when SCC returned to below 200,000 cells/mL post-challenge (< 21 days, 21-28 days, or > 28 days). Polymorphisms were identified using Illumina's BovineSNP50 v2 DNA BeadChip. Associations were tested using Plink software and identified 16 significant (p < 1.0 × 10^-4) single-nucleotide polymorphisms (SNPs) across the phenotypes. Most significant SNPs were in genes linked to cell signaling, migration, and apoptosis. Several have been recognized in relation to infectious processes (ATF7, SGK1, and PACRG), but others less so (TRIO, GLRA1, CELSR2, TIAM2, CPE). Further investigation of these genes and their roles in inflammation (e.g., SCC) can provide potential targets that influence resolution of mammary gland infection. Likewise, further investigation of the identified SNP with mastitis and other disease phenotypes can provide greater insight to the potential of these SNP as genetic markers.

Novel candidate genes important for asthma and hypertension comorbidity revealed from associative gene networks

BACKGROUND

Hypertension and bronchial asthma are a major issue for people's health. As of 2014, approximately one billion adults, or ~ 22% of the world population, have had hypertension. As of 2011, 235-330 million people globally have been affected by asthma and approximately 250,000-345,000 people have died each year from the disease. The development of the effective treatment therapies against these diseases is complicated by their comorbidity features. This is often a major problem in diagnosis and their treatment. Hence, in this study the bioinformatical methodology for the analysis of the comorbidity of these two diseases have been developed. As such, the search for candidate genes related to the comorbid conditions of asthma and hypertension can help in elucidating the molecular mechanisms underlying the comorbid condition of these two diseases, and can also be useful for genotyping and identifying new drug targets.

RESULTS

Using ANDSystem, the reconstruction and analysis of gene networks associated with asthma and hypertension was carried out. The gene network of asthma included 755 genes/proteins and 62,603 interactions, while the gene network of hypertension - 713 genes/proteins and 45,479 interactions. Two hundred and five genes/proteins and 9638 interactions were shared between asthma and hypertension. An approach for ranking genes implicated in the comorbid condition of two diseases was proposed. The approach is based on nine criteria for ranking genes by their importance, including standard methods of gene prioritization (Endeavor, ToppGene) as well as original criteria that take into account the characteristics of an associative gene network and the presence of known polymorphisms in the analysed genes. According to the proposed approach, the genes IL10, TLR4, and CAT had the highest priority in the development of comorbidity of these two diseases. Additionally, it was revealed that the list of top genes is enriched with apoptotic genes and genes involved in biological processes related to the functioning of central nervous system.

CONCLUSIONS

The application of methods of reconstruction and analysis of gene networks is a productive tool for studying the molecular mechanisms of comorbid conditions. The method put forth to rank genes by their importance to the comorbid condition of asthma and hypertension was employed that resulted in prediction of 10 genes, playing the key role in the development of the comorbid condition. The results can be utilised to plan experiments for identification of novel candidate genes along with searching for novel pharmacological targets.

FunGeneNet: a web tool to estimate enrichment of functional interactions in experimental gene sets

Background

Estimation of functional connectivity in gene sets derived from genome-wide or other biological experiments is one of the essential tasks of bioinformatics. A promising approach for solving this problem is to compare gene networks built using experimental gene sets with random networks. One of the resources that make such an analysis possible is CrossTalkZ, which uses the FunCoup database. However, existing methods, including CrossTalkZ, do not take into account individual types of interactions, such as protein/protein interactions, expression regulation, transport regulation, catalytic reactions, etc., but rather work with generalized types characterizing the existence of any connection between network members.

Results

We developed the online tool FunGeneNet, which utilizes the ANDSystem and STRING to reconstruct gene networks using experimental gene sets and to estimate their difference from random networks. To compare the reconstructed networks with random ones, the node permutation algorithm implemented in CrossTalkZ was taken as a basis. To study the FunGeneNet applicability, the functional connectivity analysis of networks constructed for gene sets involved in the Gene Ontology biological processes was conducted. We showed that the method sensitivity exceeds 0.8 at a specificity of 0.95. We found that the significance level of the difference between gene networks of biological processes and random networks is determined by the type of connections considered between objects. At the same time, the highest reliability is achieved for the generalized form of connections that takes into account all the individual types of connections. By taking examples of the thyroid cancer networks and the apoptosis network, it is demonstrated that key participants in these processes are involved in the interactions of those types by which these networks differ from random ones.

Conclusions

FunGeneNet is a web tool aimed at proving the functionality of networks in a wide range of sizes of experimental gene sets, both for different global networks and for different types of interactions. Using examples of thyroid cancer and apoptosis networks, we have shown that the links over-represented in the analyzed network in comparison with the random ones make possible a biological interpretation of the original gene/protein sets. The FunGeneNet web tool for assessment of the functional enrichment of networks is available at http://www-bionet.sscc.ru/fungenenet/.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4474-7) contains supplementary material, which is available to authorized users.