Genome and Transcriptome Sequencing of Populus × sibirica Identified Sex-Associated Allele-Specific Expression of the CLC Gene.

Transcriptome sequencing of leaves, catkin axes, and flowers from male and female trees of Populus × sibirica and genome sequencing of the same plants were performed for the first time. The availability of both genome and transcriptome sequencing data enabled the identification of allele-specific expression. Such an analysis was performed for genes from the sex-determining region (SDR). P. × sibirica is an intersectional hybrid between species from sections Aigeiros (Populus nigra) and Tacamahaca (Populus laurifolia, Populus suaveolens, or Populus × moskoviensis); therefore, a significant number of heterozygous polymorphisms were identified in the SDR that allowed us to distinguish between alleles. In the SDR, both allelic variants of the TCP (T-complex protein 1 subunit gamma), CLC (Chloride channel protein CLC-c), and MET1 (DNA-methyltransferase 1) genes were expressed in females, while in males, two allelic variants were expressed for TCP and MET1 but only one allelic variant prevailed for CLC. Targeted sequencing of TCP, CLC, and MET1 regions on a representative set of trees confirmed the sex-associated allele-specific expression of the CLC gene in generative and vegetative tissues of P. × sibirica. Our study brings new knowledge on sex-associated differences in Populus species.

Glutathione-related substances maintain cardiomyocyte contractile function in hypoxic conditions.

Severe hypoxia leads to decline in cardiac contractility and induces arrhythmic events in part due to oxidative damage to cardiomyocyte proteins including ion transporters. This results in compromised handling of Ca2+ ions that trigger heart contractile machinery. Here, we demonstrate that thiol-containing compounds such as N-acetylcysteine (NAC), glutathione ethyl ester (et-GSH), oxidized tetraethylglutathione (tet-GSSG), oxidized glutathione (GSSG) and S-nitrosoglutathione (GSNO) are capable of reducing negative effects of hypoxia on isolated rat cardiomyocytes. Preincubation of cardiomyocytes with 0.1 mM GSNO, 0.5 mM et-GSH, GSSG, tet-GSSG or with 10 mM NAC allows cells 5-times longer tolerate the hypoxic conditions and elicit regular Ca2+ transients in response to electric pacing. The shape of Ca2+ transients generated in the presence of GSNO, et-GSH and NAC was similar to that observed in normoxic control cardiomyocytes. The leader compound, GSNO, accelerated by 34% the recovery of normal contractile function of isolated rat heart subjected to ischemia-reperfusion. GSNO increased glutathionylation of Na,K-ATPase alpha-2 subunit, the principal ion-transporter of cardiac myocyte sarcolemma, which prevents irreversible oxidation of Na,K-ATPase and regulates its function to support normal Ca2+ ion handling in hypoxic cardiomyocytes. Altogether, GSNO appears effective cardioprotector in hypoxic conditions worth further studies toward its cardiovascular application.

Pan-Cancer Analysis of TCGA Data Revealed Promising Reference Genes for qPCR Normalization.

Quantitative PCR (qPCR) remains the most widely used technique for gene expression evaluation. Obtaining reliable data using this method requires reference genes (RGs) with stable mRNA level under experimental conditions. This issue is especially crucial in cancer studies because each tumor has a unique molecular portrait. The Cancer Genome Atlas (TCGA) project provides RNA-Seq data for thousands of samples corresponding to dozens of cancers and presents the basis for assessment of the suitability of genes as reference ones for qPCR data normalization. Using TCGA RNA-Seq data and previously developed CrossHub tool, we evaluated mRNA level of 32 traditionally used RGs in 12 cancer types, including those of lung, breast, prostate, kidney, and colon. We developed an 11-component scoring system for the assessment of gene expression stability. Among the 32 genes, PUM1 was one of the most stably expressed in the majority of examined cancers, whereas GAPDH, which is widely used as a RG, showed significant mRNA level alterations in more than a half of cases. For each of 12 cancer types, we suggested a pair of genes that are the most suitable for use as reference ones. These genes are characterized by high expression stability and absence of correlation between their mRNA levels. Next, the scoring system was expanded with several features of a gene: mutation rate, number of transcript isoforms and pseudogenes, participation in cancer-related processes on the basis of Gene Ontology, and mentions in PubMed-indexed articles. All the genes covered by RNA-Seq data in TCGA were analyzed using the expanded scoring system that allowed us to reveal novel promising RGs for each examined cancer type and identify several "universal" pan-cancer RG candidates, including SF3A1, CIAO1, and SFRS4. The choice of RGs is the basis for precise gene expression evaluation by qPCR. Here, we suggested optimal pairs of traditionally used RGs for 12 cancer types and identified novel promising RGs that demonstrate high expression stability and other features of reliable and convenient RGs (high expression level, low mutation rate, non-involvement in cancer-related processes, single transcript isoform, and absence of pseudogenes).

MethyMer: Design of combinations of specific primers for bisulfite sequencing of complete CpG islands.

We present MethyMer, a Python-based tool aimed at selecting primers for amplification of complete CpG islands. These regions are difficult in terms of selecting appropriate primers because of their low-complexity, high GC content. Moreover, bisulfite treatment, in fact, leads to the reduction of the 4-letter alphabet (ATGC) to 3-letter one (ATG, except for methylated cytosines), and this also reduces region complexity and increases mispriming potential. MethyMer has a flexible scoring system, which optimizes the balance between various characteristics such as nucleotide composition, thermodynamic features (melting temperature, dimers [Formula: see text]G, etc.), the presence of CpG sites and polyN tracts, and primer specificity, which is assessed with aligning primers to the bisulfite-treated genome using bowtie (up to three mismatches are allowed). Users are able to customize desired or limit ranges of various parameters as well as penalties for non-desired values. Moreover, MethyMer allows picking up the optimal combination of PCR primer pairs to perform the amplification of a large genomic locus, e.g. CpG island or other hard-to-study region, with minimal overlap of the individual amplicons. MethyMer incorporates ENCODE genome annotation records (promoter/enhancer/insulator), The Cancer Genome Atlas (TCGA) CpG methylation data derived with Illumina Infinium 450K microarrays, and records on correlations between TCGA RNA-Seq and CpG methylation data for 20 cancer types. These databases are included in the MethyMer release. Our tool is available at https://sourceforge.net/projects/methymer/ .

Cysteine residues 244 and 458-459 within the catalytic subunit of Na,K-ATPase control the enzyme's hydrolytic and signaling function under hypoxic conditions.

Our previous findings suggested that reversible thiol modifications of cysteine residues within the actuator (AD) and nucleotide binding domain (NBD) of the Na,K-ATPase may represent a powerful regulatory mechanism conveying redox- and oxygen-sensitivity of this multifunctional enzyme. S-glutathionylation of Cys244 in the AD and Cys 454-458-459 in the NBD inhibited the enzyme and protected cysteines' thiol groups from irreversible oxidation under hypoxic conditions. In this study mutagenesis approach was used to assess the role these cysteines play in regulation of the Na,K-ATPase hydrolytic and signaling functions. Several constructs of mouse α1 subunit of the Na,K-ATPase were produced in which Cys244, Cys 454-458-459 or Cys 244-454-458-459 were replaced by alanine. These constructs were expressed in human HEK293 cells. Non-transfected cells and those expressing murine α1 subunit were exposed to hypoxia or treated with oxidized glutathione (GSSG). Both conditions induced inhibition of the wild type Na,K-ATPase. Enzymes containing mutated mouse α1 lacking Cys244 or all four cysteines (Cys 244-454-458-459) were insensitive to hypoxia. Inhibitory effect of GSSG was observed for wild type murine Na,K-ATPase, but was less pronounced in Cys454-458-459Ala mutant and completely absent in the Cys244Ala and Cys 244-454-458-459Ala mutants. In cells, expressing wild type enzyme, ouabain induced activation of Src and Erk kinases under normoxic conditions, whereas under hypoxic conditions this effect was inversed. Cys454-458-459Ala substitution abolished Src kinase activation in response to ouabain treatment, uncoupled Src from Erk signaling, and interfered with O2-sensitivity of Na,K-ATPase signaling function. Moreover, modeling predicted that S-glutathionylation of Cys 458 and 459 should prevent inhibitory binding of Src to NBD. Our data indicate for the first time that cysteine residues within the AD and NBD influence hydrolytic as well as receptor function of the Na,K-ATPase and alter responses of the enzyme to hypoxia or upon treatment with cardiotonic steroids.

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes.

BACKGROUND: The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes. RESULTS: High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14). CONCLUSIONS: High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

Direct interaction of beta-amyloid with Na,K-ATPase as a putative regulator of the enzyme function.

By maintaining the Na(+) and K(+) transmembrane gradient mammalian Na,K-ATPase acts as a key regulator of neuronal electrotonic properties. Na,K-ATPase has an important role in synaptic transmission and memory formation. Accumulation of beta-amyloid (Aß) at the early stages of Alzheimer's disease is accompanied by reduction of Na,K-ATPase functional activity. The molecular mechanism behind this phenomenon is not known. Here we show that the monomeric Aß(1-42) forms a tight (Kd of 3 µM), enthalpy-driven equimolar complex with α1ß1 Na,K-ATPase. The complex formation results in dose-dependent inhibition of the enzyme hydrolytic activity. The binding site of Aß(1-42) is localized in the "gap" between the alpha- and beta-subunits of Na,K-ATPase, disrupting the enzyme functionality by preventing the subunits from shifting towards each other. Interaction of Na,K-ATPase with exogenous Aß(1-42) leads to a pronounced decrease of the enzyme transport and hydrolytic activity and Src-kinase activation in neuroblastoma cells SH-SY5Y. This interaction allows regulation of Na,K-ATPase activity by short-term increase of the Aß(1-42) level. However prolonged increase of Aß(1-42) level under pathological conditions could lead to chronical inhibition of Na,K-ATPase and disruption of neuronal function. Taken together, our data suggest the role of beta-amyloid as a novel physiological regulator of Na,K-ATPase.

Basal Glutathionylation of Na,K-ATPase α-Subunit Depends on Redox Status of Cells during the Enzyme Biosynthesis.

Many viruses induce oxidative stress and cause S-glutathionylation of Cys residues of the host and viral proteins. Changes in cell functioning during viral infection may be associated with glutathionylation of a number of key proteins including Na,K-ATPase which creates a gradient of sodium and potassium ions. It was found that Na,K-ATPase α-subunit has a basal glutathionylation which is not abrogated by reducing agent. We have shown that acute hypoxia leads to increase of total glutathionylation level of Na,K-ATPase α-subunit; however, basal glutathionylation of α-subunit increases under prolonged hypoxia only. The role of basal glutathionylation in Na,K-ATPase function remains unclear. Understanding significance of basal glutathionylation is complicated by the fact that there are no X-ray structures of Na,K-ATPase with the identified glutathione molecules. We have analyzed all X-ray structures of the Na,K-ATPase α-subunit from pig kidney and found that there are a number of isolated cavities with unresolved electron density close to the relevant cysteine residues. Analysis of the structures showed that this unresolved density in the structure can be occupied by glutathione associated with cysteine residues. Here, we discuss the role of basal glutathionylation of Na,K-ATPase α-subunit and provide evidence supporting the view that this modification is cotranslational.

Identification, Expression Analysis, and Target Prediction of Flax Genotroph MicroRNAs Under Normal and Nutrient Stress Conditions.

Cultivated flax (Linum usitatissimum L.) is an important plant valuable for industry. Some flax lines can undergo heritable phenotypic and genotypic changes (LIS-1 insertion being the most common) in response to nutrient stress and are called plastic lines. Offspring of plastic lines, which stably inherit the changes, are called genotrophs. MicroRNAs (miRNAs) are involved in a crucial regulatory mechanism of gene expression. They have previously been assumed to take part in nutrient stress response and can, therefore, participate in genotroph formation. In the present study, we performed high-throughput sequencing of small RNAs (sRNAs) extracted from flax plants grown under normal, phosphate deficient and nutrient excess conditions to identify miRNAs and evaluate their expression. Our analysis revealed expression of 96 conserved miRNAs from 21 families in flax. Moreover, 475 novel potential miRNAs were identified for the first time, and their targets were predicted. However, none of the identified miRNAs were transcribed from LIS-1. Expression of seven miRNAs (miR168, miR169, miR395, miR398, miR399, miR408, and lus-miR-N1) with up- or down-regulation under nutrient stress (on the basis of high-throughput sequencing data) was evaluated on extended sampling using qPCR. Reference gene search identified ETIF3H and ETIF3E genes as most suitable for this purpose. Down-regulation of novel potential lus-miR-N1 and up-regulation of conserved miR399 were revealed under the phosphate deficient conditions. In addition, the negative correlation of expression of lus-miR-N1 and its predicted target, ubiquitin-activating enzyme E1 gene, as well as, miR399 and its predicted target, ubiquitin-conjugating enzyme E2 gene, was observed. Thus, in our study, miRNAs expressed in flax plastic lines and genotrophs were identified and their expression and expression of their targets was evaluated using high-throughput sequencing and qPCR for the first time. These data provide new insights into nutrient stress response regulation in plastic flax cultivars.

Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L.

Effective fertilizer application is necessary to increase crop yields and reduce risk of plant overdosing. It is known that expression level of microRNAs (miRNAs) alters in plants under different nutrient concentrations in soil. The aim of our study was to identify and characterize miRNAs with expression alterations under excessive fertilizer in agriculturally important crop - flax (Linum usitatissimum L.). We have sequenced small RNAs in flax grown under normal and excessive fertilizer using Illumina GAIIx. Over 14 million raw reads was obtained for two small RNA libraries. 84 conserved miRNAs from 20 families were identified. Differential expression was revealed for several flax miRNAs under excessive fertilizer according to high-throughput sequencing data. For 6 miRNA families (miR395, miR169, miR408, miR399, miR398 and miR168) expression level alterations were evaluated on the extended sampling using qPCR. Statistically significant up-regulation was revealed for miR395 under excessive fertilizer. It is known that target genes of miR395 are involved in sulfate uptake and assimilation. However, according to our data alterations of the expression level of miR395 could be associated not only with excess sulfur application, but also with redundancy of other macro- and micronutrients. Furthermore expression level was evaluated for miRNAs and their predicted targets. The negative correlation between miR399 expression and expression of its predicted target ubiquitin-conjugating enzyme E2 gene was shown in flax for the first time. So we suggested miR399 involvement in phosphate regulation in L. usitatissimum. Revealed in our study expression alterations contribute to miRNA role in flax response to excessive fertilizer.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum.

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Targeting VDAC-bound hexokinase II: a promising approach for concomitant anti-cancer therapy.

INTRODUCTION: Enhancement of glucose metabolism and repression of oxidative phosphorylation followed by the Warburg effect is the common hallmark of cancer cells. Hexokinase II (HKII) plays a dual role - first, HKII up-regulation results in increased glycolysis rates. Second, association of VDAC and HKII contributes to inhibition of apoptosis through repression of the formation of mitochondrial permeability transition pores. AREAS COVERED: In this review, the role of HKII in evasion of apoptosis, aspects of HKII expression regulation, novel approaches targeting HKII and VDAC-HKII complexes and their application areas are discussed. EXPERT OPINION: The dual role of HKII in cancer cells makes it an attractive target for anti-cancer therapy. Several agents, either synthetic or plant-derived, that target hexokinase and induce VDAC-HK complex dissociation have been identified to date. Targeting hexokinase, HK-VDAC complexes as well as other glycolytic proteins not only improves the efficacy of commonly used drugs. The most prominent benefit of this approach is the ability to overcome drug resistance, for example, to cisplatin or sorafenib. In some cases, it could create an insurmountable challenge for selection of appropriate therapy. Future studies and trials should address the issue of how to transfer these approaches into clinical practice.