BMC Genomics最新文献

Background: Multiple sequence alignment (MSA) has proven extremely useful in computational biology, especially in inferring evolutionary relationships via phylogenetic analysis and providing insight into protein structure and function. An alternative to the standard MSA model is partial order alignment (POA), in which aligned sequences are represented as paths in a graph rather than rows in a matrix. While the POA model has proven useful in several applications (e.g. sequencing reads assembly and pangenome structure exploration), we lack efficient visualization tools that could highlight its advantages.

Results: We propose Sequence Flow - a web application designed to address the above problem. Sequence Flow presents the POA as a Sankey diagram, a kind of graph visualisation typically used for graphs representing flowcharts. Sequence Flow enables interactive alignment exploration, including fragment selection, highlighting a selected group of sequences, modification of the position of graph nodes, structure simplification etc. After adjustment, the visualization can be saved as a high-quality graphic file. Thanks to the use of SanKEY.js - a JavaScript library for creating Sankey diagrams, designed specifically to visualize POAs, Sequence Flow provides satisfactory performance even with large alignments.

Conclusions: We provide Sankey diagram-based POA visualization tools for both end users (Sequence Flow) and bioinformatic software developers (SanKEY.js). Sequence Flow webservice is available at https://sequenceflow.mimuw.edu.pl/ . The source code for SanKEY.js is available at https://github.com/Krzysiekzd/SanKEY.js and for Sequence Flow at https://github.com/Krzysiekzd/SequenceFlow .

Background: Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed.

Results: We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools.

Conclusions: This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Background: Cold is a significant limiting factor in productivity, particularly in northwestern and eastern China. Calcium-Dependent Protein Kinases (CDPKs), a primary calcium signaling sensor in plants, play an important role in their response to cold. Snow lotus (Sasussured involucrata Kar L) is a plant that thrives in harsh climates and grows in northwest China. However, there were no reports on the transcriptome of OE-SikCDPK1 transgenic tobacco in response to cold.

Results: When exposed to cold stress, OE-SikCDPK1 plants displayed a cold-tolerant phenotype compared to non-transgenic tobacco. Under cold conditions, relative water content reduced, relative conductivity increased, malondialdehyde levels rose, and cold-responsive gene expression increased. The OE-SikCDPK1 gene and non-transgenic tobacco were employed for research purposes. The transcriptome of leaves was sequenced using the HISAT2 sequencing platform, and the data were used to examine gene function annotation and differentially expressed genes (DEGs). 53,022 DEGs in tobacco leaves under cold treatment were obtained. The GO enrichment results revealed that it was enriched for biological-process, defense response and other processes under cold stress. The KEGG pathway enrichment analysis revealed that the metabolic pathways of significant enrichment of DEGs under cold stress mainly involved MAPK signaling pathway transduction. The transcription factor identification results showed that the transcription factors with the largest number of differential expressions under cold treatment were mainly from WRKY, AP2, MYB, bHLH, NAC and other transcription factor families.

Conclusion: The cold tolerance mechanism of snow lotus SikCDPK1 was comprehensively analyzed at the transcriptional level for the first time using RNA-seq technology. This study demonstrates that SikCDPK1 can respond to cold by participating in the MAPK signaling pathway and regulating the expression levels of transcription factors, including WRKY, AP2, MYB, bHLH, and NAC. These results offer valuable insights for further exploration of the cold tolerance mechanism associated with SikCDPK1.

Non-typhoidal Salmonella (NTS) is one of the top causes of diarrhea worldwide. Ceftriaxone is commonly recommended as the initial treatment option for Salmonella infections due to its antibacterial effectiveness. The objective of this study was to investigate the molecular epidemiological characteristics of NTS and to compare the phenotypic and genotypic profiles of antimicrobial resistance in multidrug-resistant Salmonella strains by sequencing 329 NTS strains collected from a county-level hospital between 2018 and 2021. Multi-locus sequence typing (MLST), antimicrobial resistance genes and plasmid types were identified by BacWGSTdb 2.0 webserver. Phylogenetic analysis of all NTS strains was carried out using Snippy and Gubbins software. The transferability of ceftriaxone resistant plasmids was confirmed through plasmid conjugation assays, and verified by S1-PFGE-Southern blot assays. The predominant serotypes among all NTS strains were Typhimurium (161/329), Enteritidis (49/329) and London (45/329). The most common sequence type observed was ST34 (86/329), followed by ST19 (72/329) and ST11 (47/329). The antimicrobial resistance of Salmonella to a wide range of antimicrobials showed an overall increase. Out of these 37 (11.24%) ceftriaxone-resistant strains, with the majority of them (33/37) being bla_CTX-M. The predominant plasmid types identified were IncHI2 (14/21) and IncI1 (6/21), ranging in size from 70 kb to 360 kb. The conjugation efficiency was calculated with the high conjugation efficiency of 1.1 × 10^- 5 to 9.3 × 10^- 2. The strains varied widely, ranging from 3 to 45,024 single nucleotide polymorphisms (SNPs). There are close linkages observed among the predominant lineage, with an average of 78 SNPs between each pair of ST34 strains. The findings contribute to our understanding of the transmission and resistance mechanisms of multidrug-resistant Salmonella, thereby facilitating the development of effective control strategies.

Prediction of chromatin interactions from DNA sequence has been a significant research challenge in the last couple of years. Several solutions have been proposed, most of which are based on encoder-decoder architecture, where 1D sequence is convoluted, encoded into the latent representation, and then decoded using 2D convolutions into the Hi-C pairwise chromatin spatial proximity matrix. Those methods, while obtaining high correlation scores and improved metrics, produce Hi-C matrices that are artificial - they are blurred due to the deep learning model architecture. In our study, we propose the HiCDiffusion, sequence-only model that addresses this problem. We first train the encoder-decoder neural network and then use it as a component of the diffusion model - where we guide the diffusion using a latent representation of the sequence, as well as the final output from the encoder-decoder. That way, we obtain the high-resolution Hi-C matrices that not only better resemble the experimental results - improving the Fréchet inception distance by an average of 11 times, with the highest improvement of 56 times - but also obtain similar classic metrics to current state-of-the-art encoder-decoder architectures used for the task.

Mustard is a commercial oilseed crop worldwide infected by a highly infectious turnip mosaic virus (TuMV). In the experimental field at ICAR-IARI, New Delhi, in 2022, a 100% incidence of TuMV infection was observed in brown, black and yellow mustard. A very low aphid population suggested the possibility of seed transmission. Earlier, the virus genome was characterized by high throughput sequencing and it was a recombinant of World-B and Asian-BR isolates. The presence of TuMV in immature seeds was confirmed in eight field-grown genotypes via RT-PCR using CP-specific primers designed from the same genome sequence. TuMV was found to be localized in embryo and cotyledon, indicating its true seed-borne nature. Presence of TuMV was also confirmed by RT-PCR in the grow out plants from seeds of field grown eight infected genotypes and 9 genotypes collected from seed stock, that were grown in an aphid-free growth chamber. Further, out of 24 seedlings of Pusa Gold (seed stock) and Pusa Karishma (seeds from field grown plants), 20 and 17 seedlings were found infected with TuMV, respectively. The internally seed-borne nature of the virus leads to its early establishment at the seedling stage, leading to stunting and leaf-puckering symptoms in the progeny plants. This study is the first evidence of seed embryo infection and seedling transmission of TuMV of all the three species of mustard plants (brown, black and yellow mustard). Seed transmission of TuMV in mustard genotypes have implications for the seed exchange programme of mustard seeds.

Background: Mesenchymal stem cells (MSCs) are multipotent stem cells that are under investigation for use in clinical trials because they are capable of self-renewal and differentiating into different cell types under defined conditions. Nonetheless, the therapeutic effects of MSCs have been constrained by low engraftment rates, cell fusion, and cell survival. Various strategies have been explored to improve the therapeutic efficacy of MSCs, with platelet-derived growth factor (PDGF)-BB emerging as a promising candidate. To enhance our comprehension of the impact of PDGF-BB on the gene expression profile and chromosomal accessibility of MSCs, RNA-sequencing and analysis of chromatin accessibility profiles were conducted on three human primary MSCs in culture, both with and without stimulation by PDGF-BB.

Results: Integrative analysis of gene expression and chromatin accessibility demonstrated that PDGF-BB treatment modified the chromatin accessibility landscape, marking regions for activation or repression through the AP-1 family transcription factors TEAD, CEBP, and RUNX2. These changes in AP-1 transcription factor expression, in turn, led to cell proliferation and differentiation potential towards osteoblasts, adipocytes, or chondrocytes. The degree of MSC differentiation varies among cells isolated from different donors. The presence of an enrichment of exosome-related genes is also noted among all the differentially expressed genes.

Conclusions: In conclusion, the observed changes in AP-1 transcription factor expression not only induced cellular proliferation and differentiation, but also revealed variations in the degree of MSC differentiation based on donor-specific differences. Moreover, the enrichment of exosome-related genes among differentially expressed genes suggests a potential significant role for PDGF-BB in facilitating intercellular communication.

Background: The blunt snout bream (Megalobrama amblycephala) is an important economic freshwater fish in China with tender flesh and high nutritional value. With the cultivation of superior new varieties and the expansion of breeding scale, it becomes imperative to employ sex-control technology to cultivate monosexual populations of M. amblycephala, thereby preventing the deterioration of desirable traits. The development of specific markers capable of accurately identifying the sex of M. amblycephala would facilitate the determination of the genetic sex of the breeding population before gonad maturation, thereby expediting the processes of sex-controlled breeding of M. amblycephala.

Results: A whole-genome re-sequencing was performed for 116 females and 141 males M. amblycephala collected from nine populations. Seven candidate male-specific sequences were identified through comparative analysis of male and female genomes, which were further compared with the sequencing data of 257 individuals, and finally three male-specific sequences were generated. These three sequences were further validated by PCR amplification in 32 males and 32 females to confirm their potential as male-specific molecular markers for M. amblycephala. One of these markers showed potential applicability in M. pellegrini as well, enabling males to be identified using this specific molecular marker.

Conclusions: The study provides a high-efficiency and cost-effective approach for the genetic sex identification in two species of Megalobrama. The developed markers in this study have great potential in facilitating sex-controlled breeding of M. amblycephala and M. pellegrini, while also contributing valuable insights into the underlying mechanisms of fish sex determination.

Salinity adaptation is an important issue in aquaculture. Understanding the immediate-early response to salinity stress helps in comprehending this process. In vitro experiments using cell lines can explain cell-independent reactions without the involvement of hormones in vivo. In this study, salinity stress experiments were conducted using cell line derived from the gills of Gymnocypris przewalskii (GPG cell line) to isolate immediate-early response-related genes and miRNAs using transcriptomics, followed by bioinformatics analysis. The results showed that intracellular free Ca²⁺ appeared to be a key factor in cell sensing and initiating downstream cell signaling in response to external salinity. Additionally, cell apoptosis was the most common feature of salinity stress, with multiple signaling pathways involved in salinity-induced cell apoptosis. Furthermore, MiRNAs played a crucial role in the rapid response to salinity stress by selectively inhibiting the expression of specific genes. Additionally, for the first time in the G. przewalskii genome, Tf2 and TY3 families of transposons were found to have responsive roles to the external salinity stress. This study contributes to a better understanding of osmotic sensing in G. przewalskii and provides theoretical assistance for improving salinity adaptation in aquaculture fish species.

Background: Cashmere, named as "soft gold", derives from the secondary hair follicles (SHFs) of cashmere goat which is vital to Northwest China's economy. The cytodifferentiation stage (E120), mirroring the complete hair follicle (HF) structure of adult goats and marking a critical phase in SHF development. Therefore, this study aims to enhance the understanding of SHF development and its impact on fiber quality, informing breeding strategies.

Results: From the scRNA-seq data analysis, the intricate processes and transcriptional dynamics of inner layer cell differentiation of HFs were unveiled in this study. we identified nine cell populations during cytodifferentiation and key structures such as the hair shaft and inner root sheath. And we discovered three main inner layer lineages and seven subpopulations, clarifying their roles in specialization and signaling. Pseudotime mapping analysis showed cell evolution from early stage to mature stages marked by unique gene expressions, and the intermediate stage on the differentiation of each lineage was revealed. The identification and spatial localization of specific transcription factors, such as GATA3, LEF1 and PRDM1, as well as keratin genes highlight regulatory pathways involved in HF development, which was further validated by immunofluorescence. These findings suggested the potential strategies to improve fiber quality, and the discovery of diverse cell types and their developmental molecular mechanisms, particularly in this species-specific context, offered a nuanced view of the regulatory mechanisms driving HF development in cashmere goats.

Conclusion: Overall, these findings provide a systematic molecular atlas of skin, defining three major branches and cell states of inner layer cells of HF, and determining how the branch-specific transcription factors, keratins, and signals coordinate HF morphogenesis during cytodifferentiation stage. This research not only advances skin tissue research in goats but also holds broader implications for the understanding of HF regeneration and development across various species.