Developing an SNP dataset for efficiently evaluating soybean germplasm resources using the genome sequencing data of 3,661 soybean accessions.

BACKGROUND: Single nucleotide polymorphism (SNP) markers play significant roles in accelerating breeding and basic crop research. Several soybean SNP panels have been developed. However, there is still a lack of SNP panels for differentiating between wild and cultivated populations, as well as for detecting polymorphisms within both wild and cultivated populations. RESULTS: This study utilized publicly available resequencing data from over 3,000 soybean accessions to identify differentiating and highly conserved SNP and insertion/deletion (InDel) markers between wild and cultivated soybean populations. Additionally, a naturally occurring mutant gene library was constructed by analyzing large-effect SNPs and InDels in the population. CONCLUSION: The markers obtained in this study are associated with numerous genes governing agronomic traits, thus facilitating the evaluation of soybean germplasms and the efficient differentiation between wild and cultivated soybeans. The natural mutant gene library permits the quick identification of individuals with natural mutations in functional genes, providing convenience for accelerating soybean breeding using reverse genetics.

Chromosome-level genome assembly and annotation of rare and endangered tropical bivalve, Tridacna crocea.

Tridacna crocea is an ecologically important marine bivalve inhabiting tropical coral reef waters. High quality and available genomic resources will help us understand the population structure and genetic diversity of giant clams. This study reports a high-quality chromosome-scale T. crocea genome sequence of 1.30 Gb, with a scaffold N50 and contig N50 of 56.38 Mb and 1.29 Mb, respectively, which was assembled by combining PacBio long reads and Hi-C sequencing data. Repetitive sequences cover 71.60% of the total length, and a total of 25,440 protein-coding genes were annotated. A total of 1,963 non-coding RNA (ncRNA) were determined in the T. crocea genome, including 62 micro RNA (miRNA), 58 small nuclear RNA (snRNA), 83 ribosomal RNA (rRNA), and 1,760 transfer RNA (tRNA). Phylogenetic analysis revealed that giant clams diverged from oyster about 505.7 Mya during the evolution of bivalves. The genome assembly presented here provides valuable genomic resources to enhance our understanding of the genetic diversity and population structure of giant clams.

MicroRNA 4407 modulates nodulation in soybean by repressing a root-specific ISOPENTENYLTRANSFERASE (GmIPT3).

MicroRNAs (miRNAs) are important regulators of plant biological processes, including soybean nodulation. One miRNA, miR4407, was identified in soybean roots and nodules. However, the function of miR4407 in soybean is still unknown. MiR4407, unique to soybean, positively regulates lateral root emergence and root structures and represses a root-specific ISOPENTENYLTRANSFERASE (GmIPT3). By altering the expression of miR4407 and GmIPT3, we investigated the role of miR4407 in lateral root and nodule development. Both miR4407 and GmIPT3 are expressed in the inner root cortex and nodule primordia. Upon rhizobial inoculation, miR4407 was downregulated while GmIPT3 was upregulated. Overexpressing miR4407 reduced the number of nodules in transgenic soybean hairy roots while overexpressing the wild-type GmIPT3 or a miR4407-resistant GmIPT3 mutant (mGmIPT3) significantly increased the nodule number. The mechanism of miR4407 and GmIPT3 functions was also linked to autoregulation of nodulation (AON), where miR4407 overexpression repressed miR172c and activated its target, GmNNC1, turning on AON. Exogenous CK mimicked the effects of GmIPT3 overexpression on miR172c, supporting the notion that GmIPT3 regulates nodulation by enhancing root-derived CK. Overall, our data revealed a new miRNA-mediated regulatory mechanism of nodulation in soybean. MiR4407 showed a dual role in lateral root and nodule development.

Comparative genomic analyses reveal the genetic basis of the yellow-seed trait in Brassica napus.

Yellow-seed trait is a desirable breeding characteristic of rapeseed (Brassica napus) that could greatly improve seed oil yield and quality. However, the underlying mechanisms controlling this phenotype in B. napus plants are difficult to discern because of their complexity. Here, we assemble high-quality genomes of yellow-seeded (GH06) and black-seeded (ZY821). Combining in-depth fine mapping of a quantitative trait locus (QTL) for seed color with other omics data reveal BnA09MYB47a, encoding an R2R3-MYB-type transcription factor, as the causal gene of a major QTL controlling the yellow-seed trait. Functional studies show that sequence variation of BnA09MYB47a underlies the functional divergence between the yellow- and black-seeded B. napus. The black-seed allele BnA09MYB47aZY821, but not the yellow-seed allele BnA09MYB47aGH06, promotes flavonoid biosynthesis by directly activating the expression of BnTT18. Our discovery suggests a possible approach to breeding B. napus for improved commercial value and facilitates flavonoid biosynthesis studies in Brassica crops.

A natural non-synonymous single nucleotide polymorphism in GmbHLH113 negates its inhibitory effect on root hair elongation in soybean.

Root hair length (RHL) is an important character that affects nutrient acquisition in plants. The regulatory network in soybean controlling RHL is yet to be fully understood. In this study, we identified a quantitative trait locus (QTL) regulating RHL. One candidate causal gene in this QTL (GmbHLH113), preferentially expressed in root hairs, was annotated as encoding a basic helix-loop-helix transcription factor. In wild soybeans, the allelic type of GmbHLH113 with a glycine in the 13th residue, which was associated with a reduction in RHL, was shown to localize in the nucleus and activate gene transcription. Another allelic type with a single nucleotide polymorphism that resulted in a glutamate in the 13th residue is fixed in cultivated soybeans, and it lost the ability to localize to the nucleus or negatively regulate RHL. The ectopic expression of GmbHLH113 from W05 in Arabidopsis root hairs resulted in shorter RHL and reduced phosphorus (P) accumulation in shoots. Hence, a loss-of-function allele in cultivated soybeans might have been selected during domestication due to its association with a longer RHL and improved nutrient acquisition.

QTL analyses of soybean root system architecture revealed genetic relationships with shoot-related traits.

KEY MESSAGE: The genetic basis of soybean root system architecture (RSA) and the genetic relationship between shoot and RSA were revealed by integrating data from recombinant inbred population grafting and QTL mapping. Variations in root system architecture (RSA) affect the functions of roots and thus play vital roles in plant adaptations and agricultural productivity. The aim of this study was to unravel the genetic relationship between RSA traits and shoot-related traits in soybean. This study characterized RSA variability at seedling stage in a recombinant inbred population, derived from a cross between cultivated soybean C08 and wild soybean W05, and performed high-resolution quantitative trait locus (QTL) mapping. In total, 34 and 41 QTLs were detected for RSA-related and shoot-related traits, respectively, constituting eight QTL clusters. Significant QTL correspondence was found between shoot biomass and RSA-related traits, consistent with significant correlations between these phenotypes. RSA-related QTLs also overlapped with selection regions in the genome, suggesting the cultivar RSA could be a partial consequence of domestication. Using reciprocal grafting, we confirmed that shoot-derived signals affected root development and the effects were controlled by multiple loci. Meanwhile, RSA-related QTLs were found to co-localize with four soybean flowering-time loci. Consistent with the phenotypes of the parental lines of our RI population, diminishing the function of flowering controlling E1 family through RNA interference (RNAi) led to reduced root growth. This implies that the flowering time-related genes within the RSA-related QTLs are actually contributing to RSA. To conclude, this study identified the QTLs that determine RSA through controlling root growth indirectly via regulating shoot functions, and discovered superior alleles from wild soybean that could be used to improve the root structure in existing soybean cultivars.

The whole-genome assembly of an endangered Salicaceae species: Chosenia arbutifolia (Pall.) A. Skv.

BACKGROUND: As a fast-growing tree species, Chosenia arbutifolia has a unique but controversial taxonomic status in the family Salicaceae. Despite its importance as an industrial material, in ecological protection, and in landscaping, C. arbutifolia is seriously endangered in Northeast China because of artificial destruction and its low reproductive capability. RESULTS: To clarify its phylogenetic relationships with other Salicaceae species, we assembled a high-quality chromosome-level genome of C. arbutifolia using PacBio High-Fidelity reads and Hi-C sequencing data, with a total size of 338.93 Mb and contig N50 of 1.68 Mb. Repetitive sequences, which accounted for 42.34% of the assembly length, were identified. In total, 33,229 protein-coding genes and 11,474 small noncoding RNAs were predicted. Phylogenetic analysis suggested that C. arbutifolia and poplars diverged approximately 15.3 million years ago, and a large interchromosomal recombination between C. arbutifolia and other Salicaceae species was discovered. CONCLUSIONS: Our study provides insights into the genome architecture and systematic evolution of C. arbutifolia, as well as comprehensive information for germplasm protection and future functional genomic studies.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs.

BACKGROUND: The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics. RESULTS: We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme (LYZ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type (TNNC1) gene might increase the proportion of slow muscle fibers and affect meat quality. CONCLUSIONS: Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

Multi-omics analysis reveals the mechanism of seed coat color formation in Brassica rapa L.

KEYMESSAGE: Multi-omics analysis of the transcriptome, metabolome and genome identified major and minor loci and candidate genes for seed coat color and explored the mechanism of flavonoid metabolites biosynthesis in Brassica rapa. Yellow seed trait is considered an agronomically desirable trait with great potential for improving seed quality of Brassica crops. Mechanisms of the yellow seed trait are complex and not well understood. In this study, we performed an integrated metabolome, transcriptome and genome-wide association study (GWAS) on different B. rapa varieties to explore the mechanisms underlying the seed coat color formation. A total of 2,499 differentially expressed genes and 116 differential metabolites between yellow and black seeds with strong association with the flavonoid biosynthesis pathway was identified. In addition, 330 hub genes involved in the seed coat color formation, and the most significantly differential flavonoids biosynthesis were detected based on weighted gene co-expression network analysis. Metabolite GWAS analysis using the contents of 42 flavonoids in developing seeds of 159 B. rapa lines resulted in the identification of 1,626 quantitative trait nucleotides (QTNs) and 37 chromosomal intervals, including one major locus on chromosome A09. A combination of QTNs detection, transcriptome and functional analyses led to the identification of 241 candidate genes that were associated with different flavonoid metabolites. The flavonoid biosynthesis pathway in B. rapa was assembled based on the identified flavonoid metabolites and candidate genes. Furthermore, BrMYB111 members (BraA09g004490.3C and BraA06g034790.3C) involved in the biosynthesis of taxifolin were functionally analyzed in vitro. Our findings lay a foundation and provide a reference for systematically investigating the mechanism of seed coat color in B. rapa and in the other plants.

The Poly-Glutamate Motif of GmMATE4 Regulates Its Isoflavone Transport Activity.

Multidrug and toxic compound extrusion (MATE) transporters in eukaryotes have been characterized to be antiporters that mediate the transport of substrates in exchange for protons. In plants, alkaloids, phytohormones, ion chelators, and flavonoids have been reported to be the substrates of MATE transporters. Structural analyses have been conducted to dissect the functional significance of various motifs of MATE proteins. However, an understanding of the functions of the N- and C-termini has been inadequate. Here, by performing phylogenetic analyses and protein sequence alignment of 14 representative plant species, we identified a distinctive N-terminal poly-glutamate motif among a cluster of MATE proteins in soybean. Amongst them, GmMATE4 has the most consecutive glutamate residues at the N-terminus. A subcellular localization study showed that GmMATE4 was localized at the vacuolar membrane-like structure. Protein charge prediction showed that the mutation of the glutamate residues to alanine would reduce the negative charge at the N-terminus. Using yeast as the model, we showed that GmMATE4 mediated the transport of daidzein, genistein, glycitein, and glycitin. In addition, the glutamate-to-alanine mutation reduced the isoflavone transport capacity of GmMATE4. Altogether, we demonstrated GmMATE4 as an isoflavone transporter and the functional significance of the N-terminal poly-glutamate motif of GmMATE4 for regulating the isoflavone transport activity.

Comparative analysis of the value of amide proton transfer-weighted imaging and diffusion kurtosis imaging in evaluating the histological grade of cervical squamous carcinoma.

BACKGROUND: Uterine cervical cancer (UCC) was the fourth leading cause of cancer death among women worldwide. The conventional MRI hardly revealing the microstructure information. This study aimed to compare the value of amide proton transfer-weighted imaging (APTWI) and diffusion kurtosis imaging (DKI) in evaluating the histological grade of cervical squamous carcinoma (CSC) in addition to routine diffusion-weighted imaging (DWI). METHODS: Forty-six patients with CSC underwent pelvic DKI and APTWI. The magnetization transfer ratio asymmetry (MTRasym), apparent diffusion coefficient (ADC), mean diffusivity (MD) and mean kurtosis (MK) were calculated and compared based on the histological grade. Correlation coefficients between each parameter and histological grade were calculated. RESULTS: The MTRasym and MK values of grade 1 (G1) were significantly lower than those of grade 2 (G2), and those parameters of G2 were significantly lower than those of grade 3 (G3). The MD and ADC values of G1 were significantly higher than those of G2, and those of G2 were significantly higher than those of G3. MTRasym and MK were both positively correlated with histological grade (r = 0.789 and 0.743, P <  0.001), while MD and ADC were both negatively correlated with histological grade (r = - 0.732 and - 0.644, P <  0.001). For the diagnosis of G1 and G2 CSCs, AUC (APTWI+DKI + DWI) > AUC (DKI + DWI) > AUC (APTWI+DKI) > AUC (APTWI+DWI) > AUC (MTRasym) > AUC (MK) > AUC (MD) > AUC (ADC), where the differences between AUC (APTWI+DKI + DWI), AUC (DKI + DWI) and AUC (ADC) were significant. For the diagnosis of G2 and G3 CSCs, AUC (APTWI+DKI + DWI) > AUC (APTWI+DWI) > AUC (APTWI+DKI) > AUC (DKI + DWI) > AUC (MTRasym) > AUC (MK) > AUC (MD > AUC (ADC), where the differences between AUC (APTWI+DKI + DWI), AUC (APTWI+DWI) and AUC (ADC) were significant. CONCLUSION: Compared with DWI and DKI, APTWI is more effective in identifying the histological grades of CSC. APTWI is recommended as a supplementary scan to routine DWI in CSCs.

Genome sequencing and population resequencing provide insights into the genetic basis of domestication and diversity of vegetable soybean.

Vegetable soybean is one of the most important vegetables in China, and the demand for this vegetable has markedly increased worldwide over the past two decades. Here, we present a high-quality de novo genome assembly of the vegetable soybean cultivar Zhenong 6 (ZN6), which is one of the most popular cultivars in China. The 20 pseudochromosomes cover 94.57% of the total 1.01 Gb assembly size, with contig N50 of 3.84 Mb and scaffold N50 of 48.41 Mb. A total of 55 517 protein-coding genes were annotated. Approximately 54.85% of the assembled genome was annotated as repetitive sequences, with the most abundant long terminal repeat transposable elements. Comparative genomic and phylogenetic analyses with grain soybean Williams 82, six other Fabaceae species and Arabidopsis thaliana genomes highlight the difference of ZN6 with other species. Furthermore, we resequenced 60 vegetable soybean accessions. Alongside 103 previously resequenced wild soybean and 155 previously resequenced grain soybean accessions, we performed analyses of population structure and selective sweep of vegetable, grain, and wild soybean. They were clearly divided into three clades. We found 1112 and 1047 genes under selection in the vegetable soybean and grain soybean populations compared with the wild soybean population, respectively. Among them, we identified 134 selected genes shared between vegetable soybean and grain soybean populations. Additionally, we report four sucrose synthase genes, one sucrose-phosphate synthase gene, and four sugar transport genes as candidate genes related to important traits such as seed sweetness and seed size in vegetable soybean. This study provides essential genomic resources to promote evolutionary and functional genomics studies and genomically informed breeding for vegetable soybean.

Erratum: The Chromosome-Scale Genome of Melon Dissects Genetic Architecture of Important Agronomic Traits.

[This corrects the article DOI: 10.1016/j.isci.2020.101422.].

The genome variation and developmental transcriptome maps reveal genetic differentiation of skeletal muscle in pigs.

Natural and artificial directional selections have resulted in significantly genetic and phenotypic differences across breeds in domestic animals. However, the molecular regulation of skeletal muscle diversity remains largely unknown. Here, we conducted transcriptome profiling of skeletal muscle across 27 time points, and performed whole-genome re-sequencing in Landrace (lean-type) and Tongcheng (obese-type) pigs. The transcription activity decreased with development, and the high-resolution transcriptome precisely captured the characterizations of skeletal muscle with distinct biological events in four developmental phases: Embryonic, Fetal, Neonatal, and Adult. A divergence in the developmental timing and asynchronous development between the two breeds was observed; Landrace showed a developmental lag and stronger abilities of myoblast proliferation and cell migration, whereas Tongcheng had higher ATP synthase activity in postnatal periods. The miR-24-3p driven network targeting insulin signaling pathway regulated glucose metabolism. Notably, integrated analysis suggested SATB2 and XLOC_036765 contributed to skeletal muscle diversity via regulating the myoblast migration and proliferation, respectively. Overall, our results provide insights into the molecular regulation of skeletal muscle development and diversity in mammals.

A high-quality genome assembly of the eggplant provides insights into the molecular basis of disease resistance and chlorogenic acid synthesis.

The eggplant (Solanum melongena L.) is one of the most important Solanaceae crops, ranking third for total production and economic value in its genus. Herein, we report a high-quality, chromosome-scale eggplant reference genome sequence of 1155.8 Mb, with an N50 of 93.9 Mb, which was assembled by combining PacBio long reads and Hi-C sequencing data. Repetitive sequences occupied 70.1% of the assembly length, and 35,018 high-confidence protein-coding genes were annotated based on multiple sources. Comparative analysis revealed 646 species-specific families and 364 positive selection genes, conferring distinguishing traits on the eggplant. We performed genome-wide comparative identification of disease resistance genes and discovered an expanded gene family of bacterial spot resistance in eggplant and pepper, but not in tomato and potato. The genes involved in chlorogenic acid synthesis were comprehensively characterized. Highly similar chromosomal distribution patterns of polyphenol oxidase genes were observed in the eggplant, tomato, and potato genomes. The eggplant reference genome sequence will not only facilitate evolutionary studies of the Solanaceae but also facilitate their breeding and improvement.

Chromosome-level genome assembly for the largemouth bass Micropterus salmoides provides insights into adaptation to fresh and brackish water.

Largemouth bass (LMB; Micropterus salmoides) has been an economically important fish in North America, Europe, and China. This study obtained a chromosome-level genome assembly of LMB using PacBio and Hi-C sequencing. The final assembled genome is 964 Mb, with contig N50 and scaffold N50 values of 1.23 Mb and 36.48 Mb, respectively. Combining with RNA sequencing data, we annotated a total of 23,701 genes. Chromosomal assembly and syntenic analysis proved that, unlike most Perciformes with the popular haploid chromosome number of 24, LMB has only 23 chromosomes (Chr), among which the Chr1 seems to be resulted from a chromosomal fusion event. LMB is phylogenetically closely related to European seabass and spotted seabass, diverging 64.1 million years ago (mya) from the two seabass species. Eight gene families comprising 294 genes associated with ionic regulation were identified through positive selection, transcriptome and genome comparisons. These genes involved in iron facilitated diffusion (such as claudin, aquaporins, sodium channel protein and so on) and others related to ion active transport (such as sodium/potassium-transporting ATPase and sodium/calcium exchanger). The claudin gene family, which is critical for regulating cell tight junctions and osmotic homeostasis, showed a significant expansion in LMB with 27 family members and 68 copies for salinity adaptation. In summary, we reported the first high-quality LMB genome, and provided insights into the molecular mechanisms of LMB adaptation to fresh and brackish water. The chromosome-level LMB genome will also be a valuable genomic resource for in-depth biological and evolutionary studies, germplasm conservation and genetic breeding of LMB.

The Chromosome-Scale Genome of Melon Dissects Genetic Architecture of Important Agronomic Traits.

Comparative and evolutionary genomics analyses are the powerful tools to provide mechanistic insights into important agronomic traits. Here, we completed a chromosome-scale assembly of the "neglected" but vital melon subspecies Cucumis melo ssp. agrestis using single-molecule real-time sequencing, Hi-C, and an ultra-dense genetic map. Comparative genomics analyses identified two targeted genes, UDP-sugar pyrophosphorylase and α-galactosidase, that were selected during evolution for specific phloem transport of oligosaccharides in Cucurbitaceae. Association analysis of transcriptome and the DNA methylation patterns revealed the epigenetic regulation of sucrose accumulation in developing fruits. We constructed the melon recombinant inbred lines to uncover Alkaline/Neutral Invertase (CINV), Sucrose-Phosphatase 2 (SPP2), α-galactosidase, and ß-galactosidase loci related to sucrose accumulation and an LRR receptor-like serine/threonine-protein kinase associated with gummy stem blight resistance. This study provides essential genomic resources enabling functional genomics studies and the genomics-informed breeding pipelines for improving the fruit quality and disease resistance traits.

The chromosome-level draft genome of Dalbergia odorifera.

BACKGROUND: Dalbergia odorifera T. Chen (Fabaceae) is an International Union for Conservation of Nature red-listed tree. This tree is of high medicinal and commercial value owing to its officinal, insect-proof, durable heartwood. However, there is a lack of genome reference, which has hindered development of studies on the heartwood formation. FINDINGS: We presented the first chromosome-scale genome assembly of D. odorifera obtained on the basis of Illumina paired-end sequencing, Pacific Biosciences single-molecule real-time sequencing, 10x Genomics linked reads, and Hi-C technology. We assembled 97.68% of the 653.45 Mb D. odorifera genome with scaffold N50 and contig sizes of 56.16 and 5.92 Mb, respectively. Ten super-scaffolds corresponding to the 10 chromosomes were assembled, with the longest scaffold reaching 79.61 Mb. Repetitive elements account for 54.17% of the genome, and 30,310 protein-coding genes were predicted from the genome, of which ∼92.6% were functionally annotated. The phylogenetic tree showed that D. odorifera diverged from the ancestor of Arabidopsis thaliana and Populus trichocarpa and then separated from Glycine max and Cajanus cajan. CONCLUSIONS: We sequence and reveal the first chromosome-level de novo genome of D. odorifera. These studies provide valuable genomic resources for the research of heartwood formation in D. odorifera and other timber trees. The high-quality assembled genome can also be used as reference for comparative genomics analysis and future population genetic studies of D. odorifera.

SMRT sequencing of the Oryza rufipogon genome reveals the genomic basis of rice adaptation.

Asian cultivated rice is believed to have been domesticated from a wild progenitor, Oryza rufipogon, offering promising sources of alleles for world rice improvement. Here we first present a high-quality chromosome-scale genome of the typical O. rufipogon. Comparative genomic analyses of O. sativa and its two wild progenitors, O. nivara and O. rufipogon, identified many dispensable genes functionally enriched in the reproductive process. We detected millions of genomic variants, of which large-effect mutations could affect agronomically relevant traits. We demonstrate how lineage-specific expansion of gene families may have contributed to the formation of reproduction isolation. We document thousands of genes with signatures of positive selection that are mainly involved in the reproduction and response to biotic- and abiotic stresses. We show that selection pressures may serve as forces to govern substantial genomic alterations that form the genetic basis of rapid evolution of mating and reproductive systems under diverse habitats.

[Glucosides of chaenomeles speciosa attenuate ischemia/reperfusion-induced brain injury by regulating NF-κB P65/TNF-α in mouse model].

OBJECTIVE: To investigate the effect and mechanism of glucosides of chaenomeles speciosa (GCS) on ischemia/reperfusion-induced brain injury in mouse model. METHODS: Fifty 8-week C57BL/C mice were randomly divided into five groups with 10 in each group:sham group, model group, GCS 30 mg/kg group, GCS 60 mg/kg group and GCS 90 mg/kg group, and the GCS was administrated by gavage (once a day) for 14 d. HE staining was performed to investigate the cell morphology; the Zea-Longa scores were measured for neurological activity; TUNEL staining was performed to investigate the cell apoptosis; ELISA was used to detected the oxidative stress and inflammation; Western Blot was performed to investigate the key pathway and neurological functional molecules. RESULTS: Compared with the sham group, the brain tissues in model group were seriously damaged, presenting severe cell apoptosis, oxidative stress and inflammation, associated with increased NF-κB P65 and TNF-α levels as well as decreased myelin associate glycoprotein (MAG) and oligodendrocyte-myelin glycoprotein (OMgp)levels (all P<0.01). Compared with the model group, the brain tissues in GCS groups were ameliorated, and cell apoptosis, oxidative stress and inflammation were inhibited, associated with decreased NF-κB P65 and TNF-α levels as well as increased MAG and OMgp levels (all P<0.01), which were more markedly in GCS 60 mg/kg group. CONCLUSIONS: GCS can inhibit the NF-κB P65 and TNF-α, reduce the oxidative stress and inflammation, decrease the cell apoptosis in mouse ischemia/reperfusion-induced brain injury model, and 60 mg/kg GCS may be the optimal dose.