Enhancing SNV identification in whole-genome sequencing data through the incorporation of known genetic variants into the minimap2 index.

MOTIVATION: Alignment of reads to a reference genome sequence is one of the key steps in the analysis of human whole-genome sequencing data obtained through Next-generation sequencing (NGS) technologies. The quality of the subsequent steps of the analysis, such as the results of clinical interpretation of genetic variants or the results of a genome-wide association study, depends on the correct identification of the position of the read as a result of its alignment. The amount of human NGS whole-genome sequencing data is constantly growing. There are a number of human genome sequencing projects worldwide that have resulted in the creation of large-scale databases of genetic variants of sequenced human genomes. Such information about known genetic variants can be used to improve the quality of alignment at the read alignment stage when analysing sequencing data obtained for a new individual, for example, by creating a genomic graph. While existing methods for aligning reads to a linear reference genome have high alignment speed, methods for aligning reads to a genomic graph have greater accuracy in variable regions of the genome. The development of a read alignment method that takes into account known genetic variants in the linear reference sequence index allows combining the advantages of both sets of methods. RESULTS: In this paper, we present the minimap2_index_modifier tool, which enables the construction of a modified index of a reference genome using known single nucleotide variants and insertions/deletions (indels) specific to a given human population. The use of the modified minimap2 index improves variant calling quality without modifying the bioinformatics pipeline and without significant additional computational overhead. Using the PrecisionFDA Truth Challenge V2 benchmark data (for HG002 short-read data aligned to the GRCh38 linear reference (GCA_000001405.15) with parameters k = 27 and w = 14) it was demonstrated that the number of false negative genetic variants decreased by more than 9500, and the number of false positives decreased by more than 7000 when modifying the index with genetic variants from the Human Pangenome Reference Consortium.

Deep Semantic Segmentation of Angiogenesis Images.

Angiogenesis is the development of new blood vessels from pre-existing ones. It is a complex multifaceted process that is essential for the adequate functioning of human organisms. The investigation of angiogenesis is conducted using various methods. One of the most popular and most serviceable of these methods in vitro is the short-term culture of endothelial cells on Matrigel. However, a significant disadvantage of this method is the manual analysis of a large number of microphotographs. In this regard, it is necessary to develop a technique for automating the annotation of images of capillary-like structures. Despite the increasing use of deep learning in biomedical image analysis, as far as we know, there still has not been a study on the application of this method to angiogenesis images. To the best of our knowledge, this article demonstrates the first tool based on a convolutional Unet++ encoder-decoder architecture for the semantic segmentation of in vitro angiogenesis simulation images followed by the resulting mask postprocessing for data analysis by experts. The first annotated dataset in this field, AngioCells, is also being made publicly available. To create this dataset, participants were recruited into a markup group, an annotation protocol was developed, and an interparticipant agreement study was carried out.

Europeans and Americans of European origin show differences between their biological pathways related to the major histocompatibility complex.

In this study, we analysed biological pathway diversity among Europeans and Northern Americans of European origin, the groups of people that share a common genetic ancestry but live in different geographic regions. We used a novel complex approach for analysing genomic data: we studied the total effects of multiple weak selection signals, accumulated from independent SNPs within a pathway. We found significant differences between immunity-related biological pathways from the two groups. All identified pathways included genes belonging to the major histocompatibility complex (MHC) system, which plays an important role in adaptive immune responses. We suggest that the ways of evolution were different for the MHC-I and MHC-II gene groups at least in Europeans and Americans of European origin. We hypothesise that the observed variability between the two populations was triggered by selection pressures due to the different pathogen landscapes and pathogen loads on the two continents. Our findings can be important for epidemic prevention and control, as well as for analysing processes related to allergies, organ transplantation, and autoimmune diseases.

Biomedical Applications of Non-Small Cell Lung Cancer Spheroids.

Lung malignancies accounted for 11% of cancers worldwide in 2020 and remained the leading cause of cancer deaths. About 80% of lung cancers belong to non-small cell lung cancer (NSCLC), which is characterized by extremely high clonal and morphological heterogeneity of tumors and development of multidrug resistance. The improvement of current therapeutic strategies includes several directions. First, increasing knowledge in cancer biology results in better understanding of the mechanisms underlying malignant transformation, alterations in signal transduction, and crosstalk between cancer cells and the tumor microenvironment, including immune cells. In turn, it leads to the discovery of important molecular targets in cancer development, which might be affected pharmaceutically. The second direction focuses on the screening of novel drug candidates, synthetic or from natural sources. Finally, "personalization" of a therapeutic strategy enables maximal damage to the tumor of a patient. The personalization of treatment can be based on the drug screening performed using patient-derived tumor xenografts or in vitro patient-derived cell models. 3D multicellular cancer spheroids, generated from cancer cell lines or tumor-isolated cells, seem to be a helpful tool for the improvement of current NSCLC therapies. Spheroids are used as a tumor-mimicking in vitro model for screening of novel drugs, analysis of intercellular interactions, and oncogenic cell signaling. Moreover, several studies with tumor-derived spheroids suggest this model for the choice of "personalized" therapy. Here we aim to give an overview of the different applications of NSCLC spheroids and discuss the potential contribution of the spheroid model to the development of anticancer strategies.

Convolutional Neural Networks Can Predict Retinal Differentiation in Retinal Organoids.

We have developed a deep learning-based computer algorithm to recognize and predict retinal differentiation in stem cell-derived organoids based on bright-field imaging. The three-dimensional "organoid" approach for the differentiation of pluripotent stem cells (PSC) into retinal and other neural tissues has become a major in vitro strategy to recapitulate development. We decided to develop a universal, robust, and non-invasive method to assess retinal differentiation that would not require chemical probes or reporter gene expression. We hypothesized that basic-contrast bright-field (BF) images contain sufficient information on tissue specification, and it is possible to extract this data using convolutional neural networks (CNNs). Retina-specific Rx-green fluorescent protein mouse embryonic reporter stem cells have been used for all of the differentiation experiments in this work. The BF images of organoids have been taken on day 5 and fluorescent on day 9. To train the CNN, we utilized a transfer learning approach: ImageNet pre-trained ResNet50v2, VGG19, Xception, and DenseNet121 CNNs had been trained on labeled BF images of the organoids, divided into two categories (retina and non-retina), based on the fluorescent reporter gene expression. The best-performing classifier with ResNet50v2 architecture showed a receiver operating characteristic-area under the curve score of 0.91 on a test dataset. A comparison of the best-performing CNN with the human-based classifier showed that the CNN algorithm performs better than the expert in predicting organoid fate (84% vs. 67 ± 6% of correct predictions, respectively), confirming our original hypothesis. Overall, we have demonstrated that the computer algorithm can successfully recognize and predict retinal differentiation in organoids before the onset of reporter gene expression. This is the first demonstration of CNN's ability to classify stem cell-derived tissue in vitro.

Comparative Analysis of the Transcriptome, Proteome, and miRNA Profile of Kupffer Cells and Monocytes.

Macrophage populations in most mammalian organs consist of cells of different origin. Resident macrophages originate from erythromyeloid precursors of the yolk sac wall; maintenance of the numbers of such macrophages in postnatal ontogenesis is practically independent of bone marrow haematopoiesis. The largest populations of the resident macrophages of embryonic origin are found in the central nervous system (microglia) and liver (Kupffer cells). In contrast, skin dermis and mucous membranes become predominantly colonized by bone marrow-derived monocytes that show pronounced functional and phenotypic plasticity. In the present study, we compared Kupffer cells and monocytes using the immunophenotype, gene expression profile, proteome, and pool of microRNA. The observed differences did not consider the resident liver macrophages as purely M2 macrophages or state that monocytes have purely M1 features. Monocytes show signs of high plasticity and sensitivity to pathogen-associated molecular patterns (e.g., high levels of transcription for Tlr 2, 4, 7, and 8). In contrast, the resident liver macrophages were clearly involved in the regulation of specific organ functions (nitrogen metabolism, complement system protein synthesis).