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.

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.

RNA and DNA G-quadruplexes bind to human dicer and inhibit its activity.

Guanine (G)-rich single-stranded nucleic acids can adopt G-quadruplex structures. Accumulating evidence indicates that G-quadruplexes serve important regulatory roles in fundamental biological processes such as DNA replication, transcription, and translation, while aberrant G-quadruplex formation is linked to genome instability and cancer. Understanding the biological functions played by G-quadruplexes requires detailed knowledge of their protein interactome. Here, we report that both RNA and DNA G-quadruplexes are bound by human Dicer in vitro. Using in vitro binding assays, mutation studies, and computational modeling we demonstrate that G-quadruplexes can interact with the Platform-PAZ-Connector helix cassette of Dicer, the region responsible for anchoring microRNA precursors (pre-miRNAs). Consequently, we show that G-quadruplexes efficiently and stably inhibit the cleavage of pre-miRNA by Dicer. Our data highlight the potential of human Dicer for binding of G-quadruplexes and allow us to propose a G-quadruplex-driven sequestration mechanism of Dicer regulation.

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.

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.

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.

Functional characterization of RNA fragments using high-throughput interactome screening.

Populations of small eukaryotic RNAs, in addition to relatively well recognized molecules such as miRNAs or siRNAs, also contain fragments derived from all classes of constitutively expressed non-coding RNAs. It has been recently demonstrated that the formation and accumulation of RNA fragments (RFs) is cell-/tissue-specific and depends on internal and external stimuli. Unfortunately, the mechanisms underlying RF biogenesis and function remain unclear. To better understand them, we employed RNA pull-down and mass spectrometry methods to characterize the interactions of seven RFs originating from tRNA, snoRNA and snRNA. By integrating our results with publicly available data on physical protein-protein interactions, we constructed an RF interactome network. We determined that the RF interactome comprises proteins generally different from those that interact with their parental full length RNAs. Proteins captured by the RFs were involved in mRNA splicing, tRNA processing, DNA recombination/replication, protein biosynthesis and carboxylic acid metabolism. Our data suggest that RFs can be endogenous aptamer-like molecules and potential players in recently revealed RNA-protein regulatory networks. SIGNIFICANCE: In the recent decade it has become evident that RNAs with well-known functions (for example tRNA, snoRNA or rRNA) can be cleaved to yield short fragments, whose role in cells remains only partially characterized. At the same time, unconventional interactions between mRNA and proteins without RNA-binding domains have been demonstrated, revealing novel layers of possible RNA-mediated regulation. Considering the above, we hypothesized that RNA fragments (RFs) can be endogenous aptamer-like molecules that unconventionally interact with proteins. In this study we identified protein partners of seven selected RFs. We found that RFs bind different set of proteins than their parental full length RNAs and identified proteins differentially bound by the particular RFs. These observations suggest biological relevance of the discovered interactions. Our data provide a novel perspective on the significance of RFs and point to this pool of molecules as to a rich collection of potential components of the recently discovered RNA-protein regulatory networks.

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.

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.

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.

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.

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.

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.

RNA-Seq-based analysis of differential gene expression associated with hepatitis C virus infection in a cell culture.

Hepatitis C virus (HCV) infection is one of the major causes of chronic liver diseases. Unfortunately, the mechanisms of HCV infection-induced liver injury and host-virus interactions are still not well recognized. To better understand these processes we determined the changes in the host gene expression that occur during HCV infection of Huh-7.5 cells. As a result, we identified genes that may contribute to the immune and metabolic cellular responses to infection. Pathway enrichment analysis indicated that HCV induced an increased expression of genes involved in mitogen-activated protein kinases signaling, adipocytokine signaling, cell cycle and nitrogen metabolism. In addition, the enrichment analyses of processes and molecular functions revealed that the up-regulated genes were mainly implicated in the negative regulation of phosphorylation. Construction of the pathway-gene-process network enabled exploration of a much more complex landscape of molecular interactions. Consequently, several essential processes altered by HCV infection were identified: negative regulation of cell cycle, response to endoplasmic reticulum stress, response to reactive oxygen species, toll-like receptor signaling and pattern recognition receptor signaling. The analyses of genes whose expression was decreased upon HCV infection showed that the latter were engaged in the metabolism of lipids and amino acids. Moreover, we observed disturbance in the cellular antiviral defense. Altogether, our results demonstrated that HCV infection elicits host response that includes a very wide range of cellular mechanisms. Our findings significantly broaden the understanding of complex processes that accompany HCV infection. Consequently, they may be used for developing new host-oriented therapeutic strategies.

Copy number variation of genes involved in the hepatitis C virus-human interactome.

Copy number variation (CNV) is a newly discovered form of intra-species genetic polymorphism that is defined as deletions or duplications of genome segments ranging from 1 kbp to several Mbp. CNV accounts for the majority of the genetic variation observed in humans (CNV regions cover more than 10% of the human genome); therefore, it may significantly influence both the phenotype and susceptibility to various diseases. Unfortunately, the impact of CNV on a number of diseases, including hepatitis C virus (HCV) infection, remains largely unexplored. Here, we analyzed 421 human genes encoding proteins that have been shown to interact with HCV proteins or genomic RNA (proteins from the HCV-human interactome). We found that 19 of the 421 candidate genes are located in putative CNV regions. For all of these genes, copy numbers were determined for European, Asiatic and African populations using the multiplex ligation-dependent amplification (MLPA) method. As a result, we identified 4 genes, IGLL1, MLLT4, PDPK1, PPP1R13L, for which the CN-genotype ranged from 1 to 6. All of these genes are involved in host-virus interaction; thus, their polymorphism has a potential impact on the development of HCV infection and/or therapy outcome.

A gene expression microarray for Nicotiana benthamiana based on de novo transcriptome sequence assembly.

BACKGROUND: Nicotiana benthamiana has been widely used in laboratories around the world for studying plant-pathogen interactions and posttranscriptional gene expression silencing. Yet the exploration of its transcriptome has lagged behind due to the lack of both adequate sequence information and genome-wide analysis tools, such as DNA microarrays. Despite the increasing use of high-throughput sequencing technologies, the DNA microarrays still remain a popular gene expression tool, because they are cheaper and less demanding regarding bioinformatics skills and computational effort. RESULTS: We designed a gene expression microarray with 103,747 60-mer probes, based on two recently published versions of N. benthamiana transcriptome (v.3 and v.5). Both versions were reconstructed from RNA-Seq data of non-strand-specific pooled-tissue libraries, so we defined the sense strand of the contigs prior to designing the probe. To accomplish this, we combined a homology search against Arabidopsis thaliana proteins and hybridization to a test 244k microarray containing pairs of probes, which represented individual contigs. We identified the sense strand in 106,684 transcriptome contigs and used this information to design an Nb-105k microarray on an Agilent eArray platform. Following hybridization of RNA samples from N. benthamiana roots and leaves we demonstrated that the new microarray had high specificity and sensitivity for detection of differentially expressed transcripts. We also showed that the data generated with the Nb-105k microarray may be used to identify incorrectly assembled contigs in the v.5 transcriptome, by detecting inconsistency in the gene expression profiles, which is indicated using multiple microarray probes that match the same v.5 primary transcripts. CONCLUSIONS: We provided a complete design of an oligonucleotide microarray that may be applied to the research of N. benthamiana transcriptome. This, in turn, will allow the N. benthamiana research community to take full advantage of microarray capabilities for studying gene expression in this plant. Additionally, by defining the sense orientation of over 106,000 contigs, we substantially improved the functional information on the N. benthamiana transcriptome. The simple hybridization-based approach for detecting the sense orientation of computationally assembled sequences can be used for updating the transcriptomes of other non-model organisms, including cases where no significant homology to known proteins exists.

Revealing a new activity of the human Dicer DUF283 domain in vitro.

The ribonuclease Dicer is a multidomain enzyme that plays a fundamental role in the biogenesis of small regulatory RNAs (srRNAs), which control gene expression by targeting complementary transcripts and inducing their cleavage or repressing their translation. Recent studies of Dicer's domains have permitted to propose their roles in srRNA biogenesis. For all of Dicer's domains except one, called DUF283 (domain of unknown function), their involvement in RNA substrate recognition, binding or cleavage has been postulated. For DUF283, the interaction with Dicer's protein partners has been the only function suggested thus far. In this report, we demonstrate that the isolated DUF283 domain from human Dicer is capable of binding single-stranded nucleic acids in vitro. We also show that DUF283 can act as a nucleic acid annealer that accelerates base-pairing between complementary RNA/DNA molecules in vitro. We further demonstrate an annealing activity of full length human Dicer. The overall results suggest that Dicer, presumably through its DUF283 domain, might facilitate hybridization between short RNAs and their targets. The presented findings reveal the complex nature of Dicer, whose functions may extend beyond the biogenesis of srRNAs.