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

Identification and Characterization of the LecRLKs Gene Family in Maize, and Its Role Under Biotic and Abiotic Stress

Plant lectin receptor-like kinases (LecRLKs) are plant membrane protein receptor kinases. Lectin-like receptor kinases play a crucial role in regulating plant growth, development, and responses to environmental stimuli. It can rapidly respond to both biotic and abiotic stresses while mediating mechanisms of plant immune responses. This study represents the first identification of the LecRLK family genes in maize. It analyzes the gene structure, chromosomal locations, phylogenetic classification, promoter homoeotropic elements, and expression patterns under both biotic and abiotic stresses. The results indicate that these genes possess kinase and transmembrane domains, along with specific L-type or G-type extracellular domains. Most ZmLecRLK gene promoters contain cis-acting elements that are responsive to known hormones and stressors. Furthermore, these genes have been identified as being sensitive to both biotic and abiotic stresses. This discovery establishes a significant theoretical foundation for the selection of corn varieties in adverse environments. Additionally, it provides a basis for further in-depth exploration of the molecular regulatory mechanisms of LecRLK family genes.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Molecular Advances to Combat Different Biotic and Abiotic Stresses in Linseed (Linum usitatissimum L.): A Comprehensive Review

Flax, or linseed, is considered a "superfood", which means that it is a food with diverse health benefits and potentially useful bioactive ingredients. It is a multi-purpose crop that is prized for its seed oil, fibre, nutraceutical, and probiotic qualities. It is suited to various habitats and agro-ecological conditions. Numerous abiotic and biotic stressors that can either have a direct or indirect impact on plant health are experienced by flax plants as a result of changing environmental circumstances. Research on the impact of various stresses and their possible ameliorators is prompted by such expectations. By inducing the loss of specific alleles and using a limited number of selected varieties, modern breeding techniques have decreased the overall genetic variability required for climate-smart agriculture. However, gene banks have well-managed collectionns of landraces, wild linseed accessions, and auxiliary Linum species that serve as an important source of novel alleles. In the past, flax-breeding techniques were prioritised, preserving high yield with other essential traits. Applications of molecular markers in modern breeding have made it easy to identify quantitative trait loci (QTLs) for various agronomic characteristics. The genetic diversity of linseed species and the evaluation of their tolerance to abiotic stresses, including drought, salinity, heavy metal tolerance, and temperature, as well as resistance to biotic stress factors, viz., rust, wilt, powdery mildew, and alternaria blight, despite addressing various morphotypes and the value of linseed as a supplement, are the primary topics of this review.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Integrated omics approaches for flax improvement under abiotic and biotic stress: Current status and future prospects

Flax (Linum usitatissimum L.) or linseed is one of the important industrial crops grown all over the world for seed oil and fiber. Besides oil and fiber, flax offers a wide range of nutritional and therapeutic applications as a feed and food source owing to high amount of α-linolenic acid (omega-3 fatty acid), lignans, protein, minerals, and vitamins. Periodic losses caused by unpredictable environmental stresses such as drought, heat, salinity-alkalinity, and diseases pose a threat to meet the rising market demand. Furthermore, these abiotic and biotic stressors have a negative impact on biological diversity and quality of oil/fiber. Therefore, understanding the interaction of genetic and environmental factors in stress tolerance mechanism and identification of underlying genes for economically important traits is critical for flax improvement and sustainability. In recent technological era, numerous omics techniques such as genomics, transcriptomics, metabolomics, proteomics, phenomics, and ionomics have evolved. The advancements in sequencing technologies accelerated development of genomic resources which facilitated finer genetic mapping, quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection in major cereal and oilseed crops including flax. Extensive studies in the area of genomics and transcriptomics have been conducted post flax genome sequencing. Interestingly, research has been focused more for abiotic stresses tolerance compared to disease resistance in flax through transcriptomics, while the other areas of omics such as metabolomics, proteomics, ionomics, and phenomics are in the initial stages in flax and several key questions remain unanswered. Little has been explored in the integration of omic-scale data to explain complex genetic, physiological and biochemical basis of stress tolerance in flax. In this review, the current status of various omics approaches for elucidation of molecular pathways underlying abiotic and biotic stress tolerance in flax have been presented and the importance of integrated omics technologies in future research and breeding have been emphasized to ensure sustainable yield in challenging environments.

Heat-responsive microRNAs and phased small interfering RNAs in reproductive development of flax

Plants will face increased heat stress due to rising global temperatures. Heat stress affects plant reproductive development and decreases productivity; however, the underlying molecular mechanisms of these processes are poorly characterized. Plant small RNAs (sRNAs) have important regulatory roles in plant reproductive development following abiotic stress responses. We generated sRNA transcriptomes of reproductive bud stages at three different time points to identify sRNA-mediated pathways responsive to heat stress in flax (Linum usitatissimum). With added sRNA transcriptomes of vegetative tissues, we comprehensively annotated miRNA and phasiRNA-encoding genes (PHAS) in flax. We identified 173 miRNA genes, of which 42 are newly annotated. Our analysis revealed that 141 miRNA genes were differentially accumulated between tissue types, while 18 miRNA genes were differentially accumulated in reproductive tissues following heat stress, including members of miR482/2118 and miR2275 families, known triggers of reproductive phasiRNAs. Furthermore, we identified 68 21-PHAS flax loci from protein-coding and noncoding regions, four 24-PHAS loci triggered by miR2275, and 658 24-PHAS-like loci with unknown triggers, derived mostly from noncoding regions. The reproductive phasiRNAs are mostly downregulated in response to heat stress. Overall, we found that several previously unreported miRNAs and phasiRNAs are responsive to heat stress in flax reproductive tissues.

Impact TMPRSS2-ERG Molecular Subtype on Prostate Cancer Recurrence

Currently, seven molecular subtypes of prostate cancer (PCa) are known, the most common of which being the subtype characterized by the presence of the TMPRSS2-ERG fusion transcript. While there is a considerable amount of work devoted to the influence of this transcript on the prognosis of the disease, data on its role in the progression and prognosis of PCa remain controversial. The present study is devoted to the analysis of the association between the TMPRSS2-ERG transcript and the biochemical recurrence of PCa. The study included two cohorts: the RNA-Seq sample of Russian patients with PCa (n = 72) and the TCGA-PRAD data (n = 203). The results of the analysis of the association between the TMPRSS2-ERG transcript and biochemical recurrence were contradictory. The differential expression analysis (biochemical recurrence cases versus biochemical recurrence-free) and the gene set enrichment analysis revealed a list of genes involved in major cellular pathways. The GNL3, QSOX2, SSPO, and SYS1 genes were selected as predictors of the potential prognostic model (AUC = 1.000 for a cohort of Russian patients with PCa and AUC = 0.779 for a TCGA-PRAD cohort).

Local Overheating of Biotissue Labeled With Upconversion Nanoparticles Under Yb3+ Resonance Excitation

Local overheating of biotissue is a critical step for biomedical applications, such as photothermal therapy, enhancement of vascular permeability, remote control of drug release, and so on. Overheating of biological tissue when exposed to light is usually realized by utilizing the materials with a high-absorption cross section (gold, silica, carbon nanoparticles, etc.). Here, we demonstrate core/shell NaYF4:Yb3+, Tm3+/NaYF4 upconversion nanoparticles (UCNPs) commonly used for bioimaging as promising near-infrared (NIR) absorbers for local overheating of biotissue. We assume that achievable temperature of tissue labeled with nanoparticles is high enough because of Yb3+ resonance absorption of NIR radiation, whereas the use of auxiliary light-absorbing materials or shells is optional for photothermal therapy. For this purpose, a computational model of tissue heating based on the energy balance equations was developed and verified with the experimentally obtained thermal-graphic maps of a mouse in response to the 975-nm laser irradiation. Labeling of biotissue with UCNPs was found to increase the local temperature up to 2°C compared to that of the non-labeled area under the laser intensity lower than 1 W/cm2. The cellular response to the UCNP-initiated hyperthermia at subcritical ablation temperatures (lower than 42°C) was demonstrated by measuring the heat shock protein overexpression. This indicates that the absorption cross section of Yb3+ in UCNPs is relatively large, and microscopic temperature of nanoparticles exceeds the integral tissue temperature. In summary, a new approach based on the use of UCNP without any additional NIR absorbers was used to demonstrate a simple approach in the development of photoluminescent probes for simultaneous bioimaging and local hyperthermia.

Identification and Characterization of microRNAs in the Developing Seed of Linseed Flax (Linum usitatissimum L.)

Seed development plays an important role during the life cycle of plants. Linseed flax is an oil crop and the seed is a key organ for fatty acids synthesis and storage. So it is important to understand the molecular mechanism of fatty acid biosynthesis during seed development. In this study, four small RNA libraries from early seeds at 5, 10, 20 and 30 days after flowering (DAF) were constructed and used for high-throughput sequencing to identify microRNAs (miRNAs). A total of 235 miRNAs including 114 known conserved miRNAs and 121 novel miRNAs were identified. The expression patterns of these miRNAs in the four libraries were investigated by bioinformatics and quantitative real-time polymerase chain reaction (qPCR) analysis. It was found that several miRNAs, including Lus-miRNA156a was significantly correlated with seed development process. In order to confirm the actual biological function of Lus-miRNA156a, over-expression vector was constructed and transformed to Arabidopsis. The phenotypes of homozygous transgenic lines showed decreasing of oil content and most of the fatty acid content in seeds as well as late flowering time. The results provided a clue that miRNA156a participating the fatty acid biosynthesis pathway and the detailed molecular mechanism of how it regulates the pathway needs to be further investigated.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

Heredity determined by the environment: Lamarckian ideas in modern molecular biology

Inheritance of acquired characteristics (IAC) is a well-documented phenomenon occurring both in eukaryotes and prokaryotes. However, it is not included in current biological theories, and the risks of IAC induction are not assessed by genetic toxicology. Furthermore, different kinds of IAC (transgenerational and intergenerational inheritance, genotrophic changes, dauermodifications, vernalization, and some others) are traditionally considered in isolation, thus impeding the development of a comprehensive view on IAC as a whole. Herein, we discuss all currently known kinds of IAC as well as their mechanisms, if unraveled. We demonstrate that IAC is a special case of genotype × environment interactions requiring certain genotypes and, as a rule, prolonged exposure to the inducing influence. Most mechanisms of IAC are epigenetic; these include but not limited to DNA methylation, histone modifications, competition of transcription factors, induction of non-coding RNAs, inhibition of plastid translation, and curing of amyloid and non-amyloid prions. In some cases, changes in DNA sequences or host-microbe interactions are involved as well. The only principal difference between IAC and other environmentally inducible hereditary changes such as the effects of radiation is the origin of the changes: in case of IAC they are definite (determined by the environment), while the others are indefinite (arise from environmentally provoked molecular stochasticity). At least some kinds of IAC are adaptive and could be regarded as the elements of natural selection, though non-canonical in their origin and molecular nature. This is a probable way towards synthesis of the Lamarckian and Darwinian evolutionary conceptions. Applied issues of IAC are also discussed.

The cinnamyl alcohol dehydrogenase gene family is involved in the response to Fusarium oxysporum in resistant and susceptible flax genotypes

Flax (Linum usitatissimum L.) is used for the production of textile, oils, pharmaceuticals, and composite materials. Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. lini, is a very harmful disease that reduces flax production. Flax cultivars that are resistant to Fusarium wilt have been developed, and the genes that are involved in the host response to F. oxysporum have been identified. However, the mechanisms underlying resistance to this pathogen remain unclear. In the present study, we used transcriptome sequencing data obtained from susceptible and resistant flax genotypes grown under control conditions or F. oxysporum infection. Approximately 250 million reads, generated with an Illumina NextSeq instrument, were analyzed. After filtering to exclude the F. oxysporum transcriptome, the remaining reads were mapped to the L. usitatissimum genome and quantified. Then, the expression levels of cinnamyl alcohol dehydrogenase (CAD) family genes, which are known to be involved in the response to F. oxysporum, were evaluated in resistant and susceptible flax genotypes. Expression alterations in response to the pathogen were detected for all 13 examined CAD genes. The most significant differences in expression between control and infected plants were observed for CAD1B, CAD4A, CAD5A, and CAD5B, with strong upregulation of CAD1B, CAD5A, and CAD5B and strong downregulation of CAD4A. When plants were grown under the same conditions, the expression levels were similar in all studied flax genotypes for most CAD genes, and statistically significant differences in expression between resistant and susceptible genotypes were only observed for CAD1A. Our study indicates the strong involvement of CAD genes in flax response to F. oxysporum but brings no evidence of their role as resistance gene candidates. These findings contribute to the understanding of the mechanisms underlying the response of flax to F. oxysporum infection and the role of CAD genes in stress resistance.

The Role of MicroRNAs in Genome Response to Plant-Lepidoptera Interaction

RNA interference is a known phenomenon of plant immune responses, involving the regulation of gene expression. The key components triggering the silencing of targeted sequences are double-stranded RNA molecules. The regulation of host-pathogen interactions is controlled by miRNA molecules, which regulate the expression of host resistance genes or the genes of the pathogen. The review focused on basic principles of RNA interference as a gene-silencing-based defense mechanism and the role of miRNA molecules in insect genomes. RNA interference as a tool for plant protection management is discussed. The review summarizes current miRNA-based biotechnology approaches for plant protection management.

The CIMP-high phenotype is associated with energy metabolism alterations in colon adenocarcinoma

BACKGROUND

CpG island methylator phenotype (CIMP) is found in 15-20% of malignant colorectal tumors and is characterized by strong CpG hypermethylation over the genome. The molecular mechanisms of this phenomenon are not still fully understood. The development of CIMP is followed by global gene expression alterations and metabolic changes. In particular, CIMP-low colon adenocarcinoma (COAD), predominantly corresponded to consensus molecular subtype 3 (CMS3, "Metabolic") subgroup according to COAD molecular classification, is associated with elevated expression of genes participating in metabolic pathways.

METHODS

We performed bioinformatics analysis of RNA-Seq data from The Cancer Genome Atlas (TCGA) project for CIMP-high and non-CIMP COAD samples with DESeq2, clusterProfiler, and topGO R packages. Obtained results were validated on a set of fourteen COAD samples with matched morphologically normal tissues using quantitative PCR (qPCR).

RESULTS

Upregulation of multiple genes involved in glycolysis and related processes (ENO2, PFKP, HK3, PKM, ENO1, HK2, PGAM1, GAPDH, ALDOA, GPI, TPI1, and HK1) was revealed in CIMP-high tumors compared to non-CIMP ones. Most remarkably, the expression of the PKLR gene, encoding for pyruvate kinase participating in gluconeogenesis, was decreased approximately 20-fold. Up to 8-fold decrease in the expression of OGDHL gene involved in tricarboxylic acid (TCA) cycle was observed in CIMP-high tumors. Using qPCR, we confirmed the increase (4-fold) in the ENO2 expression and decrease (2-fold) in the OGDHL mRNA level on a set of COAD samples.

CONCLUSIONS

We demonstrated the association between CIMP-high status and the energy metabolism changes at the transcriptomic level in colorectal adenocarcinoma against the background of immune pathway activation. Differential methylation of at least nine CpG sites in OGDHL promoter region as well as decreased OGDHL mRNA level can potentially serve as an additional biomarker of the CIMP-high status in COAD.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

Aluminum Responsive Genes in Flax (Linum usitatissimum L.)

Flax (Linum usitatissimum L.) is a multipurpose crop which is used for the production of textile, oils, composite materials, pharmaceuticals, etc. Soil acidity results in a loss of seed and fiber production of flax, and aluminum toxicity is a major factor that depresses plant growth and development in acid conditions. In the present work, we evaluated gene expression alterations in four flax genotypes with diverse tolerance to aluminum exposure. Using RNA-Seq approach, we revealed genes that are differentially expressed under aluminum stress in resistant (Hermes, TMP1919) and sensitive (Lira, Orshanskiy) cultivars and selectively confirmed the identified alterations using qPCR. To search for differences in response to aluminum between resistant and sensitive genotypes, we developed the scoring that allowed us to suggest the involvement of MADS-box and NAC transcription factors regulating plant growth and development and enzymes participating in cell wall modifications in aluminum tolerance in flax. Using Gene Ontology (GO) enrichment analysis, we revealed that glutathione metabolism, oxidoreductase, and transmembrane transporter activities are the most affected by the studied stress in flax. Thus, we identified genes that are involved in aluminum response in resistant and sensitive genotypes and suggested genes that contribute to flax tolerance to the aluminum stress.

Characterization of repeated DNA sequences in genomes of blue-flowered flax

BACKGROUND

Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method.

RESULTS

It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes.

CONCLUSIONS

The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.

Intrusive Growth of Phloem Fibers in Flax Stem: Integrated Analysis of miRNA and mRNA Expression Profiles

Phloem fibers are important elements of plant architecture and the target product of many fiber crops. A key stage in fiber development is intrusive elongation, the mechanisms of which are largely unknown. Integrated analysis of miRNA and mRNA expression profiles in intrusivelygrowing fibers obtained by laser microdissection from flax (Linum usitatissimum L.) stem revealed all 124 known flax miRNA from 23 gene families and the potential targets of differentially expressed miRNAs. A comparison of the expression between phloem fibers at different developmental stages, and parenchyma and xylem tissues demonstrated that members of miR159, miR166, miR167, miR319, miR396 families were down-regulated in intrusively growing fibers. Some putative target genes of these miRNA families, such as those putatively encoding growth-regulating factors, an argonaute family protein, and a homeobox-leucine zipper family protein were up-regulated in elongating fibers. miR160, miR169, miR390, and miR394 showed increased expression. Changes in the expression levels of miRNAs and their target genes did not match expectations for the majority of predicted target genes. Taken together, poorly understood intrusive fiber elongation, the key process of phloem fiber development, was characterized from a miRNA-target point of view, giving new insights into its regulation.

Flax (Linum usitatissimum L.) response to non-optimal soil acidity and zinc deficiency

BACKGROUND

Flax (Linum usitatissimum L.) is grown for fiber and seed production. Unfavorable environments, such as nutrient deficiency and non-optimal soil acidity, decrease the quantity and quality of yield. Cultivation of tolerant to stress varieties can significantly reduce the crop losses. Understanding the mechanisms of flax response to the stresses and identification of resistance gene candidates will help in breeding of improved cultivars. In the present work, the response of flax plants to increased pH level and zinc (Zn) deficiency was studied.

RESULTS

We performed high-throughput transcriptome sequencing of two flax cultivars with diverse tolerance to increased pH level and Zn deficiency: Norlin (tolerant) and Mogilevsky (sensitive). Sixteen cDNA libraries were created from flax plants grown under control conditions, increased pH level, Zn deficiency, and both stresses simultaneously, and about 35 million reads were obtained for each experiment type. Unfavorable pH resulted in significantly stronger gene expression alterations compared to Zn deficiency. Ion homeostasis, oxidoreductase activity, cell wall, and response to stress Gene Ontology terms were the most affected by unfavorable pH and Zn deficiency both in tolerant and sensitive flax cultivars. Upregulation of genes encoding metal transporters was identified under increased pH level, Zn deficiency, and both stresses simultaneously. Under Zn deficiency, only in tolerant cultivar Norlin, we revealed the induction of several photosynthesis-related genes and, in this way, this tolerant genotype could overcome unfavorable effects of reduced Zn content.

CONCLUSIONS

We identified genes with expression alterations in flax under non-optimal soil acidity and Zn deficiency based on high-throughput sequencing data. These genes are involved in diverse processes, including ion transport, cell wall biogenesis, and photosynthesis, and could play an important role in flax response to the studied stresses. Moreover, genes with distinct expression changes between examined tolerant and sensitive genotypes could determine the mechanisms of flax tolerance to non-optimal soil acidity and Zn deficiency.

Deep Sequencing Revealed a CpG Methylation Pattern Associated With ALDH1L1 Suppression in Breast Cancer

Hypermethylation of promoter CpG islands is generally recognized epigenetic mechanism responsible for gene silencing in cancer. However, molecular details on how this epigenetic mark triggers the process of gene downregulation are still elusive. Here, we used deep bisulfite sequencing and qPCR analysis to investigate the pattern of CpG methylation of ALDH1L1 promoter region and its association with the gene expression level in 16 paired breast cancer (BC) samples of different clinical stages. Expression of ALDH1L1 gene was suppressed in all examined BC samples up to 200-fold, and average hypermethylation level of the promoter region correlated positively with ALDH1L1 downregulation. We determined the role of every individual CpG site within the ALDH1L1 promoter, including upstream untranscribed region, first untranslated exon, and the start of the first intron, in aberrant gene expression by correlation analysis. The search revealed CpG sites which methylation has the highest impact on intensity of gene transcription. The majority of such CpG sites are located in a compact region in the first intron of the ALDH1L1 gene. These results assist in unraveling of dynamic nature of CpG promoter hypermethylation as well as demonstrate the efficiency of deep bisulfite sequencing in search for novel epigenetic markers in cancer.

HK3 overexpression associated with epithelial-mesenchymal transition in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is a common cancer worldwide. The main cause of death in CRC includes tumor progression and metastasis. At molecular level, these processes may be triggered by epithelial-mesenchymal transition (EMT) and necessitates specific alterations in cell metabolism. Although several EMT-related metabolic changes have been described in CRC, the mechanism is still poorly understood.

RESULTS

Using CrossHub software, we analyzed RNA-Seq expression profile data of CRC derived from The Cancer Genome Atlas (TCGA) project. Correlation analysis between the change in the expression of genes involved in glycolysis and EMT was performed. We obtained the set of genes with significant correlation coefficients, which included 21 EMT-related genes and a single glycolytic gene, HK3. The mRNA level of these genes was measured in 78 paired colorectal cancer samples by quantitative polymerase chain reaction (qPCR). Upregulation of HK3 and deregulation of 11 genes (COL1A1, TWIST1, NFATC1, GLIPR2, SFPR1, FLNA, GREM1, SFRP2, ZEB2, SPP1, and RARRES1) involved in EMT were found. The results of correlation study showed that the expression of HK3 demonstrated a strong correlation with 7 of the 21 examined genes (ZEB2, GREM1, TGFB3, TGFB1, SNAI2, TWIST1, and COL1A1) in CRC.

CONCLUSIONS

Upregulation of HK3 is associated with EMT in CRC and may be a crucial metabolic adaptation for rapid proliferation, survival, and metastases of CRC cells.

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.

Upregulation of NETO2 gene in colorectal cancer

BACKGROUND

Neuropilin and tolloid-like 2 (NETO2) is a single-pass transmembrane protein that has been shown primarily implicated in neuron-specific processes. Upregulation of NETO2 gene was also detected in several cancer types. In colorectal cancer (CRC), it was associated with tumor progression, invasion, and metastasis, and seems to be involved in epithelial-mesenchymal transition (EMT). However, the mechanism of NETO2 action is still poorly understood.

RESULTS

We have revealed significant increase in the expression of NETO2 gene and deregulation of eight EMT-related genes in CRC. Four of them were upregulated (TWIST1, SNAIL1, LEF1, and FOXA2); the mRNA levels of other genes (FOXA1, BMP2, BMP5, and SMAD7) were decreased. Expression of NETO2 gene was weakly correlated with that of genes involved in the EMT process.

CONCLUSIONS

We found considerable NETO2 upregulation, but no significant correlation between the expression of NETO2 and EMT-related genes in CRC. Thus, NETO2 may be involved in CRC progression, but is not directly associated with EMT.

miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.)

Acid soils limit agricultural production worldwide. Major reason of crop losses in acid soils is the toxicity of aluminum (Al). In the present work, we investigated expression alterations of microRNAs in flax (Linum usitatissimum L.) plants under Al stress. Flax seedlings of resistant (TMP1919 and G1071/4_k) and sensitive (Lira and G1071/4_o) to Al cultivars and lines were exposed to AlCl₃ solution for 4 and 24 hours. Twelve small RNA libraries were constructed and sequenced using Illumina platform. In total, 97 microRNAs from 18 conserved families were identified. miR319, miR390, and miR393 revealed expression alterations associated with Al treatment of flax plants. Moreover, for miR390 and miR393, the alterations were distinct in sensitive and resistant to Al genotypes. Expression level changes of miR319 and miR390 were confirmed using qPCR analysis. In flax, potential targets of miR319 are TCPs, miR390–TAS3 and GRF5, and miR393–AFB2-coding transcripts. TCPs, TAS3, GRF5, and AFB2 participate in regulation of plant growth and development. The involvement of miR319, miR390, and miR393 in response to Al stress in flax was shown here for the first time. We speculate that these microRNAs play an important role in Al response via regulation of growth processes in flax plants.

Involvement of Small RNAs in Phosphorus and Sulfur Sensing, Signaling and Stress: Current Update

Plants require several essential mineral nutrients for their growth and development. These nutrients are required to maintain physiological processes and structural integrity in plants. The root architecture has evolved to absorb nutrients from soil and transport them to other parts of the plant. Nutrient deficiency affects several physiological and biological processes in plants and leads to reduction in crop productivity and yield. To compensate this adversity, plants have developed adaptive mechanisms to enhance the acquisition, conservation, and mobilization of these nutrients under deficient or adverse conditions. In addition, plants have evolved an intricate nexus of complex signaling cascades, which help in nutrient sensing and uptake as well as to maintain nutrient homeostasis. In recent years, small non-coding RNAs such as micro RNAs (miRNAs) and endogenous small interfering RNAs have emerged as important component in regulating plant stress responses. A set of these small RNAs (sRNAs) have been implicated in regulating various processes involved in nutrient uptake, assimilation, and deficiency. In response to phosphorus (P) and sulphur (S) deficiencies, role of sRNAs, miR395 and miR399, have been identified to be instrumental; however, many more miRNAs might be involved in regulating the plant response to these nutrient stresses. These sRNAs modulate expression of target genes in response to P and S deficiencies and regulate their uptake and utilization for proper growth and development of the plant. This review summarizes the current understanding of uptake, sensing, and signaling of P and S and highlights the regulatory role of sRNAs in adaptive responses to these nutrient stresses in plants.

Glutathione S-transferases and UDP-glycosyltransferases Are Involved in Response to Aluminum Stress in Flax

About 30% of the world's ice-free land area is occupied by acid soils. In soils with pH below 5, aluminum (Al) releases to the soil solution, and becomes highly toxic for plants. Therefore, breeding of varieties that are resistant to Al is needed. Flax (Linum usitatissimum L.) is grown worldwide for fiber and seed production. Al toxicity in acid soils is a serious problem for flax cultivation. However, very little is known about mechanisms of flax resistance to Al and the genetics of this resistance. In the present work, we sequenced 16 transcriptomes of flax cultivars resistant (Hermes and TMP1919) and sensitive (Lira and Orshanskiy) to Al, which were exposed to control conditions and aluminum treatment for 4, 12, and 24 h. In total, 44.9-63.3 million paired-end 100-nucleotide reads were generated for each sequencing library. Based on the obtained high-throughput sequencing data, genes with differential expression under aluminum exposure were revealed in flax. The majority of the top 50 up-regulated genes were involved in transmembrane transport and transporter activity in both the Al-resistant and Al-sensitive cultivars. However, genes encoding proteins with glutathione transferase and UDP-glycosyltransferase activity were in the top 50 up-regulated genes only in the flax cultivars resistant to aluminum. For qPCR analysis in extended sampling, two UDP-glycosyltransferases (UGTs), and three glutathione S-transferases (GSTs) were selected. The general trend of alterations in the expression of the examined genes was the up-regulation under Al stress, especially after 4 h of Al exposure. Moreover, in the flax cultivars resistant to aluminum, the increase in expression was more pronounced than that in the sensitive cultivars. We speculate that the defense against the Al toxicity via GST antioxidant activity is the probable mechanism of the response of flax plants to aluminum stress. We also suggest that UGTs could be involved in cell wall modification and protection from reactive oxygen species (ROS) in response to Al stress in L. usitatissimum. Thus, GSTs and UGTs, probably, play an important role in the response of flax to Al via detoxification of ROS and cell wall modification.

Gene expression profiling of flax (Linum usitatissimum L.) under edaphic stress

BACKGROUND

Cultivated flax (Linum usitatissimum L.) is widely used for production of textile, food, chemical and pharmaceutical products. However, various stresses decrease flax production. Search for genes, which are involved in stress response, is necessary for breeding of adaptive cultivars. Imbalanced concentration of nutrient elements in soil decrease flax yields and also results in heritable changes in some flax lines. The appearance of Linum Insertion Sequence 1 (LIS-1) is the most studied modification. However, LIS-1 function is still unclear.

RESULTS

High-throughput sequencing of transcriptome of flax plants grown under normal (N), phosphate deficient (P), and nutrient excess (NPK) conditions was carried out using Illumina platform. The assembly of transcriptome was performed, and a total of 34924, 33797, and 33698 unique transcripts for N, P, and NPK sequencing libraries were identified, respectively. We have not revealed any LIS-1 derived mRNA in our sequencing data. The analysis of high-throughput sequencing data allowed us to identify genes with potentially differential expression under imbalanced nutrition. For further investigation with qPCR, 15 genes were chosen and their expression levels were evaluated in the extended sampling of 31 flax plants. Significant expression alterations were revealed for genes encoding WRKY and JAZ protein families under P and NPK conditions. Moreover, the alterations of WRKY family genes differed depending on LIS-1 presence in flax plant genome. Besides, we revealed slight and LIS-1 independent mRNA level changes of KRP2 and ING1 genes, which are adjacent to LIS-1, under nutrition stress.

CONCLUSIONS

Differentially expressed genes were identified in flax plants, which were grown under phosphate deficiency and excess nutrition, on the basis of high-throughput sequencing and qPCR data. We showed that WRKY and JAS gene families participate in flax response to imbalanced nutrient content in soil. Besides, we have not identified any mRNA, which could be derived from LIS-1, in our transcriptome sequencing data. Expression of LIS-1 flanking genes, ING1 and KRP2, was suggested not to be nutrient stress-induced. Obtained results provide new insights into edaphic stress response in flax and the role of LIS-1 in these process.