CIMB, Vol. 44, Pages 6404-6427: Comprehensive Analysis of Glutamate Receptor-like Genes in Rice (Oryza sativa L.): Genome-Wide Identification, Characteristics, Evolution, Chromatin Accessibility, gcHap Diversity, Population Variation and Expression Analysis

Conceptualization, Y.S.; Data curation, W.Z. and H.L.; Formal analysis, W.Z. and M.X.; Funding acquisition, Y.S. and W.Z.; Resources, H.L., F.Z., Z.C., G.A., S.H., Y.H., X.W., J.L., X.Z., Y.Z. and S.R.; Validation, W.Z.; Visualization, S.H.; Writing—original draft, Y.S., W.Z. and M.X.; Writing—review and editing, Y.S., G.A. and S.R. All authors have read and agreed to the published version of the manuscript.

Figure 1. Characteristics of the glutamate receptor-like (GLR) genes in rice. (a) Chromosome location of the OsGLR genes; (b) phylogenetic tree of GLR genes in rice, tomato and Arabidopsis thaliana; (c) intron–exon distribution of the OsGLR genes; (d) conserved motif prediction of the OsGLR genes. The conserved motif can provide insights into how patterns of residue conservation and divergence in a protein family relate to functional properties, and can provide useful links to more detailed information that may be helpful in understanding that those sequences have all four subfamily/structure/function relationships. There are four subfamilies, two large subfamilies and two small subfamilies. Arabidopsis thaliana has only three subfamilies (groups 1, 2 and 4).

Figure 1. Characteristics of the glutamate receptor-like (GLR) genes in rice. (a) Chromosome location of the OsGLR genes; (b) phylogenetic tree of GLR genes in rice, tomato and Arabidopsis thaliana; (c) intron–exon distribution of the OsGLR genes; (d) conserved motif prediction of the OsGLR genes. The conserved motif can provide insights into how patterns of residue conservation and divergence in a protein family relate to functional properties, and can provide useful links to more detailed information that may be helpful in understanding that those sequences have all four subfamily/structure/function relationships. There are four subfamilies, two large subfamilies and two small subfamilies. Arabidopsis thaliana has only three subfamilies (groups 1, 2 and 4).

Cimb 44 00437 g001

Figure 2. Prediction of transmembrane helices of the 26 glutamate receptor-like (GLR) genes in rice. The plot shows the posterior probabilities of inside/outside/TM helix. The plot is obtained by calculating the total probability that a residue sits in helix, inside or outside summed over all possible paths through the model.

Figure 2. Prediction of transmembrane helices of the 26 glutamate receptor-like (GLR) genes in rice. The plot shows the posterior probabilities of inside/outside/TM helix. The plot is obtained by calculating the total probability that a residue sits in helix, inside or outside summed over all possible paths through the model.

Cimb 44 00437 g002

Figure 3. Three-dimensional (3D) structure prediction of the glutamate receptor-like (GLR) genes in rice. Green helices represent α-helices, blue arrows indicate β-strands and faint lines indicate coil. Confidence is colored from high (red) to low (blue) using a rainbow spectrum.

Figure 3. Three-dimensional (3D) structure prediction of the glutamate receptor-like (GLR) genes in rice. Green helices represent α-helices, blue arrows indicate β-strands and faint lines indicate coil. Confidence is colored from high (red) to low (blue) using a rainbow spectrum.

Cimb 44 00437 g003

Figure 4. Cis-acting regulatory elements of the 26 glutamate receptor-like (GLR) genes in rice. Different colors indicate different cis-acting elements. ABRE: ABA-Responsive Element; ACE: Angiotensin-Converting Enzymes; LTR: Low-Temperature Responsiveness; MBS: Multichain Binding Site; MRE: MYB-Recognizing Elements; MSA-like: Mitosis-Specific Activator; SP1: Stress-related Protein.

Figure 4. Cis-acting regulatory elements of the 26 glutamate receptor-like (GLR) genes in rice. Different colors indicate different cis-acting elements. ABRE: ABA-Responsive Element; ACE: Angiotensin-Converting Enzymes; LTR: Low-Temperature Responsiveness; MBS: Multichain Binding Site; MRE: MYB-Recognizing Elements; MSA-like: Mitosis-Specific Activator; SP1: Stress-related Protein.

Cimb 44 00437 g004

Figure 5. The functional interaction network of the 26 glutamate receptor-like (GLR) genes in rice. The network view summarizes the network of predicted associations for a particular group of proteins. The network nodes are proteins. The edges represent the predicted functional associations. The thickness of the lines indicates the degree of confidence prediction of the interaction. The middle red protein represents the rice glutamate receptor-like gene or the homolog of rice glutamate receptor gene. Red line—indicates the presence of fusion evidence; green line—neighborhood evidence; blue line—co-occurrence evidence; purple line—experimental evidence; yellow line—text-mining evidence; light blue line—database evidence; black line—co-expression evidence.

Figure 5. The functional interaction network of the 26 glutamate receptor-like (GLR) genes in rice. The network view summarizes the network of predicted associations for a particular group of proteins. The network nodes are proteins. The edges represent the predicted functional associations. The thickness of the lines indicates the degree of confidence prediction of the interaction. The middle red protein represents the rice glutamate receptor-like gene or the homolog of rice glutamate receptor gene. Red line—indicates the presence of fusion evidence; green line—neighborhood evidence; blue line—co-occurrence evidence; purple line—experimental evidence; yellow line—text-mining evidence; light blue line—database evidence; black line—co-expression evidence.

Cimb 44 00437 g005

Figure 6. Chromatin accessibility of the glutamate receptor-like (GLR) genes in rice. Six tissues (root (RT), young leaf (YL), flag leaf (FL), young panicle (YP), lemma and palea (LP), and stamen and pistil (SP)) of Zhenshan 97 (an indica/xian variety) were collected for ATAC-seq experiment.

Figure 6. Chromatin accessibility of the glutamate receptor-like (GLR) genes in rice. Six tissues (root (RT), young leaf (YL), flag leaf (FL), young panicle (YP), lemma and palea (LP), and stamen and pistil (SP)) of Zhenshan 97 (an indica/xian variety) were collected for ATAC-seq experiment.

Cimb 44 00437 g006

Figure 7. Gene presence (a) and distribution frequency (b) of the glutamate receptor-like (GLR) genes in 3010 rice accessions. The first heatmap indicates the distribution in subspecies of this gene. The second heatmap indicates the presence frequency in subgroups of this gene. A total of 3010 rice accessions were divided into five subspecies, including JAP (japonica, 801 accessions), IND (indica, 1764 accessions), AUS (aus/boro, 221 accessions), ARO (aromatic basmati/sadri, 101 accessions) and ADM (admixed, 123 accessions). All of these were further grouped into 12 groups according to the classification of their corresponding rice accessions. These groups include four subgroups (IG1, IG2, IG3, IG4, IG5) of Indica subspecies AUSG6, four subgroups (JG7, JG8, JG9, JG10) of Japonica subspecies AROG11, and admixtures (ADM).

Figure 7. Gene presence (a) and distribution frequency (b) of the glutamate receptor-like (GLR) genes in 3010 rice accessions. The first heatmap indicates the distribution in subspecies of this gene. The second heatmap indicates the presence frequency in subgroups of this gene. A total of 3010 rice accessions were divided into five subspecies, including JAP (japonica, 801 accessions), IND (indica, 1764 accessions), AUS (aus/boro, 221 accessions), ARO (aromatic basmati/sadri, 101 accessions) and ADM (admixed, 123 accessions). All of these were further grouped into 12 groups according to the classification of their corresponding rice accessions. These groups include four subgroups (IG1, IG2, IG3, IG4, IG5) of Indica subspecies AUSG6, four subgroups (JG7, JG8, JG9, JG10) of Japonica subspecies AROG11, and admixtures (ADM).

Cimb 44 00437 g007

Figure 8. The gene-coding sequence haplotype (gcHap) diversity of the glutamate receptor-like (GLR) genes in modern and 3010 rice accessions. The 3010 rice accessions were divided into five subspecies, including Xian (indica, 1764 accessions), Geng (japonica, 801 accessions), Aus (aus/boro, 221 accessions), Aro (aromatic basmati/sadri, 101 accessions) and Adm (admixed, 123 accessions).

Figure 8. The gene-coding sequence haplotype (gcHap) diversity of the glutamate receptor-like (GLR) genes in modern and 3010 rice accessions. The 3010 rice accessions were divided into five subspecies, including Xian (indica, 1764 accessions), Geng (japonica, 801 accessions), Aus (aus/boro, 221 accessions), Aro (aromatic basmati/sadri, 101 accessions) and Adm (admixed, 123 accessions).

Cimb 44 00437 g008

Figure 9. Tissue expression patterns of the glutamate receptor-like (GLR) genes in rice. Gene expression data generated by the Affymetrix ATH1 array are normalized by MAS 5.0 and RMA methods with a trimmed mean target intensity (TGT) value.

Figure 9. Tissue expression patterns of the glutamate receptor-like (GLR) genes in rice. Gene expression data generated by the Affymetrix ATH1 array are normalized by MAS 5.0 and RMA methods with a trimmed mean target intensity (TGT) value.

Cimb 44 00437 g009

Figure 10. Expression analysis of the glutamate receptor-like (GLR) genes in rice. The shoot (a) and root (b) expression variation of OsGLR genes under a variety of hormone treatment conditions. Blue is low expression, red is high expression. Expression variation of OsGLR genes under drought (c), flood (d) and cold (e) stress conditions. Green is low expression, red is high expression.

Figure 10. Expression analysis of the glutamate receptor-like (GLR) genes in rice. The shoot (a) and root (b) expression variation of OsGLR genes under a variety of hormone treatment conditions. Blue is low expression, red is high expression. Expression variation of OsGLR genes under drought (c), flood (d) and cold (e) stress conditions. Green is low expression, red is high expression.

Cimb 44 00437 g010

Table 1. Semi-quantitative primers of the glutamate receptor-like (GLR) genes in rice.

Table 1. Semi-quantitative primers of the glutamate receptor-like (GLR) genes in rice.

Gene idGene NameThe Upstream Primer (5’-3’)Downstream Primers (5’-3’)Os09t0431100OsGLR9.7GATGGTTGCTGATGGGGCATTCCTTCAGGAACACCCACGTCCOs09t0429000OsGLR9.3ATAACAGCTTCTGGGTATGCGATATGGCTGCATCATATACCCCTAGOs06t0680500OsGLR6.12CTAAACACAATCGACGAGTACGCAATCTCTCTGGAATGCGAATCCCOs04t0585200OsGLR4.1ACCGTCAATTTGTATCAGTTGATAGATAAGCTCCGAATAGCTTGGATTCOs06t0188700OsGLR6.4TCGTGGTGGACATGACGAGCCTTGATAAGGTCTTCGGCGGCOs07t0103100OsGLR7.1ATCATCCAAGGTCTGCAGGTGATAGGAAGTAAGGGTATTGGGAGGAOs02t0787600OsGLR2.2TGTTCGACGAGGTCATGAAGATTCATCTGCCTTCTGTGACGACAOs09t0428300OsGLR9.1AGAAAGGCAGAGGAATTCCATGTCCAGCAAATCAAGAACTGCAGATOs03t0234200OsUBQ AAGAAGCTGAAGCATCCAGCCCAGGACAAGATGATCTGCC

Table 2. Genome-wide identification and characteristics of the glutamate receptor-like (GLR) genes in rice.

Table 2. Genome-wide identification and characteristics of the glutamate receptor-like (GLR) genes in rice.

SubfamilyRAP IDGene NameNo. of Amino AcidsMolecular
Weight (Da)Isoelectric Point (PI)Aliphatic
IndexStart (bp)End (bp)LocationIOs02g0787600OsGLR2.2988107,252.16.483.9433,469,13833,475,994P(9/14),Ch(3/14)IOs06g0155000OsGLR6.121722,390.236.2976.542,836,4462,837,851Ch(11/14),N(3/14)IOs06g0188400OsGLR6.248152,305.658.3491.734,459,0074,463,882Cy(7/14),Ch(3/14)IOs09g0431100OsGLR9.7955100,789.96.1886.1515,752,94915,762,817P(7/14),V(3/14)IIOs09g0428300OsGLR9.187296,824.567.9996.915,564,67515,569,470P(12/14)IIOs09g0428600OsGLR9.254160,230.776.1588.8415,574,05815,580,334E(5/14), Ch(3/14), V(3/14)IIOs09g0429000OsGLR9.3933104,346.785.9692.1815,588,26915,594,378P(13/14)IIOs09g0429200OsGLR9.4948105,520.296.3390.8815,599,14915,606,593P(8/14),ER(3/14)IIOs09g0429400OsGLR9.5934104,377.75.9791.7515,610,71515,616,401P(11/14)IIOs09g0429500OsGLR9.6946105,792.315.9492.3215,620,33515,625,678P(13/14)IIOs09g0431200OsGLR9.8950106,149.725.4190.3915,765,94515,770,717P(13/14)IIIOs06g0188600OsGLR6.367074,295.455.8780.94,472,3524,476,661P(8/14)IIIOs06g0188700OsGLR6.470276,300.59.5281.154,481,2404,487,783P(10/14)IIIOs06g0188800OsGLR6.519721,457.388.4584.674,489,5884,493,909Ch(10/14)IIIOs06g0189100OsGLR6.639743,568.456.0882.924,500,7664,502,665P(6/14),V(4/14),ER(3/14)IIIOs06g0190000OsGLR6.736539,821.644.8879.014,556,4874,558,077E(4/14),V(3/14)IIIOs06g0190500OsGLR6.853159,380.029.189.764,571,0984,574,904P(8/14),M(3/14)IIIOs06g0190700OsGLR6.954259,762.918.7382.24,587,9874,591,068P(8/14)IIIOs06g0190800OsGLR6.1082091,791.688.4192.214,593,9164,597,945P(9/14), V(3/14)IVOs02g0117500OsGLR2.1944105,303.795.6794.29918,797923,968P(13/14)IVOs04g0585200OsGLR4.154059,717.236.3695.6729,565,75429,570,387P(10/14)IVOs06g0246600OsGLR6.1114715,686.374.0590.277,607,1087,607,939Ch(4/14),E(5/14),M(3/14)IVOs06g0680500OsGLR6.1289098,881.675.3490.3928,345,36628,350,758P(11/14)IVOs07g0103100OsGLR7.1956104,440.946.2292.01191,053196,429P(10/14),ER(3/14)IVOs07g0522600OsGLR7.263769,593.235.3593.0520,221,07720,225,924P(5/14),ER(4/14),V(3/14)IVOs09g0485700OsGLR9.99510,497.734.5483.1618,740,82618,741,187P(10/14),Cy(3/14)

Table 3. KAKS analysis of the glutamate receptor-like (GLR) genes in rice.

Table 3. KAKS analysis of the glutamate receptor-like (GLR) genes in rice.

Seq_1Seq_2KaKsKa_KsMyaOsGLR6.8OsGLR6.90.170.320.5324.53OsGLR9.5OsGLR9.30.060.270.2220.93OsGLR9.5OsGLR9.60.070.280.2521.61OsGLR9.6OsGLR9.30.060.270.2221.07

留言 (0)

沒有登入
gif