IC4R018-RNA-Seq-2015-25447912 - Revision history

Guangyi: /* The Background of This Project */

2016-07-13T14:35:50Z

‎The Background of This Project

Guangyi: /* Illumina sequencing */

2016-07-13T14:34:12Z

‎Illumina sequencing

Guangyi at 14:33, 13 July 2016

2016-07-13T14:33:13Z

Guangyi at 07:58, 13 July 2016

2016-07-13T07:58:25Z

Guangyi: /* Research Findings */

2016-07-13T07:57:31Z

‎Research Findings

Guangyi: /* Research Findings */

2016-07-13T07:55:53Z

‎Research Findings

Guangyi: /* Research Findings */

2016-07-13T07:54:35Z

‎Research Findings

Guangyi: /* Research Findings */

2016-07-13T07:51:45Z

‎Research Findings

Guangyi: /* Research Findings */

2016-07-13T07:40:14Z

‎Research Findings

Guangyi: /* Research Findings */

2016-07-13T07:36:16Z

‎Research Findings

@@ Line 6: / Line 6: @@
 *The Oryza sativa ammonium transporter,OsAMT1;1, more constitutive expression in shoots and roots, was a prominent member of the OsAMT1 gene family that involved in NH4 + transport in rice plants (Sonoda et al., 2003a). Meanwhile, NH4 + uptake had been reported to be subjected to negative feedback, supposedly from nitrogen metabolites (Lee and Rudge, 1986; Morgan and Jackson,1988; Sonoda et al., 2003b); the early intermediates of nitrogen metabolism such as glutamine might control the expression of ammonium transporter genes in rice and serve as indicators of both the environmental and the cytosolic nitrogen status in rice.
 *With the rapid development of bioengineering and molecular biological technologies, transgenic approaches have broadly been applied to enhance nitrogen uptake and assimilation in rice by the overexpression of novel transgenes. Gene expression, the NH4 + permeability as well as the uptake rate in overexpressed OsAMT1;1 lines were much higher than wild type (WT), meanwhile, higher NH4 + contents also promoted higher expression levels of genes in the nitrogen assimilation pathway, resulting in greater nitrogen assimilates, chlorophyll, starch, sugars, and grain yield in transgenic lines than in the WT under both suboptimal and optimal nitrogen conditions, these results implicated that OsAMT1;1 had the potential to improve nitrogen use efficiency and grain yield in rice (Ranathunge et al., 2014).
 <br>
 ==Plant Culture & Treatment==
 *Rice seeds (O. sativa ssp. Japonica ‘Nipponbare’) were surfacesterilized with 2.5% NaClO for 30 min, rinsed thoroughly with distilled water, then soaked for 24 h in distilled water at 37 °C in the dark.Seeds were then transferred to a plastic mesh placed over the surface of distilled water (in a container of ca. 8 L capacity) and allowed to germinate for 6 days before providing light at 27 °C. Uniform seedlings were selected in batches of three and then transferred to a tank containing 8 L of IRRI nutrient solution (1.25 mM NH4NO3, 0.3 mM KH2PO4,0.35 mM K2SO4, 1 mM CaCl2·2 H2O, 1 mM MgSO4·7 H2O, 0.5 mM Na2SiO3, 20 μM NaFeEDTA, 20 μM H3BO3, 9 μM MnCl2·4 H2O, 0.32 μM CuSO4·5 H2O, 0.77 μM ZnSO4·7 H2O and 0.39 μM Na2MoO4·2 H2O, pH5.8). The outer walls and the bottom of the containers were painted black to prevent light from entering inside the tanks (to prevent algal growth). Plants were grown in a greenhouse at 27 °C/25 °C (day/night) for 10 days with a 16 h light/8 h dark regime (light intensity was 400 μmol m−2 s−1, the relative humidity was kept at 70%). The solution was exchanged every 2 days.

@@ Line 16: / Line 16: @@
 <br>
 ==Illumina sequencing==
-*To obtain a comprehensive understanding of the genes that were affected by changes in ammonium availability, separate RNA samples
+*To obtain a comprehensive understanding of the genes that were affected by changes in ammonium availability, separate RNA samples were prepared respectively from roots and shoots subjected to nitrogen-free (0 mM) and high ammonium (10 mM) for 4 h with 1 mM ammonium as the control, and sequenced using the Illumina sequencing platform.
 <br>
 ==Research Findings==
 *EggNOG (evolutionary genealogy of genes: non-supervised orthologous groups) was used to classify orthologous gene products. The orthologous group genes in this database are annotated with functional description lines (derived by identifying a common denominator for the genes based on their various annotations), with functional categories (i.e. derived from the original cluster of orthologous groups categories) (Jensen et al., 2008). The identified sequences can be aligned to the eggNOG database to predict and classify their possible functions.However, only 311 among 165,411 putative proteins from 6 libraries were classified functionally into 21 molecular families according to eggNOG database (Fig. 1).

@@ Line 14: / Line 14: @@
 *The R package DESeq (Anders and Huber, 2010) was performed to identify the differentially expressed genes. This method represents the widely accepted and accurate analysis approaches of RNA-Seq data.The researchers first mapped high-quality reads to the reference genome of Japonica variety Nipponbare (ftp://ftp.ensemblgenomes.org/pub/release-21/plants/fasta/oryza_sativa/dna/) to calculate the number of reads mapped to each gene in six samples. These raw read counts were then used as the input of DESeq to get the normalized signal for each gene,for this purpose we used the Tophat2 software, version 2.0.3 (Trapnell et al., 2009) with default settings. Aligned reads were assembled by the Cufflinks software (version 2.0.2) (Trapnell et al., 2010, 2012).Those with p-value b 0.05 were considered as significant differential expression, and a 2-fold variance was used to identify the genes differentially expressed between every two libraries (nitrogen-free compared to control and high ammonium compared to control, respectively).Then we performed GO (http://geneontology.org/) and KEGG (http://www.genome.jp/kegg/) function enrichment analysis to the differentially expressed genes (DEGs).
 *Primers designed (Primer Premier 5 software) for each gene were given in Table S1. qRT-PCR was done with a BioRad CFX96™ Real-time System C1000 Thermal Cycler using SYBR® Premix Ex TaqTM (TliRNaseH Plus, TaKaRa, Japan) according to the manufacturer's protocol.
 <br>
 ==Research Findings==

@@ Line 7: / Line 7: @@
 *With the rapid development of bioengineering and molecular biological technologies, transgenic approaches have broadly been applied to enhance nitrogen uptake and assimilation in rice by the overexpression of novel transgenes. Gene expression, the NH4 + permeability as well as the uptake rate in overexpressed OsAMT1;1 lines were much higher than wild type (WT), meanwhile, higher NH4 + contents also promoted higher expression levels of genes in the nitrogen assimilation pathway, resulting in greater nitrogen assimilates, chlorophyll, starch, sugars, and grain yield in transgenic lines than in the WT under both suboptimal and optimal nitrogen conditions, these results implicated that OsAMT1;1 had the potential to improve nitrogen use efficiency and grain yield in rice (Ranathunge et al., 2014).
 *Following the discovery of the major role of the enzyme couple glutamine synthetase (GS) and glutamate synthase (GOGAT) in ammonium assimilation in higher plants (Good et al., 2004; Miflin and Lea, 1976),several laboratories had focused on understanding the mechanisms controlling the regulation of this pathway in plants including rice (Cai et al.,2009; Harrison et al., 2000; Kusano et al., 2011).
+<br>
 ==Plant Culture & Treatment==
 *Rice seeds (O. sativa ssp. Japonica ‘Nipponbare’) were surfacesterilized with 2.5% NaClO for 30 min, rinsed thoroughly with distilled water, then soaked for 24 h in distilled water at 37 °C in the dark.Seeds were then transferred to a plastic mesh placed over the surface of distilled water (in a container of ca. 8 L capacity) and allowed to germinate for 6 days before providing light at 27 °C. Uniform seedlings were selected in batches of three and then transferred to a tank containing 8 L of IRRI nutrient solution (1.25 mM NH4NO3, 0.3 mM KH2PO4,0.35 mM K2SO4, 1 mM CaCl2·2 H2O, 1 mM MgSO4·7 H2O, 0.5 mM Na2SiO3, 20 μM NaFeEDTA, 20 μM H3BO3, 9 μM MnCl2·4 H2O, 0.32 μM CuSO4·5 H2O, 0.77 μM ZnSO4·7 H2O and 0.39 μM Na2MoO4·2 H2O, pH5.8). The outer walls and the bottom of the containers were painted black to prevent light from entering inside the tanks (to prevent algal growth). Plants were grown in a greenhouse at 27 °C/25 °C (day/night) for 10 days with a 16 h light/8 h dark regime (light intensity was 400 μmol m−2 s−1, the relative humidity was kept at 70%). The solution was exchanged every 2 days.
@@ Line 14: / Line 14: @@
 *The R package DESeq (Anders and Huber, 2010) was performed to identify the differentially expressed genes. This method represents the widely accepted and accurate analysis approaches of RNA-Seq data.The researchers first mapped high-quality reads to the reference genome of Japonica variety Nipponbare (ftp://ftp.ensemblgenomes.org/pub/release-21/plants/fasta/oryza_sativa/dna/) to calculate the number of reads mapped to each gene in six samples. These raw read counts were then used as the input of DESeq to get the normalized signal for each gene,for this purpose we used the Tophat2 software, version 2.0.3 (Trapnell et al., 2009) with default settings. Aligned reads were assembled by the Cufflinks software (version 2.0.2) (Trapnell et al., 2010, 2012).Those with p-value b 0.05 were considered as significant differential expression, and a 2-fold variance was used to identify the genes differentially expressed between every two libraries (nitrogen-free compared to control and high ammonium compared to control, respectively).Then we performed GO (http://geneontology.org/) and KEGG (http://www.genome.jp/kegg/) function enrichment analysis to the differentially expressed genes (DEGs).
 *Primers designed (Primer Premier 5 software) for each gene were given in Table S1. qRT-PCR was done with a BioRad CFX96™ Real-time System C1000 Thermal Cycler using SYBR® Premix Ex TaqTM (TliRNaseH Plus, TaKaRa, Japan) according to the manufacturer's protocol.
+<br>
 ==Research Findings==
 *EggNOG (evolutionary genealogy of genes: non-supervised orthologous groups) was used to classify orthologous gene products. The orthologous group genes in this database are annotated with functional description lines (derived by identifying a common denominator for the genes based on their various annotations), with functional categories (i.e. derived from the original cluster of orthologous groups categories) (Jensen et al., 2008). The identified sequences can be aligned to the eggNOG database to predict and classify their possible functions.However, only 311 among 165,411 putative proteins from 6 libraries were classified functionally into 21 molecular families according to eggNOG database (Fig. 1).
@@ Line 32: / Line 32: @@
 [[File:IC4R018-RNA-Seq-2015-25447912-f3b.png |center |thumb |10000px |'''Fig. 3. DEGs in rice roots and shoots under nitrogen-free (0 mM NH4 +) and high ammonium (10 mM NH4 +) conditions. b. DEGs in rice roots and shoots under high ammonium condition (10 mM NH4 +) compared to the control condition (1 mM NH4 +). Base on DESeq software,genes with p-value b 0.05were consideredas significantly differential expression,anda 2-fold variance was used to identify the genes differentially expressed between every two libraries.''']]
 <br>
 ==Labs working on this Project==
 *State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, China

@@ Line 18: / Line 18: @@
 *EggNOG (evolutionary genealogy of genes: non-supervised orthologous groups) was used to classify orthologous gene products. The orthologous group genes in this database are annotated with functional description lines (derived by identifying a common denominator for the genes based on their various annotations), with functional categories (i.e. derived from the original cluster of orthologous groups categories) (Jensen et al., 2008). The identified sequences can be aligned to the eggNOG database to predict and classify their possible functions.However, only 311 among 165,411 putative proteins from 6 libraries were classified functionally into 21 molecular families according to eggNOG database (Fig. 1).
 <br>
-[[File:IC4R018-RNA-Seq-2015-25447912-f1.png |center |thumb |10000px |'''''']]
+[[File:IC4R018-RNA-Seq-2015-25447912-f1.png |center |thumb |10000px |''''''Fig. 1. EggNOG annotations of identified proteins from six libraries under varied nitrogen supplies (0 mM NH4 +, 1 mM NH4 + and 10 mM NH4 +). Histogram presentation of the protein function category. All identified proteins were aligned to the eggNOG database and functionally classified into at least 21molecular families. The x-axis represented the function category the proteins belong to and the y-axis the protein numbers.'''''']]
 <br>
 *Based on sequence homology, 137,157 out of 165,411 genes were categorized into 53 functional groups that belong to the biological process, cellular component and molecular function clusters (Fig. 2).

@@ Line 22: / Line 22: @@
 *Based on sequence homology, 137,157 out of 165,411 genes were categorized into 53 functional groups that belong to the biological process, cellular component and molecular function clusters (Fig. 2).
 <br>
-[[File:IC4R018-RNA-Seq-2015-25447912-f2.png |center |thumb |10000px |'''''']]
+[[File:IC4R018-RNA-Seq-2015-25447912-f2.png |center |thumb |10000px |''''''Fig. 2. Gene Ontology classification of the genes identified from six libraries under varied nitrogen supplies (0 mM NH4 +, 1 mM NH4 + and 10 mM NH4 +). Histogram presentation of the Gene Ontology classification. The results were summarized in the three main GO categories: biological process, cellular component and molecular function. The x-axis represented the functional category of the genes and the y-axis the gene numbers.'''''']]
 <br>
 *The results are derived from the differential expression analysis experiment between every two libraries (the nitrogen-free transcripts compared to the control and the high ammonium transcripts compared to the control, respectively). When compared the nitrogen-free transcripts with the control ones (1 mM), of 394 DEGs in root transcripts,143 DEGs were up-regulated and 251 DEGs were down-regulated, and of 468 DEGs in shoot transcripts, 119 DEGs were up-regulated and 349 DEGs were down-regulated respectively (Fig. 3a).

@@ Line 30: / Line 30: @@
 *When transcripts from the high ammonium condition were compared with the control, 63 DEGs in roots and 115 DEGs in shoots were obtained respectively. In roots, 6 DEGs showed higher expression and 57 lower expression, and in shoots, 93 DEGs showed higher expression and 22 lower expression (Fig. 3b).
 <br>
-[[File:IC4R018-RNA-Seq-2015-25447912-f3b.png |center |thumb |10000px |'''''']]
+[[File:IC4R018-RNA-Seq-2015-25447912-f3b.png |center |thumb |10000px |'''Fig. 3. DEGs in rice roots and shoots under nitrogen-free (0 mM NH4 +) and high ammonium (10 mM NH4 +) conditions. b. DEGs in rice roots and shoots under high ammonium condition (10 mM NH4 +) compared to the control condition (1 mM NH4 +). Base on DESeq software,genes with p-value b 0.05were consideredas significantly differential expression,anda 2-fold variance was used to identify the genes differentially expressed between every two libraries.''']]
 <br>

@@ Line 26: / Line 26: @@
 *The results are derived from the differential expression analysis experiment between every two libraries (the nitrogen-free transcripts compared to the control and the high ammonium transcripts compared to the control, respectively). When compared the nitrogen-free transcripts with the control ones (1 mM), of 394 DEGs in root transcripts,143 DEGs were up-regulated and 251 DEGs were down-regulated, and of 468 DEGs in shoot transcripts, 119 DEGs were up-regulated and 349 DEGs were down-regulated respectively (Fig. 3a).
 <br>
-[[File:IC4R018-RNA-Seq-2015-25447912-f3a.png |center |thumb |10000px |'''''']]
+[[File:IC4R018-RNA-Seq-2015-25447912-f3a.png |center |thumb |10000px |''''''Fig. 3. DEGs in rice roots and shoots under nitrogen-free (0 mM NH4 +) and high ammonium (10 mM NH4 +) conditions. a. DEGs in rice roots and shoots under nitrogen-free condition (0 mM NH4 +) compared to the control condition (1 mM NH4 +). Base on DESeq software, genes with p-value b 0.05 were considered as significantly differential expression,and a 2-foldvariancewas usedto identify the genes differentially expressedbetween every two libraries.'''''']]
 <br>
 *When transcripts from the high ammonium condition were compared with the control, 63 DEGs in roots and 115 DEGs in shoots were obtained respectively. In roots, 6 DEGs showed higher expression and 57 lower expression, and in shoots, 93 DEGs showed higher expression and 22 lower expression (Fig. 3b).