Search

Background Lysine succinylation is a naturally occurring post-translational modification (PTM) that is ubiquitous in organisms. Lysine succinylation plays important roles in regulating protein structure and function as well as cellular metabolism. Global lysine succinylation at the proteomic level has been identified in a variety of species; however, limited information on lysine succinylation in plant species, especially paper mulberry, is available. Paper mulberry is not only an important plant in traditional Chinese medicine, but it is also a tree species with significant economic value. Paper mulberry is found in the temperate and tropical zones of China. The present study analyzed the effects of lysine succinylation on the growth, development, and physiology of paper mulberry. Results A total of 2097 lysine succinylation sites were identified in 935 proteins associated with the citric acid cycle (TCA cycle), glyoxylic acid and dicarboxylic acid metabolism, ribosomes and oxidative phosphorylation; these pathways play a role in carbon fixation in photosynthetic organisms and may be regulated by lysine succinylation. The modified proteins were distributed in multiple subcellular compartments and were involved in a wide variety of biological processes, such as photosynthesis and the Calvin-Benson cycle. Conclusion Lysine-succinylated proteins may play key regulatory roles in metabolism, primarily in photosynthesis and oxidative phosphorylation, as well as in many other cellular processes. In addition to the large number of succinylated proteins associated with photosynthesis and oxidative phosphorylation, some proteins associated with the TCA cycle are succinylated. Our study can serve as a reference for further proteomics studies of the downstream effects of succinylation on the physiology and biochemistry of paper mulberry.

The second ISCB-Asia conference of the International Society for Computational Biology took place December 17-19, 2012, in Shenzhen, China. The conference was co-hosted by BGI as the first Shenzhen Conference on Computational Genomics (SCCG). 45 talks were presented at ISCB-Asia/SCCG 2012. The topics covered included software tools, reproducible computing, next-generation sequencing data analysis, transcription and mRNA regulation, protein structure and function, cancer genomics and personalized medicine. Nine of the proceedings track talks are included as full papers in this supplement. In this report we first give a short overview of the conference by listing some statistics and visualizing the talk abstracts as word clouds. Then we group the talks by topic and briefly summarize each one, providing references to related publications whenever possible. Finally, we close with a few comments on the success of this conference.

Background In social groups, dominant individuals may socially inhibit reproduction of subordinates using aggressive interactions or, in the case of highly eusocial insects, pheromonal communication. It has been hypothesized these two modes of reproductive inhibition utilize conserved pathways. Here, we use a comparative framework to investigate the chemical and genomic underpinnings of reproductive dominance in the primitively eusocial wasp Polistes metricus. Our goals were to first characterize transcriptomic and chemical correlates of reproductive dominance and second, to test whether dominance-associated mechanisms in paper wasps overlapped with aggression or pheromone-related gene expression patterns in other species. To explore whether conserved molecular pathways relate to dominance, we compared wasp transcriptomic data to previous studies of gene expression associated with pheromonal communication and queen-worker differences in honey bees, and aggressive behavior in bees, Drosophila, and mice. Results By examining dominant and subordinate females from queen and worker castes in early and late season colonies, we found that cuticular hydrocarbon profiles and genome-wide patterns of brain gene expression were primarily associated with season/social environment rather than dominance status. In contrast, gene expression patterns in the ovaries were associated primarily with caste and ovary activation. Comparative analyses suggest genes identified as differentially expressed in wasp brains are not related to queen pheromonal communication or caste in bees, but were significantly more likely to be associated with aggression in other insects (bees, flies), and even a mammal (mice). Conclusions This study provides the first comprehensive chemical and molecular analysis of reproductive dominance in paper wasps. We found little evidence for a chemical basis for reproductive dominance in P. metricus, and our transcriptomic analyses suggest that different pathways regulate dominance in paper wasps and pheromone response in bees. Furthermore, there was a substantial impact of season/social environment on gene expression patterns, indicating the important role of external cues in shaping the molecular processes regulating behavior. Interestingly, genes associated with dominance in wasps were also associated with aggressive behavior in bees, solitary insects and mammals. Thus, genes involved in social regulation of reproduction in Polistes may have conserved functions associated with aggression in insects and other taxa.

Background Efficient use of feed resources for farm animals is a critical concern in animal husbandry. Numerous genetic and nutritional studies have been conducted to investigate feed efficiency during the regular laying cycle of chickens. However, by prolonging the laying period of layers, the performance of feed utilization in the late-laying period becomes increasingly important. In the present study, we measured daily feed intake (FI), residual feed intake (RFI) and feed conversion ratio (FCR) of 808 hens during 81–82 weeks of age to evaluate genetic properties and then used a genome-wide association study (GWAS) to reveal the genetic determinants. Results The heritability estimates for the investigated traits were medium and between 0.15 and 0.28 in both pedigree- and genomic-based estimates, whereas the genetic correlations among these traits were high and ranged from 0.49 to 0.90. Three genome-wide significant SNPs located on chromosome 1 (GGA1) were detected for FCR. Linkage disequilibrium (LD) and conditional GWA analysis indicated that these 3 SNPs were highly correlated with one another, located at 13.55–45.16 Kb upstream of gga-miR-15a. Results of quantitative real-time polymerase chain reaction (qRT-PCR) analysis in liver tissue showed that the expression of gga-miR-15a was significantly higher in the high FCR birds than that in the medium or low FCR birds. Bioinformatics analysis further revealed that gga-mir-15a could act on many target genes, such as forkhead box O1 (FOXO1) that is involved in the insulin-signaling pathway, which influences nutrient metabolism in many organisms. Additionally, some suggestively significant variants, located on GGA3 and GGA9, were identified to associate with FI and RFI. Conclusions This GWA analysis was conducted on feed intake and efficiency traits for chickens and was innovative for application in the late laying period. Our findings can be used as a reference in the genomic breeding programs for increasing the efficiency performance of old hens and to improve our understanding of the molecular determinants for feed efficiency. Electronic supplementary material The online version of this article (10.1186/s12864-017-4092-9) contains supplementary material, which is available to authorized users.

Background DNA copy number profiles from microarray and sequencing experiments sometimes contain wave artefacts which may be introduced during sample preparation and cannot be removed completely by existing preprocessing methods. Besides, large derivative log ratio spread (DLRS) of the probes correlating with poor DNA quality is sometimes observed in genome screening experiments and may lead to unreliable copy number profiles. Depending on the extent of these artefacts and the resulting misidentification of copy number alterations/variations (CNA/CNV), it may be desirable to exclude such samples from analyses or to adapt the downstream data analysis strategy accordingly. Results Here, we propose a method to distinguish reliable genomic copy number profiles from those containing heavy wave artefacts and/or large DLRS. We define four features that adequately summarize the copy number profiles for reliability assessment, and train a classifier on a dataset of 1522 copy number profiles from various microarray platforms. The method can be applied to predict the reliability of copy number profiles irrespective of the underlying microarray platform and may be adapted for those sequencing platforms from which copy number estimates could be computed as a piecewise constant signal. Further details can be found at https://github.com/baudisgroup/CNARA. Conclusions We have developed a method for the assessment of genomic copy number profiling data, and suggest to apply the method in addition to and after other state-of-the-art noise correction and quality control procedures. CNARA could be instrumental in improving the assessment of data used for genomic data mining experiments and support the reliable functional attribution of copy number aberrations especially in cancer research. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3074-7) contains supplementary material, which is available to authorized users.

Background Most studies on the paper mulberry are mainly focused on the medicated and pharmacology, fiber quality, leaves feed development, little is known about its mechanism of adaptability to abiotic stress. Physiological measurement, transcriptomics and proteomic analysis were employed to understand its response to cold stress in this study. Methods The second to fourth fully expanded leaves from up to down were harvested at different stress time points forthe transmission electron microscope (TEM) observation. Physiological characteristics measurement included the relative electrolyte leakage (REL), SOD activity assay, soluble sugar content, and Chlorophyll fluorescence parameter measurement. For screening of differentially expressed genes, the expression level of every transcript in each sample was calculated by quantifying the number of Illumina reads. To identify the differentially expressed protein, leaves of plants under 0, 6, 12, 24, 48 and 72 h cold stress wereharvested for proteomic analysis. Finally, real time PCR was used to verify the DEG results of the RNA-seq and the proteomics data. Results Results showed that at the beginning of cold stress, respiratory metabolism was decreased and the transportation and hydrolysis of photosynthetic products was inhibited, leading to an accumulation of starch in the chloroplasts. Total of 5800 unigenes and 38 proteins were affected, including the repressed expression of photosynthesis and the enhanced expression in signal transduction, stress defense pathway as well as secondary metabolism. Although the transcriptional level of a large number of genes has been restored after 12 h, sustained cold stress brought more serious injury to the leaf cells, including the sharp rise of the relative electrolyte leakage, the declined Fv/Fm value, swelled chloroplast and the disintegrated membrane system. Conclusion The starch accumulation and the photoinhibition might be the main adaptive mechanism of the paper mulberry responded to cold stress. Most of important, enhancing the transport and hydrolysis of photosynthetic products could be the potential targets for improving the cold tolerance of the paper mulberry. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2047-6) contains supplementary material, which is available to authorized users.

Background In biological systems, metabolomics can not only contribute to the discovery of metabolic signatures for disease diagnosis, but is very helpful to illustrate the underlying molecular disease-causing mechanism. Therefore, identification of disease-related metabolites is of great significance for comprehensively understanding the pathogenesis of diseases and improving clinical medicine. Results In the paper, we propose a disease and literature driven metabolism prediction model (DLMPM) to identify the potential associations between metabolites and diseases based on latent factor model. We build the disease glossary with disease terms from different databases and an association matrix based on the mapping between diseases and metabolites. The similarity of diseases and metabolites is used to complete the association matrix. Finally, we predict potential associations between metabolites and diseases based on the matrix decomposition method. In total, 1,406 direct associations between diseases and metabolites are found. There are 119,206 unknown associations between diseases and metabolites predicted with a coverage rate of 80.88%. Subsequently, we extract training sets and testing sets based on data increment from the database of disease-related metabolites and assess the performance of DLMPM on 19 diseases. As a result, DLMPM is proven to be successful in predicting potential metabolic signatures for human diseases with an average AUC value of 82.33%. Conclusion In this paper, a computational model is proposed for exploring metabolite-disease pairs and has good performance in predicting potential metabolites related to diseases through adequate validation. The results show that DLMPM has a better performance in prioritizing candidate diseases-related metabolites compared with the previous methods and would be helpful for researchers to reveal more information about human diseases.

Background Bacteria are an essential component of the earth`s biota and affect circulation of matters through their metabolic activity. They also play an important role in the carbon and nitrogen cycle in the deep-sea environment. In this paper, two strains from deep-sea sediments were investigated in order to understand nitrogen cycling involved in the deep-sea environment. Results In this paper, the basic genomic information of two strains was obtained by whole genome sequencing. The Cobetia amphilecti N-80 and Halomonas profundus 13 genome sizes are 4,160,095 bp with a GC content of 62.5% and 5,251,450 bp with a GC content of 54.84%. Through a comparison of functional analyses, we predicted the possible C and N metabolic pathways of the two strains and determined that Halomonas profundus 13 could use more carbon sources than Cobetia amphilecti N-80. The main genes associated with N metabolism in Halomonas profundus 13 are narG, narY, narI, nirS, norB, norC, nosZ, and nirD. On the contrast, nirD, using NH4+ for energy, plays a main role in Cobetia amphilecti N-80. Both of them have the same genes for fixing inorganic carbon: icd, ppc, fdhA, accC, accB, accD, and accA. Conclusion In this study, the whole genomes of two strains were sequenced to clarify the basic characteristics of their genomes, laying the foundation for further studying nitrogen-metabolizing bacteria. Halomonas profundus 13 can utilize more carbon sources than Cobetia amphilecti N-80, as indicated by API as well as COG and KEGG prediction results. Finally, through the analysis of the nitrification and denitrification abilities as well as the inorganic carbon fixation ability of the two strains, the related genes were identified, and the possible metabolic pathways were predicted. Together, these results provide molecular markers and theoretical support for the mechanisms of inorganic carbon fixation by deep-sea microorganisms.

Background The genome architecture mapping (GAM) technique can capture genome-wide chromatin interactions. However, besides the known systematic biases in the raw GAM data, we have found a new type of systematic bias. It is necessary to develop and evaluate effective normalization methods to remove all systematic biases in the raw GAM data. Results We have detected a new type of systematic bias, the fragment length bias, in the genome architecture mapping (GAM) data, which is significantly different from the bias of window detection frequency previously mentioned in the paper introducing the GAM method but is similar to the bias of distances between restriction sites existing in raw Hi-C data. We have found that the normalization method (a normalized variant of the linkage disequilibrium) used in the GAM paper is not able to effectively eliminate the new fragment length bias at 1 Mb resolution (slightly better at 30 kb resolution). We have developed an R package named normGAM for eliminating the new fragment length bias together with the other three biases existing in raw GAM data, which are the biases related to window detection frequency, mappability, and GC content. Five normalization methods have been implemented and included in the R package including Knight-Ruiz 2-norm (KR2, newly designed by us), normalized linkage disequilibrium (NLD), vanilla coverage (VC), sequential component normalization (SCN), and iterative correction and eigenvector decomposition (ICE). Conclusions Based on our evaluations, the five normalization methods can eliminate the four biases existing in raw GAM data, with VC and KR2 performing better than the others. We have observed that the KR2-normalized GAM data have a higher correlation with the KR-normalized Hi-C data on the same cell samples indicating that the KR-related methods are better than the others for keeping the consistency between the GAM and Hi-C experiments. Compared with the raw GAM data, the normalized GAM data are more consistent with the normalized distances from the fluorescence in situ hybridization (FISH) experiments. The source code of normGAM can be freely downloaded from http://dna.cs.miami.edu/normGAM/.

Background Our knowledge about the ecological role of bacterial antimicrobial peptides (bacteriocins) in the human gut is limited, particularly in relation to their role in the diversification of the gut microbiota during early life. The aim of this paper was therefore to address associations between bacteriocins and bacterial diversity in the human gut microbiota. To investigate this, we did an extensive screening of 2564 healthy human gut metagenomes for the presence of predicted bacteriocin-encoding genes, comparing bacteriocin gene presence to strain diversity and age. Results We found that the abundance of bacteriocin genes was significantly higher in infant-like metagenomes ( Conclusions Our findings indicate that bacteriocins may play a role in the strain diversification during the infant gut microbiota establishment.

Background Piwi-interacting RNAs (piRNAs) have been proven to be closely associated with human diseases. The identification of the potential associations between piRNA and disease is of great significance for complex diseases. Traditional “wet experiment” is time-consuming and high-priced, predicting the piRNA-disease associations by computational methods is of great significance. Methods In this paper, a method based on the embedding transformation graph convolution network is proposed to predict the piRNA-disease associations, named ETGPDA. Specifically, a heterogeneous network is constructed based on the similarity information of piRNA and disease, as well as the known piRNA-disease associations, which is applied to extract low-dimensional embeddings of piRNA and disease based on graph convolutional network with an attention mechanism. Furthermore, the embedding transformation module is developed for the problem of embedding space inconsistency, which is lightweighter, stronger learning ability and higher accuracy. Finally, the piRNA-disease association score is calculated by the similarity of the piRNA and disease embedding. Results Evaluated by fivefold cross-validation, the AUC of ETGPDA achieves 0.9603, which is better than the other five selected computational models. The case studies based on Head and neck squamous cell carcinoma and Alzheimer’s disease further prove the superior performance of ETGPDA. Conclusions Hence, the ETGPDA is an effective method for predicting the hidden piRNA-disease associations.

An amendment to this paper has been published and can be accessed via the original article.

Background Genomic selection (GS) is revolutionizing plant and animal breeding. However, still its practical implementation is challenging since it is affected by many factors that when they are not under control make this methodology not effective. Also, due to the fact that it is formulated as a regression problem in general has low sensitivity to select the best candidate individuals since a top percentage is selected according to a ranking of predicted breeding values. Results For this reason, in this paper we propose two methods to improve the prediction accuracy of this methodology. One of the methods consist in reformulating the GS (nowadays formulated as a regression problem) methodology as a binary classification problem. The other consists only in a postprocessing step that adjust the threshold used for classification of the lines predicted in its original scale (continues scale) to guarantee similar sensitivity and specificity. The postprocessing method is applied for the resulting predictions after obtaining the predictions using the conventional regression model. Both methods assume that we defined with anticipation a threshold, to divide the training data as top lines and not top lines, and this threshold can be decided in terms of a quantile (for example 80%, 90%, etc.) or as the average (or maximum) of the performance of the checks. In the reformulation method it is required to label as one those lines in the training set that are equal or larger than the specified threshold and as zero otherwise. Then we train a binary classification model with the conventional inputs, but using the binary response variable in place of the continuous response variable. The training of the binary classification should be done to guarantee a more similar sensitivity and specificity, to guarantee a reasonable probability of classification of the top lines. Conclusions We evaluated the proposed models in seven data sets and we found that the two proposed methods outperformed by large margin the conventional regression model (by 402.9% in terms of sensitivity, by 110.04% in terms of F1 score and by 70.96% in terms of Kappa coefficient, with the postprocessing methods). However, between the two proposed methods the postprocessing method was better than the reformulation as binary classification model. The simple postprocessing method to improve the accuracy of the conventional genomic regression models avoid the need to reformulate the conventional regression models as binary classification models with similar or better performance, that significantly improve the selection of the top best candidate lines. In general both proposed methods are simple and can easily be adopted for use in practical breeding programs, with the guarantee that will improve significantly the selection of the top best candidates lines.

Background Cross breeding is an important way to improve livestock performance. As an important livestock and poultry resource in Henan Province of China, Bohuai goat was formed by crossing Boer goat and Huai goat. After more than 20 years of breeding, BoHuai goats showed many advantages, such as fast growth, good reproductive performance, and high meat yield. In order to better develop and protect Bohuai goats, we sequenced the whole genomes of 30 BoHuai goats and 5 Huai goats to analyze the genetic diversity, population structure and genomic regions under selection of BoHuai goat. Furthermore, we used 126 published genomes of world-wide goat to characterize the genomic variation of BoHuai goat. Results The results showed that the nucleotide diversity of BoHuai goats was lower and the degree of linkage imbalance was higher than that of other breeds. The analysis of population structure showed that BoHuai goats have obvious differences from other goat breeds. In addition, the BoHuai goat is more closely related to the Boer goat than the Huai goat and is highly similar to the Boer goat. Group by selection signal in the BoHuai goat study, we found that one region on chromosome 7 shows a very strong selection signal, which suggests that it could well be the segment region under the intense artificial selection results. Through selective sweeps, we detected some genes related to important traits such as lipid metabolism (LDLR, STAR, ANGPTL8), fertility (STAR), and disease resistance (CD274, DHPS, PDCD1LG2). Conclusion In this paper, we elucidated the genomic variation, ancestry composition, and selective signals related to important economic traits in BoHuai goats. Our studies on the genome of BoHuai goats will not only help to understand the characteristics of the crossbred but also provide a basis for the improvement of cross-breeding programs.

Although single molecule sequencing is still improving, the lengths of the generated sequences are inevitably an advantage in genome assembly. Prior work that utilizes long reads to conduct genome assembly has mostly focused on correcting sequencing errors and improving contiguity of de novo assemblies. In this paper, we propose a disassembling-reassembling approach for both correcting structural errors in the draft assembly and scaffolding a target assembly based on error-corrected single molecule sequences. In particular, the disassembling stage identifies potential structural errors in the draft assembly and disassembles the draft assembly to remove such erroneous positions, ensuring that the resulting contigs (called validated segments in this paper) are structurally consistent with the new sequence data. The reassembling stage bridges the validated segments to build larger contigs by utilizing the remaining unaligned long reads in the disassembling stage, aiming to improve the contiguity of the final assembly. To achieve this goal, we formulate a maximum alternating path cover problem. We prove that this problem is NP-hard, and solve it by a 2-approximation algorithm. Our experimental results show that our approach can improve the structural correctness of target assemblies in the cost of some contiguity, even with smaller amounts of long reads. Finally, we show that our reassembling process can also serve as a competitive scaffolder using a well-established assembly benchmark.

Background Paeonia ludlowii (Stern & G. Taylor D.Y. Hong) belongs to the peony group of the genus Paeonia in the Paeoniaceae family and is now classified as a “critically endangered species” in China. Reproduction is important for this species, and its low fruiting rate has become a critical factor limiting both the expansion of its wild population and its domestic cultivation. Results In this study, we investigated possible causes of the low fruiting rate and ovule abortion in Paeonia ludlowii. We clarified the characteristics of ovule abortion and the specific time of abortion in Paeonia ludlowii, and used transcriptome sequencing to investigate the mechanism of abortion of ovules in Paeonia ludlowii. Conclusions In this paper, the ovule abortion characteristics of Paeonia ludlowii were systematically studied for the first time and provide a theoretical basis for the optimal breeding and future cultivation of Paeonia ludlowii.

Background Gmelina arborea Roxb is a fast-growing tree species of commercial importance for tropical countries due to multiple industrial uses of its wood. Wood is primarily composed of thick secondary cell walls of xylem cells which imparts the strength to the wood. Identification of the genes involved in the secondary cell wall biosynthesis as well as their cognate regulators is crucial to understand how the production of wood occurs and serves as a starting point for developing breeding strategies to produce varieties with improved wood quality, better paper pulping or new potential uses such as biofuel production. In order to gain knowledge on the molecular mechanisms and gene regulation related with wood development in white teak, a de novo sequencing and transcriptome assembly approach was used employing secondary cell wall synthesizing cells from young white teak trees. Results For generation of transcriptome, RNA-seq reads were assembled into 110,992 transcripts and 49,364 genes were functionally annotated using plant databases; 5071 GO terms and 25,460 SSR markers were identified within xylem transcripts and 10,256 unigenes were assigned to KEGG database in 130 pathways. Among transcription factor families, C2H2, C3H, bLHLH and MYB were the most represented in xylem. Differential gene expression analysis using leaves as a reference was carried out and a total of 20,954 differentially expressed genes were identified including monolignol biosynthetic pathway genes. The differential expression of selected genes (4CL, COMT, CCoAOMT, CCR and NST1) was validated using qPCR. Conclusions We report the very first de novo transcriptome of xylem-related genes in this tropical timber species of commercial importance and constitutes a valuable extension of the publicly available transcriptomic resource aimed at fostering both basic and breeding studies.

Background The study of protein complexes and protein functional modules has become an important method to further understand the mechanism and organization of life activities. The clustering algorithms used to analyze the information contained in protein-protein interaction network are effective ways to explore the characteristics of protein functional modules. Results This paper conducts an intensive study on the problems of low recognition efficiency and noise in the overlapping structure of protein functional modules, based on topological characteristics of PPI network. Developing a protein function module recognition method ECTG based on Topological Features and Gene expression data for Protein Complex Identification. Conclusions The algorithm can effectively remove the noise data reflected by calculating the topological structure characteristic values in the PPI network through the similarity of gene expression patterns, and also properly use the information hidden in the gene expression data. The experimental results show that the ECTG algorithm can detect protein functional modules better.

Background Both the protein domains and transcript structures influence protein functional variation. The genomic location of both protein domains and transcript structural features can be described using the genomic coordinates of their encoded sequences. However, the coordinates of protein domains and transcriptional features often differ greatly, and it is difficult to view them in combination at the genome-wide level. In this paper, we describe the development of a new tool that allows users to visualize domains and transcript features together, using either built-in or uploaded genome datasets, and export publication-ready figures. Results We developed a user-friendly, independent R package and Shiny web application named “VisProDom”. VisProDom consists of a genome-wide database containing entire annotated transcripts merged with annotated protein domains from the Pfam database. The built-in dataset includes 82 files, which merge genome general feature format (GFF) annotations with rpsblast tabular outputs from protein sequence searches in the Pfam database. Multiple genomes can be simultaneously screened for protein domains or transcript names. VisProDom includes step-by-step introductions and clickable elements for ease of use. Conclusion VisProDom can display hundreds of transcripts alongside protein domains and export figures in a publication-ready format. This makes it a valuable tool for molecular evolution and comparative genomics.

Background Crossbreeding is an important way to improve production beef cattle performance. Pinan cattle is a new hybrid cattle obtained from crossing Piedmontese bulls with Nanyang cows. After more than 30 years of cross-breeding, Pinan cattle show a variety of excellent characteristics, including fast growth, early onset of puberty, and good meat quality. In this study, we analyzed the genetic diversity, population structure, and genomic region under the selection of Pinan cattle based on whole-genome sequencing data of 30 Pinan cattle and 169 published cattle genomic data worldwide. Results Estimating ancestry composition analysis showed that the composition proportions for our Pinan cattle were mainly Piedmontese and a small amount of Nanyang cattle. The analyses of nucleotide diversity and linkage disequilibrium decay indicated that the genomic diversity of Pinan cattle was higher than that of European cattle and lower than that of Chinese indigenous cattle. De-correlated composite of multiple selection signals, which combines four different statistics including θπ, CLR, FST, and XP-EHH, was computed to detect the signatures of selection in the Pinan cattle genome. A total of 83 genes were identified, affecting many economically important traits. Functional annotation revealed that these selected genes were related to immune (BOLA-DQA2, BOLA-DQB, LSM14A, SEC13, and NAALADL2), growth traits (CYP4A11, RPL26, and MYH10), embryo development (REV3L, NT5E, CDX2, KDM6B, and ADAMTS9), hornless traits (C1H21orf62), and climate adaptation (ANTXR2). Conclusion In this paper, we elucidated the genomic characteristics, ancestry composition, and selective signals related to important economic traits in Pinan cattle. These results will provide the basis for further genetic improvement of Pinan cattle and reference for other hybrid cattle related studies.

Background Third-generation single molecule sequencing technologies can sequence long reads, which is advancing the frontiers of genomics research. However, their high error rates prohibit accurate and efficient downstream analysis. This difficulty has motivated the development of many long read error correction tools, which tackle this problem through sampling redundancy and/or leveraging accurate short reads of the same biological samples. Existing studies to asses these tools use simulated data sets, and are not sufficiently comprehensive in the range of software covered or diversity of evaluation measures used. Results In this paper, we present a categorization and review of long read error correction methods, and provide a comprehensive evaluation of the corresponding long read error correction tools. Leveraging recent real sequencing data, we establish benchmark data sets and set up evaluation criteria for a comparative assessment which includes quality of error correction as well as run-time and memory usage. We study how trimming and long read sequencing depth affect error correction in terms of length distribution and genome coverage post-correction, and the impact of error correction performance on an important application of long reads, genome assembly. We provide guidelines for practitioners for choosing among the available error correction tools and identify directions for future research. Conclusions Despite the high error rate of long reads, the state-of-the-art correction tools can achieve high correction quality. When short reads are available, the best hybrid methods outperform non-hybrid methods in terms of correction quality and computing resource usage. When choosing tools for use, practitioners are suggested to be careful with a few correction tools that discard reads, and check the effect of error correction tools on downstream analysis. Our evaluation code is available as open-source at https://github.com/haowenz/LRECE.

In this introduction article, we summarize the 2020 International Conference on Intelligent Biology and Medicine (ICIBM 2020) conference which was held on August 9–10, 2020 (virtual conference). We then briefly describe the nine research articles included in this supplement issue. ICIBM 2020 hosted four scientific sections covering current topics in bioinformatics, computational biology, genomics, biomedical informatics, among others. A total of 75 original manuscripts were submitted to ICIBM 2020. All the papers were under rigorous review (at least three reviewers), and highly ranked manuscripts were selected for oral presentation and supplement issues. This genomics supplement issue included nine manuscripts. These articles cover methods and applications for single cell RNA sequencing, multi-omics data integration for gene regulation, gene fusion detection from long-read RNA sequencing, gene co-expression analysis of metabolic pathways in cancer, integrative genome-wide association studies (GWAS) of subcortical imaging phenotype in Alzheimer’s disease, as well as deep learning methods for protein structure prediction, metabolic pathway membership inference, and horizontal gene transfer (HGT) insertion sites prediction.

Background DNA N4-methylcytosine (4mC) is a critical epigenetic modification and has various roles in the restriction-modification system. Due to the high cost of experimental laboratory detection, computational methods using sequence characteristics and machine learning algorithms have been explored to identify 4mC sites from DNA sequences. However, state-of-the-art methods have limited performance because of the lack of effective sequence features and the ad hoc choice of learning algorithms to cope with this problem. This paper is aimed to propose new sequence feature space and a machine learning algorithm with feature selection scheme to address the problem. Results The feature importance score distributions in datasets of six species are firstly reported and analyzed. Then the impact of the feature selection on model performance is evaluated by independent testing on benchmark datasets, where ACC and MCC measurements on the performance after feature selection increase by 2.3% to 9.7% and 0.05 to 0.19, respectively. The proposed method is compared with three state-of-the-art predictors using independent test and 10-fold cross-validations, and our method outperforms in all datasets, especially improving the ACC by 3.02% to 7.89% and MCC by 0.06 to 0.15 in the independent test. Two detailed case studies by the proposed method have confirmed the excellent overall performance and correctly identified 24 of 26 4mC sites from the C.elegans gene, and 126 out of 137 4mC sites from the D.melanogaster gene. Conclusions The results show that the proposed feature space and learning algorithm with feature selection can improve the performance of DNA 4mC prediction on the benchmark datasets. The two case studies prove the effectiveness of our method in practical situations.

Background Inflammation has been considered to be central to the onset, progression, and outcome of infectious diseases, especially as one of the hallmarks of cancer. Non-coding RNAs (ncRNAs), such as miRNAs and lncRNAs, have emerged as vital regulators in control of immune and inflammatory processes, and also play important roles in the inflammatory disease and immunotherapy. Results In this study, we presented a database ncRI, which documented experimentally verified ncRNAs in inflammatory diseases, from published articles. Each entry contained the detailed information about ncRNA name, inflammatory diseases, mechanism, experimental techniques (e.g., microarray, RNA-seq, qRT-PCR), experimental samples (cell line and/or tissue), expression patterns of ncRNA (up-regulated or down-regulated), reference information (PubMed ID, year of publication, title of paper) and so on. Collectively, ncRI recorded 11,166 entries that include 1976 miRNAs, 1377 lncRNAs and 107 other ncRNAs across 3 species (human, mouse, and rat) from more than 2000 articles. All these data are free for users to search, browse and download. Conclusion In summary, the presented database ncRI provides a relatively comprehensive credible repository about ncRNAs and their roles in inflammatory diseases, and will be helpful for research on immunotherapy. The ncRI is now freely available to all users at http://www.jianglab.cn/ncRI/.

MotivationRepetitive DNA is a major component of plant genomes and is thought to be a driver of evolutionary novelty. Describing variation in repeat content among individuals and between populations is key to elucidating the evolutionary significance of repetitive DNA. However, the cost of producing references genomes has limited large-scale intraspecific comparisons to a handful of model organisms where multiple reference genomes are available.ResultsWe examine repeat content variation in the genomes of 94 elite inbred maize lines using graph-based repeat clustering, a reference-free and rapid assay of repeat content. We examine population structure using genome-wide repeat profiles and demonstrate the stiff-stalk and non-stiff-stalk heterotic populations are homogenous with regard to global repeat content. In contrast and similar to previously reported results, the same individuals show clear differentiation, and aggregate into two populations, when examining population structure using genome-wide SNPs. Additionally, we develop a novel kmer based technique to examine the chromosomal distribution of repeat clusters in silico and show a cluster dependent statistically significant association with gene density.ConclusionOur results indicate that repeat content variation in the heterotic populations of maize has not diverged and is uncoupled from population stratification at SNP loci. We also show that repeat families exhibit divergent patterns with regard to chromosomal distribution, some repeat clusters accumulate in regions of high gene density, whereas others aggregate in regions of low gene density.Author’s contributionsSRB and AB conceived the study, SRB performed the bioinformatic analysis, SRB wrote the paper with input from AB. email contacts: Simon Renny-Byfield: simon.renny-byfield@corteva.com, Andy Baumgarten: andy.baumgarten@corteva.com

Background Zinc finger homeodomain (ZHD) protein is a plant-specific transcription factor and a potential regulator of phosphoenolpyruvate carboxylase (PEPCase)-coding genes, and it also participates in plant growth regulation and abiotic stress responses. To study the function of MsZF-HD genes in the alkaline stress response, this paper assessed biological information and performed transcriptome analysis of the MsZF-HD gene family by using the genomes of two different varieties of alfalfa (XinJiangDa Ye and Zhongmu No. 1). Results In total, 49 and 11 MsZF-HD genes were identified in the two different varieties respectively, including the alleles of XinJiangDa Ye. According to their phylogenetic relationships, the 60 MsZF-HD genes were divided into 5 ZHD subfamilies and 1 MIF subfamily. A total of 88.3% of MsZF-HD genes do not contain introns and are unevenly distributed among the 6 chromosomes of alfalfa. A collinearity analysis indicated that 26 genes of XinJiangDa Ye have no orthologous genes in Zhongmu No. 1, although these genes (such as ZHD-X1–2, ZHD-X3–2 and ZHD-X4–2) have homologous genes in Arabidopsis thaliana, Medicago truncatula and Glycine max. Through RNA-seq and qRT–PCR verification, it was found that MsZF-HD genes are downregulated to participate in the alkaline stress response. Conclusion The results of this study may lay the foundation for the cloning and functional study of MsZF-HD genes and provide a theoretical basis for revealing the difference between XinJiangDa Ye and Zhongmu No. 1 at the genome level.

An amendment to this paper has been published and can be accessed via the original article.

Background Calliptamus italicus is a dominant species in the desert and semi-desert grassland. It is widely distributed throughout many regions such as Asia, Europe, North Africa and the Mediterranean, and has enormous destructive potential for agriculture and animal husbandry. The C. italicus overwintering as eggs in the soil through diapause, and the cold tolerance of locust eggs is the key to their ability to survive the winter smoothly to maintain the population. Results Transcriptome analysis of C. italicus eggs was carried out in this paper in constant low temperature acclimation, natural low temperature acclimation and room temperature. The differentially expressed genes related to cold tolerance were screened out, the differences in expression patterns under different low temperature acclimation were analyzed, and the genes in the significantly up-regulated pathways may play an important role in cold tolerance. The results show that different domestication modes can induce C. italicus eggs to express a large number of genes to alleviate low temperature damage, but C. italicus eggs are more sensitive to changes in temperature. Compared with the control, there are 8689 DEGs at constant low temperature and 14,994 DEGs at natural low temperature. KEGG analysis showed that DEGs were mainly enriched in pathways related to metabolism and biological systems under constant low temperature, and were mainly enriched in pathways related to biological systems and environmental information processing under natural low temperature. In addition, RNAi technology was used to further verify the regulation of genes in the significantly enriched up-regulated pathways on C. italicus eggs, and it was confirmed that the hatching rate of C. italicus eggs at low temperature was significantly reduced after interference. Conclusions Transcriptome analysis of C. italicus eggs treated at different temperatures provided a theoretical basis for further understanding the adaptation mechanism of C. italicus eggs to low temperature. In addition, four potential RNAi target genes were verified in the eggs of C. italicus for the first time, providing new ideas for effective control of this species.

Background Northern blotting is still used as a gold standard for validation of the data obtained from high-throughput whole transcriptome-based methods. However, its disadvantages of lower sensitivity, labor-intensive operation, and higher quality of RNA required limit its utilization in a routine molecular biology laboratory to monitor gene expression at RNA level. Therefore, it is necessary to optimize the traditional Northern protocol to make the technique more applicable for standard use. Results In this paper, we report modifications and tips used to improve the traditional Northern protocol for the detection of mRNAs in total RNA. To maximize the retention of specifically bound radiolabeled probes on the blot, posthybridization washes were performed under only with moderate-stringency until the level of radioactivity retained on the filter decreased to 20~50 counts per second, rather than normally under high and low stringency sequentially for scheduled time or under only high stringent condition. Successful detection of the low-expression gene using heterologous DNA probes in 20 µg of total RNA after a two-day exposure suggested an improvement in detection sensitivity. Quantitatively controlled posthybridization washes combined with an ethidium bromide-prestaining RNA procedure to directly visualize prestained RNA bands at any time during electrophoresis or immediately after electrophoresis, which made the progress of the Northern procedure to be monitored and evaluated step by step, thereby making the experiment reliable and controllable. We also report tips used in the modified Northern protocol, including the moderate concentration of formaldehyde in the gel, the accessory capillary setup, and the staining jar placed into an enamel square tray with a lid used for hybridization. Using our modified Northern protocol, eight rounds of rehybridization could be performed on a single blot. The modification made and tips used ensured the efficient proceeding of the experiment and the resulting good performance, but without using special reagents or equipment. Conclusions The modified Northern protocol improved detection sensitivity and made the experiment easy, less expensive, reliable, and controllable, and can be employed in a routine molecular biology laboratory to detect low-expressed mRNAs with heterologous DNA probes in total RNA.

Background Semantic similarity between Gene Ontology (GO) terms is a fundamental measure for many bioinformatics applications, such as determining functional similarity between genes or proteins. Most previous research exploited information content to estimate the semantic similarity between GO terms; recently some research exploited word embeddings to learn vector representations for GO terms from a large-scale corpus. In this paper, we proposed a novel method, named GO2Vec, that exploits graph embeddings to learn vector representations for GO terms from GO graph. GO2Vec combines the information from both GO graph and GO annotations, and its learned vectors can be applied to a variety of bioinformatics applications, such as calculating functional similarity between proteins and predicting protein-protein interactions. Results We conducted two kinds of experiments to evaluate the quality of GO2Vec: (1) functional similarity between proteins on the Collaborative Evaluation of GO-based Semantic Similarity Measures (CESSM) dataset and (2) prediction of protein-protein interactions on the Yeast and Human datasets from the STRING database. Experimental results demonstrate the effectiveness of GO2Vec over the information content-based measures and the word embedding-based measures. Conclusion Our experimental results demonstrate the effectiveness of using graph embeddings to learn vector representations from undirected GO and GOA graphs. Our results also demonstrate that GO annotations provide useful information for computing the similarity between GO terms and between proteins.

Background Recently, erdafitinib (Balversa), the first targeted therapy drug for genetic alteration, was approved to metastatic urothelial carcinoma. Cancer genomics research has been greatly encouraged. Currently, a large number of gene regulatory networks between different states have been constructed, which can reveal the difference states of genes. However, they have not been applied to the subtypes of Muscle-invasive bladder cancer (MIBC). Results In this paper, we propose a method that construct gene regulatory networks under different molecular subtypes of MIBC, and analyse the regulatory differences between different molecular subtypes. Through differential expression analysis and the differential network analysis of the top 100 differential genes in the network, we find that SERPINI1, NOTUM, FGFR1 and other genes have significant differences in expression and regulatory relationship between MIBC subtypes. Conclusions Furthermore, pathway enrichment analysis and differential network analysis demonstrate that Neuroactive ligand-receptor interaction and Cytokine-cytokine receptor interaction are significantly enriched pathways, and the genes contained in them are significant diversity in the subtypes of bladder cancer.

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Background Sanghuangporus sanghuang is a well-known traditional medicinal mushroom associated with mulberry. Despite the properties of this mushroom being known for many years, the regulatory mechanisms of bioactive compound biosynthesis in this medicinal mushroom are still unclear. Lysine malonylation is a posttranslational modification that has many critical functions in various aspects of cell metabolism. However, at present we do not know its role in S. sanghuang. In this study, a global investigation of the lysine malonylome in S. sanghuang was therefore carried out. Results In total, 714 malonyl modification sites were matched to 255 different proteins. The analysis indicated that malonyl modifications were involved in a wide range of cellular functions and displayed a distinct subcellular localization. Bioinformatics analysis indicated that malonylated proteins were engaged in different metabolic pathways, including glyoxylate and dicarboxylate metabolism, glycolysis/gluconeogenesis, and the tricarboxylic acid (TCA) cycle. Notably, a total of 26 enzymes related to triterpene and polysaccharide biosynthesis were found to be malonylated, indicating an indispensable role of lysine malonylation in bioactive compound biosynthesis in S. sanghuang. Conclusions These findings suggest that malonylation is associated with many metabolic pathways, particularly the metabolism of the bioactive compounds triterpene and polysaccharide. This paper represents the first comprehensive survey of malonylation in S. sanghuang and provides important data for further study on the physiological function of lysine malonylation in S. sanghuang and other medicinal mushrooms.

Background All diseases containing genetic material undergo genetic evolution and give rise to heterogeneity including cancer and infection. Although these illnesses are biologically very different, the ability for phylogenetic retrodiction based on the genomic reads is common between them and thus tree-based principles and assumptions are shared. Just as the different frequencies of tumor genomic variants presupposes the existence of multiple tumor clones and provides a handle to computationally infer them, we postulate that the different variant frequencies in viral reads offers the means to infer multiple co-infecting sublineages. Results We present a common methodological framework to infer the phylogenomics from genomic data, be it reads of SARS-CoV-2 of multiple COVID-19 patients or bulk DNAseq of the tumor of a cancer patient. We describe the Concerti computational framework for inferring phylogenies in each of the two scenarios.To demonstrate the accuracy of the method, we reproduce some known results in both scenarios. We also make some additional discoveries. Conclusions Concerti successfully extracts and integrates information from multi-point samples, enabling the discovery of clinically plausible phylogenetic trees that capture the heterogeneity known to exist both spatially and temporally. These models can have direct therapeutic implications by highlighting “birth” of clones that may harbor resistance mechanisms to treatment, “death” of subclones with drug targets, and acquisition of functionally pertinent mutations in clones that may have seemed clinically irrelevant. Specifically in this paper we uncover new potential parallel mutations in the evolution of the SARS-CoV-2 virus. In the context of cancer, we identify new clones harboring resistant mutations to therapy.

Background The animal phylum Cnidaria consists of six classes or subphyla: Hydrozoa, Scyphozoa, Cubozoa, Staurozoa, Anthozoa, and Endocnidozoa. Cnidarians have an early evolutionary origin, diverging before the emergence of the Bilateria. Extant members from this phylum, therefore, are important resources for understanding the evolution of the nervous system. Cnidarian nervous systems are strongly peptidergic. Using genomics, we have recently shown that three neuropeptide families (the X1PRX2amides, GRFamides, and GLWamides) are wide-spread in four (Scyphozoa, Cubozoa, Staurozoa, Anthozoa) out of six cnidarian classes or subphyla, suggesting that these three neuropeptide families emerged in the common cnidarian ancestor. In the current paper, we analyze the remaining cnidarian class, Hydrozoa, and the subphylum Endocnidozoa, to make firm conclusions about the evolution of neuropeptide genes in Cnidaria. Results We analyzed sixteen hydrozoan species with a sequenced genome or transcriptome, using a recently developed software program for discovering neuropeptide genes. These species belonged to various hydrozoan subclasses and orders, among them the laboratory models Hydra, Hydractinia, and Clytia. We found that each species contained three to five neuropeptide families. A common feature for all hydrozoans was that they contained genes coding for (i) X1PRX2amide peptides, (ii) GRFamide peptides, and (iii) GLWamide peptides. These results support our previous conclusions that these three neuropeptide families evolved early in evolution. In addition to these three neuropeptide families, hydrozoans expressed up to two other neuropeptide gene families, which, however, were only occurring in certain animal groups. Endocnidozoa (Myxozoa) are microscopically small endoparasites, which are strongly reduced. For long, it was unknown to which phylum these parasites belonged, but recently they have been associated with cnidarians. We analyzed nine endocnidozoan species and found that two of them (Polypodium hydriforme and Buddenbrockia plumatellae) expressed neuropeptide genes. These genes coded for neuropeptides belonging to the GRFamide and GLWamide families with structures closely resembling them from hydrozoans. Conclusions We found X1PRX2amide, GRFamide, and GLWamide peptides in all species belonging to the Hydrozoa, confirming that these peptides originated in the common cnidarian ancestor. In addition, we discovered GRFamide and GLWamide peptide genes in some members of the Endocnidozoa, thereby linking these parasites to Hydrozoa.

Background Patient-derived xenografts (PDX) mice models play an important role in preclinical trials and personalized medicine. Sharing data on the models is highly valuable for numerous reasons – ethical, economical, research cross validation etc. The EurOPDX Consortium was established 8 years ago to share such information and avoid duplicating efforts in developing new PDX mice models and unify approaches to support preclinical research. EurOPDX Data Portal is the unified data sharing platform adopted by the Consortium. Main body In this paper we describe the main features of the EurOPDX Data Portal (https://dataportal.europdx.eu/), its architecture and possible utilization by researchers who look for PDX mice models for their research. The Portal offers a catalogue of European models accessible on a cooperative basis. The models are searchable by metadata, and a detailed view provides molecular profiles (gene expression, mutation, copy number alteration) and treatment studies. The Portal displays the data in multiple tools (PDX Finder, cBioPortal, and GenomeCruzer in future), which are populated from a common database displaying strictly mutually consistent views. (Short) Conclusion EurOPDX Data Portal is an entry point to the EurOPDX Research Infrastructure offering PDX mice models for collaborative research, (meta)data describing their features and deep molecular data analysis according to users’ interests.

Background Paris polyphylla is a herb widely used in traditional Chinese medicine to treat various diseases. Stem rot diseases seriously affected the yield of P. polyphylla in subtropical areas of China. Therefore, cost-effective, chemical-free, eco-friendly strategies to control stem rot on P. polyphylla are valuable and urgently needed. Results In this paper, we reported the biocontrol efficiency of Paenibacillus peoriae HJ-2 and its complete genome sequence. Strain HJ-2 could serve as a potential biocontrol agent against stem rot on P. polyphylla in the greenhouse and field. The genome of HJ-2 consists of a single 6,001,192 bp chromosome with an average GC content of 45% and 5,237 predicted protein coding genes, 39 rRNAs and 108 tRNAs. The phylogenetic tree indicated that HJ-2 is most closely related to P. peoriae IBSD35. Functional analysis of genome revealed numerous genes/gene clusters involved in plant colonization, biofilm formation, plant growth promotion, antibiotic and resistance inducers synthesis. Moreover, metabolic pathways that potentially contribute to biocontrol mechanisms were identified. Conclusions This study revealed that P. peoriae HJ-2 could serve as a potential BCA against stem rot on P. polyphylla. Based on genome analysis, the genome of HJ-2 contains more than 70 genes and 12 putative gene clusters related to secondary metabolites, which have previously been described as being involved in chemotaxis motility, biofilm formation, growth promotion, antifungal activity and resistance inducers biosynthesis. Compared with other strains, variation in the genes/gene clusters may lead to different antimicrobial spectra and biocontrol efficacies.

Background With the rapid increase in the amount of Protein-Protein Interaction (PPI) data, the establishment of an event-centered PPI ontology that contains temporal and spatial vocabularies is urgently needed to clarify PPI biological annotations. In this paper, we propose a precisely designed schema - PPIO (PPI Ontology) for representing the biological context of PPIs. Results Inspired by the event model and the distinct characteristics of PPI events, PPIO consists of six core aspects of the information required for reporting a PPI event, including the interactor (who), the biological process (when), the subcellular location (where), the interaction type (how), the biological function (what) and the detection method (which). PPIO is implemented through the integration of appropriate terms from the corresponding vocabularies/ontologies, e.g., Gene Ontology, Protein Ontology, PSI-MI/MOD, etc. To assess PPIO, an approach based on PPIO in developed to extract PPI biological annotations from an open standard corpus “BioCreAtIvE-PPI”. The experiment results demonstrate PPIO’s high performance, a precision of 0.69, a recall of 0.72 and an F-score of 0.70. Conclusions PPIO is a well-constructed essential ontology in the interpretation of PPI biological context. The results of the experiments conducted on the BioCreAtIvE corpus demonstrate that PPIO is able to facilitate PPI annotation extraction from biomedical literature effectively and enrich essential annotation for PPIs.

An amendment to this paper has been published and can be accessed via the original article.

In their paper entitled “Comparative analysis of Faecalibacterium prausnitzii genomes shows a high level of genome plasticity and warrants separation into new species-level taxa” (BMC Genomics (2018) 19:931), Fitzgerald et al. proposed a neotype strain for F. prausnitzii (strain A2–165 = DSM 17677 = JCM 31915) and assigned strain ATCC 27768 = NCIMB 13872 to a newly established taxon, Faecalibacterium moorei nom. Nov. These proposals contravene the Rules of the International Code of Nomenclature of Prokaryotes (ICNP). Neotype strains can only be established following a formal proposal in the International Journal of Systematic and Evolutionary Microbiology in accordance with Rule 18c and Appendix 7 of the ICNP. A proposed neotype becomes an established neotype after 2 years, provided that no objections were submitted to the Judicial Commisson of the International Committee on Systematics of Prokaryotes within the first year following publication of the request. F. moorei as proposed by Fitzgerald et al. is a later homotypic synonym of F. prausnitzii. It cannot be a ‘nom. nov.’ (nomen novum): based on Rule 34a of the ICNP: this term is only used when an author transfers a species to another genus or a subspecies to another species as a new combination, but the original specific epithet cannot be used as ‘comb. nov.’ (combinatio nova) as a result of homonymy. Moreover, ATCC 27768 and NCIMB 13872 cannot be proposed as the type strain of F. moorei as these remain permanently associated with the type strain of F. prausnitzii unless the Judicial Commission of the ICSP will decide otherwise.

Background Determination of genome-wide DNA methylation is significant for both basic research and drug development. As a key epigenetic modification, this biochemical process can modulate gene expression to influence the cell differentiation which can possibly lead to cancer. Due to the involuted biochemical mechanism of DNA methylation, obtaining a precise prediction is a considerably tough challenge. Existing approaches have yielded good predictions, but the methods either need to combine plenty of features and prerequisites or deal with only hypermethylation and hypomethylation. Results In this paper, we propose a deep learning method for prediction of the genome-wide DNA methylation, in which the Methylation Regression is implemented by Convolutional Neural Networks (MRCNN). Through minimizing the continuous loss function, experiments show that our model is convergent and more precise than the state-of-art method (DeepCpG) according to results of the evaluation. MRCNN also achieves the discovery of de novo motifs by analysis of features from the training process. Conclusions Genome-wide DNA methylation could be evaluated based on the corresponding local DNA sequences of target CpG loci. With the autonomous learning pattern of deep learning, MRCNN enables accurate predictions of genome-wide DNA methylation status without predefined features and discovers some de novo methylation-related motifs that match known motifs by extracting sequence patterns. Electronic supplementary material The online version of this article (10.1186/s12864-019-5488-5) contains supplementary material, which is available to authorized users.

Post-translational modification (PTM), which is a biological process, tends to modify proteome that leads to changes in normal cell biology and pathogenesis. In the recent times, there has been many reported PTMs. Out of the many modifications, phosphoglycerylation has become particularly the subject of interest. The experimental procedure for identification of phosphoglycerylated residues continues to be an expensive, inefficient and time-consuming effort, even with a large number of proteins that are sequenced in the post-genomic period. Computational methods are therefore being anticipated in order to effectively predict phosphoglycerylated lysines. Even though there are predictors available, the ability to detect phosphoglycerylated lysine residues still remains inadequate. We have introduced a new predictor in this paper named EvolStruct-Phogly that uses structural and evolutionary information relating to amino acids to predict phosphoglycerylated lysine residues. Benchmarked data is employed containing experimentally identified phosphoglycerylated and non-phosphoglycerylated lysines. We have then extracted the three structural information which are accessible surface area of amino acids, backbone torsion angles, amino acid’s local structure conformations and profile bigrams of position-specific scoring matrices. EvolStruct-Phogly showed a noteworthy improvement in regards to the performance when compared with the previous predictors. The performance metrics obtained are as follows: sensitivity 0.7744, specificity 0.8533, precision 0.7368, accuracy 0.8275, and Mathews correlation coefficient of 0.6242. The software package and data of this work can be obtained from https://github.com/abelavit/EvolStruct-Phogly or www.alok-ai-lab.com

Background Next-Generation Sequencing (NGS) has been widely accepted as an essential tool in molecular biology. Reduced costs and automated analysis pipelines make the use of NGS data feasible even for small labs, yet the methods for interpreting the data are not sophisticated enough to account for the amount of information. Results We propose s ·nr, a Visual Analytics tool that provides simple yet powerful visual interfaces for displaying and querying NGS data. It allows researchers to explore their own data in the context of experimental data deposited in public repositories, as well as to extract specific data sets with similar gene expression signatures. We tested s ·nr on 1543 RNA-Seq based mouse differential expression profiles derived from the public ArrayExpress platform. We provide the repository of processed data with this paper. Conclusion s ·nr, easily deployable utilizing its containerized implementation, empowers researchers to analyze and relate their own RNA-Seq as well as to provide interactive and contextual crosstalk with data from public repositories. This allows users to deduce novel and unbiased hypotheses about the underlying molecular processes. Demo Login demo/demo: snr.sf.mpg.de (Tested with Google Chrome) Electronic supplementary material The online version of this article (10.1186/s12864-018-5396-0) contains supplementary material, which is available to authorized users.

Background Single molecule measurements of DNA polymerization kinetics provide a sensitive means to detect both secondary structures in DNA and deviations from primary chemical structure as a result of modified bases. In one approach to such analysis, deviations can be inferred by monitoring the behavior of DNA polymerase using single-molecule, real-time sequencing with zero-mode waveguide. This approach uses a Single Molecule Real Time (SMRT)-sequencing measurement of time between fluorescence pulse signals from consecutive nucleosides incorporated during DNA replication, called the interpulse duration (IPD). Results In this paper we present an analysis of loci with high IPDs in two genomes, a bacterial genome (E. coli) and a eukaryotic genome (C. elegans). To distinguish the potential effects of DNA modification on DNA polymerization speed, we paired an analysis of native genomic DNA with whole-genome amplified (WGA) material in which DNA modifications were effectively removed. Adenine modification sites for E. coli are known and we observed the expected IPD shifts at these sites in the native but not WGA samples. For C. elegans, such differences were not observed. Instead, we found a number of novel sequence contexts where IPDs were raised relative to the average IPDs for each of the four nucleotides, but for which the raised IPD was present in both native and WGA samples. Conclusion The latter results argue strongly against DNA modification as the underlying driver for high IPD segments for C. elegans, and provide a framework for separating effects of DNA modification from context-dependent DNA polymerase kinetic patterns inherent in underlying DNA sequence for a complex eukaryotic genome.

An amendment to this paper has been published and can be accessed via the original article.

This correspondence responds to the critique by Butler et al. (BMC Genomics 22:241, 2021) of our recent paper on transposable element (TE) persistence. We address the three main objections raised by Butler et al. After running a series of additional simulations that were inspired by the authors’ criticisms, we are able to present a more nuanced understanding of the conditions that generate long-term persistence.

An amendment to this paper has been published and can be accessed via the original article.

Background Motif search is an important step in extracting meaningful patterns from biological data. The general problem of motif search is intractable and there is a pressing need to develop efficient, exact and approximation algorithms to solve this problem. In this paper, we present several novel, exact, sequential and parallel algorithms for solving the (l,d) Edit-distance-based Motif Search (EMS) problem: given two integers l,d and n biological strings, find all strings of length l that appear in each input string with atmost d errors of types substitution, insertion and deletion. Methods One popular technique to solve the problem is to explore for each input string the set of all possible l-mers that belong to the d-neighborhood of any substring of the input string and output those which are common for all input strings. We introduce a novel and provably efficient neighborhood exploration technique. We show that it is enough to consider the candidates in neighborhood which are at a distance exactly d. We compactly represent these candidate motifs using wildcard characters and efficiently explore them with very few repetitions. Our sequential algorithm uses a trie based data structure to efficiently store and sort the candidate motifs. Our parallel algorithm in a multi-core shared memory setting uses arrays for storing and a novel modification of radix-sort for sorting the candidate motifs. Results The algorithms for EMS are customarily evaluated on several challenging instances such as (8,1), (12,2), (16,3), (20,4), and so on. The best previously known algorithm, EMS1, is sequential and in estimated 3 days solves up to instance (16,3). Our sequential algorithms are more than 20 times faster on (16,3). On other hard instances such as (9,2), (11,3), (13,4), our algorithms are much faster. Our parallel algorithm has more than 600 % scaling performance while using 16 threads. Conclusions Our algorithms have pushed up the state-of-the-art of EMS solvers and we believe that the techniques introduced in this paper are also applicable to other motif search problems such as Planted Motif Search (PMS) and Simple Motif Search (SMS). Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2789-9) contains supplementary material, which is available to authorized users.

Background Single-cell (sc) sequencing performs unbiased profiling of individual cells and enables evaluation of less prevalent cellular populations, often missed using bulk sequencing. However, the scale and the complexity of the sc datasets poses a great challenge in its utility and this problem is further exacerbated when working with larger datasets typically generated by consortium efforts. As the scale of single cell datasets continues to increase exponentially, there is an unmet technological need to develop database platforms that can evaluate key biological hypotheses by querying extensive single-cell datasets. Large single-cell datasets like Human Cell Atlas and COVID-19 cell atlas (collection of annotated sc datasets from various human organs) are excellent resources for profiling target genes involved in human diseases and disorders ranging from oncology, auto-immunity, as well as infectious diseases like COVID-19 caused by SARS-CoV-2 virus. SARS-CoV-2 infections have led to a worldwide pandemic with massive loss of lives, infections exceeding 7 million cases. The virus uses ACE2 and TMPRSS2 as key viral entry associated proteins expressed in human cells for infections. Evaluating the expression profile of key genes in large single-cell datasets can facilitate testing for diagnostics, therapeutics, and vaccine targets, as the world struggles to cope with the on-going spread of COVID-19 infections. Main body In this manuscript we describe REVEAL: SingleCell, which enables storage, retrieval, and rapid query of single-cell datasets inclusive of millions of cells. The array native database described here enables selecting and analyzing cells across multiple studies. Cells can be selected using individual metadata tags, more complex hierarchical ontology filtering, and gene expression threshold ranges, including co-expression of multiple genes. The tags on selected cells can be further evaluated for testing biological hypotheses. One such example includes identifying the most prevalent cell type annotation tag on returned cells. We used REVEAL: SingleCell to evaluate the expression of key SARS-CoV-2 entry associated genes, and queried the current database (2.2 Million cells, 32 projects) to obtain the results in Conclusion In this paper, we introduce the REVEAL: SingleCell database that addresses immediate needs for SARS-CoV-2 research and has the potential to be used more broadly for many precision medicine applications. We used the REVEAL: SingleCell database as a reference to ask questions relevant to drug development and precision medicine regarding cell type and co-expression for genes that encode proteins necessary for SARS-CoV-2 to enter and reproduce in cells.

BackgroundGenetic improvement of root system architecture is essential to improve water and nutrient use efficiency of crops or to boost their productivity under stress or non-optimal soil conditions. One hundred ninety-two Ethiopian durum wheat accessions comprising 167 historical landraces and 25 modern cultivars were assembled for GWAS analysis to identify QTLs for root system architecture (RSA) traits and genotyped with a high-density 90 K wheat SNP array by Illumina.ResultsUsing a non-roll, paper-based root phenotyping platform, a total of 2880 seedlings and 14,947 seminal roots were measured at the three-leaf stage to collect data for total root length (TRL), total root number (TRN), root growth angle (RGA), average root length (ARL), bulk root dry weight (RDW), individual root dry weight (IRW), bulk shoot dry weight (SDW), presence of six seminal roots per seedling (RT6) and root shoot ratio (RSR). Analysis of variance revealed highly significant differences between accessions for all RSA traits. Four major (− log10P≥ 4) and 34 nominal (− log10P≥ 3) QTLs were identified and grouped in 16 RSA QTL clusters across chromosomes. A higher number of significant RSA QTL were identified on chromosome 4B particularly for root vigor traits (root length, number and/or weight).ConclusionsAfter projecting the identified QTLs on to a high-density tetraploid consensus map along with previously reported RSA QTL in both durum and bread wheat, fourteen nominal QTLs were found to be novel and could potentially be used to tailor RSA in elite lines. The major RGA QTLs on chromosome 6AL detected in the current study and reported in previous studies is a good candidate for cloning the causative underlining sequence and identifying the beneficial haplotypes able to positively affect yield under water- or nutrient-limited conditions.