Genome-wide identification of low phosphorus responsive microRNAs in two soybean genotypes by high-throughput sequencing.

MicroRNAs (miRNAs) have been reported to be correlated with various stress responses in soybean, but only a few miRNAs have been demonstrated to respond to low phosphorus (LP) stress. To unravel the response mechanisms of miRNAs to low-P stress, the roots of two representative soybean genotypes with different P efficiency, Nannong94-156 (a LP-tolerant genotype) and Bogao (a LP-sensitive genotype), were used for the construction of RNA sequencing (RNA-seq) libraries under low/normal-P treatment by high-throughput sequencing. In total, 603 existing miRNAs and 1699 novel miRNAs belonging to 248 and 1582 families in all samples were identified, respectively. Among these miRNAs, 777 miRNAs were differentially expressed (DE) across different P levels and genotypes. Furthermore, putative targets of DE miRNAs were predicted, and these miRNAs mainly targeted ERF (ethylene responsive factor), auxin response factors (ARF), zinc finger protein, MYB, and NAC domain transcription factors. Gene ontology (GO) analysis showed that targets of DE miRNAs were significantly enriched in binding, metabolic processes, biological regulation, response to stress, and phosphorus metabolic processes. In addition, the expression profiles of chosen P-responsive miRNAs and target genes were validated by quantitative real-time PCR (qRT-PCR). Our study focused on genome-wide miRNA identification in two representative soybean genotypes under low-P stress. Overall, the DE miRNAs across different P levels and genotypes and their putative target genes will provide useful information for further study of miRNAs mediating low-P response and facilitate improvements in soybean breeding.

Al-Tolerant Barley Mutant hvatr.g Shows the ATR-Regulated DNA Damage Response to Maleic Acid Hydrazide.

ATR, a DNA damage signaling kinase, is required for cell cycle checkpoint regulation and detecting DNA damage caused by genotoxic factors including Al3+ ions. We analyzed the function of the HvATR gene in response to chemical clastogen-maleic acid hydrazide (MH). For this purpose, the Al-tolerant barley TILLING mutant hvatr.g was used. We described the effects of MH on the nuclear genome of hvatr.g mutant and its WT parent cv. "Sebastian", showing that the genotoxic effect measured by TUNEL test and frequency of cells with micronuclei was much stronger in hvatr.g than in WT. MH caused a significant decrease in the mitotic activity of root cells in both genotypes, however this effect was significantly stronger in "Sebastian". The impact of MH on the roots cell cycle, analyzed using flow cytometry, showed no differences between the mutant and WT.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings.

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

Assessment of the antioxidant, cytotoxic, and genotoxic potential of the Annona muricata leaves and their influence on genomic stability.

The popular use of Annona muricata L. is based upon a range of medicinal purposes, and the plant exhibits biological activities including antihyperglycemic, antiparasitic, and antitumor activities. The objectives of this study were to examine the antioxidant, cytotoxic, and genotoxic potential of the hydroalcoholic extract of A. muricata leaves (AMEs), as well as its effects on genotoxicity induced by methyl methanesulfonate (MMS) and hydrogen peroxide (H2O2). The results using 2,2-diphenyl-1-picrylhydrazyl assay showed that AME was able to scavenge 44.71% of free radicals. The extract significantly reduced the viability of V79 cells in the clonogenic assay at concentrations ≥8 µg/ml. No significant differences in micronucleus (MN) frequency were observed between V79 cell cultures treated with different concentrations of the extract (0.125, 0.25, 0.5, and 1 µg/ml) and negative control. When AME concentrations were combined with MMS, data revealed no marked differences from mutagen alone. In contrast, significant reductions in the frequencies of MN were noted in cultures treated with AME combined with H2O2 compared to H2O2 alone. In vivo studies found no significant differences in the frequencies of micronucleated polychromatic erythrocytes (MNPCEs) between animals treated with different AME doses compared to control. Animals treated with AME doses of 125 and 250 mg/kg and MMS exhibited significantly higher frequencies of MNPCE compared to mutagen alone. In conclusion, under current experimental conditions, AME was not genotoxic and exerted a modulatory effect on DNA damage depending upon the experimental conditions. The extract did not influence markedly MMS-induced genotoxicity in in vitro test system. However, the extract increased DNA damage induced by mutagen in mice. In V79 cells, AME reduced the genotoxicity produced by H2O2, and this protective effect was attributed in part to the antioxidant activity of AME.

DNA METHYLATION ANALYSIS DURING THE OPTIMIZATION OF Agrobacterium-MEDIATED TRANSFORMATION OF SOYBEAN.

Soybean is recognized as one of the plants which are very difficult to be transformed. Considering the low transformation efficiency of soybean, we aimed to determine the effect of 6-benzylaminopurine (6-BA), shoot induction time, and infection time of Agrobacterium on the clonal propagation of Glycine max. Results showed that 1.6 mg/L 6-BA could be optimal to promote the induction of adventitious shoots. An induction time of 15 d was considered optimal for the actual experiment involving soybean shoot induction. Agrobacterium was cultured until an OD600 = 0.8 was reached for an infection time of 30 min; this infection time may be optimal to promote soybean transformation. Whole genome DNA methylation was analyzed by high-performance liquid chromatography (HPLC)-assisted quantification, and DNA methylation result is consistent with the phenotypic data of shoot development. In addition, two methylation-related genes (Decrease in DNA methylation 1 and DNA methyltransferases chromomethylase 2) were analyzed to determine expression differences by qRT-PCR in the shoots that were developed under different experimental conditions. In general, the expression values of these genes were normally downregulated under the recommended experimental conditions of soybean regeneration. This study showed the overall methylation changes in the in vitro culture of soybean, as affected by several variable parameters, which is useful to promote the transformation efficiency of soybean.

Measuring the damage of heavy metal cadmium in rice seedlings by SRAP analysis combined with physiological and biochemical parameters.

BACKGROUND: Cadmium (Cd) is one of the most poisonous pollutants, and Cd pollution has become the limiting factor of rice production and quality improvement. Therefore it is of significant importance to monitor Cd toxicity by the detection of Cd contamination in rice with biomarkers. In the present study, sequence-related amplified polymorphism (SRAP) and physiological and biochemical methods were applied to determine the toxicological effects of Cd stress on rice. RESULTS: With increasing Cd concentration and duration, the content of chlorophyll in the two rice varieties W7 and M63 decreased and that of malondialdehyde increased. This tendency was more apparent in M63. The antioxidant enzymes superoxide dismutase and peroxidase both increased significantly compared with controls. SRAP polymerase chain reaction results indicated significant differences between Cd treatments and controls in terms of SRAP profile, as well as genotypic differences. The genomic template stability (GTS) decreased with increasing Cd concentration and duration. Under the same treatment conditions, the GTS of W7 was higher than that of M63. Comparison analysis revealed that the changes in physiological and biochemical parameters of rice seedlings under Cd stress had a good correlation with the changes in SRAP profile. Furthermore, the changes in SRAP profile showed enhanced sensitivity in the roots of rice seedlings. CONCLUSION: The SRAP profile and physiological and biochemical parameters could act as appropriate biomarkers for the measurement of Cd contamination during rice production.

Detection of genome DNA methylation change in spinach induced by 5-azaC.

DNA methylation has been implicated in the regulation of gene expression, genome imprinting, and chromatin remodeling in eukaryotes. In this study, we analyzed possible alterations in levels and patterns of cytosine methylation in male and female spinach plants after treatment with demethylation agent 5-azacytidine (5-azaC) using two methods: (1) direct determination of 5-methylcytidine (5 mC) amounts in genomic DNA by high-performance liquid chromatography (HPLC) separation and quantification of nucleosides and (2) methylation-sensitive inter-simple sequence repeat (MS-ISSR) technique. HPLC analysis revealed that the DNA methylation events in male and female spinach leaves markedly decreased upon 30 µM 5-azaC treatment, and the methylation level gradually decreased with the increase in 5-azaC concentration. To study the altered DNA methylation patterns in spinach after 5-azaC treatment, untreated and 500 µM 5-azaC-treated samples were analyzed by MS-ISSR assay. A total of 385 informative profiles were resolved using 35 ISSR primer sets. MS-ISSR analysis showed various altered methylation patterns between untreated and 5-azaC-treated spinach plants. These alterations were mainly demethylation events, which were largely consistent with the HPLC results. Both HPLC and MS-ISSR analyses showed that the changes in DNA methylation levels and patterns were similar in male and female spinach leaves, which implies that sex was not the main factor influencing DNA methylation levels and patterns in the vegetative organs of spinach. This study could provide a molecular basis of the altered DNA methylation induced by 5-azaC, and lay a foundation for further investigation of the relationship between methylation and sex determination and development in this dioecious plant spinach.

Benefit-of-doubt (BOD) scoring: a sequencing-based method for SNP candidate assessment from high to medium read number data sets.

Identification of single nucleotide polymorphisms (SNPs) is a key element in sequence-based genetic analysis. Next generation sequencing offers a cost-effective basis to generate the necessary, large sequence data sets, and bioinformatic methods are being developed to process sequencing machine readouts. We were interested in detection of SNPs in a 350 kb region of an EMS-mutagenized Arabidopsis chromosome 3. The region was selectively analyzed using PCR-generated, overlapping fragments for Solexa sequencing. The ensuing reads provided a high coverage and were processed bioinformatically. In order to assess the SNP candidates obtained with a frequently used alignment program and SNP caller, we developed an additional method that allows the identification of high confidence SNP loci. The method can easily be applied to complete genome sequence data of sufficient coverage.

A suite of new genes defining salinity stress tolerance in seedlings of contrasting rice genotypes.

Salinity is one of the major constraints adversely influencing crop productivity. Saltol QTL is a major QTL associated with Na⁺-K⁺ ratio and seedling stage salinity tolerance in rice. With an aim to understand the contribution of individual genes localized within saltol towards salinity tolerance, we analysed the transcript abundance of a set of these genes in seedlings of contrasting genotypes of rice. We hypothesize that this approach may be helpful in identifying new 'candidate genes' for improving salinity tolerance in crops. For this purpose, seedlings of Oryza sativa cv. IR64 (sensitive) and the landrace Pokkali (tolerant) were subjected to short/long durations of salinity. qRT-PCR analysis clearly exhibited differential regulation of genes encoding signaling related protein (SRPs), where higher transcript abundance for most of them was observed in Pokkali than IR64 under non-stress conditions, thereby indicating towards well preparedness of the former to handle stress, in anticipation. Genes encoding proteins of unknown function (PUFs), though, constitute a considerable portion of plant genome, have so far been neglected in most studies. Time course analysis of these genes showed a continuous increase in their abundance in Pokkali, while in IR64, their abundance increased till 24 h followed by a clear decrease, thereby justifying their nomenclature as 'salinity induced factors' (SIFs). This is the first report showing possible involvement of SIFs localized within salinity related QTL towards salinity stress response. Based on the phenotypes of insertional mutants, it is proposed that these SIFs may have a putative function in vegetative growth (SIFVG), fertility (SIFF), viability (SIFV) or early flowering (SIFEF).

Genome-wide identification of Brassica napus microRNAs and their targets in response to cadmium.

MicroRNAs (miRNAs) are a distinct class of small RNAs in plants that not only regulate biological processes but also regulate response to environmental stresses. The toxic heavy metal cadmium (Cd) induces expression of several miRNAs in rapeseed (Brassica napus), but it is not known on a genome-wide scale how the expression of miRNAs and their target genes, is regulated by Cd. In this study, four small RNA libraries and four degradome libraries were constructed from Cd-treated and non-Cd-treated roots and shoots of B. napus seedlings. Using high-throughput sequencing, the study identified 84 conserved and non-conserved miRNAs (belonging to 37 miRNA families) from Cd-treated and non-treated B. napus, including 19 miRNA members that were not identified before. Some of the miRNAs were validated by RNA gel blotting. Most of the identified miRNAs were found to be differentially expressed in roots/shoots or regulated by Cd exposure. The study simultaneously identified 802 targets for the 37 (24 conserved and 13 non-conserved) miRNA families, from which there are 200, 537, and 65 targets, belonging to categories I, II, and III, respectively. In category I alone, many novel targets for miRNAs were identified and shown to be involved in plant response to Cd.

Identification of EMS-induced causal mutations in a non-reference Arabidopsis thaliana accession by whole genome sequencing.

The most frequently used method to identify mutations induced by a commonly used mutagen, EMS (ethyl methane sulfonate), in Arabidopsis thaliana has been map-based cloning. The first step of this method is crossing a mutant with a plant of another accession as it requires polymorphisms between accessions for linkage analysis. Therefore, to perform the method routinely, it is greatly preferred to use accession combinations between which enough polymorphisms are already known. Further, it requires laborious examination of a large number of F2 recombinants using many markers to detect each polymorphism. After linkage analysis narrows down the chromosomal region containing the causal mutation, sequencing candidate genes one by one within the region is necessary until the mutation is finally identified. Overall, this method is generally time-consuming and labor intensive, and it becomes harder when multiple loci are involved in phenotypes. A few recent reports showed that causal mutations induced by EMS could be identified by deep-sequencing technologies with less labor compared with the conventional method when mutants were generated in the Arabidopsis reference Columbia background whose genome organization is well known. Here we report that we succeeded in rapid identification of EMS-induced causal mutations in a non-reference accession background, whose whole genome sequence is not publicly available, using one round of whole genome sequencing. Moreover, in our case, we could monitor the causal locus and the transgenic reporter locus simultaneously, implying that this methodology could theoretically be applicable to analyzing even complex traits. We describe the pipeline of this methodology and discuss its characteristics.

Genome-wide analysis of primary auxin-responsive Aux/IAA gene family in maize (Zea mays. L.).

The phytohormone auxin is important in various aspects of organism growth and development. Aux/IAA genes encoding short-lived nuclear proteins are responsive primarily to auxin induction. Despite their physiological importance, systematic analysis of Aux/IAA genes in maize have not yet been reported. In this paper, we presented the isolation and characterization of maize Aux/IAA genes in whole-genome scale. A total of 31 maize Aux/IAA genes (ZmIAA1 to ZmIAA31) were identified. ZmIAA genes are distributed in all the maize chromosomes except chromosome 2. Aux/IAA genes expand in the maize genome partly due to tandem and segmental duplication events. Multiple alignment and motif display results revealed major maize Aux/IAA proteins share all the four conserved domains. Phylogenetic analysis indicated Aux/IAA family can be divided into seven subfamilies. Putative cis-acting regulatory DNA elements involved in auxin response, light signaling transduction and abiotic stress adaption were observed in the promoters of ZmIAA genes. Expression data mining suggested maize Aux/IAA genes have temporal and spatial expression pattern. Collectively, these results will provide molecular insights into the auxin metabolism, transport and signaling research.

Structural analysis of 83-kb genomic DNA from Thellungiella halophila: sequence features and microcolinearity between salt cress and Arabidopsis thaliana.

Salt cress (Thellungiella halophila) has become a desirable plant model for molecular analysis of the mechanisms of salt tolerance. Analysis of its physiological action and expressed EST has resulted in better understanding. However, less is known about its genomic features. Here we determined a continuous sequence approximately 83 kb from a salt cress BAC clone, providing the first insight into the genomic feature for this species. The gene density is approximately one gene per 3.6 kb in this sequence. Many types of repetitive sequences are present in this salt cress sequence, including LTR retroelements, DNA transposons and a number of simple sequence repeats. Comparison of sequence similarity indicated that salt cress shares a close relationship with Arabidopsis. Extensive conservation and high-level microcolinearity were uncovered for both genomes. Our study also indicated that genomic DNA alternations (involving chromosome inversion, sequence loss and gene translocation) contributed to the genomic discrepancies between salt cress and Arabidopsis.

Increased tolerance to organic xenobiotics following recent allopolyploidy in Spartina (Poaceae).

Genome doubling or polyploidy is a widespread phenomenon in plants where it has important evolutionary consequences affecting the species distribution and ecology. PAHs are ubiquitous organic pollutants, which represent a major environmental concern. Recent data showed that tolerance to organic xenobiotics involve specific signaling pathways, and detoxifying gene sets referred as 'the xenome'. However, no data are available about how polyploidy impacts tolerance to organic xenobiotics. In the present paper, we investigated PAH tolerance following allopolyploidization in Spartina alterniflora, S. maritima and their derived allopolyploid species S. anglica. We performed comparative analyses of cellular compartmentalization, photosynthetic indices, and oxidative stress markers under phenanthrene-induced stress, and found that S. anglica exhibit increased tolerance compared to its parents. Based on 52 genes potentially involved in phenanthrene detoxification previously identified in A. thaliana, we investigated the Spartina xenome using genomic and transcriptomic available resources. Subsequently, we focused on GSTs, a ubiquitous enzymes class involved in organic xenobiotic detoxification. We examined expression profiles of selected genes by RT-qPCR, and revealed various patterns of parental expression alteration in the allopolyploid. The impacts of allopolyploidization on phenanthrene-induced stress and their potential ecological implications are discussed. The neo-allopolyploid S. anglica appears as a potential candidate for phytoremediation in PAH-polluted marshes.

Genome-wide transcriptional profiling for elucidating the effects of brassinosteroids on Glycine max during early vegetative development.

Soybean is a widely grown grain legume and one of the most important economic crop species. Brassinosteroids play a crucial role in plant vegetative growth and reproductive development. However, it remains unclear how BRs regulate the developmental processes in soybean, and the molecular mechanism underlying soybean early development is largely unexplored. In this study, we first characterized how soybean early vegetative growth was specifically regulated by the BR biosynthesis inhibitor propiconazole; this characterization included shortened root and shoot lengths, reduced leaf area, and decreased chlorophyll content. In addition, the growth inhibition induced by Pcz could be rescued by exogenous brassinolide application. The RNA-seq technique was employed to investigate the BR regulatory networks during soybean early vegetative development. Identification and analysis of differentially expressed genes indicated that BRs orchestrate a wide range of cellular activities and biological processes in soybean under various BR concentrations. The regulatory networks between BRs and multiple hormones or stress-related pathways were investigated. The results provide a comprehensive view of the physiological functions of BRs and new insights into the molecular mechanisms at the transcriptional level of BR regulation of soybean early development.

Development of a Pedigreed Sorghum Mutant Library.

Induced mutagenesis is a powerful approach to generate variations for elucidation of gene function and to create new traits for breeding. Here, we described a procedure to develop a pedigreed mutant library through chemical mutagenesis with ethylmethane sulfonate (EMS) treated seeds in sorghum and discussed its potential to generate new traits for sorghum improvement. Unlike random mutagenesis, a pedigreed mutant library, once properly established, can serve as a powerful resource to isolate and recover mutations of both agronomical and biological importance. With the development of affordable and high-throughput next-generation sequencing technologies, identification of causal mutations from a mutant library with a uniform genetic background becomes increasingly efficient and cost-effective. Fast causal gene discovery from mutant libraries combined with precise genome editing techniques will accelerate incorporation of new traits and revolutionize crop breeding.

Rice-arsenate interactions in hydroponics: whole genome transcriptional analysis.

Rice (Oryza sativa) varieties that are arsenate-tolerant (Bala) and -sensitive (Azucena) were used to conduct a transcriptome analysis of the response of rice seedlings to sodium arsenate (AsV) in hydroponic solution. RNA extracted from the roots of three replicate experiments of plants grown for 1 week in phosphate-free nutrient with or without 13.3 muM AsV was used to challenge the Affymetrix (52K) GeneChip Rice Genome array. A total of 576 probe sets were significantly up-regulated at least 2-fold in both varieties, whereas 622 were down-regulated. Ontological classification is presented. As expected, a large number of transcription factors, stress proteins, and transporters demonstrated differential expression. Striking is the lack of response of classic oxidative stress-responsive genes or phytochelatin synthases/synthatases. However, the large number of responses from genes involved in glutathione synthesis, metabolism, and transport suggests that glutathione conjugation and arsenate methylation may be important biochemical responses to arsenate challenge. In this report, no attempt is made to dissect differences in the response of the tolerant and sensitive variety, but analysis in a companion article will link gene expression to the known tolerance loci available in the BalaxAzucena mapping population.

Analysis of aluminum toxicity in Hordeum vulgare roots with an emphasis on DNA integrity and cell cycle.

Barley is one of the cereals that are most sensitive to aluminum (Al). Al in acid soils limits barley growth and development and, as a result, its productivity. The inhibition of root growth is a widely accepted indicator of Al stress. Al toxicity is affected by many factors including the culture medium, pH, Al concentration and the duration of the treatment. However, Al can act differently in different species and still Al toxicity in barley deserves study. Since the mechanism of Al toxicity is discussed we cytogenetically describe the effects of different doses of bioavailable Al on the barley nuclear genome-mitotic activity, cell cycle profile and DNA integrity. At the same time, we tested an established deep-water culture (DWC) hydroponics system and analyzed the effects of Al on the root system parameters using WinRHIZO software. We demonstrated the cytotoxic and genotoxic effect of Al in barley root cells. We showed that Al treatment significantly reduced the mitotic activity of the root tip cells and it also induced micronuclei and damaged nuclei. The DNA-damaging effect of Al was observed using the TUNEL test. We define the inhibitory influence of Al on DNA replication in barley. Analysis with the labelling and detection of 5-ethynyl-2'-deoxyuridin (EdU) showed that the treatment with Al significantly decreased the frequency of S phase cells. We also demonstrated that Al exposure led to changes in the cell cycle profile of barley root tips. The delay of cell divisions observed as increased frequency of cells in G2/M phase after Al treatment was reported using flow cytometry.

Expression of EuNOD-ARP1 encoding auxin-repressed protein homolog is upregulated by auxin and localized to the fixation zone in root nodules of Elaeagnus umbellata.

Root nodule formation is controlled by plant hormones such as auxin. Auxin-repressed protein (ARP) genes have been identified in various plant species but their functions are not clear. We have isolated a full-length cDNA clone (EuNOD-ARP1) showing high sequence homology to previously identified ARP genes from root nodules of Elaeagnus umbellata. Genomic Southern hybridization showed that there are at least four ARP-related genes in the genome of E. umbellata. The cDNA clone encodes a polypeptide of 120 amino acid residues with no signal peptide or organelle-targeting signals, indicating that it is a cytosolic protein. Its cytosolic location was confirmed using Arabidopsis protoplasts expressing a EuNOD-ARP1:smGFP fusion protein. Northern hybridization showed that EuNOD-ARP1 expression was higher in root nodules than in leaves or uninoculated roots. Unlike the ARP genes of strawberry and black locust, which are negatively regulated by exogenous auxin, EuNOD-ARP1 expression is induced by auxin in leaf tissue of E. umbellata. In situ hybridization revealed that EuNOD-ARP1 is mainly expressed in the fixation zone of root nodules.

Mapping and Cloning of Chemical Induced Mutations by Whole-Genome Sequencing of Bulked Segregants.

Here we describe a method of gene mapping and cloning of chemical induced mutations using next-generation sequencing on their backcrossed segregants. This method utilizes polymorphisms induced by chemical mutagens for mapping and therefore does not require outcrossing to a different background or a large segregating population. It can be used to clone causal mutations generated in various accessions or mutant backgrounds for dissecting complex processes such as pattern recognition receptor pathways in plant immunity.