Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 12, Article number: 13829 (2022 ) Cite this article
Sulfur is an essential element required for plant growth and development, physiological processes and stress responses. Sulfur-encoding biosynthetic genes are involved in the primary sulfur assimilation pathway, regulating various mechanisms at the gene, cellular and system levels, and in the biosynthesis of sulfur-containing compounds (SCCs). In this study, the SCC-encoding biosynthetic genes in rice were identified using a sulfur-dependent model plant, the Arabidopsis. A total of 139 AtSCC from Arabidopsis were used as reference sequences in search of putative rice SCCs. At similarity index > 30%, the similarity search against Arabidopsis SCC query sequences identified 665 putative OsSCC genes in rice. The gene synteny analysis showed a total of 477 syntenic gene pairs comprised of 89 AtSCC and 265 OsSCC biosynthetic genes in Arabidopsis and rice, respectively. Phylogenetic tree of the collated (AtSCCs and OsSCCs) SCC-encoding biosynthetic genes were divided into 11 different clades of various sizes comprised of branches of subclades. In clade 1, nearing equal representation of OsSCC and AtSCC biosynthetic genes imply the most ancestral lineage. A total of 25 candidate Arabidopsis SCC homologs were identified in rice. The gene ontology enrichment analysis showed that the rice-Arabidopsis SCC homologs were significantly enriched in the following terms at false discovery rate (FDR) < 0.05: (i) biological process; sulfur compound metabolic process and organic acid metabolic processes, (ii) molecular function; oxidoreductase activity, acting on paired donors with incorporation or reduction of molecular oxygen and (iii) KEGG pathway; metabolic pathways and biosynthesis of secondary metabolites. At less than five duplicated blocks of separation, no tandem duplications were observed among the SCC biosynthetic genes distributed in rice chromosomes. The comprehensive rice SCC gene description entailing syntenic events with Arabidopsis, motif distribution and chromosomal mapping of the present findings offer a foundation for rice SCC gene functional studies and advanced strategic rice breeding.
Sulfur (S) is an important macronutrient for plant growth and development, immunity, and stress mitigation. In sulfur-deficient soils, plants invoke stress resistance and xenobiotic detoxification1,2. Plant S assimilation is translated into sulfur-containing compounds (SCCs), a class of important secondary metabolites. Plants utilize freely available sulfate in the soil to synthesize SCCs for growth and functional metabolisms. The primary S assimilation pathway integrates carbon, nitrogen, and S for the synthesis of various SCCs such as glutathione, S-adenosylmethionine, S-methylmethionine, sulfoquinovosyldiacylglycerol, ferredoxin and thiol-group containing plant defensins3. Both the thioredoxins and gluthathiones are redox modulators with detoxifying abilities. In view of the ecological perspective, various vital biological functions which include oxidative stress mitigation, heavy-metal detoxification4,5 and plant defense responses against biotic factors2,6 are regulated by SCCs.
Rice (Oryza sativa L.) is the second most preferred food crop consumed worldwide, after wheat. Cultivated in over 114 countries around the world, rice feeds half the world population (3 billion people) and warrants global food security7. It is predicted that rice production exceeding 800 million tonnes is required to meet the calorie demand of the expected world population in 20258. With climate change in the chart of global issues, abiotic stresses are strongly impacting rice productivity. Major limiting factors in the rice production system includes drought, heat, cold and salinity. In others, waterlogged paddy soils inherent toxic elements such as Cd, As and Fe. Rapid response to stressors regulates stress mitigation responses which include transmembrane transport, glutathione metabolism, signal transduction, and redox control9. In rice, S-associated genes, metabolites and proteins have shown involvement in abiotic stress responses and mitigation. For example, in Cd and As co-contaminated soils, the glutathione metabolism-related genes (Oso1g05367700 and Oso1g0530900) were significantly up-regulated relative to the control conditions. During rice drought stress response, the glutathione S-transferase activities were significant increased10. In another study, glutathione peroxidases and thiol-based antioxidant enzymes regulated the ABA-independent osmotic stress signalling in rice11. Although the role of SCC-encoding genes and SCCs in rice stress response have been documented by numerous studies, little is known about the SCC gene distribution and pattern, and putative functions at the rice genome scale. The SCC genome-level information is important to shed new information and knowledge in innovative rice breeding strategies.
Plant SCC distribution varies greatly with species. In the Brassicaceae family, more than two hundred different types of glucosinolates (GLSs) with potent roles in defense responses have been reported12,13. The GLS-myrosinase defense system gets activated during a pathogen attack to form unstable aglycone intermediates. Thereafter, a range of toxic volatile compounds (isothiocyanates, nitriles, and thiocyanates) is produced during hydrolysis for deterrence against the invading pathogen/pests14. In others, camalexin, an indole-type phytoalexin SCC is produced for adaptivity against abiotic stress and pathogen attack, alike2. Camalexin derived from tryptophan is converted to indole-3-acetaldoxime, which later switches into indole-3-acetonitrile upon dehydration15. Arabidopsis (Brassicaceae) and rice from the grass family (Poaceae) are S-dependent families. With about 10–30% of S expressed in the plant tissues, the first is ranked as the most S-dependent family3,16,17,18,19.
In this study, the SCC-encoding biosynthetic genes in rice are identified and characterized using Arabidopsis as the reference genome model of an S-dependent plant family. The Arabidopsis genome is an excellent reference for the identification of S-encoding biosynthetic genes in rice. There is a burst of SCC-related functional experiments and databases20,21 extensively reported in Arabidopsis; low-affinity sulphate transporters22; S dioxygenase activity in ETHE1 knockout mutant23; S deficiency responsive genes24; Arabidopsis S metabolome25; S-containing secondary metabolites from Arabidopsis2. The synteny and similarity of the Arabidopsis-rice SCC homologous sequences are visualized and the enrichment analysis along a cross-comparison of the corresponding motif sequences is provided to gain information on the extent of similarities. The findings extent to compare and capture the Arabidopsis-rice evolutionary relationship, predict the ecological functions of SCC genes in rice and provide the genetic basis for stress mitigation and defense response enhancement in rice breeding.
Arabidopsis thaliana and O. sativa genome sequences and genome annotations were obtained from the Phytozome v13.0 database (https://phytozome-next.jgi.doe.gov/)26, Arabidopsis Information Resource (TAIR) v10.0 (https://www.arabidopsis.org)27 and O. sativa Genome Annotation Project Database (RGAP) v7.0 (http://rice.uga.edu/)28. The Arabidopsis genome was set as reference sequence against the rice (query) sequences.
The SCC-encoding genes in Arabidopsis (AtSCC) were mined from AraCyc version 14.0 (https://pmn.plantcyc.org/)29 using the following keywords: (i) glucosinolate, and (ii) camalexin. The AtSCC biosynthetic protein sequences were designated as query for the identification of corresponding homologs (O. sativa SCC biosynthetic genes) in the rice genome using the BLAST program (http://blast.ncbi.nlm.nih.gov)30. Reciprocal searching was applied using BLASTP default parameters: e-value = 1e-10 and sequence similarity > 30%. The gene positions were determined by parsing the genome annotation file and the BLAST output. The genomic feature information (General Feature Format) file was concatenated as the input data for subsequent analysis.
The Multiple Collinearity Scan Toolkit X software (MCScanX) was employed for the identification of collinear blocks of homologous sequences and multiple alignment of collinear blocks to the chromosomes. Input files were executed by the MCScan function and the expected number of occurrences (E) of the collinear blocks was calculated31. The following default parameters were applied: E-value cut-off = 1e-05 and match_size = 5. The collinear blocks of interspecies were labelled as AtSCC and OsSCC, denoting A. thaliana and O. sativa, respectively. All rice-Arabidopsis collinear blocks of gene pairs (two interspecies chromosomal positions) were identified and visualized using Rcircos software32.
A multiple sequence alignment of the rice-Arabidopsis SCC-encoding biosynthetic genes was performed using the Multiple Sequence Comparison by Log-Expectation (MUSCLE) (https://www.ebi.ac.uk/Tools/msa/muscle/) with the following settings: gap open penalty = − 2.9, gap extension = 0, and hydrophobicity multiplier = 1.233. Phylogenetic analysis was performed using Molecular Evolutionary Genetics Analysis (MEGA) v7.0 (http://megasoftware.net)34. The maximum-likelihood (ML) by Tamura-Nei substitution model and phylogeny test using 1000 replicates of the bootstrap method were applied. The ML phylogenetic tree was visualized and annotated using the Interactive Tree Of Life (iTOL) v4.0 (http://itol.embl.de)35.
The exon–intron architecture of AtSCC and OsSCC biosynthetic genes was visualized using the Gene Structure Display Server 2.036. Conserved motifs were identified using the Multiple Expectation Maximization for Motif Elicitation (MEME) v4.11.3 (http://meme-suite.org/) tool with the following parameters: the number of motifs = 10, motif site distributions mode = 0/1 occurrence per sequence (zoops)37. The consensus motif sequences were annotated using Database of protein domains, families and functional sites (PROSITE) (https://prosite.expasy.org)38, Pfam, database for protein families v35.0 (http://pfam.xfam.org/)39 and Conserved Domain Database v3.19 (CDD) (https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml))40. The chromosoma gene loci were mapped using the Chromosome Map Tools available in TAIR (https://www.arabidopsis.org/jsp/ChromosomeMap/tool.jsp)41 and Oryzabase (http://viewer.shigen.info/oryzavw/maptool/MapTool.do)42 of A. thaliana and O. sativa genes, respectively. Genes separated by less than five genetic loci within 5 to 100 kb were scored as tandem duplications.
Functional enrichment analysis of the SCC-encoding biosynthetic genes (AtSCC and OsSCC) was performed using ShinyGO v.0.75 (http://bioinformatics.sdstate.edu/go75/) with p-value cut-off set at false discovery rate (FDR) = 0.05: (i) Gene ontology classification43 and (ii) KEGG pathway enrichment44. The A. thaliana and O. sativa Japonica genomes were set as reference datasets. The 20 top-most significantly enriched AtSCC and OsSCC genes were identified using the Venn webserver (https://bioinformatics.psb.ugent.be/webtools/Venn/).
A total of 139 AtSCC biosynthetic genes were obtained from a rapid search performed with the following descriptions: (i) glucosinolate activation (herbivore attack and intact plant cell) pathways, (ii) aliphatic glucosinolate (derived from homomethionine, dihomomethionine, trihomomethionine, hexahomomethionine, pentahomomethionine, and tetrahomomethionine), (iii) indolic glucosinolate (tryptophan derivative), (iv) aromatic glucosinolate (phenylalanine derivative) and (v) camalexin. The sequence homology search identified a total of 838 SCC biosynthetic genes in O. sativa. A total of 173 sequences were discarded due to low sequence similarity (< 30%) and the remaining 665 candidates were subjected to synteny analysis. Under various combinations, a total of 477 syntenic gene pairs with 89 AtSCC and 265 OsSCC biosynthetic genes were identified (Supplementary 1). The syntenic gene pairs were randomly distributed across the chromosomes with sizes, as annotated by the gene number (GN). In rice, the syntenic GN distribution were as following: OsChr1; syntenic GN = 45, OsChr2; syntenic GN = 32, OsChr3; syntenic GN = 28, OsChr11; syntenic GN = 23, OsChr4, OsChr 7 and OsChr 9; syntenic GN = 20, OsChr12; syntenic GN = 18, OsChr8; syntenic GN = 16, and OsChr5; GN = 13. In A. thaliana, the highest number of syntenic genes were distributed in AtChr1 (syntenic GN = 31), followed by AtChr5 (syntenic GN = 19), AtChr3 (syntenic GN = 15), AtChr2 (syntenic GN = 13) and AtChr4 (syntenic GN = 11) (Fig. 1).
The Arabidopsis thaliana (At)-Oryza sativa (Os) sulfur-containing compound (SCC) encoding biosynthetic synteny gene pairs identified with MCScanX. There are 477 syntenic gene pairs (represented by connecting colour lines) between 89 AtSCC and 265 OsSCC biosynthetic genes. The numbering on AT and OS labels denotes the chromosome number.
The distribution of syntenic gene pairs (SPs) was higher in AtChr1 (SPs = 267) and AtChr5 (SPs = 113) in comparison to AtChr2 (SP = 37) and AtChr3 (SP = 33). Overall, a total of 41 AtSCC and 25 OsSCC biosynthetic genes were linked with at least four synteny blocks. Ten OsSCC biosynthetic genes from the cytochrome P450 gene family with at least five or more synteny blocks were identified as following: CYP89D1 (LOC_Os01g24810), CYP706C2 (LOC_Os01g50490), CYP73A35P (LOC_Os01g60450), CYP71AA3 (LOC_Os01g72740), CYP71U3 (LOC_Os02g17760), CYP51H4 (LOC_Os02g21810), CYP73A40 (LOC_Os02g26770), CYP86E1 (LOC_Os02g38290), CYP81A6 (LOC_Os03g55240) and CYP735A4 (LOC_Os09g23820) (Supplementary 1).
The phylogenetic tree comprised of 89 AtSCC and 265 OsSCC biosynthetic genes show 11 different clades of various sizes, as annotated by the gene number (GN). Clade 8 emerged as the largest group with GN = 65, followed by clade 7 (GN = 59), clade 2 (GN = 46), clade 6 (GN = 44), clade 9 (GN = 40), clade 11 (GN = 29), clade 4 (GN = 26), and clade 1 and clade 10 with GN = 20, each. Clade 5 and clade 2 were the smallest in size, with GN = 3 and GN = 2, respectively. There were 7 clades comprised of OsSCC and AtSCC biosynthetic genes in combination: clade 1, clade 4, clade 6, clade 8, clade 9, clade 10 and clade 11. Clade 1 showed nearing an equal number of OsSCC and AtSCC biosynthetic genes. In clade 1, AtNIT2 (At3g44300), AtNIT1 (At3g44310), AtNIT4 (At5g22300) and OsNRT2 (LOC_Os02g42330) were present together (Fig. 2).
Phylogenetic analysis of collated sulfur-encoding biosynthetic genes in Arabidopsis thaliana and rice (Oryza sativa). The tree is constructed with MEGA software.
In clade 4, synteny events between the AtBGLU34 (At1g47600) biosynthetic gene and Os6BGLU24 (LOC_Os06g21570), Os4BGLU9 (LOC_Os04g39814), Os11BGLU37 (LOC_Os11g08120), Os8BGLU27 (LOC_Os08g39860) and Os9BGLU29 (LOC_Os09g31410) biosynthetic genes were identified. Likewise, clade 8 showed a collinear relationship between the OsSOT (LOC_Os09g08190) biosynthetic gene and the AtSOT18 (At1g74090) biosynthetic gene. In clade 9, AtACO9 (At5g43440) was grouped together with Os2ODD25 (LOC_Os03g32470), OsFLS1 (LOC_Os09g18450), Os2ODD16 (LOC_Os01g24980), and Os2ODD26 (LOC_Os03g63900), whilst OsHIS1 (LOC_Os02g17940) was paired with AtACO4 (At1g03400) and AtACO8 (At3g61400). There were three syntenic pairs identified in clade 10: (i) OsCOMTL4 (LOC_Os02g57760)-AtIGMT5 (At1g76790) biosynthetic genes, (ii) OsCOMTL5 (LOC_Os04g09604)-AtIGMT1 (At1g21100) biosynthetic genes and, (iii) OsCOMT (LOC_Os08g06100)-AtIGMT1 (At1g21100) biosynthetic genes. In clade 11, both OsGTF (LOC Os11g04860) and OsIAGLU (LOC Os09g11290) biosynthetic genes were identified as syntenic pairs of AtUGT74B1 (At1g24100). No syntenic evidence was present in clade 6 (Figs. 1 and 2).
A total of ten conserved motifs were identified from A. thaliana and O. sativa SCC biosynthetic genes in clade 1, clade 4, clade 6, clade 8, clade 9, clade 10 and clade 11. The detailed motif sequence information and annotations are provided in Supplementary 2. The motif distribution was similar within the clade level. All the SCC-encoding biosynthetic genes contained at least one motif, whereas a total of 14 genes displayed all 10 motifs with mosaic patterning. No apparent pattern was observed among the motifs within the different species. The following motifs were annotated as isopropyl malate dehydrogenase (IPMDH): motifs 2, 3, 4, 7, 9 and 10. All the SCC-encoding biosynthetic genes in clade 4 displayed motif 6 (annotated as glucosidase) (Supplementary 2). Motif 1 and motif 5 contain the conserved sulfotransferase domain (Fig. 3). In clade 9, at least nine different motifs were consistently present in the member genes. Motifs 1, 2, 4 and 5 are annotated with the O-methyltransferase domain. The conserved motifs 1, 2, 4, 8 and 9 were also described as UDP-glycosyltransferase (Fig. 3) (Supplementary 2).
Motif distribution structure of Arabidopsis thaliana and Oryza sativa sulfur-encoding biosynthetic genes grouped by clades. The A. thaliana (ATXXXXXXX) and O. sativa (OsXXXXXXXX) gene IDs are written in black and red, respectively. Detailed information on the motif sequence information and annotation is available in Supplementary 2.
Generally, the number of exons (EN) and introns (IN) in Arabidopsis and rice displayed no apparent trend by species. Nevertheless, similar exon–intron architecture was observed among the clades of collated AtSCC and OsSCC biosynthetic genes. The number of EN in AtSCC and OsSCC biosynthetic genes ranged from 1 to 13 (Fig. 4). The AtBGLU34 (At1g47600), Os6BGLU24 (LOC_ Os06g21570) and OsBGLU27 (LOC_Os08g39860) biosynthetic genes showed the highest exon and intron distribution with EN = 13 and IN = 12, respectively. There were eight SCC-encoding biosynthetic genes with EN = 1 and EN = 3, followed by seven SCC-encoding biosynthetic genes with EN = 2, five SCC-encoding biosynthetic genes with EN = 4–7, and three SCC-encoding biosynthetic genes with EN = 11.
Illustration of the Arabidopsis thaliana and Oryza sativa sulfur-encoding biosynthetic gene structure. Genes are grouped according to clades. The A. thaliana (ATXXXXXXX) and O. sativa (OsXXXXXXXX) gene IDs are written in black and red, respectively. Exons are indicated as yellow round-corner rectangles and introns with solid black lines.
The exon–intron architecture of AtSCC and OsSCC syntenic gene pairs are described as follows: the AtNIT2-OsNIT2 syntenic gene pair in clade 1 shared a similar number of exons (EN = 5), whereas, in OsIPMDH-AtIMD1/3 syntenic gene pair, a total of 11 exons were distributed in OsIPMDH and about 8–9 exons in AtIMD1 and AtIMD3. In sub-clade 4, the OsBGLU24 and OsBGLU27 biosynthetic genes displayed 13 exons as that of the AtBGLU34 biosynthetic gene except for OsBGLU29, OsBGLU9, and OsBGLU35 (EN = 6–7). The EN in the remaining clades displayed a similar trend; AtSOT18-OsSOT in clade 8 (EN = 1), AtACO4-OsFLS1/Os2ODD25 and AtGSL-OH—Os2ODD25 in clade 9 (EN = 3) (Fig. 4). Four non-syntenic genes with the same exon number are present in clades 6 and 8. The following syntenic gene pairs displayed dissimilarities in the EN: (i) AtACO4-OsHIS1, (ii) AtACO9-Os2ODD16, (iii) AtACO8-OsHIS1, (iv) AtACO8-Os2ODD16, (v) AtIGMT5-OsCOMTL4 and (vi) AtIGMT1-OsCOMTL5. The rice OsHIS1, Os2ODD16, OsCOMTL4 and OsCOMTL5 biosynthetic genes gained one exon, while their syntenic pairs AtACO4, AtACO8, AtACO9, AtIGMT1 and AtIGMT5 lost one exon. Two exon gains were observed in OsIPMDH, OsIPMS1, and OsIPMS2 biosynthetic genes, in contrast to two exon losses in each AtIMD1 and AtMAM1 biosynthetic genes. A total of 25 rice SCC-encoding biosynthetic genes in established synteny and similarity against motifs distributions and exon–intron structure of 18 AtSCC biosynthesis genes (Table 1).
The GO and pathway enrichment analysis of rice and Arabidopsis SCC-encoding biosynthetic genes revealed a total of 206, 149 and 37 hits (terms) in biological process (BP), molecular function (MF) and KEGG pathway, respectively. The number of hit terms, commonly enriched among the rice and Arabidopsis SCC-encoding biosynthetic genes are as follows: BP; 30, MF; 34 and KEGG pathway; 9. In BP, the most significantly enriched terms among the rice SCC-encoding biosynthetic genes are sulfation, hormone biosynthetic process and hormone metabolic process whereas, in Arabidopsis SCC biosynthetic genes, the following terms were significantly enriched: (i) S-glycoside metabolic process, glycosinolate metabolite process and glucosinolate metabolic processes. In both the rice and Arabidopsis SCC biosynthetic genes, sulfur compound metabolic process and organic acid metabolic process were commonly present.
In MF, oxidoreductase activity, acting on paired donors with incorporation or reduction of molecular oxygen was the most significantly enriched (with more than 180 hits) term in both rice and Arabidopsis SCC biosynthetic genes. Other terms enriched at a relatively high extent are as follow: (i) oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen, NAD(P)H as one donor and incorporation, (ii) monooxygenase activity, (iii) iron ion binding, (iv) heme binding, (v) tetrapyrrole binding, (vi) metal ion binding and (vii) N, N-dimethylalanine monooxygenase activity.
The KEGG pathway enrichment showed involvement of the rice-Arabidopsis homologous genes in 10 different signalling pathways. The highest number of genes were significantly enriched in the metabolic pathways and biosynthesis of secondary metabolites with a total number of genes of 80 and 67, respectively. The tryptophan metabolism and 2-oxocarboxylic acid metabolism were fairly high at 25 and 21, respectively (Fig. 6).
Highly conserved SCC biosynthetic genes were physically mapped on the Arabidopsis and rice genomes. The SCC biosynthetic gene distribution in Arabidopsis and rice chromosomes are unequal (Fig. 5). In Arabidopsis, chromosome 1 showed the highest gene number (GN) = 10, followed by chromosome 5 (GN = 8), chromosome 2, (GN = 5), chromosome 3 (GN = 3) and chromosome 4 (GN = 1). The rice SCC biosynthetic genes are distributed in all the 12 chromosomes except chromosomes 5 and 7. Chromosomes 1, 3, 4, 6, 8, 10, 11 and 12 contain one to three OsSCC biosynthetic genes, and the highest number of OsSCC biosynthetic genes (GN = 4) are distributed on chromosomes 2 and 9. No tandem duplications are observed among the SCC biosynthetic genes; no two gene loci are arranged in close proximity and genes are separated by more than five duplicated blocks (Fig. 5). The OsIPMS1 and OsIPMS2 encoding proteins have the longest protein length (635 aa) in rice and AtCYP79B3 (543 aa) in Arabidopsis. OsIAGLU and AtIPMI2 are the shortest protein-encoding gene in rice (113 aa) and Arabidopsis (256 aa), respectively. More than half of the proteins encoded by the SCC biosynthetic genes are acidic, with a theoretical pI (isoelectric point) ranging from 4.63 to 6.24 (Arabidopsis) and 5.1 to 6.8 (rice). The average molecular weight (MW) of AtSCC biosynthetic genes is 45.94 kDa and 43.24 kDa for the OsSCC biosynthetic genes (Table 2).
Sulfur-containing compound (SCC) encoding biosynthetic gene distribution in A. thaliana and O. sativa chromosomes. Grey bars represent the physical maps. The chromosomes are numbered accordingly: A. thaliana;1–5 and O. sativa;1–12. Short lines on grey bars represent the locations of SCCs biosynthetic genes (labelled in red) on each physical map. The different colour boxes expressed adjacent to the gene ID represent the clades.
Sulfur (S) is a secondary macronutrient that regulates plant physiology, growth and developmental processes such as photosynthesis, biosynthesis of sulfur-containing compounds (SCCs) and hormone biosynthesis. It is the 4th major nutrient for crop production after nitrogen, phosphorus and potassium. In higher plants, the S acquisition and assimilation consumes high energy. The S element is taken up by plants as sulphate ions mainly via roots and a small amount can be absorbed through leaves. In rice, the S element, S-containing genes and associated SCCs are critically involved in stress-responsive mechanisms45.
For example, the glutathione S-transferase (GST), a detoxification enzyme ubiquitously present in vertebrates and invertebrates plays an important role in xenobiotic compound detoxification. GST activity is associated with oxidative stress protection as it acts as a mediating substrate in various biochemical reactions, interacts with phytohormones and redox metabolites, and coordinates stress-induced signalling events10. Glutathione (GSH) mediates abiotic and biotic stress resistance using the ROS-scavenging mechanism of the first defense line system in crop plants46. Extensive studies have evident GSH-mediated tolerance mechanisms against salinity, drought, heavy metal toxicity, chilling and herbicides in rice, wheat, barley, soybean and canola47. The effect of S amendment on plant defense response had contributed to similar evidence. As such, the soil amendment of S-containing fertilizer on wheat varieties increased resistance against brown rust and improved the overall productivity48.
Rice yield-impeding factors include pest and pathogen, climate, weather, soil infertility, heavy metal contamination and others. Presently, rice yield enhancement strategies are vigorously carried out by tapping into various aspects of rice biology. Genetic studies, molecular breeding, genetic engineering, heterosis breeding and population improvement are amongst the most sought-after tools utilized in modern rice breeding49,50,51. Since a large number of studies on rice S and SCCs have been linked to stress mechanisms and defense responses, a comprehensive annotation of SCC-encoding genes in the rice genome is important to necessitate enhanced manipulation strategies in breeding approaches52,53,54,55,56.
In this study, a total of 665 OsSCC biosynthetic genes were identified as the homologs of AtSCC query sequences. A total of 477 syntenic gene pairs (Arabidopsis-rice) and 25 rice SCC biosynthetic genes (AtSCC homologs) were obtained using a comprehensive analysis entailing synteny, phylogenetic, conserved motif distribution and gene structure. The synteny analysis identified the gene order and compared the genomic structural changes of the target genes. Shared synteny assumes a common ancestor/evolutionary origin and a syntenic fragment shares a similar function57,58. A small number of genes identified as Arabidopsis-rice syntenies , suggests the early Angiosperm divergence of monophyletic monocot from its eudicot relatives59. The monocot rice genome with 5 chromosomes typically diverged from the eudicot Arabidopsis genome (7 chromosomes) of a higher chromosome number60. The synteny analysis of Arabidopsis-rice SCC biosynthetic genes implies the ancient existence of SCC biosynthetic genes, even before the divergence of the Arabidopsis-rice (eudicot-monocot).
The SCC biosynthetic gene distribution pattern suggests the occurrence of an expansion event during evolution which could have possibly gone through gene co-localization or inter-chromosomal translocation61. The phylogenetic and gene structure pattern of the SCC-encoding biosynthetic genes suggest exon loss and gain events during Arabidopsis-rice (eudicot-monocot) evolution. The exon–intron arrangement pattern in 25 AtSCC and 18 OsSCC suggests that the species-specific genome features are conserved62. The mosaic patterning of the SCC gene exon–intron regions could be associated with evolutionary forces that shaped the SCC biosynthetic gene structure dynamics.
Motifs are frequently occurring (conserved) regions within a DNA sequence. Found within the regulatory regions such as promoters and 3î UTRs, the 4–10 base pair motifs carry significant genome regulatory functions. Two species are likely to be close relatives if they share a high content of common motifs63. During speciation, mutations lead to either an accumulation or loss of motifs (motif turnover) and thus, a motif content analysis is often regarded as more advantageous than the counterpart sequence similarity search analysis. Our results showed that at least 10 different motifs identified in the Arabidopsis and rice SCC-encoding biosynthetic genes have similar distribution patterns by clades.
For instance, in clade 1, six motifs were annotated as 3-isopropylmalate dehydrogenase despite differences in the DNA and protein sequences. Likewise in Clade 4, about 7 different motifs are annotated as glycosyl hydrolase family 1 whereas, in Clade 10, there are 4 motifs corresponding to O-methyltransferase domain (Supplementary 2). The OsSCC biosynthetic genes identified in this study showed potential functional roles in plant defense response. In clade 1, LOC_Os02g42330 (nitrilase 1), the syntenic pair of At3g44300 (nitrilase 2) was reported to participate in the tryptophan-dependent pathway of auxin biosynthesis in rice64. Three OsSCC biosynthetic genes from clade 10 were characterized as O-methyltransferase, a key gene in Arabidopsis indolic glucosinolate modification. As shown in Table 1, five -glucosidase genes from clade 4 showed syntenies with glucosidase 34 (AtBGLU34). AtBGLU34 plays a major role in response to salt stress65 and indolic glucosinolate biosynthesis66 in Arabidopsis.
The SCC biosynthetic genes distributed among the unique phylogenetic clades, carrying similar motif pattern are possibly sharing a similar function. The unique motifs in each clade could be associated with specific functional roles of the SCC biosynthetic genes. The current findings shed insights on the potential functional roles of SCC biosynthetic genes in rice as more than half of the genes were putatively involved in the biosynthesis of aliphatic glucosinolate and indolic glucosinolate. Based on the gene ontology and pathway enrichment analysis, the Arabidopsis-rice homologous SCC-encoding genes were significantly enriched in the sulfur compound metabolic process (BP), oxidoreductase activity, acting on paired donors with incorporation or reduction of molecular oxygen (MF) and biosynthesis of secondary metabolites (KEGG pathway) (Fig. 6). This may suggest the role of the SCC-encoding genes in S assimilation, whereby the reduction of sulphate ion to sulphide and subsequent S-containing amino acids (methionine and cysteine) via the adenosine phosphosulphate pyrophosphate (APS) and phosphoadenosine phosphosulphate (PAPS) is catalyzed by the participating enzyme activities.
Gene ontology (GO) and pathway enrichment analysis. The bubble plot represents the top 20 significantly enriched terms of the Arabidopsis-rice homologous SCC-encoding geens. The GO terms are presented in (i-ii) biological process and (iii-iv) molecular functions whereas the KEGG pathways are presented in (v-vii). Red arrows represent the terms shared among the Arabidopsis-rice orthologous genes. The results are visualized at P < 0.05 using ShinyGO v0.75 (http://bioinformatics.sdstate.edu/go75/).
In plant breeding strategies, exploiting the naturally occurring genetic variation is of utmost fundamental in controlling genes of agronomic importance. Physical maps of rice SCC biosynthetic genes provided in this study could be harnessed for chromosomal region manipulated breeding techniques such as the target chromosome-segment substitution67 and hotspot chromosomal regional positioning of desirable candidate genes68. The findings enable the selection of desirable target rice genes which are tightly linked to S and SCC-encoding genes with a putative functional role in stress response mechanisms.
Rice SCCs biosynthetic genes show syntenic associations with Arabidopsis homologs (AtSCCs). The high degree of conservation between the AtSCC and OsSCC genes suggests long conservation history which could be implicated in SCC gene functions in plant defense response. The present findings not only identified the rice SCC-encoding genes (OsSCC) but also stretch further to include chromosomal level-mapping to better inform new directions in rice functional research and breeding manipulation strategies.
All open-source genomic datasets analysed in this study are available in the Phytozome v13.0 database (https://phytozome-next.jgi.doe.gov/), Arabidopsis Information Resource v10.0 (TAIR) (https://www.arabidopsis.org) and O. sativa Genome Annotation Project Database v7.0 (RGAP) (http://rice.uga.edu/).
Zhang, N. et al. Sulfur deficiency exacerbates phytotoxicity and residues of imidacloprid through suppression of thiol-dependent detoxification in lettuce seedlings. Environ. Pollut. 291, 118221 (2021).
CAS PubMed Article Google Scholar
Bednarek, P. Sulfur-containing secondary metabolites from arabidopsis thaliana and other Brassicaceae with function in plant immunity. ChemBioChem 13(13), 1846–1859 (2012).
CAS PubMed Article Google Scholar
Aarabi, F., Naake, T., Fernie, A. R. & Hoefgen, R. Coordinating sulfur pools under sulfate deprivation. Trends Plant Sci. 25, 12 (2020).
Yadav, P. & Srivastava, S. Effect of thioura application on root, old leaf and young leaf of two contrasting rice varities (Oryza sativa L.) grown in arsenic contaminated soil. Environ. Technol. Innov. 21, 101368 (2021).
Kopriva, S., Calderwood, A., Weckoppp, S. C. & Koprivova, A. Plant sulfur and big data. Plant Sci. 241, 1–10 (2015).
CAS PubMed Article Google Scholar
Capaldi, F. R. et al. Sulfur metabolism and stress defense responses in plants. Trop. Plant Biol. 8, 60–73 (2015).
ADS CAS Article Google Scholar
Chen, D. et al. N6-methyladenosine methylation analysis reveals transcriptome-wide expression response to salt stress in rice roots. Environ. Exp. Bot. 201, 104945 (2022).
Zhang, H., Li, Y. Y. & Zhu, J. K. Developing naturally stress-resistant crops for a sustainable agriculture. Nat. Plants 4(12), 989–996 (2018).
Huang, Y. et al. Transcriptomic (RNA-seq) analysis of genes responsive to both cadmium and arsenic stress in rice root. Sci. Total Environ. 666, 445–460 (2019).
ADS CAS PubMed Article Google Scholar
Hasanuzzaman, M. et al. Exogenous proline and glycine betaine mediated upregulation of antioxidant defense and glyoxalase systems provides better protection against salt-induced oxidative stress in two rice (Oryza sativa L.) varieties. Biomed. Res. Int. 2014, 757219 (2014).
PubMed PubMed Central Google Scholar
Zhou, H. et al. Rice glutathione peroxidase1-mediated oxidation of bZIP68 positively regulates ABA-independent osmotic stress. Mol. Plant 15, 651–670 (2022).
CAS PubMed Article Google Scholar
Agerbirk, N. & Olsen, C. E. Glucosinolate structures in evolution. Phytochemistry 77, 16–45 (2012).
CAS PubMed Article Google Scholar
Kunstler, A., Gullner, G., Adam, A. L., Nagy, J. K. & Kiraly, L. The versatile roles of sulfur-containing biomolecules in plant defense-a road to disease resistance. Plants 9, 1705 (2020).
CAS PubMed Central Article Google Scholar
Textor, S. & Gershenzon, J. Herbivore induction of the glucosinolate-myrosinase defense system: Major trends, biochemical bases and ecological significance. Phytochem. Rev. 8, 149–170 (2009).
Pedras, M. S. C., Yaya, E. E., Glawischnig, E. & Links, D. A. The phytoalexins from cultivated and wild crucifers: Chemistry and biology. Nat. Prod. Rep. 28(8), 1381–1405 (2011).
CAS PubMed Article Google Scholar
Bell, L. et al. Taste and flavor perceptions of glucosinolates, isothiocyanates, and related compounds. Mol. Nutr. Food Res. 62(18), e1700990 (2018).
PubMed Article CAS Google Scholar
Beck, T. K., Jensen, S., Bjoern, G. K. & Kidmose, U. The masking effect of sucrose on perception of bitter compounds in brassica vegetables. J. Sens. Stud. 29(3), 190–200 (2014).
Hawkesford, M. et al. Functions of macronutrients. In Marschnerís Mineral Nutrition of Higher Plants 3rd edn (ed. Marschner, P.) 135–189 (Academic Press, 2012).
Burow, M., Wittstock, U. & Gershenzon, J. Sulfur-containing secondary metabolites and their role in plant defense. In Sulfur Metabolism in Phototrophic Organisms (eds Hell, R. et al.) 201–222 (Springer, 2008).
Harun, S., Abdullah-Zawawi, M. R., Goh, H. H. & Mohamed-Hussein, Z. A. A comprehensive gene inventory for glucosinolate biosynthetic pathway in Arabidopsis thaliana. J. Agric. Food Chem. 68(28), 7281–7297 (2020).
CAS PubMed Article Google Scholar
Ashari, K. S., Abdullah-Zawawi, M. R. & Harun, S. Reconstruction of the transcriptional regulatory network in arabidopsis thaliana aliphatic glucosinolate biosynthetic pathway. Sains Malays. 47(12), 2993–3002 (2018).
Zuber, H. et al. Sultr4 mutant seeds of Arabidopsis have an enhanced sulphate content and modified proteome suggesting metabolic adaptations to altered sulphate compartmentalization. BMC Plant Biol. 10, 78 (2010).
PubMed PubMed Central Article CAS Google Scholar
Meghan, M. et al. Arabidopsis ETHE1 encodes a sulfur dioxygenase that is essential for embryo and endosperm development. Plant Physiol. 160(1), 226–236 (2012).
Wawrzynska, A. & Sirko, A. To control and to be controlled: Understanding the Arabidopsis SLIM1 function in sulfur deficiency through comprehensive investigation of the EIL protein family. Front Plant Sci. 5, 575 (2014).
PubMed PubMed Central Google Scholar
Glaser, K. et al. Exploring the Arabidopsis sulfur metabolome. Plant J. 77(1), 31–45 (2013).
PubMed Article CAS Google Scholar
Goodstein, D. M. et al. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 40, D1178–D1186 (2012).
CAS PubMed Article Google Scholar
Huala, E. et al. The Arabidopsis information resource (TAIR): A comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant. Nucleic Acids Res. 29(1), 102–105 (2001).
CAS PubMed PubMed Central Article Google Scholar
Ouyang, S. et al. The TIGR rice genome annotation resource: Improvements and new features. Nucleic Acids Res. 35, D883–D887 (2007).
CAS PubMed Article Google Scholar
Mueller, L. A., Zhang, P. & Rhee, S. Y. AraCyc: A biochemical pathway database for Arabidopsis. Plant Physiol. 13(2), 453–460 (2015).
Boratyn, G. M. et al. BLAST: A more efficient report with usability improvements. Nucleic Acids Res. 41, W29–W33 (2013).
PubMed PubMed Central Article Google Scholar
Wang, Y. et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40(7), e49 (2012).
ADS CAS PubMed PubMed Central Article Google Scholar
Zhang, H., Meltzer, P. & Davis, S. RCircos: An R package for circos 2D track plots. BMC Bioinform. 14, 244 (2013).
Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32(5), 1792–1797 (2004).
CAS PubMed PubMed Central Article Google Scholar
Kumar, S., Stecher, G. & Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33(7), 1870–1874 (2016).
CAS PubMed PubMed Central Article Google Scholar
Letunic, I. & Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245 (2016).
CAS PubMed PubMed Central Article Google Scholar
Hu, B. et al. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 31(8), 1296–1297 (2014).
PubMed PubMed Central Article Google Scholar
Bailey, T. L. et al. MEME suite: Tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208 (2009).
CAS PubMed PubMed Central Article Google Scholar
Sigrist, C. J. A. et al. PROSITE, a protein domain database for functional characterization and annotation. Nucleic Acids Res. 38, D161–D166 (2010).
CAS PubMed Article Google Scholar
Finn, R. D. et al. Pfam: The protein families database. Nucleic Acids Res. 42, D222–D230 (2014).
CAS PubMed Article Google Scholar
Marchler-Bauer, A. et al. CDD: NCBIís conserved domain database. Nucleic Acids Res. 43, D222–D226 (2015).
CAS PubMed Article Google Scholar
Rhee, S. Y. et al. The Arabidopsis information resource (TAIR): A model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Nucleic Acids Res. 31(1), 224–228 (2003).
CAS PubMed Article Google Scholar
Kurata, N. & Yamazaki, Y. Oryzabase. An integrated biological and genome information database for rice. Plant Physiol. 140(1), 12–17 (2006).
CAS PubMed PubMed Central Article Google Scholar
Ge, S. X., Jung, D. & Yao, R. ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 36(8), 2628–2629 (2020).
CAS PubMed Article Google Scholar
Kanehisa, M. & Goto, S. KEGG: Kyoto encyclopedia of genes and Genomes. Nucleic Acids Res. 28, 27–30 (2000).
CAS PubMed PubMed Central Article Google Scholar
Mridha, D. & Roychowdhury, T. Effect of sulfate application on inhibition of arsenic bioaccumulation in rice (Oryza sativa L.) with consequent health risk assessment of cooked rice arsenic on human: A pot to plate study. Environ. Pollut. 293, 118561 (2022).
CAS PubMed Article Google Scholar
Cao, F. et al. Exogenous glutathione-mediated abiotic stress tolerance in plants. In Glutathione in Plant Growth, Development, Vol. 18 (Springer International Publishing AG, 2017).
Koh, Y. S. et al. Mitigation of environmental stress-impacts in plants: Role of sole and combinatory exogenous application of glutathione. Front. Plant Sci. 12, 791205 (2021).
PubMed PubMed Central Article Google Scholar
Kurmanbayeva, M. et al. Influence of new sulfur-containig fertilizers on performance of wheat yield. Saudi J. Biol. Sci. 28(8), 4644–4655 (2021).
CAS PubMed PubMed Central Article Google Scholar
Asmuni, M. I., Ismail, A. & Abd Aziz, S. N. Morpho-physiological responses of rice towards submergence tolerance. Int. J. Agric. Biol. 22(1), 35–42 (2019).
Nadarajah, K., Abdul Hamid, N. W. & Abdul Rahman, N. S. SA-mediated regulation and control of abiotic stress tolerance in rice. Int. J. Mol. Sci. 22(11), 5591 (2021).
CAS PubMed PubMed Central Article Google Scholar
Mishra, A. et al. Broad spectrum blast resistance alleles in newly developed Malaysian rice (Oryza sativa L.) genotypes. Euphytica 217, 8 (2021).
Dwarika, B., Sandeep, K. & Thandiwe, N. Nitrogen and sulfur fertilizers effects on growth and yield of Brassica carinata in South Dakota. Agron. J. 113, 1945–1960 (2021).
Huang, L. et al. An interventional study of rice for reducing cadmium exposure in a Chinese industrial town. Environ. Int. 122, 301–309 (2019).
CAS PubMed Article Google Scholar
Jia, L. et al. Rice bioinformatics in the genomic era: Status and perspectives. Crop J. 9(3), 609–621 (2021).
Harada, E., Kusano, T. & Sano, H. Differential expression of genes encoding enzymes involved in sulfur assimilation pathways in response to wounding and jasmonate in Arabidopsis thaliana. J. Plant Physiol. 156(2), 272–276 (2000).
Wermter, N. S., Rohn, S. & Hanchen, F. S. Seasonal variation of glucosinolate hydrolysis products in commercial white and red cabbages (Brassica oleracea var. capitate). Foods. 9, 1682 (2020).
CAS PubMed Central Article Google Scholar
Tang, H. et al. Synteny and collinearity in plant genomes. Science 320, 486–488 (2008).
ADS CAS PubMed Article Google Scholar
Lyons, E. et al. Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids. Plant Physiol. 148, 1772–1781 (2008).
CAS PubMed PubMed Central Article Google Scholar
Hertweck, K. L. et al. Phylogenetics, divergence times and diversification from three genomic partitions in monocots. Bot. J. Linn. Soc. 178, 375–393 (2015).
Abrouk, M. et al. Palaeogenomics of plants: Synteny-based modelling of extinct ancestors. Trends Plant Sci. 15(9), 479–487 (2010).
CAS PubMed Article Google Scholar
Moreno-Hagelsieb, G. et al. Transcription unit conservation in the three domains of life: A perspective from Escherichia coli. Trends Genet. 17(4), 175–177 (2001).
CAS PubMed Article Google Scholar
Roy, S. W. & Penny, D. Patterns of intron loss and gain in plants: intron loss-dominated evolution and genome-wide comparison of O. sativa and A. thaliana. Mol. Biol. Evol. 24(1), 171–181 (2007).
CAS PubMed Article Google Scholar
Cserhati, M. Motif content comparison between monocot and dicot species. Genomics Data 3, 128–136 (2015).
PubMed PubMed Central Article Google Scholar
Abu-Zaitoon, Y. M. Phylogenetic analysis of putative genes involved in the tryptophan-dependent pathway of auxin biosynthesis in rice. Appl. Biochem. Biotechnol. 172(5), 2480–2495 (2014).
CAS PubMed Article Google Scholar
Xu, Z. Y. et al. A vacuolar-glucosidase homolog that possesses glucose-conjugated abscisic acid hydrolyzing activity plays an important role in osmotic stress responses in Arabidopsis. Plant Cell 112, 2184–2199 (2012).
Nakano, R. T. et al. PYK10 myrosinase reveals a functional coordination between endoplasmic reticulum bodies and glucosinolates in Arabidopsis thaliana. Plant J. 89(2), 204–220 (2017).
CAS PubMed Article Google Scholar
Zhang, G. et al. Target chromosome-segment substitution: A way to breeding by design in rice. Crop J. 9, 658–668 (2021).
Raza, Q. et al. Meta-analysis of grain iron and zinc associated QTLs identified hotspot chromosomal regions and positional candidate genes for breeding biofortified rice. Plant Sci. 288, 110214 (2019).
CAS PubMed Article Google Scholar
This work was supported by ERGS/1/2013/STG07/UKM/02/3 awarded to ZAMH by the Malaysian Ministry of Higher Education and Ministry of Science, Technology and Innovation. The Ph.D. scholarship obtained by the first author (MRAZ) is funded by Skim Zamalah Universiti Penyelidikan, Universiti Kebangsaan Malaysia.
Institute of Systems Biology (INBIOSIS), Universiti Kebangsaan Malaysia, 43600, Bangi, Selangor, Malaysia
Muhammad-Redha Abdullah-Zawawi, Nisha Govender, Nor Azlan Nor Muhammad, Zamri Zainal & Zeti-Azura Mohamed-Hussein
UKM Medical Molecular Biology Institute (UMBI), UKM Medical Centre, Jalan Yaíacob Latiff, Bandar Tun Razak, 56000, Cheras, Kuala Lumpur, Malaysia
Faculty of Science and Technology, Universiti Kebangsaan Malaysia, 43600, Bangi, Selangor, Malaysia
Norfarhan Mohd-Assaad, Zamri Zainal & Zeti-Azura Mohamed-Hussein
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Conceptualization, M.R.A.Z. and Z.A.M.H.; methodology, investigation, writing-original draft preparation, M.R.A.Z.; validation, writing-review and editing, N.G., N.A.N.M., N.M.A., Z.Z. and Z.A.M.H.; visualization, supervision, writing-revision and improved the draft manuscript, Z.A.M.H.; supervision and funding acquisition. All authors reviewed the manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Abdullah-Zawawi, MR., Govender, N., Muhammad, N.A.N. et al. Genome-wide analysis of sulfur-encoding biosynthetic genes in rice (Oryza sativa L.) with Arabidopsis as the sulfur-dependent model plant. Sci Rep 12, 13829 (2022). https://doi.org/10.1038/s41598-022-18068-0
DOI: https://doi.org/10.1038/s41598-022-18068-0
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.