Genome-wide analysis of rice cis-natural antisense transcription under cadmium exposure using strand-specific RNA-Seq.

BACKGROUND: The elucidation of novel transcripts and their expression in response to various stress conditions is necessary to understand the transcriptional network of plants as an adaptation to biotic and abiotic stresses. We performed strand-specific RNA-Seq (ssRNA-Seq) on rice exposed to cadmium (Cd) for 24 h and investigated the expression of cis-natural antisense transcripts (cis-NATs), a class of endogenous coding or non-protein-coding RNAs with sequence complementarity to the opposite strands of RAP transcripts. RESULTS: Many RAP transcripts possessed cis-NATs and these cis-NATs were responsive to some extent. Cis-NATs were upregulated from 26, 266 and 409 RAP gene loci, while 2054, 2501 and 2825 RAP transcripts were upregulated from 38,123 RAP loci under high Cd exposure in roots at 1, 12 and 24 h, respectively. In addition, most of the upregulated cis-NATs showed little upregulation under ABA or cold treatment. A number of cis-NATs were upregulated from less than 35 RAP gene loci in different tissue and time-point combinations under low Cd exposure, suggesting that cis-NATs respond to environmental stress. Furthermore, 409 RAP transcripts with upregulated cis-NATs were classified into three groups based on the expression of the RAP transcripts from the opposite DNA strand, including 138 upregulated, 128 invariable, and 143 downregulated transcripts, although the responses of cis-NATs and RAP transcripts were not always correlated. CONCLUSIONS: We have shown that the cis-NATs identified by ssRNA-Seq analysis are novel genes and that some of them are stress-specific and show different responses depending on the degree of stress and tissue. These results improve our understanding of the complete molecular mechanism of plant adaptation to Cd exposure.

A novel open-type biosensor for the in-situ monitoring of biochemical oxygen demand in an aerobic environment.

Biochemical oxygen demand (BOD) is a widely used index of water-quality assessment. Since bioelectrochemical BOD biosensors require anaerobic conditions for anodic reactions, they are not directly used in aerobic environments such as aeration tanks. Normally, the BOD biosensors are closed-type, where the anode is packed inside a closed chamber to avoid exposure to oxygen. In this study, a novel bioelectrochemical open-type biosensor was designed for in-situ monitoring of BOD during intermittent aeration. The open-type anode, without any protection against exposure to oxygen, was directly inserted into an intermittently aerated tank filled with livestock wastewater. Anodic potential was controlled using a potentiostat. Interestingly, this novel biosensor generated similar levels of current under both aerating and non-aerating conditions, and showed a logarithmic correlation (R2 > 0.9) of current with BOD concentrations up to 250 mg/L. Suspended solids in the wastewater attached to and covered the whole anode, presumably leading to the production of anaerobic conditions inside the covered anode via biological oxygen removal. Exoelectrogenic anaerobes (Geobacter spp.) were detected inside the covered anode using the 16S-rRNA gene. This biosensor will have various practical applications, such as the automatic control of aeration intensity and the in-situ monitoring of natural water environments.

Genomic dissection of a 'Fuji' apple cultivar: re-sequencing, SNP marker development, definition of haplotypes, and QTL detection.

'Fuji' is one of the most popular and highly-produced apple cultivars worldwide, and has been frequently used in breeding programs. The development of genotypic markers for the preferable phenotypes of 'Fuji' is required. Here, we aimed to define the haplotypes of 'Fuji' and find associations between haplotypes and phenotypes of five traits (harvest day, fruit weight, acidity, degree of watercore, and flesh mealiness) by using 115 accessions related to 'Fuji'. Through the re-sequencing of 'Fuji' genome, total of 2,820,759 variants, including single nucleotide polymorphisms (SNPs) and insertions or deletions (indels) were detected between 'Fuji' and 'Golden Delicious' reference genome. We selected mapping-validated 1,014 SNPs, most of which were heterozygous in 'Fuji' and capable of distinguishing alleles inherited from the parents of 'Fuji' (i.e., 'Ralls Janet' and 'Delicious'). We used these SNPs to define the haplotypes of 'Fuji' and trace their inheritance in relatives, which were shown to have an average of 27% of 'Fuji' genome. Analysis of variance (ANOVA) based on 'Fuji' haplotypes identified one quantitative trait loci (QTL) each for harvest time, acidity, degree of watercore, and mealiness. A haplotype from 'Delicious' chr14 was considered to dominantly cause watercore, and one from 'Ralls Janet' chr1 was related to low-mealiness.

Expression of a putative dioxygenase gene adjacent to an insertion mutation is involved in the short internodes of columnar apples (Malus × domestica).

Determining the molecular mechanism of fruit tree architecture is important for tree management and fruit production. An apple mutant 'McIntosh Wijcik', which was discovered as a bud mutation from 'McIntosh', exhibits a columnar growth phenotype that is controlled by a single dominant gene, Co. In this study, the mutation and the Co gene were analyzed. Fine mapping narrowed the Co region to a 101 kb region. Sequence analysis of the Co region and the original wild-type co region identified an insertion mutation of an 8202 bp long terminal repeat (LTR) retroposon in the Co region. Segregation analysis using a DNA marker based on the insertion polymorphism showed that the LTR retroposon was closely associated with the columnar growth phenotype. RNA-seq and RT-PCR analysis identified a promising Co candidate gene (91071-gene) within the Co region that is specifically expressed in 'McIntosh Wijcik' but not in 'McIntosh'. The 91071-gene was located approximately 16 kb downstream of the insertion mutation and is predicted to encode a 2-oxoglutarate-dependent dioxygenase involved in an unknown reaction. Overexpression of the 91071-gene in transgenic tobaccos and apples resulted in phenotypes with short internodes, like columnar apples. These data suggested that the 8202 bp retroposon insertion in 'McIntosh Wijcik' is associated with the short internodes of the columnar growth phenotype via upregulated expression of the adjacent 91071-gene. Furthermore, the DNA marker based on the insertion polymorphism could be useful for the marker-assisted selection of columnar apples.

Detection of Novel QTLs Regulating Grain Size in Extra-Large Grain Rice (Oryza sativa L.) Lines.

BACKGROUND: Grain size is an important trait that affects rice yield. Although many genes that contribute to grain size have been cloned from mutants or by quantitative trait locus (QTL) analysis based on bi-parental mapping, the molecular mechanisms underlying grain-size determination remain poorly understood. In this study, we identified the lines with the largest grain size and detected novel QTLs affecting the grain size. RESULTS: We screened the National Institute for Agrobiological Sciences Genebank database and identified two rice lines, BG23 with the widest grain and LG10 with the longest grain. Using these two lines, we performed QTL analysis for grain size. Eight QTLs were detected during the QTL analyses using F2 populations derived from crosses between the large-grain lines BG23 or LG10 and the middle-size grain cultivars Nipponbare and Kasalath. Both BG23 and LG10 possessed large-grain alleles of four major QTLs: GW2, GS3, qSW5/GW5, and GW8. Other three minor QTLs were derived from BG23. However, these QTLs did not explain the differences in grain size between these two lines. Additionally, four QTLs for grain length or width were detected in an F2 population derived from a cross between BG23 and LG10; this population lacked the strong effects of the four major QTLs shared by both parent plants. Of these newly detected QTLs, the effects of two QTLs, GL3b and GL6, were confirmed by progeny testing. Comparison of the length of inner epidermal cells in plants homozygous for BG23 and LG10 alleles indicated that GL3b and GL6 genes regulate cell elongation and cell division, respectively. CONCLUSIONS: In this study, we detected 12 loci including 14 QTLs regulating grain size from two lines with largest grains available in Japanese stock. Of these loci, we confirmed the effect of two gene loci and mapped their candidate region. Identification of novel genes regulating grain size will contribute to our understanding of the molecular mechanisms controlling grain size.

Genome-Wide Transcriptome Analysis of Cadmium Stress in Rice.

Rice growth is severely affected by toxic concentrations of the nonessential heavy metal cadmium (Cd). To elucidate the molecular basis of the response to Cd stress, we performed mRNA sequencing of rice following our previous study on exposure to high concentrations of Cd (Oono et al., 2014). In this study, rice plants were hydroponically treated with low concentrations of Cd and approximately 211 million sequence reads were mapped onto the IRGSP-1.0 reference rice genome sequence. Many genes, including some identified under high Cd concentration exposure in our previous study, were found to be responsive to low Cd exposure, with an average of about 11,000 transcripts from each condition. However, genes expressed constitutively across the developmental course responded only slightly to low Cd concentrations, in contrast to their clear response to high Cd concentration, which causes fatal damage to rice seedlings according to phenotypic changes. The expression of metal ion transporter genes tended to correlate with Cd concentration, suggesting the potential of the RNA-Seq strategy to reveal novel Cd-responsive transporters by analyzing gene expression under different Cd concentrations. This study could help to develop novel strategies for improving tolerance to Cd exposure in rice and other cereal crops.

Enhanced electrical power generation using flame-oxidized stainless steel anode in microbial fuel cells and the anodic community structure.

BACKGROUND: Carbon-based materials are commonly used as anodes in microbial fuel cells (MFCs), whereas metal and metal-oxide-based materials are not used frequently because of low electrical output. Stainless steel is a low-cost material with high conductivity and physical strength. In this study, we investigated the power generation using flame-oxidized (FO) stainless steel anodes (SSAs) in single-chambered air-cathode MFCs. The FO-SSA performance was compared to the performance of untreated SSA and carbon cloth anode (CCA), a common carbonaceous electrode. The difference in the anodic community structures was analyzed using high-throughput sequencing of the V4 region in 16S rRNA gene. RESULTS: Flame oxidation of SSA produced raised node-like sites, predominantly consisting of hematite (Fe2O3), on the surface, as determined by X-ray diffraction spectroscopy. The flame oxidation enhanced the maximum power density (1063 mW/m(2)) in MFCs, which was 184 and 24 % higher than those for untreated SSA and CCA, respectively. The FO-SSA exhibited 8.75 and 2.71 times higher current production than SSA and CCA, respectively, under potentiostatic testing conditions. Bacteria from the genus Geobacter were detected at a remarkably higher frequency in the biofilm formed on the FO-SSA (8.8-9.2 %) than in the biofilms formed on the SSA and CCA (0.7-1.4 %). Bacterial species closely related to Geobacter metallireducens (>99 % identity in the gene sequence) were predominant (93-96 %) among the genus Geobacter in the FO-SSA biofilm, whereas bacteria with a 100 % identity to G. anodireducens were abundant (>55 %) in the SSA and CCA biofilms. CONCLUSIONS: This is the first demonstration of power generation using an FO-SSA in MFCs. Flame oxidation of the SSA enhances electricity production in MFCs, which is higher than that with the common carbonaceous electrode, CCA. The FO-SSA is not only inexpensive but also can be prepared using a simple method. To our knowledge, this study reveals, for the first time, that the predominant Geobacter species in the biofilm depends on the anode material. The high performance of the FO-SSA could result from the particularly high population of bacteria closely related to G. metallireducens in the biofilm.

Comprehensive detection of bacterial carbohydrate-active enzyme coding genes expressed in cow rumen.

To find the abundant and characteristic fibrolytic enzyme-coding gene expressed in fiber-associating microbiota, a metatranscriptomic data set was obtained from fiber-associating microbiota, and it was compared with that of rumen fluid-floating microbiota and two metagenomic data sets. Fibrolytic rumen bacteria associate with plant polysaccharide and hydrolyze it in the rumen. We obtained a metatranscriptomic assembly from fiber-associating microbiota in three ruminally fistulated Holstein cows fed timothy (Phleum pratense) hay. Each metatranscriptomic data set involved over a thousand of the glycoside hydrolase (GH) gene transcripts that accounted for about 1% of total protein coding gene transcripts. Three-quarters of the total GH gene transcripts were dominated by non-structural oligosaccharide-acting hydrolase gene transcripts. In the fiber-associating microbiota, endo-cellulase coding gene families, especially GHs 9 and 5, were abundantly detected, and GHs 9, 11, 30 and 43, carbohydrate esterase 8 and carbohydrate-binding module 6 were characteristically detected. Most fibrolytic gene transcripts assigned to Fibrobacter succinogenes were detected in fiber-associating sections, and GHs 45, 44, 74, 11, 30 and 16 were Fibrobacter-characteristically detected. The metatranscriptomic assembly highlighted the characteristic fibrolytic enzymes expressed in the fiber-associated rumen microbiota and offered access to the fibrolytic activities in each fibrolytic bacteria.

TENOR: Database for Comprehensive mRNA-Seq Experiments in Rice.

Here we present TENOR (Transcriptome ENcyclopedia Of Rice, http://tenor.dna.affrc.go.jp), a database that encompasses large-scale mRNA sequencing (mRNA-Seq) data obtained from rice under a wide variety of conditions. Since the elucidation of the ability of plants to adapt to various growing conditions is a key issue in plant sciences, it is of great interest to understand the regulatory networks of genes responsible for environmental changes. We used mRNA-Seq and performed a time-course transcriptome analysis of rice, Oryza sativa L. (cv. Nipponbare), under 10 abiotic stress conditions (high salinity; high and low phosphate; high, low and extremely low cadmium; drought; osmotic; cold; and flood) and two plant hormone treatment conditions (ABA and jasmonic acid). A large number of genes that were responsive to abiotic stresses and plant hormones were detected by differential expression analysis. Furthermore, several responsive genes were found to encode transcription factors that could control the transcriptional network of stress responses, but the timing of the induction of these genes was not uniform across conditions. A significant number of cis-regulatory elements were enriched in the promoter regions of the responsive genes and were shared among conditions. These data suggest that some key components of gene regulation networks are shared between different stress signaling pathways. All the resources (novel genes identified from mRNA-Seq data, expression profiles, co-expressed genes and cis-regulatory elements) can be searched for and are available in TENOR.

Construction of a high-density mutant library in soybean and development of a mutant retrieval method using amplicon sequencing.

BACKGROUND: Functions of most genes predicted in the soybean genome have not been clarified. A mutant library with a high mutation density would be helpful for functional studies and for identification of novel alleles useful for breeding. Development of cost-effective and high-throughput protocols using next generation sequencing (NGS) technologies is expected to simplify the retrieval of mutants with mutations in genes of interest. RESULTS: To increase the mutation density, seeds of the Japanese elite soybean cultivar Enrei were treated with the chemical mutagen ethyl methanesulfonate (EMS); M2 seeds produced by M1 plants were treated with EMS once again. The resultant library, which consisted of DNA and seeds from 1536 plants, revealed large morphological and physiological variations. Based on whole-genome re-sequencing analysis of 12 mutant lines, the average number of base changes was 12,796 per line. On average, 691 and 35 per line were missense and nonsense mutations, respectively. Two screening strategies for high resolution melting (HRM) analysis and indexed amplicon sequencing were designed to retrieve the mutants; the mutations were confirmed by Sanger sequencing as the final step. In comparison with HRM screening of several genes, indexed amplicon sequencing allows one to scan a longer sequence range and skip screening steps and to know the sequence information of mutation because it uses systematic DNA pooling and the index of NGS reads, which simplifies the discovery of mutants with amino acid substitutions. CONCLUSIONS: A soybean mutant library with a high mutation density was developed. A high mutation density (1 mutation/74 kb) was achieved by repeating the EMS treatment. The mutation density of our library is sufficiently high to obtain a plant in which a gene is nonsense mutated. Thus, our mutant library and the indexed amplicon sequencing will be useful for functional studies of soybean genes and have a potential to yield useful mutant alleles for soybean breeding.

Genome-wide transcriptome analysis reveals that cadmium stress signaling controls the expression of genes in drought stress signal pathways in rice.

Plant growth is severely affected by toxic concentrations of the non-essential heavy metal cadmium (Cd). Comprehensive transcriptome analysis by RNA-Seq following cadmium exposure is required to further understand plant responses to Cd and facilitate future systems-based analyses of the underlying regulatory networks. In this study, rice plants were hydroponically treated with 50 µM Cd for 24 hours and ∼60,000 expressed transcripts, including transcripts that could not be characterized by microarray-based approaches, were evaluated. Upregulation of various ROS-scavenging enzymes, chelators and metal transporters demonstrated the appropriate expression profiles to Cd exposure. Gene Ontology enrichment analysis of the responsive transcripts indicated the upregulation of many drought stress-related genes under Cd exposure. Further investigation into the expression of drought stress marker genes such as DREB suggested that expression of genes in several drought stress signal pathways was activated under Cd exposure. Furthermore, qRT-PCR analyses of randomly selected Cd-responsive metal transporter transcripts under various metal ion stresses suggested that the expression of Cd-responsive transcripts might be easily affected by other ions. Our transcriptome analysis demonstrated a new transcriptional network linking Cd and drought stresses in rice. Considering our data and that Cd is a non-essential metal, the network underlying Cd stress responses and tolerance, which plants have developed to adapt to other stresses, could help to acclimate to Cd exposure. Our examination of this transcriptional network provides useful information for further studies of the molecular mechanisms of plant adaptation to Cd exposure and the improvement of tolerance in crop species.

A BAC physical map of aus rice cultivar 'Kasalath', and the map-based genomic sequence of 'Kasalath' chromosome 1.

Comparative analysis using available genomic resources within closely related species is an effective way to investigate genomic sequence and structural diversity. Rice (Oryza sativa L.) has undergone significant physiological and morphological changes during its domestication and local adaptation. We present a complete bacterial artificial chromosome (BAC) physical map for the aus rice cultivar 'Kasalath', which covers 90% of the sequence of temperate japonica rice cultivar 'Nipponbare'. Examination of physical distances between computational and experimental measurements of 'Kasalath' BAC insert size revealed the presence of more than 500 genomic regions that appear to have significant chromosome structural changes between the two cultivars. In particular, a genomic region on the long arm of 'Kasalath' chromosome 11 carrying a disease-resistance gene cluster was greatly expanded relative to the 'Nipponbare' genome. We also decoded 41.37 Mb of high-quality genomic sequence from 'Kasalath' chromosome 1. Extensive comparisons of chromosome 1 between 'Kasalath' and 'Nipponbare' led to the discovery of 317,843 single-nucleotide polymorphisms (SNPs) and 66,331 insertion/deletion (indel) sites. Nearly two-thirds of the expressed genes on rice chromosome 1 carried natural variations involving SNPs and/or indels that resulted in substitutions, insertions or deletions of amino acids in one cultivar relative to the other. We also observed gain and loss of genes caused by large indels. This study provides an important framework and an invaluable dataset for further understanding of the molecular mechanisms underlying the evolution and functions of the rice genome.