Evolutionary Bioinformatics最新文献

A common task in bioinformatics is to compare DNA sequences to identify similarities between organisms at the sequence level. An approach to such comparison is the dot-plots, a 2-dimensional graphical representation to analyze DNA or protein alignments. Dot-plots alignment software existed before the sequencing revolution, and now there is an ongoing limitation when dealing with large-size sequences, resulting in very long execution times. High-Performance Computing (HPC) techniques have been successfully used in many applications to reduce computing times, but so far, very few applications for graphical sequence alignment using HPC have been reported. Here, we present G-SAIP (Graphical Sequence Alignment in Parallel), a software capable of spawning multiple distributed processes on CPUs, over a supercomputing infrastructure to speed up the execution time for dot-plot generation up to 1.68× compared with other current fastest tools, improve the efficiency for comparative structural genomic analysis, phylogenetics because the benefits of pairwise alignments for comparison between genomes, repetitive structure identification, and assembly quality checking.

Autosomal dominant hyper-IgE syndrome (AD-HIES) is linked to dominant negative mutations of the STAT3 protein whose molecular basis for dysfunction is unclear and presenting with a variety of clinical manifestations with only supportive treatment. To establish the relationship between the impact of STAT3 mutations in different domains and the severity of the clinical manifestations, 105 STAT3 mutations were analyzed for their impact on protein stability, flexibility, function, and binding affinity using in Silico approaches. Our results showed that 73% of the studied mutations have an impact on the physicochemical properties of the protein, altering the stability, flexibility and function to varying degrees. In particular, mutations affecting the DNA binding domain (DBD) and the Src Homology 2 (SH2) have a significant impact on the protein structure and disrupt its interaction either with DNA or other STAT3 to form a heterodomain complex, leading to severe clinical phenotypes. Collectively, this study suggests that there is a close relationship between the domain involving the mutation, the degree of variation in the properties of the protein and the degree of loss of function ranging from partial loss to complete loss, explaining the variability of clinical manifestations between mild and severe.

MicroRNAs (miRNAs) are a class of non-coding RNAs that act as regulators of disease. An evolutionary approach to the disease could reveal factors such as diagnosis, treatment, and prognosis prediction. The expression patterns of transposable element (TEs)-derived miRNAs could help elucidate diseases, and their evolutionary patterns are also valuable. The 34 miRNAs were compared in terms of stage survival and tumor status in 33 carcinomas from TCGA. Expression levels were compared using a t-test and presented as differentially expressed miRNAs (DEMs). For DEMs showing statistically specific expression patterns for 3 conditions (normal and cancer, early and advanced stage, and survival), interactions with related genes in 10 species, including humans, were compared. The enrichment term was discovered for the gene-miRNA interactions. In 18 out of the 33 carcinomas, at least one miRNA was retrieved with P < .05 and |fold change| >.05. A total of 128 DEMs for the 9 miRNAs were identified. Based on the TargetScan database, interactions between miRNAs and genes in 10 species, including humans, were confirmed. The evolutionarily conserved miR-130a was observed in all 10 species, whereas miR-151a was only observed in humans. GO terms of related genes were selected for the miRNAs commonly found in each species. Evolutionary analysis of TE-derived disease-associated miRNAs was performed, and the evolutionarily conserved miR-130a-related carcinomas included renal and thyroid cancers. Human and rhesus monkey-specific miR-625 is associated with various carcinomas.

Based on the important role of antibiotic treatment in the research of the interaction between Wolbachia and insect hosts, this study aimed to identify the most suitable antibiotic and concentration for Wolbachia elimination in the P. xylostella, and to investigate the effect of Wolbachia and antibiotic treatment on the bacterial community of P. xylostella. Our results showed that the Wolbachia-infected strain was plutWB1 of supergroup B in the P. xylostella population collected in Nepal in this study; 1 mg/mL rifampicin could remove Wolbachia infection in P. xylostella after 1 generation of feeding treatment and the toxic effect was relatively low; among the 29 samples of adult P. xylostella in our study (10 WU samples, 10 WA samples, and 9 WI samples), 52.5% of the sequences were of Firmicutes and 47.5% were of Proteobacteria, with the dominant genera being mainly Carnobacterium (46.2%), Enterobacter (10.1%), and Enterococcus (6.2%); Moreover, antibiotic removal of Wolbachia infection in P. xylostella and transfer to normal conditions for 10 generations no longer significantly affected the bacterial community of P. xylostella. This study provides a theoretical basis for the elimination method of Wolbachia in the P. xylostella, as well as a reference for the elimination method of Wolbachia in other Wolbachia-infected insect species, and a basis for the study of the extent and duration of the effect of antibiotic treatment on the bacterial community of the P. xylostella.

Hepatocellular carcinoma (HCC) is the most common primary malignancy of the liver. Although the RNA modification N6-methyladenine (m6A) has been reported to be involved in HCC carcinogenesis, early diagnostic markers and promising personalized therapeutic targets are still lacking. In this study, we identified that 19 m6A regulators and 34 co-expressed lncRNAs were significantly upregulated in HCC samples; based on these factors, we established a prognostic signal of HCC associated with 9 lncRNAs and 19 m6A regulators using LASSO Cox regression analysis. Kaplan-Meier survival estimate revealed correlations between the risk scores and patients' OS in the training and validation dataset. The ROC curve demonstrated that the risk score-based curve has satisfactory prediction efficiency for both training and validation datasets. Multivariate Cox's proportional hazard regression analysis indicated that the risk score was an independent risk factor within the training and validation dataset. In addition, the risk score could distinguish HCC patients from normal non-cancerous samples and HCC samples of different pathological grades. Eventually, 232 mRNAs were co-expressed with these 9 lncRNAs according to GSE101685 and GSE112790; these mRNAs were enriched in cell cycle and cell metabolic activities, drug metabolism, liver disease-related pathways, and some important cancer related pathways such as p53, MAPK, Wnt, RAS and so forth. The expression of the 9 lncRNAs was significantly higher in HCC samples than that in the neighboring non-cancerous samples. Altogether, by using the Consensus Clustering, PCA, ESTIMATE algorithm, LASSO regression model, Kaplan-Meier survival assessment, ROC curve analysis, and multivariate Cox's proportional hazard regression model analysis, we established a prognostic marker consisting of 9 m6A regulator-related lncRNAs that markers may have prognostic and diagnostic potential for HCC.

B12D-Like is a member of the B12D domain-containing protein family, which includes several transmembrane proteins in plants. In this study, the cDNA of PgB12D-Like from Pennisetum glaucum subsp. monodii (Maire) Brunken was sequenced and characterized. The 446-bp cDNA for PgB12D-Like encodes for a deduced protein of 95 amino acids. The PgB12D-Like protein contains a B12D domain and a transmembrane helix embedded in the mitochondrial membrane. Cis-regulatory elements analysis reveals binding sites for various transcription factors involved in responses to stress, light, and plant hormones in the putative promoter sequence for PgB12D-Like. Several proteins involved in floral organ development were also found to have binding sites in the PgB12D-Like promoter, such as agamous-like proteins and squamosa promoter binding proteins. Real-time PCR reveals high expression of PgB12D-Like in flowers during heading, whereas its expression in a 4-day-old seedling shoot was the lowest. Moreover, cold, drought, and heat stress were found to upregulate PgB12D-Like, whereas gibberellic acid downregulated its expression in seedlings. The present study helps to uncover the function of the B12D-Like in response to plant hormones and abiotic stress during P. glaucum development.