The landscape of the COVID-19 pandemic in Poland emerging from epidemiological and genomic data.

The COVID-19 pandemic has profoundly affected all aspects of our lives. Through real-time monitoring and rapid vaccine implementation, we succeeded in suppressing the spread of the disease and mitigating its consequences. Finally, conclusions can be summarized and drawn. Here, we use the example of Poland, which was seriously affected by the pandemic. Compared to other countries, Poland has not achieved impressive results in either testing or vaccination, which may explain its high mortality (case fatality rate, CFR 1.94%). Through retrospective analysis of data collected by the COVID-19 Data Portal Poland, we found significant regional differences in the number of tests performed, number of cases detected, number of COVID-19-related deaths, and vaccination rates. The Masovian, Greater Poland, and Pomeranian voivodeships, the country's leaders in vaccination, reported high case numbers but low death rates. In contrast, the voivodeships in the eastern and southern parts of Poland (Subcarpathian, Podlaskie, Lublin, Opole), which documented low vaccination levels and low case numbers, had higher COVID-19-related mortality rates. The strong negative correlation between the CFR and the percentage of the population that was vaccinated in Poland supports the validity of vaccination. To gain insight into virus evolution, we sequenced more than 500 genomes and analyzed nearly 80 thousand SARS-CoV-2 genome sequences deposited in GISAID by Polish diagnostic centers. We showed that the SARS-CoV-2 variant distribution over time in Poland reflected that in Europe. Haplotype network analysis allowed us to follow the virus transmission routes and identify potential superspreaders in each pandemic wave.

Anti-SARS-CoV-2 Antibodies Level and COVID-19 Vaccine Boosters among Healthcare Workers with the Highest SARS-CoV-2 Infection Risk-Follow Up Study.

During the COVID-19 pandemic, several vaccines were developed to limit the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, due to SARS-CoV-2 mutations and uneven vaccination coverage among populations, a series of COVID-19 waves have been caused by different variants of concern (VOCs). Despite the updated vaccine formulations for the new VOC, the benefits of additional COVID-19 vaccine doses have raised many doubts, even among high-risk groups such as healthcare workers (HCWs). We examined the factors underlying hesitancy to receive COVID-19 booster vaccine doses and analysed the anti-SARS-CoV-2 IgG antibody response after booster vaccination among HCWs. Our study found that 42% of the HCWs were hesitant about the second booster dose, while 7% reported no intent to get vaccinated with any additional doses. As reasons for not vaccinating, participants most frequently highlighted lack of time, negative experiences with previous vaccinations, and immunity conferred by past infections. In addition, we found the lowest post-vaccination antibody titres among HCWs who did not receive any vaccine booster dose and the highest among HCWs vaccinated with two booster doses.

Gut Mycobiota Dysbiosis Is Associated with Melanoma and Response to Anti-PD-1 Therapy.

Recent research indicates that gut microbiota may be vital in the advancement of melanoma. In this study, we found that melanoma patients exhibited a distinct gut mycobiota structure compared with healthy participants. Candida albicans, Candida dubliniensis, and Neurospora crassa were more abundant in samples from patients with melanoma, whereas Saccharomyces cerevisiae and Debaryomyces hansenii were less abundant. During anti-PD-1 treatment, the relative amount of Malassezia restricta and C. albicans increased. A higher level of Saccharomyces paradoxus was associated with a positive response to anti-PD-1 treatment, whereas a higher level of Tetrapisispora blattae was associated with a lack of clinical benefits. High levels of M. restricta and C. albicans, elevated serum lactate dehydrogenase, and being overweight were linked to increased risk of melanoma progression and poorer response to anti-PD-1 treatment. Thus, this study has revealed melanoma-associated mycobiome dysbiosis, characterized by altered fungal composition and fungi species associated with a higher risk of melanoma progression, identifying a role for the gut mycobiome in melanoma progression.

Genetic history of East-Central Europe in the first millennium CE.

BACKGROUND: The appearance of Slavs in East-Central Europe has been the subject of an over 200-year debate driven by two conflicting hypotheses. The first assumes that Slavs came to the territory of contemporary Poland no earlier than the sixth century CE; the second postulates that they already inhabited this region in the Iron Age (IA). Testing either hypothesis is not trivial given that cremation of the dead was the prevailing custom in Central Europe from the late Bronze Age until the Middle Ages (MA). RESULTS: To address this problem, we determined the genetic makeup of representatives of the IA Wielbark- and MA Slav-associated cultures from the territory of present-day Poland. The study involved 474 individuals buried in 27 cemeteries. For 197 of them, genome-wide data were obtained. We found close genetic affinities between the IA Wielbark culture-associated individuals and contemporary to them and older northern European populations. Further, we observed that the IA individuals had genetic components which were indispensable to model the MA population. CONCLUSIONS: The collected data suggest that the Wielbark culture-associated IA population was formed by immigrants from the north who entered the region of contemporary Poland most likely at the beginning of the first millennium CE and mixed with autochthons. The presented results are in line with the hypothesis that assumes the genetic continuation between IA and MA periods in East-Central Europe.

Host Factors Associated with Gut Mycobiome Structure.

Recent studies revealed a significant role of the gut fungal community in human health. Here, we investigated the content and variation of gut mycobiota among subjects from the European population. We explored the interplay between gut fungi and various host-related sociodemographic, lifestyle, health, and dietary factors. The study included 923 participants. Fecal DNA samples were analyzed by whole-metagenome high-throughput sequencing. Subsequently, fungi taxonomic profiles were determined and accompanied by computational and statistical analyses of the association with 53 host-related factors. Fungal communities were characterized by a high prevalence of Saccharomyces, Candida, and Sporisorium. Ten factors were found to correlate significantly with the overall mycobiota variation. Most were diet related, including the consumption of chips, meat, sodas, sweetening, processed food, and alcohol, followed by age and marital status. Differences in α- and/or ß-diversity were also reported for other factors such as body mass index (BMI), job type, autoimmunological diseases, and probiotics. Differential abundance analysis revealed fungal species that exhibited different patterns of changes under specific conditions. The human gut mycobiota is dominated by yeast, including Saccharomyces, Malassezia, and Candida. Although intervolunteer variability was high, several fungal species persisted across most samples, which may be evidence that a core gut mycobiota exists. Moreover, we showed that host-related factors such as diet, age, and marital status influence the variability of gut mycobiota. To our knowledge, this is the first large and comprehensive study of the European cohort in terms of gut mycobiota associations with such an extensive and differentiated host-related set of factors. IMPORTANCE The human gut is inhabited by many organisms, including bacteria and fungi, that may affect human health. However, research on human gut mycobiome is still rare. Moreover, the large European-based cohort study is missing. Here, we analyzed the first large European cohort in terms of gut mycobiota associations with a differentiated host-related set of factors. Our results showed that chips, meat, sodas, sweetening, processed food, beer, alcohol consumption, age, and marital status were associated with the variability of gut mycobiota. Moreover, our analysis revealed changes in abundances at the fungal species level for many investigated factors. Our results can suggest potentially valuable paths for further, narrowly focused research on gut mycobiome and its impact on human health. In the coming era of gut microbiome-based precision medicine, further research into the relationship between different mycobial structures and host-related factors may result in new preventive approaches or therapeutic procedures.

Essential thrombocythaemia progression to the fibrotic phase is associated with a decrease in JAK2 and PDL1 levels.

It has been postulated that the changes in the molecular characteristics of the malignant clone(s) and the abnormal activation of JAK-STAT signaling are responsible for myeloproliferative neoplasm progression to more advanced disease phases and the immune escape of the malignant clone. The continuous JAK-STAT pathway activation leads to enhanced activity of the promoter of CD274 coding programmed death-1 receptor ligand (PD-L1), increased PD-L1 level, and the immune escape of MPN cells. The aim of study was to evaluate the PDL1 mRNA and JAK2 mRNA level in molecularly defined essential thrombocythaemia (ET) patients (pts) during disease progression to post-ET- myelofibrosis (post-ET-MF). The study group consisted of 162 ET pts, including 30 pts diagnosed with post-ET-MF. The JAK2V617F, CALR, and MPL mutations were found in 59.3%, 19.1%, and 1.2% of pts, respectively. No copy-number alternations of the JAK2, PDL1, and PDCDL1G2 (PDL2) genes were found. The level of PD-L1 was significantly higher in the JAK2V617F than in the JAK2WT, CALR mutation-positive, and triple-negative pts. The PD-L1 mRNA level was weakly correlated with both the JAK2V617F variant allele frequency (VAF), and with the JAK2V617F allele mRNA level. The total JAK2 level in post-ET-MF pts was lower than in ET pts, despite the lack of differences in the JAK2V617F VAF. In addition, the PD-L1 level was lower in post-ET-MF. A detailed analysis has shown that the decrease in JAK2 and PDL1 mRNA levels depended on the bone marrow fibrosis grade. The PDL1 expression showed no differences in relation to the genotype of the JAK2 haplotypeGGCC_46/1, hemoglobin concentration, hematocrit value, leukocyte, and platelet counts. The observed drop of the total JAK2 and PDL1 levels during the ET progression to the post-ET-MF may reflect the changes in the JAK2V617F positive clone proliferative potential and the PD-L1 level-related immunosuppressive effect. The above-mentioned hypothesis is supported by The Cancer Genome Atlas (TCGA) data, confirming a strong positive association between CD274 (encoding PD-L1), CXCR3 (encoding CXCR3), and CSF1 (encoding M-CSF) expression levels, and recently published results documenting a drop in the CXCR3 level and circulating M-CSF in patients with post-ET-MF.

Mutations in the miR-142 gene are not common in myeloproliferative neoplasms.

Recent data indicate that MIR142 is the most frequently mutated miRNA gene and one of the most frequently mutated noncoding elements in all cancers, with mutations occurring predominantly in blood cancers, especially diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma. Functional analyses show that the MIR142 alterations have profound consequences for lympho- and myelopoiesis. Furthermore, one of the targets downregulated by miR-142-5p is CD274, which encodes PD-L1 that is elevated in many cancer types, including myeloproliferative neoplasms (MPNs). To extend knowledge about the occurrence of MIR142 mutations, we sequenced the gene in a large panel of MPNs [~ 700 samples, including polycythemia vera, essential thrombocythemia, primary myelofibrosis (PMF), and chronic myeloid leukemia], neoplasm types in which such mutations have never been tested, and in panels of acute myeloid leukemia (AML), and chronic lymphocytic leukemia (CLL). We identified 3 mutations (one in a PMF sample and two others in one CLL sample), indicating that MIR142 mutations are rare in MPNs. In summary, mutations in MIR142 are rare in MPNs; however, in specific subtypes, such as PMF, their frequency may be comparable to that observed in CLL or AML.

The standardisation of the approach to metagenomic human gut analysis: from sample collection to microbiome profiling.

In recent years, the number of metagenomic studies increased significantly. Wide range of factors, including the tremendous community complexity and variability, is contributing to the challenge in reliable microbiome community profiling. Many approaches have been proposed to overcome these problems making hardly possible to compare results of different studies. The significant differences between procedures used in metagenomic research are reflected in a variation of the obtained results. This calls for the need for standardisation of the procedure, to reduce the confounding factors originating from DNA isolation, sequencing and bioinformatics analyses in order to ensure that the differences in microbiome composition are of a true biological origin. Although the best practices for metagenomics studies have been the topic of several publications and the main aim of the International Human Microbiome Standard (IHMS) project, standardisation of the procedure for generating and analysing metagenomic data is still far from being achieved. To highlight the difficulties in the standardisation of metagenomics methods, we thoroughly examined each step of the analysis of the human gut microbiome. We tested the DNA isolation procedure, preparation of NGS libraries for next-generation sequencing, and bioinformatics analysis, aimed at identifying microbial taxa. We showed that the homogenisation time is the leading factor impacting sample diversity, with the recommendation for a shorter homogenisation time (10 min). Ten minutes of homogenisation allows for better reflection of the bacteria gram-positive/gram-negative ratio, and the obtained results are the least heterogenous in terms of beta-diversity of samples microbial composition. Besides increasing the homogenisation time, we observed further potential impact of the library preparation kit on the gut microbiome profiling. Moreover, our analysis revealed that the choice of the library preparation kit influences the reproducibility of the results, which is an important factor that has to be taken into account in every experiment. In this study, a tagmentation-based kit allowed for obtaining the most reproducible results. We also considered the choice of the computational tool for determining the composition of intestinal microbiota, with Kraken2/Bracken pipeline outperforming MetaPhlAn2 in our in silico experiments. The design of an experiment and a detailed establishment of an experimental protocol may have a serious impact on determining the taxonomic profile of the intestinal microbiome community. Results of our experiment can be helpful for a wide range of studies that aim to better understand the role of the gut microbiome, as well as for clinical purposes.

Genetic tests based on the RT-PCR reaction in the diagnostics of SARS-CoV-2 infection.

INTRODUCTION: Since the SARS-CoV-2 emergence in 2019/2020, at least 158 million infections with this pathogen have been recorded, of which 3.29 million infected people have died. Due to the non-specific symptoms of SARS-CoV-2 infection, laboratory tests based on RT-PCR (reverse transcription and polymerase chain reaction) are mainly used in the diagnosis of COVID-19 disease. AIM: The aim of this study is to compare the molecular tests available on the Polish market for the diagnosis of SARS-CoV2 infection. RESULTS: Based on the data provided by the manufacturers and the performed laboratory analyses, we have shown that the available diagnostic kits differ mainly in the sensitivity and duration of the reaction. CONCLUSION: Due to the ongoing COVID-19 pandemic, the indicated parameters are key to effective control of the spread of SARS-CoV2, and therefore should be mainly taken into account when choosing and purchasing by diagnostic centres.

Transcript-Level Dysregulation of BCL2 Family Genes in Acute Myeloblastic Leukemia.

The expression of apoptosis-related BCL2 family genes, fine-tuned in normal cells, is dysregulated in many neoplasms. In acute myeloid leukemia (AML), this problem has not been studied comprehensively. To address this issue, RNA-seq data were used to analyze the expression of 26 BCL2 family members in 27 AML FAB M1 and M2 patients, divided into subgroups differently responding to chemotherapy. A correlation analysis, analysis of variance, and Kaplan-Meier analysis were applied to associate the expression of particular genes with other gene expression, clinical features, and the presence of mutations detected by exome sequencing. The expression of BCL2 family genes was dysregulated in AML, as compared to healthy controls. An upregulation of anti-apoptotic and downregulation of pro-apoptotic genes was observed, though only a decrease in BMF, BNIP1, and HRK was statistically significant. In a group of patients resistant to chemotherapy, overexpression of BCL2L1 was manifested. In agreement with the literature data, our results reveal that BCL2L1 is one of the key players in apoptosis regulation in different types of tumors. An exome sequencing data analysis indicates that BCL2 family genes are not mutated in AML, but their expression is correlated with the mutational status of other genes, including those recurrently mutated in AML and splicing-related. High levels of some BCL2 family members, in particular BIK and BCL2L13, were associated with poor outcome.

Proteomic Profiling of Leukocytes Reveals Dysregulation of Adhesion and Integrin Proteins in Chronic Kidney Disease-Related Atherosclerosis.

A progressive loss of functional nephrons defines chronic kidney disease (CKD). Complications related to cardiovascular disease (CVD) are the principal causes of mortality in CKD; however, the acceleration of CVD in CKD remains unresolved. Our study used a complementary proteomic approach to assess mild and advanced CKD patients with different atherosclerosis stages and two groups of patients with different classical CVD progression but without renal dysfunction. We utilized a label-free approach based on LC-MS/MS and functional bioinformatic analyses to profile CKD and CVD leukocyte proteins. We revealed dysregulation of proteins involved in different phases of leukocytes' diapedesis process that is very pronounced in CKD's advanced stage. We also showed an upregulation of apoptosis-related proteins in CKD as compared to CVD. The differential abundance of selected proteins was validated by multiple reaction monitoring, ELISA, Western blotting, and at the mRNA level by ddPCR. An increased rate of apoptosis was then functionally confirmed on the cellular level. Hence, we suggest that the disturbances in leukocyte extravasation proteins may alter cell integrity and trigger cell death, as demonstrated by flow cytometry and microscopy analyses. Our proteomics data set has been deposited to the ProteomeXchange Consortium via the PRIDE repository with the data set identifier PXD018596.

Transcriptome Remodeling in Gradual Development of Inverse Resistance between Paclitaxel and Cisplatin in Ovarian Cancer Cells.

Resistance to anti-cancer drugs is the main challenge in oncology. In pre-clinical studies, established cancer cell lines are primary tools in deciphering molecular mechanisms of this phenomenon. In this study, we proposed a new, transcriptome-focused approach, utilizing a model of isogenic cancer cell lines with gradually changing resistance. We analyzed trends in gene expression in the aim to find out a scaffold of resistance development process. The ovarian cancer cell line A2780 was treated with stepwise increased concentrations of paclitaxel (PTX) to generate a series of drug resistant sublines. To monitor transcriptome changes we submitted them to mRNA-sequencing, followed by the identification of differentially expressed genes (DEGs), principal component analysis (PCA), and hierarchical clustering. Functional interactions of proteins, encoded by DEGs, were analyzed by building protein-protein interaction (PPI) networks. We obtained human ovarian cancer cell lines with gradually developed resistance to PTX and collateral sensitivity to cisplatin (CDDP) (inverse resistance). In their transcriptomes, we identified two groups of DEGs: (1) With fluctuations in expression in the course of resistance acquiring; and (2) with a consistently changed expression at each stage of resistance development, constituting a scaffold of the process. In the scaffold PPI network, the cell cycle regulator-polo-like kinase 2 (PLK2); proteins belonging to the tumor necrosis factor (TNF) ligand and receptor family, as well as to the ephrin receptor family were found, and moreover, proteins linked to osteo- and chondrogenesis and the nervous system development. Our cellular model of drug resistance allowed for keeping track of trends in gene expression and studying this phenomenon as a process of evolution, reflected by global transcriptome remodeling. This approach enabled us to explore novel candidate genes and surmise that abrogation of the osteomimic phenotype in ovarian cancer cells might occur during the development of inverse resistance between PTX and CDDP.

Expression Landscape of circRNAs in Arabidopsis thaliana Seedlings and Adult Tissues.

RNA-seq is currently the only method that can provide a comprehensive landscape of circular RNA (circRNAs) in the whole organism and its particular organs. Recent years have brought an increasing number of RNA-seq-based reports on plant circRNAs. Notably, the picture they revealed is questionable and depends on the applied circRNA identification and quantification techniques. In consequence, little is known about the biogenesis and functions of circRNAs in plants. In this work, we tested two experimental and six bioinformatics procedures of circRNA analysis to determine the optimal approach for studying the profiles of circRNAs in Arabidopsis thaliana. Then using the optimized strategy, we determined the accumulation of circular and corresponding linear transcripts in plant seedlings and organs. We observed that only a small fraction of circRNAs was reproducibly generated. Among them, two groups of circRNAs were discovered: ubiquitous and organ-specific. The highest number of circRNAs with significantly increased accumulation in comparison to other organs/seedlings was found in roots. The circRNAs in seedlings, leaves and flowers originated mainly from genes involved in photosynthesis and the response to stimulus. The levels of circular and linear transcripts were not correlated. Although RNase R treatment enriches the analyzed RNA samples in circular transcripts, it may also have a negative impact on the stability of some of the circRNAs. We also showed that the normalization of NGS data by the library size is not proper for circRNAs quantification. Alternatively, we proposed four other normalization types whose accuracy was confirmed by ddPCR. Moreover, we provided a comprehensive characterization of circRNAs in A. thaliana organs and in seedlings. Our analyses revealed that plant circRNAs are formed in both stochastic and controlled processes. The latter are less frequent and likely engage circRNA-specific mechanisms. Only a few circRNAs were organ-specific. The lack of correlation between the accumulation of linear and circular transcripts indicated that their biogenesis depends on different mechanisms.

Arabidopsis thaliana cbp80, c2h2, and flk Knockout Mutants Accumulate Increased Amounts of Circular RNAs.

Circular RNAs (circRNAs) are the products of the non-canonical splicing of pre-mRNAs. In contrast to humans and animals, our knowledge of the biogenesis and function of circRNAs in plants is very scarce. To identify proteins involved in plant circRNA generation, we characterized the transcriptomes of 18 Arabidopsis thaliana knockout mutants for genes related to splicing. The vast majority (>90%) of circRNAs were formed in more than one variant; only a small fraction of circRNAs was mutant-specific. Five times more circRNA types were identified in cbp80 and three times more in c2h2 mutants than in the wild-type. We also discovered that in cbp80, c2h2 and flk mutants, the accumulation of circRNAs was significantly increased. The increased accumulation of circular transcripts was not accompanied by corresponding changes in the accumulation of linear transcripts. Our results indicate that one of the roles of CBP80, C2H2 and FLK in splicing is to ensure the proper order of the exons. In the absence of one of the above-mentioned factors, the process might be altered, leading to the production of circular transcripts. This suggests that the transition toward circRNA production can be triggered by factors sequestering these proteins. Consequently, the expression of linear transcripts might be regulated through circRNA production.

Analysis of oral microbiome from fossil human remains revealed the significant differences in virulence factors of modern and ancient Tannerella forsythia.

BACKGROUND: Recent advances in the next-generation sequencing (NGS) allowed the metagenomic analyses of DNA from many different environments and sources, including thousands of years old skeletal remains. It has been shown that most of the DNA extracted from ancient samples is microbial. There are several reports demonstrating that the considerable fraction of extracted DNA belonged to the bacteria accompanying the studied individuals before their death. RESULTS: In this study we scanned 344 microbiomes from 1000- and 2000- year-old human teeth. The datasets originated from our previous studies on human ancient DNA (aDNA) and on microbial DNA accompanying human remains. We previously noticed that in many samples infection-related species have been identified, among them Tannerella forsythia, one of the most prevalent oral human pathogens. Samples containing sufficient amount of T. forsythia aDNA for a complete genome assembly were selected for thorough analyses. We confirmed that the T. forsythia-containing samples have higher amounts of the periodontitis-associated species than the control samples. Despites, other pathogens-derived aDNA was found in the tested samples it was too fragmented and damaged to allow any reasonable reconstruction of these bacteria genomes. The anthropological examination of ancient skulls from which the T. forsythia-containing samples were obtained revealed the pathogenic alveolar bone loss in tooth areas characteristic for advanced periodontitis. Finally, we analyzed the genetic material of ancient T. forsythia strains. As a result, we assembled four ancient T. forsythia genomes - one 2000- and three 1000- year-old. Their comparison with contemporary T. forsythia genomes revealed a lower genetic diversity within the four ancient strains than within contemporary strains. We also investigated the genes of T. forsythia virulence factors and found that several of them (KLIKK protease and bspA genes) differ significantly between ancient and modern bacteria. CONCLUSIONS: In summary, we showed that NGS screening of the ancient human microbiome is a valid approach for the identification of disease-associated microbes. Following this protocol, we provided a new set of information on the emergence, evolution and virulence factors of T. forsythia, the member of the oral dysbiotic microbiome.

Characterization of RNP Networks of PUM1 and PUM2 Post-Transcriptional Regulators in TCam-2 Cells, a Human Male Germ Cell Model.

Mammalian Pumilio (PUM) proteins are sequence-specific, RNA-binding proteins (RBPs) with wide-ranging roles. They are involved in germ cell development, which has functional implications in development and fertility. Although human PUM1 and PUM2 are closely related to each other and recognize the same RNA binding motif, there is some evidence for functional diversity. To address that problem, first we used RIP-Seq and RNA-Seq approaches, and identified mRNA pools regulated by PUM1 and PUM2 proteins in the TCam-2 cell line, a human male germ cell model. Second, applying global mass spectrometry-based profiling, we identified distinct PUM1- and PUM2-interacting putative protein cofactors, most of them involved in RNA processing. Third, combinatorial analysis of RIP and RNA-Seq, mass spectrometry, and RNA motif enrichment analysis revealed that PUM1 and PUM2 form partially varied RNP-regulatory networks (RNA regulons), which indicate different roles in human reproduction and testicular tumorigenesis. Altogether, this work proposes that protein paralogues with very similar and evolutionary highly conserved functional domains may play divergent roles in the cell by combining with different sets of protein cofactors. Our findings highlight the versatility of PUM paralogue-based post-transcriptional regulation, offering insight into the mechanisms underlying their diverse biological roles and diseases resulting from their dysfunction.

Not Only Mutations Matter: Molecular Picture of Acute Myeloid Leukemia Emerging from Transcriptome Studies.

The last two decades of genome-scale research revealed a complex molecular picture of acute myeloid leukemia (AML). On the one hand, a number of mutations were discovered and associated with AML diagnosis and prognosis; some of them were introduced into diagnostic tests. On the other hand, transcriptome studies, which preceded AML exome and genome sequencing, remained poorly translated into clinics. Nevertheless, gene expression studies significantly contributed to the elucidation of AML pathogenesis and indicated potential therapeutic directions. The power of transcriptomic approach lies in its comprehensiveness; we can observe how genome manifests its function in a particular type of cells and follow many genes in one test. Moreover, gene expression measurement can be combined with mutation detection, as high-impact mutations are often present in transcripts. This review sums up 20 years of transcriptome research devoted to AML. Gene expression profiling (GEP) revealed signatures distinctive for selected AML subtypes and uncovered the additional within-subtype heterogeneity. The results were particularly valuable in the case of AML with normal karyotype which concerns up to 50% of AML cases. With the use of GEP, new classes of the disease were identified and prognostic predictors were proposed. A plenty of genes were detected as overexpressed in AML when compared to healthy control, including KIT, BAALC, ERG, MN1, CDX2, WT1, PRAME, and HOX genes. High expression of these genes constitutes usually an unfavorable prognostic factor. Upregulation of FLT3 and NPM1 genes, independent on their mutation status, was also reported in AML and correlated with poor outcome. However, transcriptome is not limited to the protein-coding genes; other types of RNA molecules exist in a cell and regulate genome function. It was shown that microRNA (miRNA) profiles differentiated AML groups and predicted outcome not worse than protein-coding gene profiles. For example, upregulation of miR-10a, miR-10b, and miR-196b and downregulation of miR-192 were found as typical of AML with NPM1 mutation whereas overexpression of miR-155 was associated with FLT3-internal tandem duplication (FLT3-ITD). Development of high-throughput technologies and microarray replacement by next generation sequencing (RNA-seq) enabled uncovering a real variety of leukemic cell transcriptomes, reflected by gene fusions, chimeric RNAs, alternatively spliced transcripts, miRNAs, piRNAs, long noncoding RNAs (lncRNAs), and their special type, circular RNAs. Many of them can be considered as AML biomarkers and potential therapeutic targets. The relations between particular RNA puzzles and other components of leukemic cells and their microenvironment, such as exosomes, are now under investigation. Hopefully, the results of this research will shed the light on these aspects of AML pathogenesis which are still not completely understood.

Altered Levels of Proteins and Phosphoproteins, in the Absence of Early Causative Transcriptional Changes, Shape the Molecular Pathogenesis in the Brain of Young Presymptomatic Ki91 SCA3/MJD Mouse.

Spinocerebellar ataxia type 3 (SCA3/MJD) is a polyQ neurodegenerative disease where the presymptomatic phase of pathogenesis is unknown. Therefore, we investigated the molecular network of transcriptomic and proteomic triggers in young presymptomatic SCA3/MJD brain from Ki91 knock-in mouse. We found that transcriptional dysregulations resulting from mutant ataxin-3 are not occurring in young Ki91 mice, while old Ki91 mice and also postmitotic patient SCA3 neurons demonstrate the late transcriptomic changes. Unlike the lack of early mRNA changes, we have identified numerous early changes of total proteins and phosphoproteins in 2-month-old Ki91 mouse cortex and cerebellum. We discovered the network of processes in presymptomatic SCA3 with three main groups of disturbed processes comprising altered proteins: (I) modulation of protein levels and DNA damage (Pabpc1, Ddb1, Nedd8), (II) formation of neuronal cellular structures (Tubb3, Nefh, p-Tau), and (III) neuronal function affected by processes following perturbed cytoskeletal formation (Mt-Co3, Stx1b, p-Syn1). Phosphoproteins downregulate in the young Ki91 mouse brain and their phosphosites are associated with kinases that interact with ATXN3 such as casein kinase, Camk2, and kinases controlled by another Atxn3 interactor p21 such as Gsk3, Pka, and Cdk kinases. We conclude that the onset of SCA3 pathology occurs without altered transcript level and is characterized by changed levels of proteins responsible for termination of translation, DNA damage, spliceosome, and protein phosphorylation. This disturbs global cellular processes such as cytoskeleton and transport of vesicles and mitochondria along axons causing energy deficit and neurodegeneration also manifesting in an altered level of transcripts at later ages.

NPM1 alternative transcripts are upregulated in acute myeloid and lymphoblastic leukemia and their expression level affects patient outcome.

BACKGROUND: Expression of the NPM1 gene, encoding nucleophosmin, is upregulated in cancers. Although more than ten NPM1 transcripts are known, the reports were usually limited to one predominant transcript. In leukemia, the NPM1 expression has not been widely studied so far. In acute myeloid leukemia (AML), the mutational status of the gene seems to play a pivotal role in carcinogenesis. Therefore, the aim of the study was to quantify alternative NPM1 transcripts in two types of acute leukemia, AML and ALL (acute lymphoblastic leukemia). METHODS: Using droplet digital PCR, we analyzed the levels of three protein-coding NPM1 transcripts in 66 samples collected from AML and ALL patients and 16 control samples. Using RNA-seq, we detected 8 additional NPM1 transcripts, including non-coding splice variants with retained introns. For data analysis, Welch two sample t-test, Pearson's correlation and Kaplan-Meier analysis were applied. RESULTS: The levels of the particular NPM1 transcripts were significantly different but highly correlated with each other in both leukemia and control samples. Transcript NPM1.1, encoding the longest protein (294 aa), had the highest level of accumulation and was one of the most abundant transcripts in the cell. Comparing to NPM1.1, the levels of the NPM1.2 and NPM1.3 transcripts, encoding a 265-aa and 259-aa proteins, were 30 and 3 times lower, respectively. All three NPM1 transcripts were proportionally upregulated in both types of leukemia compared to control samples. In AML, the levels of NPM1 transcripts decreased in complete remission and increased again with relapse of the disease. Low levels of NPM1.1 and NPM1.3 were associated with better prognosis. The contribution of non-coding transcripts to the total level of NPM1 gene seemed to be marginal, except for one short 5-end transcript accumulated at high levels in AML and control cells. Aberrant proportions of particular NPM1 splice variants could be linked to abnormal expression of genes encoding alternative splicing factors. CONCLUSIONS: The levels of the studied NPM1 transcripts were different but highly correlated with each other. Their upregulation in AML and ALL, decrease after therapy and association with patient outcome suggests the involvement of elevated NPM1 expression in the acute leukemia pathogenesis.

A mosaic genetic structure of the human population living in the South Baltic region during the Iron Age.

Despite the increase in our knowledge about the factors that shaped the genetic structure of the human population in Europe, the demographic processes that occurred during and after the Early Bronze Age (EBA) in Central-East Europe remain unclear. To fill the gap, we isolated and sequenced DNAs of 60 individuals from Kowalewko, a bi-ritual cemetery of the Iron Age (IA) Wielbark culture, located between the Oder and Vistula rivers (Kow-OVIA population). The collected data revealed high genetic diversity of Kow-OVIA, suggesting that it was not a small isolated population. Analyses of mtDNA haplogroup frequencies and genetic distances performed for Kow-OVIA and other ancient European populations showed that Kow-OVIA was most closely linked to the Jutland Iron Age (JIA) population. However, the relationship of both populations to the preceding Late Neolithic (LN) and EBA populations were different. We found that this phenomenon is most likely the consequence of the distinct genetic history observed for Kow-OVIA women and men. Females were related to the Early-Middle Neolithic farmers, whereas males were related to JIA and LN Bell Beakers. In general, our findings disclose the mechanisms that could underlie the formation of the local genetic substructures in the South Baltic region during the IA.