Genomics and Informatics最新文献

Diabetes and its related complications are associated with long term damage and failure of various organ systems. The microvascular complications of diabetes considered in this study are diabetic retinopathy, diabetic neuropathy, and diabetic nephropathy. The aim is to identify the weighted co-expressed and differentially expressed genes (DEGs), major pathways, and their miRNA, transcription factors (TFs) and drugs interacting in all the three conditions. The primary goal is to identify vital DEGs in all the three conditions. The overlapped five genes (AKT1, NFKB1, MAPK3, PDPK1, and TNF) from the DEGs and the co-expressed genes were defined as key genes, which differentially expressed in all the three cases. Then the protein-protein interaction network and gene set linkage analysis (GSLA) of key genes was performed. GSLA, gene ontology, and pathway enrichment analysis of the key genes elucidates nine major pathways in diabetes. Subsequently, we constructed the miRNA-gene and transcription factorgene regulatory network of the five gene of interest in the nine major pathways were studied. hsa-mir-34a-5p, a major miRNA that interacted with all the five genes. RELA, FOXO3, PDX1 and SREBF1 were the TFs interacting with the major five gene of interest. Finally, drug-gene interaction network elucidates five potential drugs to treat the genes of interest. This research reveals biomarker genes, miRNA, TFs, and therapeutic drugs in the key signaling pathways, which may help us, understand the processes of all three secondary microvascular problems and aid in disease detection and management.

Alfalfa (Medicago sativa) is an important food and feed crop which rich in mineral sources. The WUSCHEL-related homeobox (WOX) gene family plays important roles in plant development and identification of putative gene families, their structure, and potential functions is a primary step for not only understanding the genetic mechanisms behind various biological process but also for genetic improvement. A variety of computational tools, including MAFFT, HMMER, hidden Markov models, Pfam, SMART, MEGA, ProtTest, BLASTn, and BRAD, among others, were used. We identified 34 MsWOX genes based on a systematic analysis of the alfalfa plant genome spread in eight chromosomes. This is an expansion of the gene family which we attribute to observed chromosomal duplications. Sequence alignment analysis revealed 61 conserved proteins containing a homeodomain. Phylogenetic study sung reveal five evolutionary clades with 15 motif distributions. Gene structure analysis reveals various exon, intron, and untranslated structures which are consistent in genes from similar clades. Functional analysis prediction of promoter regions reveals various transcription binding sites containing key growth, development, and stress-responsive transcription factor families such as MYB, ERF, AP2, and NAC which are spread across the genes. Most of the genes are predicted to be in the nucleus. Also, there are duplication events in some genes which explain the expansion of the family. The present research provides a clue on the potential roles of MsWOX family genes that will be useful for further understanding their functional roles in alfalfa plants.

Recent studies have focused on the early detection of ovarian cancer (OC) using tumor materials by liquid biopsy. The mechanisms of microRNAs (miRNAs) to impact OC and signaling pathways are still unknown. This study aims to reliably perform functional analysis of previously validated circulating miRNAs' target genes by using pathfindR. Also, overall survival and pathological stage analyses were evaluated with miRNAs' target genes which are common in the The Cancer Genome Atlas and GTEx datasets. Our previous studies have validated three downregulated miRNAs (hsa-miR-885-5p, hsa-miR-1909-5p, and hsalet7d-3p) having a diagnostic value in OC patients' sera, with high-throughput techniques. The predicted target genes of these miRNAs were retrieved from the miRDB database (v6.0). Active-subnetwork-oriented Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted by pathfindR using the target genes. Enrichment of KEGG pathways assessed by the analysis of pathfindR indicated that 24 pathways were related to the target genes. Ubiquitin-mediated proteolysis, spliceosome and Notch signaling pathway were the top three pathways with the lowest p-values (p < 0.001). Ninety-three common genes were found to be differentially expressed (p < 0.05) in the datasets. No significant genes were found to be significant in the analysis of overall survival analyses, but 24 genes were found to be significant with pathological stages analysis (p < 0.05). The findings of our study provide in-silico evidence that validated circulating miRNAs' target genes and enriched pathways are related to OC and have potential roles in theranostics applications. Further experimental investigations are required to validate our results which will ultimately provide a new perspective for translational applications in OC management.

A survival prediction model has recently been developed to evaluate the prognosis of resected nonmetastatic pancreatic ductal adenocarcinoma based on a Cox model using two nationwide databases: Surveillance, Epidemiology and End Results (SEER) and Korea Tumor Registry System-Biliary Pancreas (KOTUS-BP). In this study, we applied two machine learning methods-random survival forests (RSF) and support vector machines (SVM)-for survival analysis and compared their prediction performance using the SEER and KOTUS-BP datasets. Three schemes were used for model development and evaluation. First, we utilized data from SEER for model development and used data from KOTUS-BP for external evaluation. Second, these two datasets were swapped by taking data from KOTUS-BP for model development and data from SEER for external evaluation. Finally, we mixed these two datasets half and half and utilized the mixed datasets for model development and validation. We used 9,624 patients from SEER and 3,281 patients from KOTUS-BP to construct a prediction model with seven covariates: age, sex, histologic differentiation, adjuvant treatment, resection margin status, and the American Joint Committee on Cancer 8th edition T-stage and N-stage. Comparing the three schemes, the performance of the Cox model, RSF, and SVM was better when using the mixed datasets than when using the unmixed datasets. When using the mixed datasets, the C-index, 1-year, 2-year, and 3-year time-dependent areas under the curve for the Cox model were 0.644, 0.698, 0.680, and 0.687, respectively. The Cox model performed slightly better than RSF and SVM.

Hypomyelinating leukodystrophy type 2 (HLD2), is an inherited genetic disease of the central nervous system caused by recessive mutations in the gap junction protein gamma 2 (GJC2/GJA12). HLD2 is characterized by nystagmus, developmental delay, motor impairments, ataxia, severe speech problem, and hypomyelination in the brain. The GJC2 sequence encodes connexin 47 protein (Cx47). Connexins are a group of membrane proteins that oligomerize to construct gap junctions protein. In the present study, a novel missense mutation gene c.760G>A (p.Val254Met) was identified in a patient with HLD2 by performing whole exome sequencing. Following the discovery of the new mutation in the proband, we used Sanger sequencing to analyze his affected sibling and parents. Sanger sequencing verified homozygosity of the mutation in the proband and his affected sibling. The autosomal recessive inheritance pattern was confirmed since Sanger sequencing revealed both healthy parents were heterozygous for the mutation. PolyPhen2, SIFT, PROVEAN, and CADD were used to evaluate the function prediction scores of detected mutations. Cx47 is essential for oligodendrocyte function, including adequate myelination and myelin maintenance in humans. Novel mutation p.Val254Met is located in the second extracellular domain of Cx47, both extracellular loops are highly conserved and probably induce intramolecular disulfide interactions. This novel mutation in the Cx47 gene causes oligodendrocyte dysfunction and HLD2 disorder.

Various methodologies for the genetic analysis of longitudinal data have been proposed and applied to data from large-scale genome-wide association studies (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with traits of interest and to detect SNP-time interactions. We recently proposed a grid-based Bayesian mixed model for longitudinal genetic data and showed that our Bayesian method increased the statistical power compared to the corresponding univariate method and well detected SNP-time interactions. In this paper, we further analyze longitudinal obesity-related traits such as body mass index, hip circumference, waist circumference, and waist-hip ratio from Korea Association Resource data to evaluate the proposed Bayesian method. We first conducted GWAS analyses of cross-sectional traits and combined the results of GWAS analyses through a meta-analysis based on a trajectory model and a random-effects model. We then applied our Bayesian method to a subset of SNPs selected by meta-analysis to further discover SNPs associated with traits of interest and SNP-time interactions. The proposed Bayesian method identified several novel SNPs associated with longitudinal obesity-related traits, and almost 25% of the identified SNPs had significant p-values for SNP-time interactions.

According to long-term projections, by 2030, the world's population is predicted to reach 7.5 billion individuals, and there will be roughly 27 million new cancer cases diagnosed. The global burden of breast cancer (BC) is expected to rise. According to the Ministry of Health-Iraqi Cancer Registry, cancer is the second largest cause of death after cardiovascular disease. This study investigated the interleukin-18 (IL18) single-nucleotide polymorphisms (SNPs) -607C/A rs1946518 and -137G/C rs187238 using the sequence-specific amplification-polymerase chain reaction approach. Regarding the position -607C/A, there was a highly significant difference between the observed and expected frequencies in patients and controls (χ 2 = 3.16 and χ 2 = 16.5), respectively. The AA and CA genotypes were associated with significantly increased BC risk (odds ratio [OR], 3.68; p = 0.004 and OR, 2.83; p = 0.04, respectively). Women with the A allele had a 5.03-fold increased susceptibility to BC. The C allele may be a protective allele against BC (OR, 0.19). Although position -137G/C showed no significant differences in the CC genotype distribution (p = 0.18), the frequency of the CC genotype was significantly higher in patients than in controls. In contrast, patients had a significantly higher frequency of GC genotypes than controls (p = 0.04), which was associated with an increased risk of developing BC (OR, 2.63). The G allele frequency was significantly lower in patients than in controls (55.0% vs. 76.2%, respectively). This SNP may be considered a common genotype in the Iraqi population, with the wild-type G allele having a protective function (OR, 0.19) and the mutant C allele having an environmental effect (OR, 2.63).

The influenza A viruses have high mutation rates and cause a serious health problem worldwide. Therefore, this study focused on genome characterization of the viruses isolated from Thai patients based on the next-generation sequencing technology. The nasal swabs were collected from patients with influenza-like illness in Thailand during 2017-2018. Then, the influenza A viruses were detected by reverse transcription-quantitative polymerase chain reaction and isolated by MDCK cells. The viral genomes were amplified and sequenced by Illumina MiSeq platform. Whole genome sequences were used for characterization, phylogenetic construction, mutation analysis and nucleotide diversity of the viruses. The result revealed that 90 samples were positive for the viruses including 44 of A/ H1N1 and 46 of A/H3N2. Among these, 43 samples were successfully isolated and then the viral genomes of 25 samples were completely amplified. Finally, 17 whole genomes of the viruses (A/H1N1, n=12 and A/H3N2, n=5) were successfully sequenced with an average of 232,578 mapped reads and 1,720 genome coverage per sample. Phylogenetic analysis demonstrated that the A/H1N1 viruses were distinguishable from the recommended vaccine strains. However, the A/H3N2 viruses from this study were closely related to the recommended vaccine strains. The nonsynonymous mutations were found in all genes of both viruses, especially in hemagglutinin (HA) and neuraminidase (NA) genes. The nucleotide diversity analysis revealed negative selection in the PB1, PA, HA, and NA genes of the A/H1N1 viruses. High-throughput data in this study allow for genetic characterization of circulating influenza viruses which would be crucial for preparation against pandemic and epidemic outbreaks in the future.

Genetic associations have been quantified using a number of statistical measures. Entropy-based mutual information may be one of the more direct ways of estimating the association, in the sense that it does not depend on the parametrization. For this purpose, both the entropy and conditional entropy of the phenotype distribution should be obtained. Quantitative traits, however, do not usually allow an exact evaluation of entropy. The estimation of entropy needs a probability density function, which can be approximated by kernel density estimation. We have investigated the proper sequence of procedures for combining the kernel density estimation and entropy estimation with a probability density function in order to calculate mutual information. Genotypes and their interactions were constructed to set the conditions for conditional entropy. Extensive simulation data created using three types of generating functions were analyzed using two different kernels as well as two types of multifactor dimensionality reduction and another probability density approximation method called m-spacing. The statistical power in terms of correct detection rates was compared. Using kernels was found to be most useful when the trait distributions were more complex than simple normal or gamma distributions. A full-scale genomic dataset was explored to identify associations using the 2-h oral glucose tolerance test results and γ-glutamyl transpeptidase levels as phenotypes. Clearly distinguishable single-nucleotide polymorphisms (SNPs) and interacting SNP pairs associated with these phenotypes were found and listed with empirical p-values.

The rise of newer coronavirus disease 2019 (COVID-19) variants has brought a challenge to ending the spread of COVID-19. The variants have a different fatality, morbidity, and transmission rates and affect vaccine efficacy differently. Therefore, the impact of each new variant on the spread of COVID-19 is of interest to governments and scientists. Here, we proposed mathematical SEIQRDVP and SEIQRDV3P models to predict the impact of the Omicron variant on the spread of the COVID-19 situation in South Korea. SEIQEDVP considers one vaccine level at a time while SEIQRDV3P considers three vaccination levels (only one dose received, full doses received, and full doses + booster shots received) simultaneously. The omicron variant's effect was contemplated as a weighted sum of the delta and omicron variants' transmission rate and tuned using a hyperparameter k. Our models' performances were compared with common models like SEIR, SEIQR, and SEIQRDVUP using the root mean square error (RMSE). SEIQRDV3P performed better than the SEIQRDVP model. Without consideration of the variant effect, we don't see a rapid rise in COVID-19 cases and high RMSE values. But, with consideration of the omicron variant, we predicted a continuous rapid rise in COVID-19 cases until maybe herd immunity is developed in the population. Also, the RMSE value for the SEIQRDV3P model decreased by 27.4%. Therefore, modeling the impact of any new risen variant is crucial in determining the trajectory of the spread of COVID-19 and determining policies to be implemented.