Search

The intensive application of inorganic nitrogen underlies marked increases in crop production, but imposes detrimental effects on ecosystems1,2: it is therefore crucial for future sustainable agriculture to improve the nitrogen-use efficiency of crop plants. Here we report the genetic basis of nitrogen-use efficiency associated with adaptation to local soils in rice (Oryza sativa L.). Using a panel of diverse rice germplasm collected from different ecogeographical regions, we performed a genome-wide association study on the tillering response to nitrogen-the trait that is most closely correlated with nitrogen-use efficiency in rice-and identified OsTCP19 as a modulator of this tillering response through its transcriptional response to nitrogen and its targeting to the tiller-promoting gene DWARF AND LOW-TILLERING (DLT)3,4. A 29-bp insertion and/or deletion in the OsTCP19 promoter confers a differential transcriptional response and variation in the tillering response to nitrogen among rice varieties. The allele of OsTCP19 associated with a high tillering response to nitrogen is prevalent in wild rice populations, but has largely been lost in modern cultivars: this loss correlates with increased local soil nitrogen content, which suggests that it might have contributed to geographical adaptation in rice. Introgression of the allele associated with a high tillering response into modern rice cultivars boosts grain yield and nitrogen-use efficiency under low or moderate levels of nitrogen, which demonstrates substantial potential for rice breeding and the amelioration of negative environment effects by reducing the application of nitrogen to crops.

Sepsis-induced acute kidney injury (AKI) continues to be associated with poor outcomes in critical care patients. Previous research has revealed that tetrahydrocurcumin (THC) exerts renoprotective effects in multiple nephritic disorders by modulating inflammation and oxidative stress. However, the effects of THC on sepsis-induced AKI and the underlying mechanisms remain unclear. In this study, a mouse model of sepsis-induced AKI, generated by cecal ligation and puncture operation, was used to investigate the protective effects of THC and the role of SIRT1. Histological manifestation and TUNEL analysis were observed to determine the severity of kidney damage. Levels of BUN, SCr, KIM-1, and UAlb/Cr were calculated to assess the renal function. Expressions of IL-1β, IL-6, and TNF-α were measured to evaluate the inflammatory response. MDA content, SOD, GSH, CAT, and GPx activities and DHE staining were analyzed to estimate the degree of oxidative stress. Protein expressions of SIRT1, Ac-p65, and Ac-foxo1 were detected to explore the underlying mechanisms. We observed that THC not only increased the survival rate, improved the kidney function and ameliorated the renal histological damage of septic mice, but also inhibited inflammatory response, prohibited oxidative stress, and prevented cell apoptosis in renal tissues in septic mice. Mechanistically, THC remarkably increased the expression of SIRT1, accompanied by decreased expressions of downstream molecules Ac-p65 and Ac-foxo1. Meanwhile, the beneficial effects of THC were clearly abolished by the SIRT1-specific inhibitor EX527. These results delineate that THC prevents sepsis-induced AKI by suppressing inflammation and oxidative stress through activating the SIRT1 signaling. Abbreviation: Ac-p65: acetylated p65; Ac-foxo 1: acetylated forkhead box O1; AKI: acute kidney injury; BUN: blood urea nitrogen; CAT: catalase; DHE: dihydroethidium; GPx: glutathione peroxidase; GSH: reduced glutathione; IL-1β: Interleukin-1 beta; IL-6: Interleukin-6; KIM-1: kidney injury molecule 1; MDA: malondialdehyde; SCr: serum creatinine; SIRT1: silent information regulator 1; SOD: superoxide dismutase; THC: tetrahydrocurcumin; TNF-α: tumor necrosis factor-alpha; TUNEL: TdT-mediated dUTP Nick-End Labeling; UAlb/Cr: urine micro albumin/creatinine.

DNA studies of endangered or extinct species often rely on ancient or degraded remains. The majority of ancient DNA (aDNA) extraction protocols focus on skeletal elements, with skin and hair samples rarely explored. Similar to that found in bones and teeth, DNA extracted from historical or ancient skin and fur samples is also extremely fragmented with low endogenous content due to natural degradation processes. Thus, the development of effective DNA extraction methods is required for these materials. Here, we compared the performance of two DNA extraction protocols (commercial and custom laboratory aDNA methods) on hair and skin samples from decades-old museum specimens to Iron Age archaeological material. We found that apart from the impact sample-specific taphonomic and handling history has on the quantity and quality of DNA preservation, skin yielded more endogenous DNA than hair of the samples and protocols tested. While both methods recovered DNA from ancient soft tissue, the laboratory method performed better overall in terms of DNA yield and quality, which was primarily due to the poorer performance of the commercial binding buffer in recovering aDNA.

Significance Early maturing rice cultivars usually have a low yield, suffering from a shortening maturity duration in comparison with late-maturing ones. Here we demonstrate that shortening of maturity duration with no yield penalty of rice can be achieved by the quantitative trait locus Early flowering-completely dominant (Ef-cd). Ef-cd is a long noncoding RNA transcribed from the antisense strand of the flowering activator OsSOC1 locus, which can positively regulate the expression of OsSOC1. Physiological analysis demonstrated that Ef-cd could facilitate nitrogen utilization and also improve the photosynthesis rate. This study provides breeders with a valuable genetic resource that is useful for balancing grain yield with maturity duration for further elevating food production in the world., The contradiction between “high yielding” and “early maturing” hampers further improvement of annual rice yield. Here we report the positional cloning of a major maturity duration regulatory gene, Early flowering-completely dominant (Ef-cd), and demonstrate that natural variation in Ef-cd could be used to overcome the above contradictory. The Ef-cd locus gives rise to a long noncoding RNA (lncRNA) antisense transcript overlapping the OsSOC1 gene. Ef-cd lncRNA expression positively correlates with the expression of OsSOC1 and H3K36me3 deposition. Field test comparisons of early maturing Ef-cd near-isogenic lines with their wild types as well as of the derivative early maturing hybrids with their wild-type hybrids conducted under different latitudes determined that the early maturing Ef-cd allele shortens maturity duration (ranging from 7 to 20 d) without a concomitant yield penalty. Ef-cd facilitates nitrogen utilization and also improves the photosynthesis rate. Analysis of 1,439 elite hybrid rice varieties revealed that the 16 homozygotes and 299 heterozygotes possessing Ef-cd matured significantly earlier. Therefore, Ef-cd could be a vital contributor of elite early maturing hybrid varieties in balancing grain yield with maturity duration.

The maximum quantum yield of photosystem II (as reflected by variable to maximum chlorophyll a fluorescence, Fv /Fm ) is regarded as one of the most important photosynthetic parameters. The genetic basis underlying natural variation in Fv /Fm , which shows low level of variations in plants under non-stress conditions, is not easy to be exploited using the conventional gene cloning approaches. Thus, in order to answer this question, we have followed another strategy: we used genome-wide association study (GWAS) and transgenic analysis in a rice mini-core collection. We report here that four single-nucleotide polymorphisms, located in the promoter region of β-glucosidase 5 (BGlu-5), are associated with observed variation in Fv /Fm . Indeed, our transgenic analysis showed a good correlation between BGlu-5 and Fv /Fm . Thus, our work demonstrates the feasibility of using GWAS to study natural variation in Fv /Fm , suggesting that cis-element polymorphism, affecting the BGlu-5 expression level, may, indirectly, contribute to Fv /Fm variation in rice through the gibberellin signaling pathway. Further research is needed to understand the mechanism of our novel observation.

The rooting of the SARS-CoV-2 phylogeny is important for understanding the origin and early spread of the virus. Previously published phylogenies have used different rootings that do not always provide consistent results. We investigate several different strategies for rooting the SARS-CoV-2 tree and provide measures of statistical uncertainty for all methods. We show that methods based on the molecular clock tend to place the root in the B clade, whereas methods based on outgroup rooting tend to place the root in the A clade. The results from the two approaches are statistically incompatible, possibly as a consequence of deviations from a molecular clock or excess back-mutations. We also show that none of the methods provide strong statistical support for the placement of the root in any particular edge of the tree. These results suggest that phylogenetic evidence alone is unlikely to identify the origin of the SARS-CoV-2 virus and we caution against strong inferences regarding the early spread of the virus based solely on such evidence.

Human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is most closely related, by average genetic distance, to two coronaviruses isolated from bats, RaTG13 and RmYN02. However, there is a segment of high amino acid similarity between human SARS-CoV-2 and a pangolin isolated strain, GD410721, in the receptor binding domain (RBD) of the spike protein, a pattern that can be caused by either recombination or by convergent amino acid evolution driven by natural selection. We perform a detailed analysis of the synonymous divergence, which is less likely to be affected by selection than amino acid divergence, between human SARS-CoV-2 and related strains. We show that the synonymous divergence between the bat derived viruses and SARS-CoV-2 is larger than between GD410721 and SARS-CoV-2 in the RBD, providing strong additional support for the recombination hypothesis. However, the synonymous divergence between pangolin strain and SARS-CoV-2 is also relatively high, which is not consistent with a recent recombination between them, instead it suggests a recombination into RaTG13. We also find a 14-fold increase in thedN/dSratio from the lineage leading to SARS-CoV-2 to the strains of the current pandemic, suggesting that the vast majority of non-synonymous mutations currently segregating within the human strains have a negative impact on viral fitness. Finally, we estimate that the time to the most recent common ancestor of SARS-CoV-2 and RaTG13 or RmYN02 based on synonymous divergence, is 51.71 years (95% C.I., 28.11-75.31) and 37.02 years (95% C.I., 18.19-55.85), respectively.

The acceleration of stomatal closure upon high to low light transition could improve plant water use efficiency and drought tolerance. Herein, using genome-wide association study, we showed that the genetic variation in OsNHX1 was strongly associated with the changes in τcl , the time constant of stomatal closure, in 206 rice accessions. OsNHX1 overexpression in rice resulted in a decrease in τcl , and an increase in biomass, grain yield under drought. Conversely, OsNHX1 knockout by CRISPR/CAS9 shows opposite trends for these traits. We further found three haplotypes spanning the OsNHX1 promoter and CDS regions. Two among them, HapII and HapIII, were found to be associated with a high and low τcl , respectively. A near-isogenic line (NIL, S464) was developed through replacing the genomic region harboring HapII (~10 kb) from MH63 (recipient) rice cultivar by the same sized genomic region containing Hap III from 02428 (donor). Compared with MH63, S464 shows a reduction by 35% in τcl and an increase by 40% in the grain yield under drought. However, under normal conditions, S464 maintains closely similar grain yield as MH63. The global distribution of the two OsNHX1 haplotypes is associated with the local precipitation. Taken together, the natural variation in OsNHX1 could be utilized to manipulate the stomatal dynamics for an improved rice drought tolerance.

BACKGROUND Ischemia-reperfusion injury is caused by a blood reperfusion injury in ischemic brain tissue, and usually occurs in the treatment stage of ischemic disease, which can aggravate brain tissue injury. MiR-122 is closely related to ischemia-reperfusion injury in the myocardium, kidney, and liver; however, the role in cerebral ischemia-reperfusion injury has not been established. MATERIAL AND METHODS In this study, cerebral ischemia-reperfusion injury was established in a rat model, and the control group was a sham-operated group. After ischemia-reperfusion injury for 6, 12, and 24 hours, brain tissue specimens were collected and the expression of miR-122 and DJ-1 were determined using quantitative real-time polymerase chain reaction. Flow cytometry was used to determine the reactive oxygen species (ROS) content. The modified Neurological Severity Score (mNSS) scale was used to evaluate the sensory and motor function defects of the rats. The malondialdehyde (MDA), superoxide dismutase (SOD), and enzyme activity were determined. The rats in the cerebral ischemia-reperfusion injury model were divided into 2 groups (antagomir-NC group and antagomir miR-122 group). Brain neuron RN-c cells were divided into the following 4 groups: antagomir-NC, antagomir miR-122, pIRES2-blank, and pIRES2-DJ-1. Seventy-two hours after transfection, ischemia-reperfusion treatment was carried out and conventional cultured RN-c cells were used as the control group. Flow cytometry was used to detect apoptosis and western blot was used to detect the expression of DJ-1, PTEN, AKT, and p-AKT. RESULTS The expression of miR-122 increased significantly in the process of ischemia-reperfusion damage after cerebral infarction, while the expression of DJ-1 decreased significantly. Downregulation of miR-122 significantly increased the expression of DJ-1, enhanced the activity of the PTEN/PI3K/AKT pathway, reduced cell apoptosis, and alleviated cerebral ischemia-reperfusion injury. CONCLUSIONS Inhibition of miR-122 can decrease cerebral ischemia-reperfusion injury by upregulating DJ-1-PTEN/PI3K/AKT pathway.

Domesticated species often exhibit convergent phenotypic evolution, termed the domestication syndrome, of which loss of seed dormancy is a component. To date, dormancy genes that contribute to parallel domestication across different families have not been reported. Here, we cloned the classical stay-green G gene from soybean and found that it controls seed dormancy and showed evidence of selection during soybean domestication. Moreover, orthologs in rice and tomato also showed evidence of selection during domestication. Analysis of transgenic plants confirmed that orthologs of G had conserved functions in controlling seed dormancy in soybean, rice, and Arabidopsis. Functional investigation demonstrated that G affected seed dormancy through interactions with NCED3 and PSY and in turn modulated abscisic acid synthesis. Therefore, we identified a gene responsible for seed dormancy that has been subject to parallel selection in multiple crop families. This may help facilitate the domestication of new crops.

Many studies have reported the recovery ability of umbilical cord-derived mesenchymal stem cells (UC-MSCs) for neural diseases. In this study, the authors explored the roles of UC-MSCs to treat the traumatic brain injury. Umbilical cord-derived mesenchymal stem cells were isolated from healthy neonatal rat umbilical cord immediately after delivery. The traumatic brain injury (TBI) model was formed by the classical gravity method. The authors detected the behavior changes and measured the levels of inflammatory factors, such as interleukin-lβ and tumor necrosis factor-α by enzyme linked immunosorbent assay (ELISA) at 1, 2, 3, 4 weeks after transplantation between TBI treated and untreated with UC-MSCs. Simultaneously, the expression of glial cell line-derived neurotrophic factor (GDNF) and brain derived neurotrophic factor (BDNF) were measured by real-time–polymerase chain reaction and ELISA. The authors found that the group of transplantation UC-MSCs has a significant improvement than other group treated by phosphate buffered saline. In the behavioral test, the Neurological Severity Scores of UC-MSCs + TBI group were lower than TBI group (P

Cold stress is a major factor limiting production and geographic distribution of rice (Oryza sativa). Although the growth range of japonica subspecies has expanded northward compared to modern wild rice (O. rufipogon), the molecular basis of the adaptation remains unclear. Here we report bZIP73, a bZIP transcription factor-coding gene with only one functional polymorphism (+511 G>A) between the two subspecies japonica and indica, may have facilitated japonica adaptation to cold climates. We show the japonica version of bZIP73 (bZIP73Jap) interacts with bZIP71 and modulates ABA levels and ROS homeostasis. Evolutionary and population genetic analyses suggest bZIP73 has undergone balancing selection; the bZIP73Jap allele has firstly selected from standing variations in wild rice and likely facilitated cold climate adaptation during initial japonica domestication, while the indica allele bZIP73Ind was subsequently selected for reasons that remain unclear. Our findings reveal early selection of bZIP73Jap may have facilitated climate adaptation of primitive rice germplasms., Japonica rice can grow further north than wild or indica rice and is more tolerant of cold climates. Here, the authors show that bZIP73 likely underwent selection in the early phase of rice domestication to facilitate cold tolerance in japonica by modulating ABA and ROS homeostasis.

Pre-harvest sprouting (PHS) is an unfavorable trait in cereal crops that could seriously decrease grain yield and quality. Although some PHS-associated quantitative trait loci or genes in cereals have been reported, the molecular mechanism underlying PHS remains largely elusive. Here, we characterized a rice mutant, phs8, which exhibits PHS phenotype accompanied by sugary endosperm. Map-based cloning revealed that PHS8 encodes a starch debranching enzyme named isoamylase1. Mutation in PHS8 resulted in the phytoglycogen breakdown and sugar accumulation in the endosperm. Intriguingly, with increase of sugar contents, decreased expression of OsABI3 and OsABI5 as well as reduced sensitivity to abscisic acid (ABA) were found in the phs8 mutant. Using rice suspension cell system, we confirmed that exogenous sugar is sufficient to suppress the expression of both OsABI3 and OsABI5. Furthermore, overexpression of OsABI3 or OsABI5 could partially rescue the PHS phenotype of phs8. Therefore, our study presents important evidence supporting that endosperm sugar not only acts as an essential energy source for seed germination but also determines seed dormancy and germination by affecting ABA signaling.

Nitrogen (N) is a major driving force for crop yield improvement, but application of high levels of N delays flowering, prolonging maturation and thus increasing the risk of yield losses. Therefore, traits that enable utilization of high levels of N without delaying maturation will be highly desirable for crop breeding. Here, we show that OsNRT1.1A (OsNPF6.3), a member of the rice (Oryza sativa) nitrate transporter 1/peptide transporter family, is involved in regulating N utilization and flowering, providing a target to produce high yield and early maturation simultaneously. OsNRT.1A has functionally diverged from previously reported NRT1.1 genes in plants and functions in upregulating the expression of N utilization-related genes not only for nitrate but also for ammonium, as well as flowering-related genes. Relative to the wild type, osnrt1.1a mutants exhibited reduced N utilization and late flowering. By contrast, overexpression of OsNRT1.1A in rice greatly improved N utilization and grain yield, and maturation time was also significantly shortened. These effects were further confirmed in different rice backgrounds and also in Arabidopsis thaliana Our study paves a path for the use of a single gene to dramatically increase yield and shorten maturation time for crops, outcomes that promise to substantially increase world food security.

In the present study, the expression of Beclin 1 in osteoarthritis (OA) cartilage tissue was investigated, and also its role in proliferation, apoptosis and expression of matrix metalloproteinases (MMPs) in chondrocytes obtained from patients with OA. Beclin 1 expression in cartilage tissue from OA patients, and in the age- and sex-matched controls, was detected by immunohistochemistry, semi-quantitative polymerase chain reaction and western blotting. Chondrocytes were divided into control and Beclin 1-overexpressed groups. After transfection for 48, 72 and 96 h, cell viability, apoptosis, the phosphatidylinositol-3-kinase (PI3K)/protein kinase B (Akt)/mammalian target of rapamycin (mTOR) signaling pathway and MMPs were examined. The mRNA and protein expression levels of Beclin 1 were significantly decreased in cartilage tissue from OA patients compared with the sex- and age-matched controls (P

Integrative analysis of multi-omics layers at single cell level is critical for accurate dissection of cell-to-cell variation within certain cell populations. Here we report scCAT-seq, a technique for simultaneously assaying chromatin accessibility and the transcriptome within the same single cell. We show that the combined single cell signatures enable accurate construction of regulatory relationships between cis-regulatory elements and the target genes at single-cell resolution, providing a new dimension of features that helps direct discovery of regulatory patterns specific to distinct cell identities. Moreover, we generate the first single cell integrated map of chromatin accessibility and transcriptome in early embryos and demonstrate the robustness of scCAT-seq in the precise dissection of master transcription factors in cells of distinct states. The ability to obtain these two layers of omics data will help provide more accurate definitions of “single cell state” and enable the deconvolution of regulatory heterogeneity from complex cell populations., Heterogeneity in gene expression and epigenetic states exists across individual cells. Here, the authors develop scCAT-seq, a technique for simultaneously performing ATAC-seq and RNA-seq within the same single cell.

The clarification of the genetic origins of present-day Tibetans requires an understanding of their past relationships with the ancient populations of the Tibetan Plateau. Here we successfully sequenced 67 complete mitochondrial DNA genomes of 5200 to 300-year-old humans from the plateau. Apart from identifying two ancient plateau lineages (haplogroups D4j1b and M9a1a1c1b1a) that suggest some ancestors of Tibetans came from low-altitude areas 4750 to 2775 years ago and that some were involved in an expansion of people moving between high-altitude areas 2125 to 1100 years ago, we found limited evidence of recent matrilineal continuity on the plateau. Furthermore, deep learning of the ancient data incorporated into simulation models with an accuracy of 97% supports that present-day Tibetan matrilineal ancestry received partial contribution rather than complete continuity from the plateau populations of the last 5200 years.

Melatonin, an antioxidant in both animals and plants, has been reported to have beneficial effects on the aging process. It was also suggested to play a role in extending longevity and enhancing abiotic stress resistance in plant. In this study, we demonstrate that melatonin acts as a potent agent to delay leaf senescence and cell death in rice. Treatments with melatonin significantly reduced chlorophyll degradation, suppressed the transcripts of senescence-associated genes, delayed the leaf senescence, and enhanced salt stress tolerance. Genome-wide expression profiling by RNA sequencing reveals that melatonin is a potent free radical scavenger, and its exogenous application results in enhanced antioxidant protection. Leaf cell death in noe1, a mutant with over-produced H2O2, can be relieved by exogenous application of melatonin. These data demonstrate that melatonin delays the leaf senescence and cell death and also enhances abiotic stress tolerance via directly or indirectly counteracting the cellular accumulation of H2O2.

Grain size is one of the most important determinants of crop yield in cereals. Here, we identified a dominant mutant, big grain2 (bg2-D) from our enhancer-trapping population. Genetic analysis and SiteFinding PCR (polymerase chain reaction) revealed that BG2 encodes a cytochrome P450, OsCYP78A13. Sequence search revealed that CYP78A13 has a paralogue Grain Length 3.2 (GL3.2, LOC_Os03g30420) in rice with distinct expression patterns, analysis of transgenic plants harbouring either CYP78A13 or GL3.2 showed that both can promote grain growth. Sequence polymorphism analysis with 1529 rice varieties showed that the nucleotide diversity at CYP78A13 gene body and the 20 kb flanking region in the indica varieties were markedly higher than those in japonica varieties. Further, comparison of the genomic sequence of CYP78A13 in the japonica cultivar Nipponbare and the indica cultivar 9311 showed that there were three InDels in the promoter region and eight SNPs (single nucleotide polymorphism) in its coding sequence. Detailed examination of the transgenic plants with chimaeric constructs suggested that variation in CYP78A13 coding region is responsible for the variation of grain yield. Taken together, our results suggest that the variations in CYP78A13 in the indica varieties hold potential in rice breeding for application of grain yield improvement.

Mining natural variations is a major approach to identify new options to improve crop light use efficiency. So far, successes in identifying photosynthetic parameters positively related to crop biomass accumulation through this approach are scarce, possibly due to the earlier emphasis on properties related to leaf instead of canopy photosynthetic efficiency. This study aims to uncover rice (Oryza sativa) natural variations to identify leaf physiological parameters that are highly correlated with biomass accumulation, a surrogate of canopy photosynthesis. To do this, we systematically investigated 14 photosynthetic parameters and four morphological traits in a rice population, which consists of 204 U.S. Department of Agriculture-curated minicore accessions collected globally and 11 elite Chinese rice cultivars in both Beijing and Shanghai. To identify key components responsible for the variance of biomass accumulation, we applied a stepwise feature-selection approach based on linear regression models. Although there are large variations in photosynthetic parameters measured in different environments, we observed that photosynthetic rate under low light (Alow) was highly related to biomass accumulation and also exhibited high genomic inheritability in both environments, suggesting its great potential to be used as a target for future rice breeding programs. Large variations in Alow among modern rice cultivars further suggest the great potential of using this parameter in contemporary rice breeding for the improvement of biomass and, hence, yield potential.

The domestication history of rice remains controversial, with multiple studies reaching different conclusions regarding its origin(s). These studies have generally assumed that populations of living wild rice, O. rufipogon, are descendants of the ancestral population that gave rise to domesticated rice, but relatively little attention has been paid to the origins and history of wild rice itself. Here, we investigate the genetic ancestry of wild rice by analyzing a diverse panel of rice genomes consisting of 203 domesticated and 435 wild rice accessions. We show that most modern wild rice is heavily admixed with domesticated rice through both pollen- and seed-mediated gene flow. In fact, much presumed wild rice may simply represent different stages of feralized domesticated rice. In line with this hypothesis, many presumed wild rice varieties show remnants of the effects of selective sweeps in previously identified domestication genes, as well as evidence of recent selection in flowering genes possibly associated with the feralization process. Furthermore, there is a distinct geographical pattern of gene flow from aus, indica, and japonica varieties into colocated wild rice. We also show that admixture from aus and indica is more recent than gene flow from japonica, possibly consistent with an earlier spread of japonica varieties. We argue that wild rice populations should be considered a hybrid swarm, connected to domesticated rice by continuous and extensive gene flow.

Grain size and shape are important determinants of grain weight and yield in rice. Here, we report a new major quantitative trait locus (QTL), qTGW3, that controls grain size and weight in rice. This locus, qTGW3, encodes OsSK41 (also known as OsGSK5), a member of the GLYCOGEN SYNTHASE KINASE 3/SHAGGY-like family. Rice near-isogenic lines carrying the loss-of-function allele of OsSK41 have increased grain length and weight. We demonstrate that OsSK41 interacts with and phosphorylates AUXIN RESPONSE FACTOR 4 (OsARF4). Co-expression of OsSK41 with OsARF4 increases the accumulation of OsARF4 in rice protoplasts. Loss of function of OsARF4 results in larger rice grains. RNA-sequencing analysis suggests that OsARF4 and OsSK41 repress the expression of a common set of downstream genes, including some auxin-responsive genes, during rice grain development. The loss-of-function form of OsSK41 at qTGW3 represents a rare allele that has not been extensively utilized in rice breeding. Suppression of OsSK41 function by either targeted gene editing or QTL pyramiding enhances rice grain size and weight. Thus, our study reveals the important role of OsSK41 in rice grain development and provides new candidate genes for genetic improvement of grain yield in rice and perhaps in other cereal crops.

Rice (Oryza sativa L.) is an excellent model plant in elucidation of cereal domestication. Loss of seed shattering, weakened dormancy, and changes in plant architecture were thought to be three key events in the rice domestication and creating the high-yield, uniform-germinating, and densely-planting modern rice. Loss of shattering is considered to be the direct morphological evidence for identifying domesticated rice. Two major shattering QTLs, Sh4 and qSH1, have displayed different domestication histories. Weakened seed dormancy is essential for synchronous germination in agricultural production. Genes Sdr4, qSD7-1, and qSD12 impose a global and complementary adaptation strategies in controlling seed dormancy. The prostate growth habit of wild rice is an adaptation to disturbed habitats, while the erect growth habit of rice cultivars meet the needs of compact planting, and such a plant architecture is mainly controlled by PROG1. The outcrossing habit of wild rice promotes propagation of domestication genes among different populations, while the self-pollinating habit of cultivated rice facilitates fixation of domestication genes. Currently, the researches on rice domestication mainly focus on individual genes or multiple neutral markers, and much less attention has been paid to the evolution of network controlling domestication traits. With the progress in functional genomics research, the molecular mechanism of domestication traits is emerging. Rice domestication researches based on network will be more comprehensive and better reflect rice domestica-tion process. Here, we reviewed most progresses in molecular mechanisms of rice domestication traits, in order to provide the new insights for rice domestication and molecular breeding.

Low-coverage whole-genome sequencing is an effective strategy for genome-wide association studies in humans, due to the availability of large reference panels for genotype imputation. However, it is unclear whether this strategy can be utilized in other species without reference panels. Using simulations, we show that this approach is even more relevant in inbred species such as rice (Oryza sativa L.), which are effectively haploid, allowing easy haplotype construction and imputation-based genotype calling, even without the availability of large reference panels. We sequenced 203 rice varieties with well-characterized phenotypes from the United States Department of Agriculture Rice Mini-Core Collection at an average depth of 1.5× and used the data for mapping three traits. For the first two traits, amylose content and seed length, our approach leads to direct identification of the previously identified causal SNPs in the major-effect loci. For the third trait, pericarp color, an important trait underwent selection during domestication, we identified a new major-effect locus. Although known loci can explain color variation in the varieties of two main subspecies of Asian domesticated rice, japonica and indica, the new locus identified is unique to another domesticated rice subgroup, aus, and together with existing loci, can fully explain the major variation in pericarp color in aus. Our discovery of a unique genetic basis of white pericarp in aus provides an example of convergent evolution during rice domestication and suggests that aus may have a domestication history independent of japonica and indica.

Given the continuously growing population and decreasing arable land, food shortage is becoming one of the most serious global problems in this century1. Grain size is one of the determining factors for grain yield and thus is a prime target for genetic breeding2,3. Although a number of quantitative trait loci (QTLs) associated with rice grain size have been identified in the past decade, mechanisms underlying their functions remain largely unknown4,5. Here we show that a grain-length-associated QTL, GL2, has the potential to improve grain weight and grain yield up to 27.1% and 16.6%, respectively. We also show that GL2 is allelic to OsGRF4 and that it contains mutations in the miR396 targeting sequence. Because of the mutation, GL2 has a moderately increased expression level, which consequently activates brassinosteroid responses by upregulating a large number of brassinosteroid-induced genes to promote grain development. Furthermore, we found that GSK2, the central negative regulator of rice brassinosteroid signalling, directly interacts with OsGRF4 and inhibits its transcription activation activity to mediate the specific regulation of grain length by the hormone. Thus, this work demonstrates the feasibility of modulating specific brassinosteroid responses to improve plant productivity. To understand the mechanisms of grain size control, researchers experimentally reveal that the molecular module miR396/GRF4 regulates rice grain size by activating brassinosteroid signalling. Modulating miR396/GRF4 or brassinosteroid responses can thus be used to improve crop yield.

Asian cultivated rice (Oryza sativa L.) consists of two main subspecies, indica and japonica. Indica has higher nitrate-absorption activity than japonica, but the molecular mechanisms underlying that activity remain elusive. Here we show that variation in a nitrate-transporter gene, NRT1.1B (OsNPF6.5), may contribute to this divergence in nitrate use. Phylogenetic analysis revealed that NRT1.1B diverges between indica and japonica. NRT1.1B-indica variation was associated with enhanced nitrate uptake and root-to-shoot transport and upregulated expression of nitrate-responsive genes. The selection signature of NRT1.1B-indica suggests that nitrate-use divergence occurred during rice domestication. Notably, field tests with near-isogenic and transgenic lines confirmed that the japonica variety carrying the NRT1.1B-indica allele had significantly improved grain yield and nitrogen-use efficiency (NUE) compared to the variety without that allele. Our results show that variation in NRT1.1B largely explains nitrate-use divergence between indica and japonica and that NRT1.1B-indica can potentially improve the NUE of japonica.

SiteFinding-PCR is a method for isolating flanking sequence tags (FSTs) of T-DNA insertion lines, but the efficiency needs to be improved. Here we report a computation-assisted design for the random primers used in SiteFinding- PCR. A short sequence, GCATG, was screened from the rice genome and used as the 3′ end of the random primer. When applying the optimized primer for isolating FSTs from 168 transgenic rice lines, we obtained 107 specific products, including 64 FSTs. The efficiency of obtaining FSTs using the modified version of SiteFinding-PCR increased by 73.0% compared with the method previously reported (P < 0.01, µ test). We also provide computational results for several other plant species such as maize, sorghum, Arabidopsis, foxtail millet, and Brachypodium based on the available genome data, so that the modified method could be easily adapted to other species.