Evolutionary Bioinformatics最新文献

Microbial community is ubiquitous in nature, which has a great impact on the living environment and human health. All these effects of microbial communities on the environment and their hosts are often referred to as the functions of these communities, which depend largely on the composition of the communities. The study of microbial higher-order module can help us understand the dynamic development and evolution process of microbial community and explore community function. Considering that traditional clustering methods depend on the number of clusters or the influence of data that does not belong to any cluster, this paper proposes a hypergraph clustering algorithm based on game theory to mine the microbial high-order interaction module (HCGI), and the hypergraph clustering problem naturally turns into a clustering game problem, the partition of network modules is transformed into finding the critical point of evolutionary stability strategy (ESS). The experimental results show HCGI does not depend on the number of classes, and can get more conservative and better quality microbial clustering module, which provides reference for researchers and saves time and cost. The source code of HCGI in this paper can be downloaded from https://github.com/ylm0505/HCGI.

Childhood obesity is a serious public health problem worldwide. Perturbations in the gut microbiota composition have been associated with the development of obesity in both children and adults. Probiotics, on the other hand, are proven to restore the composition of the gut microbiome which helps reduce the development of obesity. However, data on the effect of probiotics on gut microbiota and its association with childhood obesity is limited. This study aims to determine the effect of probiotics supplement intervention on gut microbiota profiles in obese and normal-weight children. A total of 37 children, 17 normal weight, and 20 overweight school children from a government school in Selangor were selected to participate in this study. Participants were further divided into intervention and control groups. The intervention groups received daily probiotic drinks while the control groups continued eating their typical diet. Fecal samples were collected from the participants for DNA extraction. The hypervariable V3 and V4 regions of 16S rRNA gene were amplified and sequenced using the Illumina MiSeq platform. No significant differences in alpha diversity were observed between normal weight and obese children in terms of the Shannon Index for evenness or species richness. However, a higher intervention effect on alpha diversity was observed among normal-weight participants compared to obese. The participants' microbiome was found to fluctuate throughout the study. Analysis of the taxa at species level showed an increase in Bacteroides ovatus among the normal weight cohort. Genus-level comparison revealed a rise in genus Lachnospira and Ruminococcus in the overweight participants after intervention, compared to the normal-weight participants. The probiotics intervention causes an alteration in gut microbiota composition in both normal and overweight children. Though the association could not be defined statistically, this study has provided an improved understanding of the intervention effect of probiotics on gut microbiome dysbiosis in an underrepresented population.

The massive worldwide spread of the SARS-CoV-2 virus is fueling the COVID-19 pandemic. Since the first whole-genome sequence was published in January 2020, a growing database of tens of thousands of viral genomes has been constructed. This offers opportunities to study pathways of molecular change in the expanding viral population that can help identify molecular culprits of virulence and virus spread. Here we investigate the genomic accumulation of mutations at various time points of the early pandemic to identify changes in mutationally highly active genomic regions that are occurring worldwide. We used the Wuhan NC_045512.2 sequence as a reference and sampled 15 342 indexed sequences from GISAID, translating them into proteins and grouping them by month of deposition. The per-position amino acid frequencies and Shannon entropies of the coding sequences were calculated for each month, and a map of intrinsic disorder regions and binding sites was generated. The analysis revealed dominant variants, most of which were located in loop regions and on the surface of the proteins. Mutation entropy decreased between March and April of 2020 after steady increases at several sites, including the D614G mutation site of the spike (S) protein that was previously found associated with higher case fatality rates and at sites of the NSP12 polymerase and the NSP13 helicase proteins. Notable expanding mutations include R203K and G204R of the nucleocapsid (N) protein inter-domain linker region and G251V of the viroporin encoded by ORF3a between March and April. The regions spanning these mutations exhibited significant intrinsic disorder, which was enhanced and decreased by the N-protein and viroporin 3a protein mutations, respectively. These results predict an ongoing mutational shift from the spike and replication complex to other regions, especially to encoded molecules known to represent major β-interferon antagonists. The study provides valuable information for therapeutics and vaccine design, as well as insight into mutation tendencies that could facilitate preventive control.

Communication systems within and between plant cells involve the transfer of ions and molecules between compartments, and are essential for development and responses to biotic and abiotic stresses. This in turn requires the regulated movement and fusion of membrane systems with their associated cargo. Recent advances in genomics has provided new resources with which to investigate the evolutionary relationships between membrane proteins across plant species. Members of the soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) are known to play important roles in vesicle trafficking across plant, animal and microbial species. Using recent public expression and transcriptomic data from 9 representative green plants, we investigated the evolution of the SNARE classes and linked protein changes to functional specialization (expression patterns). We identified an additional 3 putative SNARE genes in the model plant Arabidopsis. We found that all SNARE classes have expanded in number to a greater or lesser degree alongside the evolution of multicellularity, and that within-species expansions are also common. These gene expansions appear to be associated with the accumulation of amino acid changes and with sub-functionalization of SNARE family members to different tissues. These results provide an insight into SNARE protein evolution and functional specialization. The work provides a platform for hypothesis-building and future research into the precise functions of these proteins in plant development and responses to the environment.

Parasitic cysteine proteases are involved in parasite stage transition, invasion of host tissues, nutrient uptake, and immune evasion. The cysteine protease cathepsin F is the most abundant protein produced by fourth-stage larvae (L4) of the nematode Teladorsagia circumcincta, while its transcript is only detectable in L4 and adults. T. circumcincta cathepsin F is a recently evolved cysteine protease that does not fall clearly into either of the cathepsin L or F subfamilies. This protein exhibits characteristics of both cathepsins F and L, and its phylogenetic relationship to its closest homologs is distant, including proteins of closely related nematodes of the same subfamily.

Hepatitis B virus (HBV) infection is a major cause of acute liver failure (ALF) in China, and mortality rates are high among patients who do not receive a matched liver transplant. This study aimed to determine potential mechanisms involved in HBV-ALF pathogenesis. Gene expression profiles under access numbers GSE38941 and GSE14668 were downloaded from the Gene Expression Omnibus database, including cohorts of HBV-ALF liver tissue and normal samples. Differentially expressed genes (DEGs) with false discovery rates (FDR) <0.05 and |log₂(fold change)| >1 as thresholds were screened using the Limma package. Gene modules associated with stable disease were mined using weighed gene co-expression network analysis (WGCNA). A co-expression network was constructed and DEGs were analyzed using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. A gene-based network was constructed to explore major factors associated with disease progression. We identified 2238 overlapping DEGs as crucial gene cohorts in ALF development. Based on a WGCNA algorithm, 10 modules (modules 1-10) were obtained that ranged from 75 to 1078 genes per module. Cyclin-dependent kinase 1 (CDK1), cyclin B1 (CCNB1), and cell-division cycle protein 20 (CDC20) hub genes were screened using the co-expression network. Furthermore, 17 GO terms and 6 KEGG pathways were identified, such as cell division, immune response process, and antigen processing and presentation. Two overlapping signaling pathways that are crucial factors in HBV-ALF were screened using the Comprehensive Toxicogenomics Database (CTD). Several candidate genes including HLA-E, B2M, HLA-DPA1, and SYK were associated with HBV-ALF progression. Natural killer cell-mediated cytotoxicity and antigen presentation contributed to the progression of HBV-ALF. The HLA-E, B2M, HLA-DPA1, and SYK genes play critical roles in the pathogenesis and development of HBV-ALF.

The dysbiosis of the gut microbiome associated with ulcerative colitis (UC) has been extensively studied in recent years. However, the question of whether UC influences the spatial heterogeneity of the human gut mucosal microbiome has not been addressed. Spatial heterogeneity (specifically, the inter-individual heterogeneity in microbial species abundances) is one of the most important characterizations at both population and community scales, and can be assessed and interpreted by Taylor's power law (TPL) and its community-scale extensions (TPLEs). Due to the high mobility of microbes, it is difficult to investigate their spatial heterogeneity explicitly; however, TPLE offers an effective approach to implicitly analyze the microbial communities. Here, we investigated the influence of UC on the spatial heterogeneity of the gut microbiome with intestinal mucosal microbiome samples collected from 28 UC patients and healthy controls. Specifically, we applied Type-I TPLE for measuring community spatial heterogeneity and Type-III TPLE for measuring mixed-species population heterogeneity to evaluate the heterogeneity changes of the mucosal microbiome induced by UC at both the community and species scales. We further used permutation test to determine the possible differences between UC patients and healthy controls in heterogeneity scaling parameters. Results showed that UC did not significantly influence gut mucosal microbiome heterogeneity at either the community or mixed-species levels. These findings demonstrated significant resilience of the human gut microbiome and confirmed a prediction of TPLE: that the inter-subject heterogeneity scaling parameter of the gut microbiome is an intrinsic property to humans, invariant with UC disease.

Monitoring the mutation and evolution of the virus is important for tracing its ongoing transmission and facilitating effective vaccine development. A total of 342 complete genomic sequences of SARS-CoV-2 were analyzed in this study. Compared to the reference genome reported in December 2019, 465 mutations were found, among which, 347 occurred in only 1 sequence, while 26 occurred in more than 5 sequences. For these 26 further identified as SNPs, 14 were closely linked and were grouped into 5 profiles. Phylogenetic analysis revealed the sequences formed 2 major groups. Most of the sequences in late period (March and April) constituted the Cluster II, while the sequences before March in this study and the reported S/L and A/B/C types in previous studies were all in Cluster I. The distributions of some mutations were specific geographically or temporally, the potential effect of which on the transmission and pathogenicity of SARS-CoV-2 deserves further evaluation and monitoring. Two mutations were found in the receptor-binding domain (RBD) but outside the receptor-binding motif (RBM), indicating that mutations may only have marginal biological effects but merit further attention. The observed novel sequence divergence is of great significance to the study of the transmission, pathogenicity, and development of an effective vaccine for SARS-CoV-2.

The detection of copy number variations (CNVs) on whole-exome sequencing (WES) represents a cost-effective technique for the study of genetic variants. This approach, however, has encountered an obstacle with high false-positive rates due to biases from exome sequencing capture kits and GC contents. Although plenty of CNV detection tools have been developed, they do not perform well with all types of CNVs. In addition, most tools lack features of genetic annotation, CNV visualization, and flexible installation, requiring users to put much effort into CNV interpretation. Here, we present "inCNV," a web-based application that can accept multiple CNV-tool results, then integrate and prioritize them with user-friendly interfaces. This application helps users analyze the importance of called CNVs by generating CNV annotations from Ensembl, Database of Genomic Variants (DGV), ClinVar, and Online Mendelian Inheritance in Man (OMIM). Moreover, users can select and export CNVs of interest including their flanking sequences for primer design and experimental verification. We demonstrated how inCNV could help users filter and narrow down the called CNVs to a potentially novel CNV, a common CNV within a group of samples of the same disease, or a de novo CNV of a sample within the same family. Besides, we have provided in CNV as a docker image for ease of installation (https://github.com/saowwapark/inCNV).

Hepatocellular carcinoma (HCC) is a common malignant tumor representing more than 90% of primary liver cancer. This study aimed to identify metabolism-related biomarkers with prognostic value by developing the novel prognostic score (PS) model. Transcriptomic profiles derived from TCGA and EBIArray databases were analyzed to identify differentially expressed genes (DEGs) in HCC tumor samples compared with normal samples. The overlapped genes between DEGs and metabolism-related genes (crucial genes) were screened and functionally analyzed. A novel PS model was constructed to identify optimal signature genes. Cox regression analysis was performed to identify independent clinical factors related to prognosis. Nomogram model was constructed to estimate the predictability of clinical factors. Finally, protein expression of crucial genes was explored in different cancer tissues and cell types from the Human Protein Atlas (HPA). We screened a total of 305 overlapped genes (differentially expressed metabolism-related genes). These genes were mainly involved in "oxidation reduction," "steroid hormone biosynthesis," "fatty acid metabolic process," and "linoleic acid metabolism." Furthermore, we screened ten optimal DEGs (CYP2C9, CYP3A4, and TKT, among others) by using the PS model. Two clinical factors of pathologic stage (P < .001, HR: 1.512 [1.219-1.875]) and PS status (P <.001, HR: 2.259 [1.522-3.354]) were independent prognostic predictors by cox regression analysis. Nomogram model showed a high predicted probability of overall survival time, and the AUC value was 0.837. The expression status of 7 proteins was frequently altered in normal or differential tumor tissues, such as liver cancer and stomach cancer samples.We have identified several metabolism-related biomarkers for prognosis prediction of HCC based on the PS model. Two clinical factors were independent prognostic predictors of pathologic stage and PS status (high/low risk). The prognosis prediction model described in this study is a useful and stable method for novel biomarker identification.