Evolutionary Bioinformatics最新文献

The RNA G-quadruplex (rG4) is a kind of non-canonical high-order secondary structure with important biological functions and is enriched in untranslated regions (UTRs) of protein-coding genes. However, how rG4 structures evolve is largely unknown. Here, we systematically investigated the evolution of RNA sequences around UTR rG4 structures in 5 eukaryotic organisms. We found universal selection on UTR sequences, which facilitated rG4 formation in all the organisms that we analyzed. While G-rich sequences were preferred in the rG4 structural region, C-rich sequences were selectively not preferred. The selective pressure acting on rG4 structures in the UTRs of genes with higher G content was significantly smaller. Furthermore, we found that rG4 structures experienced smaller evolutionary selection near the translation initiation region in the 5' UTR, near the polyadenylation signals in the 3' UTR, and in regions flanking the miRNA targets in the 3' UTR. These results suggest universal selection for rG4 formation in the UTRs of eukaryotic genomes and the selection may be related to the biological functions of rG4s.

Gefitinib resistance is a serious threat in the treatment of patients with non-small cell lung cancer (NSCLC). Elucidating the underlying mechanisms and developing effective therapies to overcome gefitinib resistance is urgently needed. The differentially expressed genes (DEGs) were screened from the gene expression profile GSE122005 between gefitinib-sensitive and resistant samples. GO and KEGG analyses were performed with DAVID. The protein-protein interaction (PPI) network was established to visualize DEGs and screen hub genes. The functional roles of CCL20 in lung adenocarcinoma (LUAD) were examined using gene set enrichment analysis (GSEA). Functional analysis revealed that the DEGs were mainly concentrated in inflammatory, cell chemotaxis, and PI3K signal regulation. Ten hub genes were identified based on the PPI network. The survival analysis of the hub genes showed that CCL20 had a significant effect on the prognosis of LUAD patients. GSEA analysis showed that CCL20 high expression group was mainly enriched in cytokine-related signaling pathways. In conclusion, our analysis suggests that changes in inflammation and cytokine-related signaling pathways are closely related to gefitinib resistance in patients with lung cancer. The CCL20 gene may promote the formation of gefitinib resistance, which may serve as a new biomarker for predicting gefitinib resistance in patients with lung cancer.

Soil contamination by hydrocarbons due to oil spills has become a global concern and it has more implications in oil producing regions. Biostimulation is considered as one of the promising remediation techniques that can be adopted to enhance the rate of degradation of crude oil. The soil microbial consortia play a critical role in governing the biodegradation of total petroleum hydrocarbons (TPHs), in particular polycyclic aromatic hydrocarbons (PAHs). In this study, the degradation pattern of TPHs and PAHs of Kuwait soil biopiles was measured at three-month intervals. Then, the microbial consortium associated with oil degradation at each interval was revealed through 16S rRNA based next generation sequencing. Rapid degradation of TPHs and most of the PAHs was noticed at the first 3 months of biostimulation with a degradation rate of pyrene significantly higher compared to other PAHs counterparts. The taxonomic profiling of individual stages of remediation revealed that, biostimulation of the investigated soil favored the growth of Proteobacteria, Alphaprotobacteria, Chloroflexi, Chlorobi, and Acidobacteria groups. These findings provide a key step towards the restoration of oil-contaminated lands in the arid environment.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two sets of SARS-CoV-2 genomic sequences to identify emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, 20 variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5' untranslated region (UTR), including a group of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-ss molecular switch in the 3'UTR. Finally, 5 variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M (stem loop 2 m) selfish genetic element, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. Our analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences.

Antibiotic resistance is a major global health issue that has seen alarming rates of increase in all parts of the world over the past two decades. The surge in antibiotic resistance has resulted in longer hospital stays, higher medical costs, and elevated mortality rates. Constant attempts have been made to discover newer and more effective antimicrobials to reduce the severity of antibiotic resistance. Plant secondary metabolites, such as essential oils, have been the major focus due to their complexity and bioactive nature. However, the underlying mechanism of their antimicrobial effect remains largely unknown. Understanding the antimicrobial mode of action of essential oils is crucial in developing potential strategies for the use of essential oils in a clinical setting. Recent advances in genomics and proteomics have enhanced our understanding of the antimicrobial mode of action of essential oils. We might well be at the dawn of completing a mystery on how essential oils carry out their antimicrobial activities. Therefore, an overview of essential oils with regard to their antimicrobial activities and mode of action is discussed in this review. Recent approaches used in identifying the antimicrobial mode of action of essential oils, specifically from the perspective of genomics and proteomics, are also synthesized. Based on the information gathered from this review, we offer recommendations for future strategies and prospects for the study of essential oils and their function as antimicrobials.

The phytohormone auxin are important in all aspects of plant growth and development. The Auxin/Indole-3-Acetic Acid (Aux/IAA) gene responds to auxin induction as auxin early response gene family. Despite the physiological importance of the Aux/IAA gene, a systematic analysis of the Aux/IAA gene in Acer rubrum has not been reported. This paper describes the characterization of Acer rubrum Aux/IAA genes at the transcriptomic level and Acer yangbiense Aux/IAA genes at the genomic level, with 17 Acer rubrum AUX/IAA genes (ArAux/IAA) and 23 Acer yangbiense Aux/IAA (AyAux/IAA) genes identified. Phylogenetic analysis shows that AyAux/IAA and ArAux/IAA family genes can be subdivided into 4 groups and show strong evolutionary conservatism. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to test the expression profile of ArAux/IAA genes in different tissues under indole-3-acetic acid (IAA) treatment. Most ArAux/IAA genes are responsive to exogenous auxin and have tissue-specific expression. Overall, these results will provide molecular-level insights into auxin metabolism, transport, and signaling in Acer species.

Lung adenocarcinoma (LUAD) is a tumor with high incidence. This study aimed to identify the central genes of LUAD. LUAD were analyzed by weighted gene co-expression network (WGCNA), and differentially expressed genes (DEGs) were identified. Samples were obtained from The Cancer Genome Atlas (TCGA) and Genotype Tissue Expression (GTEx) databases and included 515 LUAD samples and 347 normal samples. The WGCNA algorithm generated a total of 10 modules. The top 2 modules (MEturquoise and MEblue) with the highest correlation to LUAD were selected. Ten Hub genes (IL6, CDH1, PECAM1, SPP1, THBS1, HGF, SNCA, CDH5, CAV1, and DLC1) were screened in the intersecting genes of DEGs and WGCNA (MEturquoise and MEblue). Only SPP1 was correlated with LUAD poor survival, indicating that SPP1 may be a key Hub gene for LUAD. The competing endogenous RNA (ceRNA) network was constructed to analyze the regulatory relationship of Hub genes, and SPP1 may be directly regulated by 4 microRNAs (miRNAs) and indirectly regulated by 49 long noncoding RNAs (lncRNAs).

A high level of mutation enables the influenza A virus to resist antibiotics previously effective against the influenza A virus. A portion of the structure of hemagglutinin HA is assumed to be well-conserved to maintain its role in cellular fusion, and the structure tends to be more conserved than sequence. We designed peptide inhibitors to target the conserved residues on the HA surface, which were identified based on structural alignment. Most of the conserved and strongly similar residues are located in the receptor-binding and esterase regions on the HA1 domain In a later step, fragments of anti-HA antibodies were gathered and screened for the binding ability to the found conserved residues. As a result, Methionine amino acid got the best docking score within the -2.8 Å radius of Van der Waals when it is interacting with Tyrosine, Arginine, and Glutamic acid. Then, the binding affinity and spectrum of the fragments were enhanced by grafting hotspot amino acid into the fragments to form peptide inhibitors. Our peptide inhibitor was able to form in silico contact with a structurally conserved region across H1, H2, and H3 HA, with the binding site at the boundary between HA1 and HA2 domains, spreading across different monomers, suggesting a new target for designing broad-spectrum antibody and vaccine. This research presents an affordable method to design broad-spectrum peptide inhibitors using fragments of an antibody as a scaffold.

ORF8 is a highly variable genomic region of SARS-CoV-2. Although non-essential and the precise functions are unknown, it has been suggested that this protein assists in SARS-CoV-2 replication in the early secretory pathway and in immune evasion. We utilized the binding partners of SARS-CoV-2 proteins in human HEK293T cells and performed genome-wide phylogenetic profiling and clustering analyses in 446 eukaryotic species to predict and discover ORF8 binding partners that share associated functional mechanisms based on co-evolution. Results classified 47 ORF8 binding partner proteins into 3 clusters (groups 1-3), which were conserved in vertebrates (group 1), metazoan (group 2), and eukaryotes (group 3). Gene ontology analysis indicated that group 1 had no significant associated biological processes, while groups 2 and 3 were associated with glycoprotein biosynthesis process and ubiquitin-dependent endoplasmic reticulum-associated degradation pathways, respectively. Collectively, our results classified potential genes that might be associated with SARS-CoV-2 viral pathogenesis, specifically related to acute respiratory distress syndrome, and the secretory pathway. Here, we discuss the possible role of ORF8 in viral pathogenesis and in assisting viral replication and immune evasion via secretory pathway, as well as the possible factors associated with the rapid evolution of ORF8.