Short 2'-O-methyl/LNA oligomers as highly-selective inhibitors of miRNA production in vitro and in vivo

MicroRNAs (miRNAs) that share identical or near-identical sequences constitute miRNA families and are predicted to act redundantly. Yet recent evidence suggests that members of the same miRNA family with high sequence similarity might have different roles and that this functional divergence might be rooted in their precursors' sequence. Current knock-down strategies such as antisense oligonucleotides (ASOs) or miRNA sponges cannot distinguish between identical or near identical miRNAs originating from different precursors to allow exploring unique functions of these miRNAs. We here develop a novel strategy based on short 2'-OMe/LNA-modified oligonucleotides to selectively target specific precursor molecules and ablate the production of individual members of miRNA families in vitro and in vivo. Leveraging the highly conserved Xenopus miR-181a family as proof-of-concept, we demonstrate that 2'-OMe/LNA-ASOs targeting the apical region of pre-miRNAs achieve precursor-selective inhibition of mature miRNA-5p production. Furthermore, we extend the applicability of our approach to the human miR-16 family, illustrating its universality in targeting precursors generating identical miRNAs. Overall, our strategy enables efficient manipulation of miRNA expression, offering a powerful tool to dissect the functions of identical or highly similar miRNAs derived from different precursors within miRNA families.

Secreted novel AID/APOBEC-like deaminase 1 (SNAD1) - a new important player in fish immunology

The AID/APOBECs are a group of zinc-dependent cytidine deaminases that catalyse the deamination of bases in nucleic acids, resulting in a cytidine to uridine transition. Secreted novel AID/APOBEC-like deaminases (SNADs), characterized by the presence of a signal peptide are unique among all of intracellular classical AID/APOBECs, which are the central part of antibody diversity and antiviral defense. To date, there is no available knowledge on SNADs including protein characterization, biochemical characteristics and catalytic activity. We used various in silico approaches to define the phylogeny of SNADs, their common structural features, and their potential structural variations in fish species. Our analysis provides strong evidence of the universal presence of SNAD1 proteins/transcripts in fish, in which expression commences after hatching and is highest in anatomical organs linked to the immune system. Moreover, we searched published fish data and identified previously, "uncharacterized proteins" and transcripts as SNAD1 sequences. Our review into immunological research suggests SNAD1 role in immune response to infection or immunization, and interactions with the intestinal microbiota. We also noted SNAD1 association with temperature acclimation, environmental pollution and sex-based expression differences, with females showing higher level. To validate in silico predictions we performed expression studies of several SNAD1 gene variants in carp, which revealed distinct patterns of responses under different conditions. Dual sensitivity to environmental and pathogenic stress highlights its importance in the fish and potentially enhancing thermotolerance and immune defense. Revealing the biological roles of SNADs represents an exciting new area of research related to the role of DNA and/or RNA editing in fish biology.

Engineered deaminases as a key component of DNA and RNA editing tools

Over recent years, zinc-dependent deaminases have attracted increasing interest as key components of nucleic acid editing tools that can generate point mutations at specific sites in either DNA or RNA by combining a targeting module (such as a catalytically impaired CRISPR-Cas component) and an effector module (most often a deaminase). Deaminase-based molecular tools are already being utilized in a wide spectrum of therapeutic and research applications; however, their medical and biotechnological potential seems to be much greater. Recent reports indicate that the further development of nucleic acid editing systems depends largely on our ability to engineer the substrate specificity and catalytic activity of the editors themselves. In this review, we summarize the current trends and achievements in deaminase engineering. The presented data indicate that the potential of these enzymes has not yet been fully revealed or understood. Several examples show that even relatively minor changes in the structure of deaminases can give them completely new and unique properties.

Cancer Stem Cell Markers-Clinical Relevance and Prognostic Value in High-Grade Serous Ovarian Cancer (HGSOC) Based on The Cancer Genome Atlas Analysis

Cancer stem cells (CSCs) may contribute to an increased risk of recurrence in ovarian cancer (OC). Further research is needed to identify associations between CSC markers and OC patients' clinical outcomes with greater certainty. If they prove to be correct, in the future, the CSC markers can be used to help predict survival and indicate new therapeutic targets. This study aimed to determine the CSC markers at mRNA and protein levels and their association with clinical presentation, outcome, and risk of recurrence in HGSOC (High-Grade Serous Ovarian Cancer). TCGA (The Cancer Genome Atlas) database with 558 ovarian cancer tumor samples was used for the evaluation of 13 CSC markers (ALDH1A1, CD44, EPCAM, KIT, LGR5, NES, NOTCH3, POU5F1, PROM1, PTTG1, ROR1, SOX9, and THY1). Data on mRNA and protein levels assessed by microarray and mass spectrometry were retrieved from TCGA. Models to predict chemotherapy response and survival were built using multiple variables, including epidemiological data, expression levels, and machine learning methodology. ALDH1A1 and LGR5 mRNA expressions indicated a higher platinum sensitivity (p = 3.50 × 10-3; p = 0.01, respectively). POU5F1 mRNA expression marked platinum-resistant tumors (p = 9.43 × 10-3). CD44 and EPCAM mRNA expression correlated with longer overall survival (OS) (p = 0.043; p = 0.039, respectively). THY1 mRNA and protein levels were associated with worse OS (p = 0.019; p = 0.015, respectively). Disease-free survival (DFS) was positively affected by EPCAM (p = 0.004), LGR5 (p = 0.018), and CD44 (p = 0.012). In the multivariate model based on CSC marker expression, the high-risk group had 9.1 months longer median overall survival than the low-risk group (p < 0.001). ALDH1A1, CD44, EPCAM, LGR5, POU5F1, and THY1 levels in OC may be used as prognostic factors for the primary outcome and help predict the treatment response.

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.

Human AGEs: an interactive spatio-temporal visualization and database of human archeogenomics

Archeogenomics is a rapidly growing interdisciplinary research field driven by the development of techniques that enable the acquisition and analysis of ancient DNA (aDNA). Recent advances in aDNA studies have contributed significantly to increasing our understanding of the natural history of humans. One of the most significant challenges facing archeogenomics is the integration of highly heterogeneous genomic, archeological, and anthropological data and their comprehensive analysis, considering changes that occur in time and space. Only this complex approach can explain the relationship between past populations in the context of migration or cultural development. To address these challenges, we developed a Human AGEs web server. It focuses on creating comprehensive spatiotemporal visualizations of genomic, archeogenomic, and archeological information, which can be provided by the user or loaded from a graph database. The interactive map application at the center of Human AGEs can display multiple layers of data in various forms, such as bubble charts, pie charts, heatmaps, or tag clouds. These visualizations can be modified using various clustering, filtering, and styling options, and the map state can be exported to a high-resolution image or saved as a session file for later use. Human AGEs, along with their tutorial, are accessible at https://archeogenomics.eu/.

Differences in the intraspecies copy number variation of Arabidopsis thaliana conserved and nonconserved miRNA genes

MicroRNAs (miRNAs) regulate gene expression by RNA interference mechanism. In plants, miRNA genes (MIRs) which are grouped into conserved families, i.e. they are present among the different plant taxa, are involved in the regulation of many developmental and physiological processes. The roles of the nonconserved MIRs-which are MIRs restricted to one plant family, genus, or even species-are less recognized; however, many of them participate in the responses to biotic and abiotic stresses. Both over- and underproduction of miRNAs may influence various biological processes. Consequently, maintaining intracellular miRNA homeostasis seems to be crucial for the organism. Deletions and duplications in the genomic sequence may alter gene dosage and/or activity. We evaluated the extent of copy number variations (CNVs) among Arabidopsis thaliana (Arabidopsis) MIRs in over 1000 natural accessions, using population-based analysis of the short-read sequencing data. We showed that the conserved MIRs were unlikely to display CNVs and their deletions were extremely rare, whereas nonconserved MIRs presented moderate variation. Transposon-derived MIRs displayed exceptionally high diversity. Conversely, MIRs involved in the epigenetic control of transposons reactivated during development were mostly invariable. MIR overlap with the protein-coding genes also limited their variability. At the expression level, a higher rate of nonvariable, nonconserved miRNAs was detectable in Col-0 leaves, inflorescence, and siliques compared to nonconserved variable miRNAs, although the expression of both groups was much lower than that of the conserved MIRs. Our data indicate that CNV rate of Arabidopsis MIRs is related with their age, function, and genomic localization.

Analysis of Arabidopsis non-reference accessions reveals high diversity of metabolic gene clusters and discovers new candidate cluster members

Metabolic gene clusters (MGCs) are groups of genes involved in a common biosynthetic pathway. They are frequently formed in dynamic chromosomal regions, which may lead to intraspecies variation and cause phenotypic diversity. We examined copy number variations (CNVs) in four Arabidopsis thaliana MGCs in over one thousand accessions with experimental and bioinformatic approaches. Tirucalladienol and marneral gene clusters showed little variation, and the latter was fixed in the population. Thalianol and especially arabidiol/baruol gene clusters displayed substantial diversity. The compact version of the thalianol gene cluster was predominant and more conserved than the noncontiguous version. In the arabidiol/baruol cluster, we found a large genomic insertion containing divergent duplicates of the CYP705A2 and BARS1 genes. The BARS1 paralog, which we named BARS2, encoded a novel oxidosqualene synthase. The expression of the entire arabidiol/baruol gene cluster was altered in the accessions with the duplication. Moreover, they presented different root growth dynamics and were associated with warmer climates compared to the reference-like accessions. In the entire genome, paired genes encoding terpene synthases and cytochrome P450 oxidases were more variable than their nonpaired counterparts. Our study highlights the role of dynamically evolving MGCs in plant adaptation and phenotypic diversity.

The emerging role of the gut mycobiome in liver diseases

In recent years, it has become clear that gut microbiota plays a major role in the human body, both in health and disease. Because of that, the gut microbiome and its impact on human well-being are getting wider and wider attention. Studies focused on the liver are not an exception. However, the majority of the analyses are concentrated on the bacterial part of the gut microbiota, while the fungi living in the human intestines are often omitted or underappreciated. This review is focused on the gut mycobiome as an important factor that should be taken into consideration regarding liver homeostasis and its perturbations. We have collected the findings in this field and we discuss their importance. We aim to emphasize the fungal compositional changes related to liver diseases and, by that, provide novel insights into the directions of liver research and gut microbiota as a therapeutic target for liver diseases.

Dimensionality reduction by UMAP for visualizing and aiding in classification of imaging flow cytometry data

Recent advances in imaging flow cytometry (IFC) have revolutionized high-throughput multiparameter analyses at single-cell resolution. Although enabling the discovery of population heterogeneities and the detection of rare events, IFC generates hyperdimensional datasets that demand innovative analytical approaches. Current methods work in a supervised manner, utilize only limited information content, or require large annotated reference datasets. Dimensionality reduction algorithms, including uniform manifold approximation and projection (UMAP), have been successfully applied to analyze the large number of parameters generated in various high-throughput techniques.

Here, we apply a workflow incorporating UMAP to analyze different IFC datasets. We demonstrate that it out-competes other popular dimensionality reduction methods in speed and accuracy. Moreover, it enables fast visualization, clustering, and tagging of unannotated objects in large-scale experiments. We anticipate that our workflow will be a robust method to address complex IFC datasets, either alone or as an upstream addition to the deep learning approaches.

•

UMAP dimensionality reduction provides fast and accurate method of IFC data analysis

•

UMAP yields improved object clustering and tagging of the multispectral IFC data

•

PCA decomposition allows multispectral signals merging for direct image embedding

From Alpha to Delta-Genetic Epidemiology of SARS-CoV-2 (hCoV-19) in Southern Poland

In Poland, the first case of SARS-CoV-2 infection was confirmed in March 2020. Since then, many circulating virus lineages fueled rapid pandemic waves which inflicted a severe burden on the Polish healthcare system. Some of these lineages were associated with increased transmissibility and immune escape. Mutations in the viral spike protein, which is responsible for host cell recognition and serves as the primary target for neutralizing antibodies, are of particular importance. We investigated the molecular epidemiology of the SARS-CoV-2 clades circulating in Southern Poland from February 2021 to August 2021. The 921 whole-genome sequences were used for variant identification, spike mutation, and phylogenetic analyses. The Pango B.1.1.7 was the dominant variant (n = 730, 89.68%) from March 2021 to July 2021. In July 2021, the B.1.1.7 was displaced by the B.1.617.2 lineage with 66.66% in July 2021 and 92.3% in August 2021 frequencies, respectively. Moreover, our results were compared with the sequencing available on the GISAID platform for other regions of Poland, the Czech Republic, and Slovakia. The analysis showed that the dominant variant in the analyzed period was B.1.1.7 in all countries and Southern Poland (Silesia). Interestingly, B.1.1.7 was replaced by B.1.617.2 earlier in Southern Poland than in the rest of the country. Moreover, in the Czech Republic and Slovakia, AY lineages were predominant at that time, contrary to the Silesia region.

Arrested in Glass: Actin within Sophisticated Architectures of Biosilica in Sponges

Actin is a fundamental member of an ancient superfamily of structural intracellular proteins and plays a crucial role in cytoskeleton dynamics, ciliogenesis, phagocytosis, and force generation in both prokaryotes and eukaryotes. It is shown that actin has another function in metazoans: patterning biosilica deposition, a role that has spanned over 500 million years. Species of glass sponges (Hexactinellida) and demosponges (Demospongiae), representatives of the first metazoans, with a broad diversity of skeletal structures with hierarchical architecture unchanged since the late Precambrian, are studied. By etching their skeletons, organic templates dominated by individual F-actin filaments, including branched fibers and the longest, thickest actin fiber bundles ever reported, are isolated. It is proposed that these actin-rich filaments are not the primary site of biosilicification, but this highly sophisticated and multi-scale form of biomineralization in metazoans is ptterned.

Cryo-EM reconstructions of BMV-derived virus-like particles reveal assembly defects in the icosahedral lattice structure

The increasing interest in virus-like particles (VLPs) has been reflected by the growing number of studies on their assembly and application. However, the formation of complete VLPs is a complex phenomenon, making it difficult to rationally design VLPs with desired features de novo. In this paper, we describe VLPs assembled in vitro from the recombinant capsid protein of brome mosaic virus (BMV). The analysis of VLPs was performed by Cryo-EM reconstructions and allowed us to visualize a few classes of VLPs, giving insight into the VLP self-assembly process. Apart from the mature icosahedral VLP practically identical with native virions, we describe putative VLP intermediates displaying non-icosahedral arrangements of capsomers, proposed to occur before the final disorder-order transition stage of icosahedral VLP assembly. Some of the described VLP classes show a lack of protein shell continuity, possibly resulting from too strong interaction with the cargo (in this case tRNA) with the capsid protein. We believe that our results are a useful prerequisite for the rational design of VLPs in the future and lead the way to the effective production of modified VLPs.

At-C-RNA database, a one-stop source for information on circRNAs in Arabidopsis thaliana in a unified format

Circular RNAs (circRNAs) are a large class of noncoding RNAs with functions that, in most cases, remain unknown. Recent genome-wide analysis of circRNAs using RNA-Seq has revealed that circRNAs are abundant and some of them conserved in plants. Furthermore, it has been shown that the expression of circRNAs in plants is regulated in a tissue-specific manner. Arabidopsis thaliana circular RNA database is a new resource designed to integrate and standardize the data available for circRNAs in a model plant A. thaliana, which is currently the best-characterized plant in terms of circRNAs. The resource integrates all applicable publicly available RNA-seq datasets. These datasets were subjected to extensive reanalysis and curation, yielding results in a unified format. Moreover, all data were normalized according to our optimized approach developed for circRNA identification in plants. As a result, the database accommodates circRNAs identified across organs and seedlings of wild-type A. thaliana and its single-gene knockout mutants for genes related to splicing. The database provides free access to unified data and search functionalities, thus enabling comparative analyses of A. thaliana circRNAs between organs, variants and studies for the first time.

Database URLhttps://plantcircrna.ibch.poznan.pl/

Polymorphism rs368234815 of interferon lambda 4 gene and spontaneous clearance of hepatitis C virus in haemodialysis patients: a case-control study

BACKGROUND

In non-uremic subjects, IFNL4 rs368234815 predicts HCV clearance. We investigated whether rs368234815 is associated with spontaneous HCV clearance in haemodialysis patients and whether it is a stronger predictor of HCV resolution than the IFNL polymorphisms already associated with HCV clearance in dialysis subjects. We also evaluated an association of rs368234815 with patients` survival and alterations in transcription factor binding sites (TFBS) caused by IFNL polymorphisms.

METHODS

Among 161 haemodialysis patients with positive anti-HCV antibodies, 68 (42.2%) spontaneously resolved HCV infection, whereas 93 remained HCV RNA positive. Patients were tested for near IFNL3 rs12980275, IFNL3 rs4803217, IFNL4 rs12979860, IFNL4 rs368234815, and near IFNL4 rs8099917. IFNL4 rs368234815 polymorphism (TT/TT, ΔG/TT, ΔG/ΔG) was genotyped by restriction fragment length polymorphism analysis; other IFNL polymorphisms - by high resolution melting curve analysis. We used the Kaplan-Meier method with the log-rank test for survival analysis. In silico analysis included the use of ENCODE TFBS ChIP-seq data, HOCOMOCO, JASPAR CORE, and CIS-BP databases, and FIMO software.

RESULTS

The probability (OR, 95%CI, P) of spontaneous HCV clearance for rs368234815 TT/TT patients was higher than for the ΔG allele carriers (2.63, 1.38-5.04, 0.003). This probability for other major homozygotes varied between 2.80, 1.45-5.43, 0.002 for rs12980275 and 2.44, 1.27-4.69, 0.007 for rs12979860. In the additive model, rs368234815 TT/TT was the strongest predictor of HCV clearance (6.38, 1.69-24.2, 0.003). Survival analysis suggested an association of the ΔG allele with mortality due to neoplasms (log-rank P = 0.005). The rs368234815 ∆G allele caused TFBS removal for PLAGL1.

CONCLUSIONS

In haemodialysis patients, the association of rs368234815 with the spontaneous HCV clearance is better than that documented for other IFNL3/IFNL4 polymorphisms only in the additive mode of inheritance. However, identifying the homozygosity in the variant ∆G allele of rs368234815 means a more potent prediction of persistent HCV infection in haemodialysis subjects that we observe in the case of the variant homozygosity of other tested IFNL3/IFNL4 polymorphisms. Removal of PLAGL1 TFBS in subjects harbouring the rs368234815 ∆G allele may contribute to cancer susceptibility. The association of rs368234815 with cancer-related mortality needs further studies in HCV-exposed subjects.

LifeTime and improving European healthcare through cell-based interceptive medicine

Here we describe the LifeTime Initiative, which aims to track, understand and target human cells during the onset and progression of complex diseases, and to analyse their response to therapy at single-cell resolution. This mission will be implemented through the development, integration and application of single-cell multi-omics and imaging, artificial intelligence and patient-derived experimental disease models during the progression from health to disease. The analysis of large molecular and clinical datasets will identify molecular mechanisms, create predictive computational models of disease progression, and reveal new drug targets and therapies. The timely detection and interception of disease embedded in an ethical and patient-centred vision will be achieved through interactions across academia, hospitals, patient associations, health data management systems and industry. The application of this strategy to key medical challenges in cancer, neurological and neuropsychiatric disorders, and infectious, chronic inflammatory and cardiovascular diseases at the single-cell level will usher in cell-based interceptive medicine in Europe over the next decade.

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.

Cancer Stem Cell-Inducing Media Activates Senescence Reprogramming in Fibroblasts

Cellular senescence is a tumor-suppressive mechanism blocking cell proliferation in response to stress. However, recent evidence suggests that senescent tumor cells can re-enter the cell cycle to become cancer stem cells, leading to relapse after cancer chemotherapy treatment. Understanding how the senescence reprogramming process is a precursor to cancer stem cell formation is of great medical importance. To study the interplay between senescence, stemness, and cancer, we applied a stem cell medium (SCM) to human embryonic fibroblasts (MRC5 and WI-38) and cancer cell lines (A549 and 293T). MRC5 and WI-38 cells treated with SCM showed symptoms of oxidative stress and became senescent. Transcriptome analysis over a time course of SCM-induced senescence, revealed a developmental process overlapping with the upregulation of genes for growth arrest and the senescence-associated secretory phenotype (SASP). We demonstrate that histone demethylases jumonji domain-containing protein D3 (Jmjd3) and ubiquitously transcribed tetratricopeptide repeat, X chromosome (Utx), which operate by remodeling chromatin structure, are implicated in the senescence reprogramming process to block stem cell formation in fibroblasts. In contrast, A549 and 293T cells cultured in SCM were converted to cancer stem cells that displayed the phenotype of senescence uncoupled from growth arrest. The direct overexpression of DNA methyltransferases (Dnmt1 and Dnmt3A), ten-eleven translocation methylcytosine dioxygenases (Tet1 and Tet3), Jmjd3, and Utx proteins could activate senescence-associated beta-galactosidase (SA-β-gal) activity in 293T cells, suggesting that epigenetic alteration and chromatin remodeling factors trigger the senescence response. Overall, our study suggests that chromatin machinery controlling senescence reprogramming is significant in cancer stem cell formation.

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.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

The RNA-RNA base pairing potential of human Dicer and Ago2 proteins

The ribonuclease Dicer produces microRNAs (miRNAs) and small interfering RNAs that are handed over to Ago proteins to control gene expression by targeting complementary sequences within transcripts. Interestingly, a growing number of reports have demonstrated that the activity of Dicer may extend beyond the biogenesis of small regulatory RNAs. Among them, a report from our latest studies revealed that human Dicer facilitates base pairing of complementary sequences present in two nucleic acids, thus acting as a nucleic acid annealer. Accordingly, in this manuscript, we address how RNA structure influences the annealing activity of human Dicer. We show that Dicer supports hybridization between a small RNA and a complementary sequence of a longer RNA in vitro, even when both complementary sequences are trapped within secondary structures. Moreover, we show that under applied conditions, human Ago2, a core component of RNA-induced silencing complex, displays very limited annealing activity. Based on the available data from new-generation sequencing experiments regarding the RNA pool bound to Dicer in vivo, we show that multiple Dicer-binding sites within mRNAs also contain miRNA targets. Subsequently, we demonstrate in vitro that Dicer but not Ago2 can anneal miRNA to its target present within mRNA. We hypothesize that not all miRNA duplexes are handed over to Ago proteins. Instead, miRNA-Dicer complexes could target specific sequences within transcripts and either compete or cooperate for binding sites with miRNA-Ago complexes. Thus, not only Ago but also Dicer might be directly involved in the posttranscriptional control of gene expression.

Spider Chitin: An Ultrafast Microwave-Assisted Method for Chitin Isolation from Caribena versicolor Spider Molt Cuticle

Chitin, as a fundamental polysaccharide in invertebrate skeletons, continues to be actively investigated, especially with respect to new sources and the development of effective methods for its extraction. Recent attention has been focused on marine crustaceans and sponges; however, the potential of spiders (order Araneae) as an alternative source of tubular chitin has been overlooked. In this work, we focused our attention on chitin from up to 12 cm-large Theraphosidae spiders, popularly known as tarantulas or bird-eating spiders. These organisms "lose" large quantities of cuticles during their molting cycle. Here, we present for the first time a highly effective method for the isolation of chitin from Caribena versicolor spider molt cuticle, as well as its identification and characterization using modern analytical methods. We suggest that the tube-like molt cuticle of this spider can serve as a naturally prefabricated and renewable source of tubular chitin with high potential for application in technology and biomedicine.

Human ribonuclease Dicer – structure and functions

Endoribonuclease III Dicer plays a crucial role in the biogenesis of small regulatory RNAs, such as microRNAs (miRNAs) and small inter­fering RNAs (siRNAs). However, this is not the only role that Dicer plays in cells. For example, it has been shown that Dicer is involved in processing of diverse classes of RNA, including tRNA and snoRNA, cleavage of repeat-element-derived RNAs, and maintenance of genome integrity. Dicer has also been found to participate in the chromosome fragmentation during apoptosis or in the inflammatory processes. More­over, a recent discovery of Dicer-binding passive sites in mRNAs and long non-coding RNAs, and its putative nucleic acid chaperone activity, has pointed out a novel regulatory role of the enzyme. Here we focus on human Dicer and review its structure and function including recent findings on miRNA-independent roles and their impact on cell biology.

Spider Chitin. The biomimetic potential and applications of Caribena versicolor tubular chitin

Diverse fields of modern technology and biomedicine can benefit from the application of ready-to-use chitin-based scaffolds. In this work we show for the first time the applicability of tubular and porous chitin from Caribena versicolor spiders as a scaffold for the development of an effective CuO/Cu(OH)₂ catalyst for the reduction of 4-nitrophenol (4-NP) to 4-aminophenol (4-AM), and as a scaffold for the tissue engineering of selected cells. The formation of CuO/Cu(OH)₂ phases on and within the chitinous tubes leads to a hybrid material with excellent catalytic performance with respect to the reduction of p-nitrophenol. On the other hand, experimental results provide for the first time strong evidence for the biocompatibility of spider chitin with different cell types, a human progenitor cell line (hPheo1), as well as cardiomyocytes differentiated from induced pluripotent stem cells (iPSC-CMs) that were cultured on a tube-like scaffold.

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.

Identification of Altered Developmental Pathways in Human Juvenile HD iPSC With 71Q and 109Q Using Transcriptome Profiling

In Huntington disease (HD) subtle symptoms in patients may occur years or even decades prior to diagnosis. HD changes at a molecular level may begin as early as in cells that are non-lineage committed such as stem cells or HD patients induced pluripotent stem cells (iPSCs) offering opportunity to enhance the understanding of the HD pathogenesis. In addition, juvenile HD non-linage committed cells were previously not directly investigated in detail by RNA-seq. In the present manuscript, we define the early HD and juvenile HD transcriptional alterations using 6 human HD iPS cell lines from two patients, one with 71 CAGs and one with 109 CAG repeats. We identified 107 (6 HD lines), 198 (3 HD71Q lines) and 217 (3 HD109Q lines) significantly dysregulated mRNAs in each comparison group. The analyses showed that many of dysregulated transcripts in HD109Q iPSC lines are involved in DNA damage response and apoptosis, such as CCND1, CDKN1A, TP53, BAX, TNFRSF10B, TNFRSF10C, TNFRSF10D, DDB2, PLCB1, PRKCQ, HSH2D, ZMAT3, PLK2, and RPS27L. Most of them were identified as downregulated and their proteins are direct interactors with TP53. HTT probably alters the level of several TP53 interactors influencing apoptosis. This may lead to accumulation of an excessive number of progenitor cells and potential disruption of cell differentiation and production of mature neurons. In addition, HTT effects on cell polarization also demonstrated in the analysis may result in a generation of incorrect progenitors. Bioinformatics analysis of transcripts dysregulated in HD71Q iPSC lines showed that several of them act as transcription regulators during the early multicellular stages of development, such as ZFP57, PIWIL2, HIST1H3C, and HIST1H2BB. Significant upregulation of most of these transcripts may lead to a global increase in expression level of genes involved in pathways critical for embryogenesis and early neural development. In addition, MS analysis revealed altered levels of TP53 and ZFP30 proteins reflecting the functional significance of dysregulated mRNA levels of these proteins which were associated with apoptosis and DNA binding. Our finding very well corresponds to the fact that mutation in the HTT gene may cause precocious neurogenesis and identifies pathways likely disrupted during development.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

ENHO, RXRA, and LXRA polymorphisms and dyslipidaemia, related comorbidities and survival in haemodialysis patients

BACKGROUND

The energy homeostasis-associated gene (ENHO), retinoid X receptor alpha gene (RXRA), and liver X receptor alpha gene (LXRA) are involved in adipogenic/lipogenic regulation. We investigated whether single-nucleotide polymorphisms in these genes (ENHO rs2281997, rs72735260; RXRA rs749759, rs10776909, rs10881578; LXRA rs2279238, rs7120118, rs11039155) are associated with dyslipidaemia, related comorbidities and survival of haemodialysis (HD) patients also tested for T-helper (Th) cell interleukin genes (IL).

METHODS

The study was carried out in 873 HD patients. Dyslipidaemia was diagnosed by the recommendations of the Kidney Disease Outcomes Quality Initiative (K/DOQI) guidelines (2003); atherogenic dyslipidaemia was referred to if the TG/HDL cholesterol ratio was equal to or higher than 3.8. Genotyping of ENHO SNPs, LXRA SNPs, and IL12A rs568408 was carried out using HRM analysis. RXRA SNPs, IL12B rs3212227, and IL18 rs360719 were genotyped using PCR-RFLP analysis. The circulating adropin concentration was determined in 126 patients by enzyme-linked immunosorbent assay. Survival probability was analysed using the Kaplan-Meier method in 440 patients followed through 7.5 years.

RESULTS

Dyslipidaemia by K/DOQI was diagnosed in 459 patients (91% revealed hyper-LDL- cholesterolaemia), atherogenic dyslipidaemia was diagnosed in 454 patients, and 231 patients were free of dyslipidaemia by both criteria. The variant allele (T) of ENHO rs2281997 was associated with the hyper-LDL cholesterolaemic pattern of dyslipidaemia by K/DOQI. The frequency of atherogenic dyslipidaemia was lower in T-allele bearers than in CC-genotype patients. The rs2281997 T allele was associated with lower cardiovascular mortality in HD patients showing atherogenic dyslipidaemia. ENHO, RXRA, and LXRA showed epistatic interactions in dyslipidaemia. Circulating adropin was lower in atherogenic dyslipidaemia than in non-atherogenic conditions. RXRA rs10776909 was associated with myocardial infarction. Bearers of LXRA rs2279238, rs7120118 or rs11039155 minor alleles showed higher mortality. ENHO SNP positions fell within the same DNase 1 hypersensitivity site expressed in the Th1 cell line. Epistatic interactions occurred between rs2281997 and Th1 IL SNPs (rs360719, rs568408).

CONCLUSIONS

Atherogenic dyslipidaemia occurs in HD patients in whom ENHO encodes less adropin. ENHO, RXRA, and LXRA SNPs, separately or jointly, are associated with dyslipidaemia, myocardial infarction, and survival in HD patients. Differences in the availability of transcription binding sites may contribute to these associations.

The Defense Response of Nicotiana benthamiana to Peanut Stunt Virus Infection in the Presence of Symptom Exacerbating Satellite RNA

Peanut stunt virus (PSV) is a widespread disease infecting legumes. The PSV strains are classified into four subgroups and some are defined by the association of satellite RNAs (satRNAs). In the case of PSV, the presence of satRNAs alters the symptoms of disease in infected plants. In this study, we elucidated the plant response to PSV-G strain, which occurs in natural conditions without satRNA. However, it was found that it might easily acquire satRNA, which exacerbated pathogenesis in Nicotiana benthamiana. To explain the mechanisms underlying PSV infection and symptoms exacerbation caused by satRNA, we carried out transcriptome profiling of N. benthamiana challenged by PSV-G and satRNA using species-specific microarrays. Co-infection of plants with PSV-G + satRNA increased the number of identified differentially expressed genes (DEGs) compared with the number identified in PSV-G-infected plants. In both treatments, the majority of up-regulated DEGs were engaged in translation, ribosome biogenesis, RNA metabolism, and response to stimuli, while the down-regulated DEGs were required for photosynthesis. The presence of satRNA in PSV-G-infected plants caused different trends in expression of DEGs associated with phosphorylation, ATP binding, and plasma membrane.

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.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1608-2) contains supplementary material, which is available to authorized users.

Sorghum CCoAOMT and CCoAOMT-like gene evolution, structure, expression and the role of conserved amino acids in protein activity

Sorghum is a crop plant that is grown for seeds, sucrose, forage and biofuel production. In all these applications, lignin is a superfluous component that decreases the efficiency of technological processes. Caffeoyl-coenzyme A O-methyltransferase (CCoAOMT) is an enzyme involved in monolignol synthesis that affects the efficiency of lignification and lignin composition. The sorghum genome harbors one CCoAOMT gene and six closely related CCoAOMT-like genes. The structures of four sorghum CCoAOMT-like enzymes suggest that these proteins might methylate caffeoyl coenzyme A and contribute to monolignol synthesis. In this study, two sorghum genes, CCoAOMT and one CCoAOMT-like, were found to be highly expressed in leaves, stems and immature seeds. The promoters of these genes possess clusters of transcription factor-binding sites specific for lignification, and this suggests that they are important for lignification. Phylogenetic analysis revealed that one sorghum CCoAOMT-like enzyme is closely related to ancestral cyanobacterial CCoAOMT-like proteins. The remaining CCoAOMT-like enzymes, including the one highly expressed in the leaves and stem, are closely related to CCoAOMT. Genes from these two groups possess different, evolutionarily conserved gene structures. The structure of the sorghum CCoAOMT-like protein from the ancestral clade was modeled and differences between enzymes from the two clades were analyzed. These results facilitate a better understanding of the evolution of genes involved in lignification, and provide valuable data for sorghum improvement through traditional breeding or molecular genetic techniques. The findings suggest that CCoAOMT-like genes might be recruited in lignification and raise questions of the frequency of such functional shifts.

MTTE: an innovative strategy for the evaluation of targeted/exome enrichment efficiency

Although currently available strategies for the preparation of exome-enriched libraries are well established, a final validation of the libraries in terms of exome enrichment efficiency prior to the sequencing step is of considerable importance. Here, we present a strategy for the evaluation of exome enrichment, i.e., the Multipoint Test for Targeted-enrichment Efficiency (MTTE), PCR-based approach utilizing multiplex ligation-dependent probe amplification with capillary electrophoresis separation. We used MTTE for the analysis of subsequent steps of the Illumina TruSeq Exome Enrichment procedure. The calculated values of enrichment-associated parameters (i.e., relative enrichment, relative clearance, overall clearance, and fold enrichment) and the comparison of MTTE results with the actual enrichment revealed the high reliability of our assay. Additionally, the MTTE assay enabled the determination of the sequence-associated features that may confer bias in the enrichment of different targets. Importantly, the MTTE is low cost method that can be easily adapted to the region of interest important for a particular project. Thus, the MTTE strategy is attractive for post-capture validation in a variety of targeted/exome enrichment NGS projects.

Gene expression profiling of acute myeloid leukemia samples from adult patients with AML-M1 and -M2 through boutique microarrays, real-time PCR and droplet digital PCR

Acute myeloid leukemia (AML) is the most common and severe form of acute leukemia diagnosed in adults. Owing to its heterogeneity, AML is divided into classes associated with different treatment outcomes and specific gene expression profiles. Based on previous studies on AML, in this study, we designed and generated an AML-array containing 900 oligonucleotide probes complementary to human genes implicated in hematopoietic cell differentiation and maturation, proliferation, apoptosis and leukemic transformation. The AML-array was used to hybridize 118 samples from 33 patients with AML of the M1 and M2 subtypes of the French-American-British (FAB) classification and 15 healthy volunteers (HV). Rigorous analysis of the microarray data revealed that 83 genes were differentially expressed between the patients with AML and the HV, including genes not yet discussed in the context of AML pathogenesis. The most overexpressed genes in AML were STMN1, KITLG, CDK6, MCM5, KRAS, CEBPA, MYC, ANGPT1, SRGN, RPLP0, ENO1 and SET, whereas the most underexpressed genes were IFITM1, LTB, FCN1, BIRC3, LYZ, ADD3, S100A9, FCER1G, PTRPE, CD74 and TMSB4X. The overexpression of the CPA3 gene was specific for AML with mutated NPM1 and FLT3. Although the microarray-based method was insufficient to differentiate between any other AML subgroups, quantitative PCR approaches enabled us to identify 3 genes (ANXA3, S100A9 and WT1) whose expression can be used to discriminate between the 2 studied AML FAB subtypes. The expression levels of the ANXA3 and S100A9 genes were increased, whereas those of WT1 were decreased in the AML-M2 compared to the AML-M1 group. We also examined the association between the STMN1, CAT and ABL1 genes, and the FLT3 and NPM1 mutation status. FLT3⁺/NPM1⁻ AML was associated with the highest expression of STMN1, and ABL1 was upregulated in FLT3⁺ AML and CAT in FLT3⁻ AML, irrespectively of the NPM1 mutation status. Moreover, our results indicated that CAT and WT1 gene expression levels correlated with the response to therapy. CAT expression was highest in patients who remained longer under complete remission, whereas WT1 expression increased with treatment resistance. On the whole, this study demonstrates that the AML-array can potentially serve as a first-line screening tool, and may be helpful for the diagnosis of AML, whereas the differentiation between AML subgroups can be more successfully performed with PCR-based analysis of a few marker genes.

Effects of G-quadruplex topology on translational inhibition by tRNA fragments in mammalian and plant systems in vitro

The folding of tRNA fragments (tRFs) into G-quadruplex structures and the implications of G-quadruplexes in translational inhibition have been studied mainly in mammalian systems. To increase our knowledge of these phenomena, we determined the influence of human and plant tRFs and model G-quadruplexes on translation in rabbit reticulocyte lysate and wheat germ extract. The efficiency of translational inhibition in the mammalian system was strongly associated with the type of G-quadruplex topology. In the plant system, the ability of a small RNA to adopt the G-quadruplex conformation was not sufficient to repress translation, indicating the importance of other structural determinants.

Corrigendum: The Generation of Mouse and Human Huntington Disease iPS Cells Suitable for In vitro Studies on Huntingtin Function

[This corrects the article on p. 253 in vol. 10, PMID: 28848389.].

The Generation of Mouse and Human Huntington Disease iPS Cells Suitable for In vitro Studies on Huntingtin Function

Huntington disease (HD) is an incurable neurodegenerative disorder caused by expansion of CAG repeats in huntingtin (HTT) gene, resulting in expanded polyglutamine tract in HTT protein. Although, HD has its common onset in adulthood, subtle symptoms in patients may occur decades before diagnosis, and molecular and cellular changes begin much earlier, even in cells that are not yet lineage committed such as stem cells. Studies in induced pluripotent stem cell (iPSC) HD models have demonstrated that multiple molecular processes are altered by the mutant HTT protein and suggested its silencing as a promising therapeutic strategy. Therefore, we aimed to generate HD iPS cells with stable silencing of HTT and further to investigate the effects of HTT knock-down on deregulations of signaling pathways e.g., p53 downregulation, present in cells already in pluripotent state. We designed a gene silencing strategy based on RNAi cassette in piggyBAC vector for constant shRNA expression. Using such system we delivered and tested several shRNA targeting huntingtin in mouse HD YAC128 iPSC and human HD109, HD71, and Control iPSC. The most effective shRNA (shHTT2) reagent stably silenced HTT in all HD iPS cells and remained active upon differentiation to neural stem cells (NSC). When investigating the effects of HTT silencing on signaling pathways, we found that in mouse HD iPSC lines expressing shRNA the level of mutant HTT inversely correlated with p53 levels, resulting in p53 level normalization upon silencing of mutant HTT. We also found that p53 deregulation continues into the NSC developmental stage and it was reversed upon HTT silencing. In addition, we observed subtle effects of silencing on proteins of Wnt/β-catenin and ERK1/2 signaling pathways. In summary, we successfully created the first mouse and human shRNA-expressing HD iPS cells with stable and continuous HTT silencing. Moreover, we demonstrated reversal of HD p53 phenotype in mouse HD iPSC, therefore, the stable knockdown of HTT is well-suited for investigation on HD cellular pathways, and is potentially useful as a stand-alone therapy or component of cell therapy. In addition, the total HTT knock-down in our human cells has further implications for mutant allele selective approach in iPSC.

Comprehensive analysis of microorganisms accompanying human archaeological remains

Metagenome analysis has become a common source of information about microbial communities that occupy a wide range of niches, including archaeological specimens. It has been shown that the vast majority of DNA extracted from ancient samples come from bacteria (presumably modern contaminants). However, characterization of microbial DNA accompanying human remains has never been done systematically for a wide range of different samples. We used metagenomic approaches to perform comparative analyses of microorganism communities present in 161 archaeological human remains. DNA samples were isolated from the teeth of human skeletons dated from 100 AD to 1200 AD. The skeletons were collected from 7 archaeological sites in Central Europe and stored under different conditions. The majority of identified microbes were ubiquitous environmental bacteria that most likely contaminated the host remains not long ago. We observed that the composition of microbial communities was sample-specific and not correlated with its temporal or geographical origin. Additionally, traces of bacteria and archaea typical for human oral/gut flora, as well as potential pathogens, were identified in two-thirds of the samples. The genetic material of human-related species, in contrast to the environmental species that accounted for the majority of identified bacteria, displayed DNA damage patterns comparable with endogenous human ancient DNA, which suggested that these microbes might have accompanied the individual before death. Our study showed that the microbiome observed in an individual sample is not reliant on the method or duration of sample storage. Moreover, shallow sequencing of DNA extracted from ancient specimens and subsequent bioinformatics analysis allowed both the identification of ancient microbial species, including potential pathogens, and their differentiation from contemporary species that colonized human remains more recently.

Small RNA fragments derived from multiple RNA classes - the missing element of multi-omics characteristics of the hepatitis C virus cell culture model

BACKGROUND

A pool of small RNA fragments (RFs) derived from diverse cellular RNAs has recently emerged as a rich source of functionally relevant molecules. Although their formation and accumulation has been connected to various stress conditions, the knowledge on RFs produced upon viral infections is very limited. Here, we applied the next generation sequencing (NGS) to characterize RFs generated in the hepatitis C virus (HCV) cell culture model (HCV-permissive Huh-7.5 cell line).

RESULTS

We found that both infected and non-infected cells contained a wide spectrum of RFs derived from virtually all RNA classes. A significant fraction of identified RFs accumulated to similar levels as miRNAs. Our analysis, focused on RFs originating from constitutively expressed non-coding RNAs, revealed three major patterns of parental RNA cleavage. We found that HCV infection induced significant changes in the accumulation of low copy number RFs, while subtly altered the levels of high copy number ones. Finally, the candidate RFs potentially relevant for host-virus interactions were identified.

CONCLUSIONS

Our results indicate that RFs should be considered an important component of the Huh-7.5 transcriptome and suggest that the main factors influencing the RF biogenesis are the RNA structure and RNA protection by interacting proteins. The data presented here significantly complement the existing transcriptomic, miRnomic, proteomic and metabolomic characteristics of the HCV cell culture model.

Mutations in human AID differentially affect its ability to deaminate cytidine and 5-methylcytidine in ssDNA substrates in vitro

Activation-induced cytidine deaminase (AID) is known for its established role in antibody production. AID induces the diversification of antibodies by deaminating deoxycytidine (C) within immunoglobulin genes. The capacity of AID to deaminate 5-methyldeoxycytidine (5 mC) and/or 5-hydroxymethyldeoxycytidine (5 hmC), and consequently AID involvement in active DNA demethylation, is not fully resolved. For instance, structural determinants of AID activity on different substrates remain to be identified. To better understand the latter issue, we tested how mutations in human AID (hAID) influence its ability to deaminate C, 5 mC, and 5 hmC in vitro. We showed that each of the selected mutations differentially affects hAID’s ability to deaminate C and 5 mC. At the same time, we did not observe hAID activity on 5 hmC. Surprisingly, we found that the N51A hAID mutant, with no detectable activity on C, efficiently deaminated 5 mC, which may suggest different requirements for C and 5 mC deamination. Homology modeling and molecular dynamics simulations revealed that the pattern of enzyme-substrate recognition is one of the important factors determining enzyme activity on C and 5 mC. Consequently, we have proposed mechanisms that explain why wild type hAID more efficiently deaminates C than 5 mC in vitro and why 5 hmC is not deaminated.

Mutations in proteasome-related genes are associated with thyroid hemiagenesis

PURPOSE

Human thyroid development is a complex and still unexplained process. Thyroid hemiagenesis is a congenital anomaly, where one of the thyroid lobes fails to develop. In the majority of patients with thyroid hemiagenesis, the genetic background remains unknown. The aim of the study was to search for novel genetic contributors to the etiology of thyroid hemiagenesis.

METHODS

A cohort of 34 sporadic patients diagnosed with thyroid hemiagenesis and one three-generation family were subjected to comprehensive genomic examination. Initially, targeted screening of associated transcription factors, known to be linked to thyroid development, was performed. As a next step, genomic examinations were applied using high-resolution microarrays, whereas for the thyroid hemiagenesis family, additionally the whole exome sequencing was performed.

RESULTS

Screening of transcription factors revealed no causative mutations in the studied cohort. Genomic examinations revealed the presence of four recurrent defects (three deletions and one duplication) affecting highly conservative proteasome genes PSMA1, PSMA3, and PSMD3. In a thyroid hemiagenesis family a splice site mutation in a proteasome gene PSMD2 (c.612T > C cDNA.1170T > C, g.3271T > C) was found in both affected mother and daughter.

CONCLUSIONS

Our results shed a new light on etiology of thyroid hemiagenesis, so far suspected to be linked only to mutations in the genes directly involved in the thyroid development. We demonstrated, for the first time, that genomic alterations in proteasome-associated genes co-occur in patients presenting this developmental anomaly.

EvOligo: A Novel Software to Design and Group Libraries of Oligonucleotides Applicable for Nucleic Acid-Based Experiments

Experimental methods based on DNA and RNA hybridization, such as multiplex polymerase chain reaction, multiplex ligation-dependent probe amplification, or microarray analysis, require the use of mixtures of multiple oligonucleotides (primers or probes) in a single test tube. To provide an optimal reaction environment, minimal self- and cross-hybridization must be achieved among these oligonucleotides. To address this problem, we developed EvOligo, which is a software package that provides the means to design and group DNA and RNA molecules with defined lengths. EvOligo combines two modules. The first module performs oligonucleotide design, and the second module performs oligonucleotide grouping. The software applies a nearest-neighbor model of nucleic acid interactions coupled with a parallel evolutionary algorithm to construct individual oligonucleotides, and to group the molecules that are characterized by the weakest possible cross-interactions. To provide optimal solutions, the evolutionary algorithm sorts oligonucleotides into sets, preserves preselected parts of the oligonucleotides, and shapes their remaining parts. In addition, the oligonucleotide sets can be designed and grouped based on their melting temperatures. For the user's convenience, EvOligo is provided with a user-friendly graphical interface. EvOligo was used to design individual oligonucleotides, oligonucleotide pairs, and groups of oligonucleotide pairs that are characterized by the following parameters: (1) weaker cross-interactions between the non-complementary oligonucleotides and (2) more uniform ranges of the oligonucleotide pair melting temperatures than other available software products. In addition, in contrast to other grouping algorithms, EvOligo offers time-efficient sorting of paired and unpaired oligonucleotides based on various parameters defined by the user.

MLPA-Based Analysis of Copy Number Variation in Plant Populations

Copy number variants (CNVs) are intraspecies duplications/deletions of large DNA segments (>1 kb). A growing number of reports highlight the functional and evolutionary impact of CNV in plants, increasing the need for appropriate tools that enable locus-specific CNV genotyping on a population scale. Multiplex ligation-dependent probe amplification (MLPA) is considered a gold standard in genotyping CNV in humans. Consequently, numerous commercial MLPA assays for CNV-related human diseases have been created. We routinely genotype complex multiallelic CNVs in human and plant genomes using the modified MLPA procedure based on fully synthesized oligonucleotide probes (90-200 nt), which greatly simplifies the design process and allows for the development of custom assays. Here, we present a step-by-step protocol for gene-specific MLPA probe design, multiplexed assay setup and data analysis in a copy number genotyping experiment in plants. As a case study, we present the results of a custom assay designed to genotype the copy number status of 12 protein coding genes in a population of 80 Arabidopsis accessions. The genes were pre-selected based on whole genome sequencing data and are localized in the genomic regions that display different levels of population-scale variation (non-variable, biallelic, or multiallelic, as well as CNVs overlapping whole genes or their fragments). The presented approach is suitable for population-scale validation of the CNV regions inferred from whole genome sequencing data analysis and for focused analysis of selected genes of interest. It can also be very easily adopted for any plant species, following optimization of the template amount and design of the appropriate control probes, according to the general guidelines presented in this paper.

Modeling of the catalytic core of Arabidopsis thaliana Dicer-like 4 protein and its complex with double-stranded RNA

Plant Dicer-like proteins (DCLs) belong to the Ribonuclease III (RNase III) enzyme family. They are involved in the regulation of gene expression and antiviral defense through RNA interference pathways. A model plant, Arabidopsis thaliana encodes four DCL proteins (AtDCL1-4) that produce different classes of small regulatory RNAs. Our studies focus on AtDCL4 that processes double-stranded RNAs (dsRNAs) into 21 nucleotide trans-acting small interfering RNAs. So far, little is known about the structures of plant DCLs and the complexes they form with dsRNA. In this work, we present models of the catalytic core of AtDCL4 and AtDCL4-dsRNA complex constructed by computational methods. We built a homology model of the catalytic core of AtDCL4 comprising Platform, PAZ, Connector helix and two RNase III domains. To assemble the AtDCL4-dsRNA complex two modeling approaches were used. In the first method, to establish conformations that allow building a consistent model of the complex, we used Normal Mode Analysis for both dsRNA and AtDCL4. The second strategy involved template-based approach for positioning of the PAZ domain and manual arrangement of the Connector helix. Our results suggest that the spatial orientation of the Connector helix, Platform and PAZ relative to the RNase III domains is crucial for measuring dsRNA of defined length. The modeled complexes provide information about interactions that may contribute to the relative orientations of these domains and to dsRNA binding. All these information can be helpful for understanding the mechanism of AtDCL4-mediated dsRNA recognition and binding, to produce small RNA of specific size.

Arabidopsis thaliana population analysis reveals high plasticity of the genomic region spanning MSH2, AT3G18530 and AT3G18535 genes and provides evidence for NAHR-driven recurrent CNV events occurring in this location

Background

Intraspecies copy number variations (CNVs), defined as unbalanced structural variations of specific genomic loci, ≥1 kb in size, are present in the genomes of animals and plants. A growing number of examples indicate that CNVs may have functional significance and contribute to phenotypic diversity. In the model plant Arabidopsis thaliana at least several hundred protein-coding genes might display CNV; however, locus-specific genotyping studies in this plant have not been conducted.

Results

We analyzed the natural CNVs in the region overlapping MSH2 gene that encodes the DNA mismatch repair protein, and AT3G18530 and AT3G18535 genes that encode poorly characterized proteins. By applying multiplex ligation-dependent probe amplification and droplet digital PCR we genotyped those genes in 189 A. thaliana accessions. We found that AT3G18530 and AT3G18535 were duplicated (2–14 times) in 20 and deleted in 101 accessions. MSH2 was duplicated in 12 accessions (up to 12-14 copies) but never deleted. In all but one case, the MSH2 duplications were associated with those of AT3G18530 and AT3G18535. Considering the structure of the CNVs, we distinguished 5 genotypes for this region, determined their frequency and geographical distribution. We defined the CNV breakpoints in 35 accessions with AT3G18530 and AT3G18535 deletions and tandem duplications and showed that they were reciprocal events, resulting from non-allelic homologous recombination between 99 %-identical sequences flanking these genes. The widespread geographical distribution of the deletions supported by the SNP and linkage disequilibrium analyses of the genomic sequence confirmed the recurrent nature of this CNV.

Conclusions

We characterized in detail for the first time the complex multiallelic CNV in Arabidopsis genome. The region encoding MSH2, AT3G18530 and AT3G18535 genes shows enormous variation of copy numbers among natural ecotypes, being a remarkable example of high Arabidopsis genome plasticity. We provided the molecular insight into the mechanism underlying the recurrent nature of AT3G18530-AT3G18535 duplications/deletions. We also performed the first direct comparison of the two leading experimental methods, suitable for assessing the DNA copy number status. Our comprehensive case study provides foundation information for further analyses of CNV evolution in Arabidopsis and other plants, and their possible use in plant breeding.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3221-1) contains supplementary material, which is available to authorized users.