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.

Background: Pairwise sequence alignment methods are widely used in biological research. The increasing number of sequences is perceived as one of the upcoming challenges for sequence alignment methods in the nearest future. To overcome this challenge several GPU (Graphics Processing Unit) computing approaches have been proposed lately. These solutions show a great potential of a GPU platform but in most cases address the problem of sequence database scanning and computing only the alignment score whereas the alignment itself is omitted. Thus, the need arose to implement the global and semiglobal Needleman-Wunsch, and Smith-Waterman algorithms with a backtracking procedure which is needed to construct the alignment. Results: In this paper we present the solution that performs the alignment of every given sequence pair, which is a required step for progressive multiple sequence alignment methods, as well as for DNA recognition at the DNA assembly stage. Performed tests show that the implementation, with performance up to 6.3 GCUPS on a single GPU for affine gap penalties, is very efficient in comparison to other CPU and GPU-based solutions. Moreover, multiple GPUs support with load balancing makes the application very scalable. Conclusions: The article shows that the backtracking procedure of the sequence alignment algorithms may be designed to fit in with the GPU architecture. Therefore, our algorithm, apart from scores, is able to compute pairwise alignments. This opens a wide range of new possibilities, allowing other methods from the area of molecular biology to take advantage of the new computational architecture. Performed tests show that the efficiency of the implementation is excellent. Moreover, the speed of our GPU-based algorithms can be almost linearly increased when using more than one graphics card.

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.

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.

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.