Evolutionary Bioinformatics最新文献

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.

Sub-Saharan Africa is courting the risk of artemisinin resistance (ARTr) emerging in Plasmodium falciparum malaria parasites. Current molecular surveillance efforts for ARTr have been built on the utility of P. falciparum kelch13 (pfk13) validated molecular markers. However, whether these molecular markers will serve the purpose of early detection of artemisinin-resistant parasites in Ghana is hinged on a pfk13 dependent evolution. Here, we tested the hypothesis that the background pfk13 genome may be present before the pfk13 ARTr-conferring variant(s) is selected and that signatures of balancing selection on these genomic loci may serve as an early warning signal of ARTr. We analyzed 12 198 single nucleotide polymorphisms (SNPs) in Ghanaian clinical isolates in the Pf3K MalariaGEN dataset that passed a stringent filtering regimen. We identified signatures of balancing selection in 2 genes (phosphatidylinositol 4-kinase and chloroquine resistance transporter) previously reported as background loci for ARTr. These genes showed statistically significant and high positive values for Tajima's D, Fu and Li's F, and Fu and Li's D. This indicates that the biodiversity required to establish a pfk13 background genome may have been primed in clinical isolates of P. falciparum from Ghana as of 2010. Despite the absence of ARTr in Ghana to date, our finding supports the current use of pfk13 for molecular surveillance of ARTr in Ghana and highlights the potential utility of monitoring malaria parasite populations for balancing selection in ARTr precursor background genes as early warning molecular signatures for the emergence of ARTr.

Wild-caught animals must cope with drastic lifestyle and dietary changes after being induced to captivity. How the gut microbiome structure of these animals will change in response receives increasing attention. The plateau zokor (Eospalax baileyi), a typic subterranean rodent endemic to the Qinghai-Tibet plateau, spends almost the whole life underground and is well adapted to the environmental pressures of both plateau and underground. However, how the gut microbiome of the plateau zokor will change in response to captivity has not been reported to date. This study compared the microbial community structure and functions of 22 plateau zokors before (the WS group) and after being kept in captivity for 15 days (the LS group, fed on carrots) using the 16S rRNA gene via high-throughput sequencing technology. The results showed that the LS group retained 973 of the 977 operational taxonomic units (OTUs) in the WS group, and no new OTUs were found in the LS group. The dominant bacterial phyla were Bacteroides and Firmicutes in both groups. In alpha diversity analysis, the Shannon, Sobs, and ACE indexes of the LS group were significantly lower than those of the WS group. A remarkable difference (P < 0.01) between groups was also detected in beta diversity analysis. The UPGMA clustering, NMDS, PCoA, and Anosim results all showed that the intergroup difference was significantly greater than the intragroup difference. And compared with the WS group, the intragroup difference of the gut microbiota in the LS group was much larger, which failed to support the assumption that similar diets should drive convergence of gut microbial communities. PICRUSt revealed that although some functional categories displayed significant differences between groups, the relative abundances of these categories were very close in both groups. Based on all the results, we conclude that as plateau zokors enter captivity for a short time, although the relative abundances of different gut microbiota categories shifted significantly, they can maintain almost all the OTUs and the functions of the gut microbiota in the wild. So, the use of wild-caught plateau zokors in gut microbial studies is acceptable if the time in captivity is short.

Hepatocellular carcinoma (HCC) is one of the common cancers with a high incidence and mortality. The human replication factor C (RFC) family contains 5 subunits that play an important role in DNA replication and DNA damage repair. RFCs are abnormally expressed in a variety of cancers; some of them are differentially expressed in HCC tissues and related to tumor growth. However, the expression, prognostic value, and effect targets of the whole RFC family in HCC are still unclear. To address these issues, we performed a multidimensional analysis of RFCs in HCC patients by Oncomine, UALCAN, GEPIA, Human protein atlas, Kaplan-Meier plotter, cBioPortal, GeneMANIA, String, and LinkedOmics. mRNA expression of RFCs was significantly increased in HCC tissues. There was a significant correlation between the expression of RFC2/3/4/5 and tumor stage of HCC patients. Besides, high mRNA expression of RFC2/4 was associated with worse overall survival (OS). Moreover, genetic alterations of RFCs were associated with worse OS in HCC patients. We found that genes co-expressed with RFC2/4 were mainly involved in biological processes, such as chromosome segregation, mitotic cell cycle phase transition, and telomere organization and they activated the cell cycle and spliceosome pathways. The gene set is mainly enriched in cancer-related kinases AURKA, ATR, CDK1, PLK1, and CHEK1. E2F family members were the key transcription factors for RFCs. Our results suggest that differentially expressed RFC2 and RFC4 are potential prognostic biomarkers in HCC and may act on E2F transcription factors and some kinase targets to dysregulate the cell cycle pathway. These efforts may provide new research directions for prognostic biomarkers and therapeutic targets in HCC.

Cycads have developed a complex root system categorized either as normal or coralloid roots. Past literatures revealed that a great diversity of key microbes is associated with these roots. This recent study aims to comprehensively determine the diversity and community structure of bacteria and fungi associated with the roots of two Cycas spp. endemic to China, Cycas debaoensis Zhong & Chen and Cycas fairylakea D.Y. Wang using high-throughput amplicon sequencing of the full-length 16S rRNA (V1-V9 hypervariable) and short fragment ITS region. The total DNA from 12 root samples were extracted, amplified, sequenced, and analyzed. Resulting sequences were clustered into 61 bacteria and 2128 fungal OTUs. Analysis of community structure revealed that the coralloid roots were dominated mostly by the nitrogen-fixer Nostocaceae but also contain other non-diazotrophic bacteria. The sequencing of entire 16S rRNA gene identified four different strains of cyanobacteria under the heterocystous genera Nostoc and Desmonostoc. Meanwhile, the top bacterial families in normal roots were Xanthobacteraceae, Burkholderiaceae, and Bacillaceae. Moreover, a diverse fungal community was also found in the roots of cycads and the predominating families were Ophiocordycipitaceae, Nectriaceae, Bionectriaceae, and Trichocomaceae. Our results demonstrated that bacterial diversity in normal roots of C. fairylakea is higher in richness and abundance than C. debaoensis. On the other hand, a slight difference, albeit insignificant, was noted for the diversity of fungi among root types and host species as the number of shared taxa is relatively high (67%). Our results suggested that diverse microbes are present in roots of cycads which potentially interact together to support cycads survival. Our study provided additional knowledge on the microbial diversity and composition in cycads and thus expanding our current knowledge on cycad-microbe association. Our study also considered the possible impact of ex situ conservation on cyanobiont community of cycads.