A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase

Application of a Novel Quantitative Trait Locus Combination to Improve Grain Shape without Yield Loss in Rice (Oryza sativa L. spp. japonica)

Grain shape is one of the key factors deciding the yield product and the market value as appearance quality in rice (Oryza sativa L.). The grain shape of japonica cultivars in Korea is quite monotonous because the selection pressure of rice breeding programs works in consideration of consumer preference. In this study, we identified QTLs associated with grain shape to improve the variety of grain shapes in Korean cultivars. QTL analysis revealed that eight QTLs related to five tested traits were detected on chromosomes 2, 5, and 10. Among them, three QTLs-qGL2 (33.9% of PEV for grain length), qGW5 (64.42% for grain width), and qGT10 (49.2% for grain thickness)-were regarded as the main effect QTLs. Using the three QTLs, an ideal QTL combination (qGL2^P + qGW5^P + qGT10^B) could be constructed on the basis of the accumulated QTL effect without yield loss caused by the change in grain shape in the population. In addition, three promising lines with a slender grain type were selected as a breeding resource with a japonica genetic background based on the QTL combination. The application of QTLs detected in this study could improve the grain shape of japonica cultivars without any linkage drag or yield loss.

Can the Wild Perennial, Rhizomatous Rice Species Oryza longistaminata be a Candidate for De Novo Domestication?

As climate change intensifies, the development of resilient rice that can tolerate abiotic stresses is urgently needed. In nature, many wild plants have evolved a variety of mechanisms to protect themselves from environmental stresses. Wild relatives of rice may have abundant and virtually untapped genetic diversity and are an essential source of germplasm for the improvement of abiotic stress tolerance in cultivated rice. Unfortunately, the barriers of traditional breeding approaches, such as backcrossing and transgenesis, make it challenging and complex to transfer the underlying resilience traits between plants. However, de novo domestication via genome editing is a quick approach to produce rice with high yields from orphans or wild relatives. African wild rice, Oryza longistaminata, which is part of the AA-genome Oryza species has two types of propagation strategies viz. vegetative propagation via rhizome and seed propagation. It also shows tolerance to multiple types of abiotic stress, and therefore O. longistaminata is considered a key candidate of wild rice for heat, drought, and salinity tolerance, and it is also resistant to lodging. Importantly, O. longistaminata is perennial and propagates also via rhizomes both of which are traits that are highly valuable for the sustainable production of rice. Therefore, O. longistaminata may be a good candidate for de novo domestication through genome editing to obtain rice that is more climate resilient than modern elite cultivars of O. sativa.

An endoplasmic reticulum-associated degradation-related E2-E3 enzyme pair controls grain size and weight through the brassinosteroid signaling pathway in rice

Grain size is an important agronomic trait, but our knowledge about grain size determination in crops is still limited. Endoplasmic reticulum (ER)-associated degradation (ERAD) is a special ubiquitin proteasome system that is involved in degrading misfolded or incompletely folded proteins in the ER. Here, we report that SMALL GRAIN 3 (SMG3) and DECREASED GRAIN SIZE 1 (DGS1), an ERAD-related E2-E3 enzyme pair, regulate grain size and weight through the brassinosteroid (BR) signaling pathway in rice (Oryza sativa). SMG3 encodes a homolog of Arabidopsis (Arabidopsis thaliana) UBIQUITIN CONJUGATING ENZYME 32, which is a conserved ERAD-associated E2 ubiquitin conjugating enzyme. SMG3 interacts with another grain size regulator, DGS1. Loss of function of SMG3 or DGS1 results in small grains, while overexpression of SMG3 or DGS1 leads to long grains. Further analyses showed that DGS1 is an active E3 ubiquitin ligase and colocates with SMG3 in the ER. SMG3 and DGS1 are involved in BR signaling. DGS1 ubiquitinates the BR receptor BRASSINOSTEROID INSENSITIVE 1 (BRI1) and affects its accumulation. Genetic analysis suggests that SMG3, DGS1, and BRI1 act together to regulate grain size and weight. In summary, our findings identify an ERAD-related E2-E3 pair that regulates grain size and weight, which gives insight into the function of ERAD in grain size control and BR signaling.

Genetic Localization and Homologous Genes Mining for Barley Grain Size

Grain size is an important agronomic trait determining barley yield and quality. An increasing number of QTLs (quantitative trait loci) for grain size have been reported due to the improvement in genome sequencing and mapping. Elucidating the molecular mechanisms underpinning barley grain size is vital for producing elite cultivars and accelerating breeding processes. In this review, we summarize the achievements in the molecular mapping of barley grain size over the past two decades, highlighting the results of QTL linkage analysis and genome-wide association studies. We discuss the QTL hotspots and predict candidate genes in detail. Moreover, reported homologs that determine the seed size clustered into several signaling pathways in model plants are also listed, providing the theoretical basis for mining genetic resources and regulatory networks of barley grain size.

A chromosome-level genome assembly of an early matured aromatic Japonica rice variety Qigeng10 to accelerate rice breeding for high grain quality in Northeast China

Early-matured aromatic japonica rice from the Northeast is the most popular rice commodity in the Chinese market. The Qigeng10 (QG10) was one of the varieties with the largest planting area in this region in recent years. It was an early-matured japonica rice variety with a lot of superior traits such as semi-dwarf, lodging resistance, long grain, aromatic and good quality. Therefore, a high-quality assembly of Qigeng10 genome is critical and useful for japonica research and breeding. In this study, we produced a high-precision QG10 chromosome-level genome by using a combination of Nanopore and Hi-C platforms. Finally, we assembled the QG10 genome into 77 contigs with an N50 length of 11.80 Mb in 27 scaffolds with an N50 length of 30.55 Mb. The assembled genome size was 378.31Mb with 65 contigs and constituted approximately 99.59% of the 12 chromosomes. We identified a total of 1,080,819 SNPs and 682,392 InDels between QG10 and Nipponbare. We also annotated 57,599 genes by the Ab initio method, homology-based technique, and RNA-seq. Based on the assembled genome sequence, we detected the sequence variation in a total of 63 cloned genes involved in grain yield, grain size, disease tolerance, lodging resistance, fragrance, and many other important traits. Finally, we identified five elite alleles (qTGW2^Nipponbare , qTGW3^Nanyangzhan , GW5^IR24 , GW6^Suyunuo , and qGW8^Basmati385 ) controlling long grain size, four elite alleles (COLD1^Nipponbare , bZIP73^Nipponbare , CTB4a^{Kunmingxiaobaigu} , and CTB2^{Kunmingxiaobaigu} ) controlling cold tolerance, three non-functional alleles (DTH7^Kitaake , Ghd7^Hejiang19 , and Hd1^Longgeng31 ) for early heading, two resistant alleles (Pia^Akihikari and Pid4^Digu ) for rice blast, a resistant allele STV11^Kasalath for rice stripe virus, an NRT1.1B^IR24 allele for higher nitrate absorption activity, an elite allele SCM3^Chugoku117 for stronger culms, and the typical aromatic gene badh2-E2 for fragrance in QG10. These results not only help us to better elucidate the genetic mechanisms underlying excellent agronomic traits in QG10 but also have wide-ranging implications for genomics-assisted breeding in early-matured fragrant japonica rice.

Identification, Fine Mapping and Application of Quantitative Trait Loci for Grain Shape Using Single-Segment Substitution Lines in Rice (Oryza sativa L.)

Quantitative trait loci (QTLs) and HQTL (heterosis QTLs) for grain shape are two major genetic factors of grain yield and quality in rice (Oryza sativa L.). Although many QTLs for grain shape have been reported, only a few are applied in production. In this study, 54 QTLs for grain shape were detected on 10 chromosomes using 33 SSSLs (single-segment substitution lines) and methods of statistical genetics. Among these, 23 exhibited significant positive additive genetic effects, including some novel QTLs, among which qTGW4-(1,2), qTGW10-2, and qTGW10-3 were three QTLs newly found in this study and should be paid more attention. Moreover, 26 HQTLs for grain shape were probed. Eighteen of these exhibited significant positive dominant genetic effects. Thirty-three QTLs for grain shape were further mapped using linkage analysis. Most of the QTLs for grain shape produced pleiotropic effects, which simultaneously controlled multiple appearance traits of grain shape. Linkage mapping of the F₂ population derived from sub-single-segment substitution lines further narrowed the interval harbouring qTGW10-3 to 75.124 kb between PSM169 and RM25753. The candidate gene was identified and could be applied to breeding applications by molecular marker-assisted selection. These identified QTLs for grain shape will offer additional insights for improving grain yield and quality in rice breeding.

Molecular bases of rice grain size and quality for optimized productivity

The accomplishment of further optimization of crop productivity in grain yield and quality is a great challenge. Grain size is one of the crucial determinants of rice yield and quality; all of these traits are typical quantitative traits controlled by multiple genes. Research advances have revealed several molecular and developmental pathways that govern these traits of agronomical importance. This review provides a comprehensive summary of these pathways, including those mediated by G-protein, the ubiquitin-proteasome system, mitogen-activated protein kinase, phytohormone, transcriptional regulators, and storage product biosynthesis and accumulation. We also generalize the excellent precedents for rice variety improvement of grain size and quality, which utilize newly developed gene editing and conventional gene pyramiding capabilities. In addition, we discuss the rational and accurate breeding strategies, with the aim of better applying molecular design to breed high-yield and superior-quality varieties.

Targeted manipulation of grain shape genes effectively improves outcrossing rate and hybrid seed production in rice

Stigma exsertion rate (SER) of the male sterile line is a key limiting factor for hybrid seed production in rice. Although a large number of quantitative trait loci associated with SER have been reported, few genes have been molecularly cloned and functionally characterized, severely hindering the genetic improvement of SER of the male sterile line and the breeding efficiency of hybrid rice. In this study, we identified three grain shape regulatory genes, GS3, GW8 and GS9, as potential candidate genes for targeted manipulation of grain shape and SER. We show that simultaneously knocking out these three genes could effectively increase SER by increasing the ratio of spikelet length/spikelet width and length of stigma and style, without negative impacts on other agronomic traits. Cellular examination and transcriptomic analyses revealed a role of these genes in coordinated regulation of transverse and longitudinal cell division in the pistils. Moreover, we demonstrate that targeted manipulation of these grain shape genes could significantly improve the outcrossing rate in both the ZH11 (a japonica variety) and Zhu6S (an indica male sterile line) backgrounds. Our results provide new insights into the mechanisms of rice SER regulation and develop an effective strategy to improve SER and out-crossing rate in rice, thus facilitating hybrid rice production.

Understanding the regulation of cereal grain filling: The way forward

During grain filling, starch and other nutrients accumulate in the endosperm; this directly determines grain yield and grain quality in crops such as rice (Oryza sativa), maize (Zea mays), and wheat (Triticum aestivum). Grain filling is a complex trait affected by both intrinsic and environmental factors, making it difficult to explore the underlying genetics, molecular regulation, and the application of these genes for breeding. With the development of powerful genetic and molecular techniques, much has been learned about the genes and molecular networks related to grain filling over the past decades. In this review, we highlight the key factors affecting grain filling, including both biological and abiotic factors. We then summarize the key genes controlling grain filling and their roles in this event, including regulators of sugar translocation and starch biosynthesis, phytohormone-related regulators, and other factors. Finally, we discuss how the current knowledge of valuable grain filling genes could be integrated with strategies for breeding cereal varieties with improved grain yield and quality.

Analysis of Domestication Loci in Wild Rice Populations

The domestication syndrome is defined as a collection of domestication-related traits that have undergone permanent genetic changes during the domestication of cereals. Australian wild rice populations have not been exposed to gene flow from domesticated rice populations. A high level of natural variation of the sequences at domestication loci (e.g., seed shattering, awn development, and grain size) was found in Australian AA genome wild rice from the primary gene pool of rice. This natural variation is much higher than that found in Asian cultivated rice and wild Asian rice populations. The Australian Oryza meridionalis populations exhibit a high level of homozygous polymorphisms relative to domesticated rice, inferring the fixation of distinct wild and domesticated alleles. Alleles of the seed shattering genes (SH4/SHA1 and OsSh1/SH1) present in the shattering-prone O. meridionalis populations are likely to be functional, while the dysfunctional alleles of these seed shattering genes are found in domesticated rice. This confirms that unlike Asian wild rice populations, Australian wild rice populations have remained genetically isolated from domesticated rice, retaining pre-domestication alleles in their wild populations that uniquely allow the impact of domestication on the rice genome to be characterized. This study also provides key information about the domestication loci in Australian wild rice populations that will be valuable in the utilization of these genetic resources in crop improvement and de novo domestication.

Genome-Wide Association Study of Rice Grain Shape and Chalkiness in a Worldwide Collection of Xian Accessions

Rice (Oryza sativa L.) appearance quality, which is mainly defined by grain shape and chalkiness, is an important target in rice breeding. In this study, we first re-sequenced 137 indica accessions and then conducted a genome-wide association study (GWAS) for six agronomic traits with the 2,998,034 derived single nucleotide polymorphisms (SNPs) by using the best linear unbiased prediction (BLUP) values for each trait. The results revealed that 195 SNPs had significant associations with the six agronomic traits. Based on the genome-wide linkage disequilibrium (LD) blocks, candidate genes for the target traits were detected within 100 kb upstream and downstream of the relevant SNP loci. Results indicate that six quantitative trait loci (QTLs) significantly associated with six traits (qTGW4.1, qTGW4.2, qGL4.1, qGL12.1, qGL12.2, qGW2.1, qGW4.1, qGW6.1, qGW8.1, qGW8.2, qGW9.1, qGW11.1, qGLWR2.1, qGLWR2.2, qGLWR4.2, qPGWC5.1 and qDEC6.1) were identified for haplotype analysis. Among these QTLs, two (qTGW4.2 and qGW6.1), were overlapped with FLO19 and OsbZIP47, respectively, and the remaining four were novel QTLs. These candidate genes were further validated by haplotype block construction.

Artificially Selected Grain Shape Gene Combinations in Guangdong Simiao Varieties of Rice (Oryza sativa L.)

BACKGROUND

Grain shape is a key trait in rice breeding. Although many QTLs and genes of grain shape have been identified, how different combinations of alleles of these genes affect grain shape is largely unknown. It is important to understand the effects of grain shape gene combinations for breeding by design. In the present study, we performed genetic dissection of the grain shapes in Guangdong Simiao varieties, a popular kind of rice in South China, to identify the effective alleles and their combination for breeding.

RESULTS

We selected two hundred nineteen indica accessions with diverse grain shapes and fifty-two Guangdong Simiao varieties with long and slender grain shapes for genome-wide selection analysis. The results showed that four (GS3, GS5, GW5 and GL7) of the twenty grain shape genes fall into the regions selected for in Guangdong Simiao varieties. Allele analysis and frequency distribution of these four genes showed that GS3allele3 and GW5allele2 accounted for 96.2%, and GL7allele2 and GS5allele2 accounted for 76.9% and 74.5% of the Simiao varieties, respectively. Further analysis of the allelic combinations showed that 30 allelic combinations were identified in the whole panel, with 28 allelic combinations found in the international indica accessions and 6 allelic combinations found in Guangdong Simiao varieties. There were mainly three combinations (combinations 17, 18 and 19) in the Guangdong Simiao varieties, with combination 19 (GS3allele3 + GW5allele2 + GL7allele2 + GS5allele2) having the highest percentage (51.9%). All three combinations carried GS3allele3 + GW5allele2, while combinations 17 (GL7allele1) and 19 (GL7allele2) showed significant differences in both grain length and length/width ratio due to differences in GL7 alleles. Pedigree analysis of Guang8B, the maintainer of the first released Simiao male sterile line Guang8A, showed that the parent lines and Guang8B carried GS3allele3 + GW5allele2 + GS5allele2, while the GL7 allele differed, resulting in significant differences in grain size.

CONCLUSION

The results suggest that specific alleles of GS3, GS5, GW5 and GL7 are the key grain shape genes used in the Guangdong Simiao varieties and selected for grain shape improvement. Combination 19 is the predominant allelic combination in the Guangdong Simiao varieties. Our current study is the first to dissect the genetics of grain shape in Guangdong Simiao varieties, and the results will facilitate molecular breeding of Guangdong Simiao varieties.

Identification of major QTLs for soybean seed size and seed weight traits using a RIL population in different environments

INTRODUCTION

The seed weight of soybean [Glycine max (L.) Merr.] is one of the major traits that determine soybean yield and is closely related to seed size. However, the genetic basis of the synergistic regulation of traits related to soybean yield is unclear.

METHODS

To understand the molecular genetic basis for the formation of soybean yield traits, the present study focused on QTLs mapping for seed size and weight traits in different environments and target genes mining.

RESULTS

A total of 85 QTLs associated with seed size and weight traits were identified using a recombinant inbred line (RIL) population developed from Guizao1×B13 (GB13). We also detected 18 environmentally stable QTLs. Of these, qSL-3-1 was a novel QTL with a stable main effect associated with seed length. It was detected in all environments, three of which explained more than 10% of phenotypic variance (PV), with a maximum of 15.91%. In addition, qSW-20-3 was a novel QTL with a stable main effect associated with seed width, which was identified in four environments. And the amount of phenotypic variance explained (PVE) varied from 9.22 to 21.93%. Five QTL clusters associated with both seed size and seed weight were summarized by QTL cluster identification. Fifteen candidate genes that may be involved in regulating soybean seed size and weight were also screened based on gene function annotation and GO enrichment analysis.

DISCUSSION

The results provide a biologically basic reference for understanding the formation of soybean seed size and weight traits.

Co-Overexpression of Two Key Source Genes, OsBMY4 and OsISA3, Improves Multiple Key Traits of Rice Seeds

Optimized source-sink interactions are determinants of both rice yield and quality. However, most source genes have not been well studied in rice, a major grain crop. In this study, OsBMY4 and OsISA3, the key β-amylase and debranching enzymes that control transient starch degradation in rice leaves, were co-overexpressed in rice in order to accelerate starch degradation efficiency and increase the sugar supply for sink organs. Systematic analyses of the transgenic rice indicated that co-overexpression of OsBMY4 and OsISA3 not only promoted rice yield and quality, but also improved seed germination and stress tolerance. Moreover, since the OsBMY4 gene has not been characterized, we generated osbmy4 mutants using CRIPSR/Cas9 gene editing, which helped to reveal the roles of β-amylase in rice yield and quality. This study demonstrated that specific modulation of the expression of some key source genes improves the source-sink balance and leads to improvements in multiple key traits of rice seeds.

Fine Mapping and Cloning of a qRA2 Affect the Ratooning Ability in Rice (Oryza sativa L.)

Ratooning ability is a key factor that influences the ratoon rice yield in areas where light and temperature are not sufficient for second-season rice. Near-isogenic lines (NILs) are the most powerful tools for the detection and precise mapping of quantitative trait loci (QTLs). In this study, using 176 NILs, we identified a novel QTL for ratooning ability in NIL128. First, we mapped the QTL between the markers Indel12-29 and Indel12-31, which encompass a region of 233 kb. The rice genome annotation indicated the existence of three candidate genes in this region that may be related to ratooning ability. Through gene prediction and cDNA sequencing, we speculated that the target gene of ratooning ability is LOC_Os02g51930 which encodes cytokinin glucosyl transferases (CGTs), hereafter named qRA2. Further analysis showed that qra2 was a 1-bp substitution in the first exon in NIL128, which resulted in the premature termination of qRA2. The results of the knockdown experiment showed that the Jiafuzhan knockdown mutants exhibited the ratooning ability phenotype of NIL128. Interestingly, the qRA2 gene was found to improve ratooning ability without affecting major agronomic traits. These results will help us better understand the genetic basis of rice ratooning ability and provide a valuable gene resource for breeding strong ratoon rice varieties.

Identification of QTNs, QTN-by-environment interactions, and their candidate genes for grain size traits in main crop and ratoon rice

Although grain size is an important quantitative trait affecting rice yield and quality, there are few studies on gene-by-environment interactions (GEIs) in genome-wide association studies, especially, in main crop (MC) and ratoon rice (RR). To address these issues, the phenotypes for grain width (GW), grain length (GL), and thousand grain weight (TGW) of 159 accessions of MC and RR in two environments were used to associate with 2,017,495 SNPs for detecting quantitative trait nucleotides (QTNs) and QTN-by-environment interactions (QEIs) using 3VmrMLM. As a result, 64, 71, 67, 72, 63, and 56 QTNs, and 0, 1, 2, 2, 2, and 1 QEIs were found to be significantly associated with GW in MC (GW-MC), GL-MC, TGW-MC, GW-RR, GL-RR, and TGW-RR, respectively. 3, 4, 7, 2, 2, and 4 genes were found to be truly associated with the above traits, respectively, while 2 genes around the above QEIs were found to be truly associated with GL-RR, and one of the two known genes was differentially expressed under two soil moisture conditions. 10, 7, 1, 8, 4, and 3 candidate genes were found by differential expression and GO annotation analysis to be around the QTNs for the above traits, respectively, in which 6, 3, 1, 2, 0, and 2 candidate genes were found to be significant in haplotype analysis. The gene Os03g0737000 around one QEI for GL-MC was annotated as salt stress related gene and found to be differentially expressed in two cultivars with different grain sizes. Among all the candidate genes around the QTNs in this study, four were key, in which two were reported to be truly associated with seed development, and two (Os02g0626100 for GL-MC and Os02g0538000 for GW-MC) were new. Moreover, 1, 2, and 1 known genes, along with 8 additional candidate genes and 2 candidate GEIs, were found to be around QTNs and QEIs for GW, GL, and TGW, respectively in MC and RR joint analysis, in which 3 additional candidate genes were key and new. Our results provided a solid foundation for genetic improvement and molecular breeding in MC and RR.

Variations in Grain Traits among Local Rice Varieties Collected More Than Half-Century Ago in Indochinese Countries

More than half-century ago, local rice varieties were collected from Indochinese countries (Cambodia, Thailand, Laos, and Vietnam). Of these, 162 local varieties were examined for 7 grain-size traits: seed length/width/thickness, brown rice length/width/thickness, and 100-seed weight. Since these traits varied considerably, a survey of functional mutations was performed in the genes related to these traits. In total, 19 markers (12 InDel and 7 dCAPS markers) were used to investigate the mutations at 14 grain-size loci of GW2, GS2, qLGY3, GS3, GL3.1, TGW3, GS5, GW5, GS6, TGW6, GW6a, GLW7, GL7, and GW8. Significant allele effects were observed with six markers detecting base substitution mutations at GW2 and GS3 and insertion/deletion mutations at GS5, GW5, and GW6a, suggesting that these mutations might have affected the grain trait and caused variation among local varieties in the Indochinese countries. In addition to grain size, the hull color, grain color, and glutinosity were also examined using a survey of loss-of-function mutations at major responsible loci. Most phenotypes were reflected based on functional mutations at these loci. Since the local varieties have wide genetic variation, they are important genetic resources for future rice breeding.

Detecting and pyramiding target QTL for plant- and grain-related traits via chromosomal segment substitution line of rice

INTRODUCTION

Plant height and grain length are important agronomic traits in rice, exhibiting a strong effect on plant architecture and grain quality of rice varieties.

METHODS

Methods: A novel rice chromosomal segment substitution line (CSSL), i.e., CSSL-Z1357, with significantly increased plant height (PH) and grain length (GL) was identified from CSSLs constructed by using Nipponbare as a receptor and a restorer line Xihui 18 as a donor. Seven agronomic traits of PH, PL, GL, GW, GPP, SPP, and TGW were phenotyped, and REML implemented in HPMIXED of SAS were used to detect the QTL for these traits. Secondary CSSLs were screened out via marker-assisted selection (MAS) to estimate the additive and epistatic effects of detected QTLs, evaluating the potential utilization of pyramiding the target QTLs for yield and quality improvement of rice varieties.

RESULTS AND DISCUSSION

Results and Discussion: CSSL-Z1357 carried nine segments from Xihui 18 with an average segment length of 4.13 Mb. The results show that the long grain of CSSL-Z1357 was caused by the increased number of surface cells and the length of the inner glume. Thirteen quantitative trait loci were identified via the F2 population of Nipponbare/CSSL-Z1357, including three each for GL (qGL-3, qGL-6, and qGL-7) and PH (qPH-1, qPH-7, and qPH-12I), among which qGL-3 increased GL by 0.23 mm with synergistic allele from CSSL-Z1357. Additionally, three single (S1 to S3), two double (D1, D2), and one triple segment (T1) substitution lines were developed in F3 via MAS. Results show that pyramiding the segments from Chr.3 (qGL-3 and qPH-3), Chr.6 (qGL-6 and qPH-6), and Chr.7 (Null and qPH-7) tended to result in better phenotype of increased GL and PH and decreased grain width, providing a potential basis for enhancing grain yield and quality in rice breeding.

The double round-robin population unravels the genetic architecture of grain size in barley

Grain number, size and weight primarily determine the yield of barley. Although the genes regulating grain number are well studied in barley, the genetic loci and the causal gene for sink capacity are poorly understood. Therefore, the primary objective of our work was to dissect the genetic architecture of grain size and weight in barley. We used a multi-parent population developed from a genetic cross between 23 diverse barley inbreds in a double round-robin design. Seed size-related parameters such as grain length, grain width, grain area and thousand-grain weight were evaluated in the HvDRR population comprising 45 recombinant inbred line sub-populations. We found significant genotypic variation for all seed size characteristics, and observed 84% or higher heritability across four environments. The quantitative trait locus (QTL) detection results indicate that the genetic architecture of grain size is more complex than previously reported. In addition, both cultivars and landraces contributed positive alleles at grain size QTLs. Candidate genes identified using genome-wide variant calling data for all parental inbred lines indicated overlapping and potential novel regulators of grain size in cereals. Furthermore, our results indicated that sink capacity was the primary determinant of grain weight in barley.

Multiplexed promoter and gene editing in wheat using a virus-based guide RNA delivery system

The low efficiency of genetic transformation and gene editing across diverse cultivars hinder the broad application of CRISPR technology for crop improvement. The development of virus-based methods of CRISPR-Cas system delivery into the plant cells holds great promise to overcome these limitations. Here, we perform direct inoculation of wheat leaves with the barley stripe mosaic virus (BSMV) transcripts to deliver guide RNAs (sgRNA) into the Cas9-expressing wheat. We demonstrate that wheat inoculation with the pool of BSMV-sgRNAs could be used to generate heritable precise deletions in the promoter region of a transcription factor and to perform multiplexed editing of agronomic genes. We transfer the high-expressing locus of Cas9 into adapted spring and winter cultivars by marker-assisted introgression and use of the BSMV-sgRNAs to edit two agronomic genes. A strategy presented in our study could be applied to any adapted cultivar for creating new cis-regulatory diversity or large-scale editing of multiple genes in biological pathways or QTL regions, opening possibilities for the effective engineering of crop genomes, and accelerating gene discovery and trait improvement efforts.

gw2.1, a new allele of GW2, improves grain weight and grain yield in rice

Grain weight is an important characteristic of grain shape and a key contributing factor to the grain yield in rice. Here, we report that gw2.1, a new allele of the Grain Width and Weight 2 (GW2) gene, regulates grain size and grain weight. A single nucleotide substitution in the coding sequence (CDS) of gw2.1 resulted in the change of glutamate to lysine (E128K) in GW2.1 protein. Complementation tests and GW2 overexpression experiments demonstrated that the missense mutation in gw2.1 was responsible for the phenotype of enlarged grain size in the mutant line jf42. The large grain trait of the near-isogenic line NIL-gw2.1 was found to result from increased cell proliferation during flower development. Meanwhile, NIL-gw2.1 was shown to increase grain yield without compromising the grain quality. The GW2 protein was localized to the cell nucleus and membrane, and interacted with CHB705, a subunit of the chromatin remodeling complex. Finally, the F₁ hybrids from crosses of NIL-gw2.1 with 7 cytoplasmic male-sterile lines exhibited large grains and desirable grain appearance. Thus, gw2.1 is a promising allele that could be applied to improve grain yield and grain appearance in rice. AVAILABILITY OF DATA AND MATERIALS: The datasets generated and/or analyzed in the study are available from the corresponding author on reasonable request.

Control of Grain Weight and Size in Rice (Oryza sativa L.) by OsPUB3 Encoding a U-Box E3 Ubiquitin Ligase

Grain weight and size, mostly determined by grain length, width and thickness, are crucial traits affecting grain quality and yield in rice. A quantitative trait locus controlling grain length and width in rice, qGS1-35.2, was previously fine-mapped in a 57.7-kb region on the long arm of chromosome 1. In this study, OsPUB3, a gene encoding a U-box E3 ubiquitin ligase, was validated as the causal gene for qGS1-35.2. The effects were confirmed firstly by using CRISPR/Cas9-based mutagenesis and then through transgenic complementation of a Cas9-free knock-out (KO) mutant. Two homozygous KO lines were produced, each having a 1-bp insertion in OsPUB3 which caused frameshift mutation and premature termination. Compared with the recipient and a transgenic-negative control, both mutants showed significant decreases in grain weight and size. In transgenic complementation populations derived from four independent T₀ plants, grain weight of transgenic-positive plants was significantly higher than transgenic-negative plants, coming with increased grain length and a less significant decrease in grain width. Based on data documented in RiceVarMap V2.0, eight haplotypes were classified according to six single-nucleotide polymorphisms (SNPs) found in the OsPUB3 coding region of 4695 rice accessions. Significant differences on grain size traits were detected between the three major haplotypes, Hap1, Hap2 and Hap3 that jointly occupy 98.6% of the accessions. Hap3 having the largest grain weight and grain length but intermediate grain width exhibits a potential for simultaneously improving grain yield and quality. In another set of 257 indica rice cultivars tested in our study, Hap1 and Hap2 remained to be the two largest groups. Their differences on grain weight and size were significant in the background of non-functional gse5, but non-significant in the background of functional GSE5, indicating a genetic interaction between OsPUB3 and GSE5. Cloning of OsPUB3 provides a new gene resource for investigating the regulation of grain weight and size.

Genetic background- and environment-independent QTL and candidate gene identification of appearance quality in three MAGIC populations of rice

Many QTL have been identified for grain appearance quality by linkage analysis (LA) in bi-parental mapping populations and by genome-wide association study (GWAS) in natural populations in rice. However, few of the well characterized genes/QTL have been successfully applied in molecular rice breeding due to genetic background (GB) and environment effects on QTL expression and deficiency of favorable alleles. In this study, GWAS and LA were performed to identify QTL for five grain appearance quality-related traits using three multi-parent advanced generation inter-cross (MAGIC) populations. A total of 22 QTL on chromosomes 1-3, 5-8 were identified by GWAS for five traits in DC1, DC2 and 8way, and four combined populations DC12 (DC1+DC2), DC18 (DC1+8way), DC28 (DC2+8way) and DC128 (DC1+DC2+8way). And a total of 42 QTL were identified on all 12 chromosomes except 10 by LA in the three single populations. Among 20 QTL identified by GWAS in DC1, DC2 and 8way, 10, four and three QTL were commonly detected in DC18, DC28, and DC128, respectively. Similarly, among 42 QTL detected by LA in the three populations, four, one and two QTL were commonly detected in DC18, DC28, and DC128, respectively. There was no QTL mapped together in DC12 by both two mapping methods, indicating that GB could greatly affect the mapping results, and it was easier to map the common QTL among populations with similar GB. The 8way population was more powerful for QTL mapping than the DC1, DC2 and various combined populations. Compared with GWAS, LA can not only identify large-effect QTL, but also identify minor-effect ones. Among 11 QTL simultaneously detected by the two methods in different GBs and environments, eight QTL corresponded to known genes, including AqGL3b and AqGLWR3a for GL and GLWR, AqGW5a, AqGLWR5, AqDEC5 and AqPGWC5 for GW, GLWR, DEC and PGWC, and AqDEC6b and AqPGWC6b for DEC and PGWC, respectively. AqGL7, AqGL3c/AqGLWR3b, AqDEC6a/AqPGWC6a, and AqPGWC7 were newly identified and their candidate genes were analyzed and inferred. It was discussed to further improve grain appearance quality through designed QTL pyramiding strategy based on the stable QTL identified in the MAGIC populations.

Phenotypic Characterization and Fine Mapping of a Major-Effect Fruit Shape QTL FS5.2 in Cucumber, Cucumis sativus L., with Near-Isogenic Line-Derived Segregating Populations

Cucumber (Cucumis sativus L.) fruit size/shape (FS) is an important yield and quality trait that is quantitatively inherited. Many quantitative trait loci (QTLs) for fruit size/shape have been identified, but very few have been fine-mapped or cloned. In this study, through marker-assisted foreground and background selections, we developed near-isogenic lines (NILs) for a major-effect fruit size/shape QTL FS5.2 in cucumber. Morphological and microscopic characterization of NILs suggests that the allele of fs5.2 from the semi-wild Xishuangbanna (XIS) cucumber (C. s. var. xishuangbannesis) reduces fruit elongation but promotes radial growth resulting in shorter but wider fruit, which seems to be due to reduced cell length, but increased cellular layers. Consistent with this, the NIL carrying the homozygous XIS allele (fs5.2) had lower auxin/IAA contents in both the ovary and the developing fruit. Fine genetic mapping with NIL-derived segregating populations placed FS5.2 into a 95.5 kb region with 15 predicted genes, and a homolog of the Arabidopsis CRABS CLAW (CsCRC) appeared to be the most possible candidate for FS5.2. Transcriptome profiling of NIL fruits at anthesis identified differentially expressed genes enriched in the auxin biosynthesis and signaling pathways, as well as genes involved in cell cycle, division, and cell wall processes. We conclude that the major-effect QTL FS5.2 controls cucumber fruit size/shape through regulating auxin-mediated cell division and expansion for the lateral and longitudinal fruit growth, respectively. The gibberellic acid (GA) signaling pathway also plays a role in FS5.2-mediated fruit elongation.

Fine mapping of qTGW2b and qGL9, two minor QTL conferring grain size and weight in rice

Rice grain size is a key determinant of both grain yield and quality. In this study, we conducted QTL mapping on grain size using a recombinant inbred line (RIL) population derived from a cross between japonica variety Beilu130 (BL130) and indica variety Jin23B (J23B). A total of twenty-two QTL related to grain length (GL), grain width (GW), grain length-to-width ratio (LWR), grain thickness (GT), and thousand grain weight (TGW) were detected under two environments, and 14 of them were repeatedly detected. Two minor QTL, qTGW2b and qGL9, were validated and further delimited to regions of 631 kb and 272 kb, respectively. Parental sequence comparison of genes expressed in inflorescence in corresponding candidate regions identified frameshifts in the exons of LOC_Os02g38690 and LOC_Os02g38780, both of which encode protein phosphatase 2C-containing protein, and LOC_Os09g29930, which encodes a BIM2 protein. Scanning electron microscopy (SEM) analysis revealed that the increase of cell size rather than cell number caused the differences in grain size between NILs of qTGW2b and qGL9. Quantitative RT-PCR analysis showed that the expression levels of EXPA4, EXPA5, EXPA6, EXPB3, EXPB4, and EXPB7 were significantly different in both qTGW2b NILs and qGL9 NILs. Our results lay the foundation for the cloning of qTGW2b and qGL9, and provide genetic materials for the improvement of rice yield and quality.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01328-2.

Developing Genetic Engineering Techniques for Control of Seed Size and Yield

Many signaling pathways regulate seed size through the development of endosperm and maternal tissues, which ultimately results in a range of variations in seed size or weight. Seed size can be determined through the development of zygotic tissues (endosperm and embryo) and maternal ovules. In addition, in some species such as rice, seed size is largely determined by husk growth. Transcription regulator factors are responsible for enhancing cell growth in the maternal ovule, resulting in seed growth. Phytohormones induce significant effects on entire features of growth and development of plants and also regulate seed size. Moreover, the vegetative parts are the major source of nutrients, including the majority of carbon and nitrogen-containing molecules for the reproductive part to control seed size. There is a need to increase the size of seeds without affecting the number of seeds in plants through conventional breeding programs to improve grain yield. In the past decades, many important genetic factors affecting seed size and yield have been identified and studied. These important factors constitute dynamic regulatory networks governing the seed size in response to environmental stimuli. In this review, we summarized recent advances regarding the molecular factors regulating seed size in Arabidopsis and other crops, followed by discussions on strategies to comprehend crops' genetic and molecular aspects in balancing seed size and yield.

Identification and Investigation of the Genetic Variations and Candidate Genes Responsible for Seed Weight via GWAS in Paper Mulberry

Seeds directly determine the survival and population size of woody plants, but the genetic basis of seed weight in woody plants remain poorly explored. To identify genetic variations and candidate genes responsible for seed weight in natural woody populations, we investigated the hundred-seed weight of 198 paper mulberry individuals from different areas. Our results showed that the hundred-seed weight of paper mulberry was significantly associated with the bioclimatic variables of sampling sites, which increased from south to north along the latitudinal-temperature gradient. Using 2,414,978 high-quality SNPs from re-sequencing data, the genome-wide association analysis of the hundred-seed weight was performed under three models, which identified 148, 19 and 12 associated genes, respectively. Among them, 25 candidate genes were directly hit by the significant SNPs, including the WRKY transcription factor, fatty acid desaturase, F-box protein, etc. Most importantly, we identified three crucial genetic variations in the coding regions of candidate genes (Bp02g2123, Bp01g3291 and Bp10g1642), and significant differences in the hundred-seed weight were detected among the individuals carrying different genotypes. Further analysis revealed that Bp02g2123 encoding a fatty acid desaturase (FAD) might be a key factor affecting the seed weight and local climate adaptation of woody plants. Furthermore, the genome-wide investigation and expression analysis of FAD genes were performed, and the results suggested that BpFADs widely expressed in various tissues and responded to multiple phytohormone and stress treatments. Overall, our study identifies valuable genetic variations and candidate genes, and provides a better understanding of the genetic basis of seed weight in woody plants.

Genome-wide association studies provide genetic insights into natural variation of seed-size-related traits in mungbean

Although mungbean (Vigna radiata (L.) R. Wilczek) is an important legume crop, its seed yield is relatively low. To address this issue, here 196 accessions with 3,607,508 SNP markers were used to identify quantitative trait nucleotides (QTNs), QTN-by-environment interactions (QEIs), and their candidate genes for seed length (SL), seed width, and 100-seed weight (HSW) in two environments. As a result, 98 QTNs and 20 QEIs were identified using 3VmrMLM, while 95, >10,000, and 15 QTNs were identified using EMMAX, GEMMA, and CMLM, respectively. Among 809 genes around these QTNs, 12 were homologous to known seed-development genes in rice and Arabidopsis thaliana, in which 10, 2, 1, and 0 genes were found, respectively, by the above four methods to be associated with the three traits, such as VrEmp24/25 for SL and VrKIX8 for HSW. Eight of the 12 genes were significantly differentially expressed between two large-seed and two small-seed accessions, and VrKIX8, VrPAT14, VrEmp24/25, VrIAR1, VrBEE3, VrSUC4, and Vrflo2 were further verified by RT-qPCR. Among 65 genes around these QEIs, VrFATB, VrGSO1, VrLACS2, and VrPAT14 were homologous to known seed-development genes in A. thaliana, although new experiments are necessary to explore these novel GEI-trait associations. In addition, 54 genes were identified in comparative genomics analysis to be associated with seed development pathway, in which VrKIX8, VrABA2, VrABI5, VrSHB1, and VrIKU2 were also identified in genome-wide association studies. This result provided a reliable approach for identifying seed-size-related genes in mungbean and a solid foundation for further molecular biology research on seed-size-related genes.

Map-based cloning and transcriptome analysis of the more-tiller and small-grain mutant in rice

A G to T nucleotide substitution of OsTSG2 led to more tillers and smaller grains in rice by participating in phytohormone signal transduction and starch and sucrose metabolism. Rice is one of the most important food crops worldwide. Grain size and tiller number are the most important factors determining rice yield. The more-tiller and small-grain 2 (tsg2) mutant in rice, developed by ethyl methanesulfonate (EMS) mutagenesis, has smaller grains, more tillers, and a higher yield per plant relative to the wild-type (WT). Based on the genetic analysis, the tsg2 traits were conferred by a single recessive nuclear gene located on the long arm of chromosome 2. After fine-mapping the OsTSG2 locus, a G to T nucleotide substitution was identified, which resulted in an A to S mutation in a highly conserved domain of the growth-regulation factor protein. The single-strand conformation polymorphism (SSCP) marker was developed based on the SNP associated with the phenotypic segregation of traits. The functional complementation of OsTSG2 from the tsg2 mutant to the WT led to an increase in grain size and weight. The differentially expressed genes (DEGs) identified by RNA sequencing were involved in phytohormone signal transduction and starch and sucrose metabolism. Enzyme-linked immunosorbent assay (ELISA) analysis detected variation in the indole acetic acid (IAA) and jasmonic acid (JA) content in the tsg2 inflorescence, while the cellular organization, degree of chalkiness, gel consistency, amylose content, and alkaline spreading value were affected in the tsg2 grains. The findings elucidated the regulatory mechanisms of the tsg2 traits. This mutant could be used in marker-assisted breeding for high-yield and good-quality rice.

Role of cytokinins in seed development in pulses and oilseed crops: Current status and future perspective

Cytokinins constitutes a vital group of plant hormones regulating several developmental processes, including growth and cell division, and have a strong influence on grain yield. Chemically, they are the derivatives of adenine and are the most complex and diverse group of hormones affecting plant physiology. In this review, we have provided a molecular understanding of the role of cytokinins in developing seeds, with special emphasis on pulses and oilseed crops. The importance of cytokinin-responsive genes including cytokinin oxidases and dehydrogenases (CKX), isopentenyl transferase (IPT), and cytokinin-mediated genetic regulation of seed size are described in detail. In addition, cytokinin expression in germinating seeds, its biosynthesis, source-sink dynamics, cytokinin signaling, and spatial expression of cytokinin family genes in oilseeds and pulses have been discussed in context to its impact on increasing economy yields. Recently, it has been shown that manipulation of the cytokinin-responsive genes by mutation, RNA interference, or genome editing has a significant effect on seed number and/or weight in several crops. Nevertheless, the usage of cytokinins in improving crop quality and yield remains significantly underutilized. This is primarily due to the multigene control of cytokinin expression. The information summarized in this review will help the researchers in innovating newer and more efficient ways of manipulating cytokinin expression including CKX genes with the aim to improve crop production, specifically of pulses and oilseed crops.

Identification and fine-mapping of a major QTL (PH1.1) conferring plant height in broomcorn millet (Panicum miliaceum)

The plant height of broomcorn millet (Panicum miliaceum) is a significant agronomic trait that is closely related to its plant architecture, lodging resistance, and final yield. However, the genes underlying the regulation of plant height in broomcorn millet are rarely reported. Here, an F2 population derived from a cross between a normal variety, "Longmi12," and a dwarf mutant, "Zhang778," was constructed. Genetic analysis for the F2 and F2:3 populations revealed that the plant height was controlled by more than one locus. A major quantitative trait locus (QTL), PH1.1, was preliminarily identified in chromosome 1 using bulked segregant analysis sequencing (BSA-seq). PH1.1 was fine-mapped to a 109-kb genomic region with 15 genes using a high-density map. Among them, longmi011482 and longmi011489, containing nonsynonymous variations in their coding regions, and longmi011496, covering multiple insertion/deletion sequences in the promoter regions, may be possible candidate genes for PH1.1. Three diagnostic markers closely linked to PH1.1 were developed to validate the PH1.1 region in broomcorn millet germplasm. These findings laid the foundation for further understanding of the molecular mechanism of plant height regulation in broomcorn millet and are also beneficial to the breeding program for developing new varieties with optimal height.

Genetic dissection of grain traits and their corresponding heterosis in an elite hybrid

Rice productivity has considerably improved due to the effective employment of heterosis, but the genetic basis of heterosis for grain shape and weight remains uncertain. For studying the genetic dissection of heterosis for grain shape/weight and their relationship with grain yield in rice, quantitative trait locus (QTL) mapping was performed on 1,061 recombinant inbred lines (RILs), which was developed by crossing xian/indica rice Quan9311B (Q9311B) and Wu-shan-si-miao (WSSM). Whereas, BC₁F₁ (a backcross F₁) was developed by crossing RILs with Quan9311A (Q9311A) combined with phenotyping in Hefei (HF) and Nanning (NN) environments. Overall, 114 (main-effect, mQTL) and 359 (epistatic QTL, eQTL) were identified in all populations (RIL, BC₁F₁, and mid-parent heterosis, H_MPs) for 1000-grain weight (TGW), grain yield per plant (GYP) and grain shape traits including grain length (GL), grain width (GW), and grain length to width ratio (GLWR). Differential QTL detection revealed that all additive loci in RILs population do not show heterotic effects, and few of them affect the performance of BC₁F₁. However, 25 mQTL not only contributed to BC₁F₁'s performance but also contributed to heterosis. A total of seven QTL regions was identified, which simultaneously affected multiple grain traits (grain yield, weight, shape) in the same environment, including five regions with opposite directions and two regions with same directions of favorable allele effects, indicating that partial genetic overlaps are existed between different grain traits. This study suggested different approaches for obtaining good grain quality with high yield by pyramiding or introgressing favorable alleles (FA) with the same direction of gene effect at the QTL regions affecting grain shape/weight and grain yield distributing on different chromosomes, or introgressing or pyramiding FA in the parents instead of fixing additive effects in hybrid as well as pyramiding the polymorphic overdominant/dominant loci between the parents and eliminating underdominant loci from the parents. These outcomes offer valuable information and strategy to develop hybrid rice with suitable grain type and weight.

Ribonuclease H-like gene SMALL GRAIN2 regulates grain size in rice through brassinosteroid signaling pathway

Grain size is a key agronomic trait that determines the yield in plants. Regulation of grain size by brassinosteroids (BRs) in rice has been widely reported. However, the relationship between the BR signaling pathway and grain size still requires further study. Here, we isolated a rice mutant, named small grain2 (sg2), which displayed smaller grain and a semi-dwarf phenotype. The decreased grain size was caused by repressed cell expansion in spikelet hulls of the sg2 mutant. Using map-based cloning combined with a MutMap approach, we cloned SG2, which encodes a plant-specific protein with a ribonuclease H-like domain. SG2 is a positive regulator downstream of GLYCOGEN SYNTHASE KINASE2 (GSK2) in response to BR signaling, and its mutation causes insensitivity to exogenous BR treatment. Genetical and biochemical analysis showed that GSK2 interacts with and phosphorylates SG2. We further found that BRs enhance the accumulation of SG2 in the nucleus, and subcellular distribution of SG2 is regulated by GSK2 kinase activity. In addition, Oryza sativa OVATE family protein 19 (OsOFP19), a negative regulator of grain shape, interacts with SG2 and plays an antagonistic role with SG2 in controlling gene expression and grain size. Our results indicated that SG2 is a new component of GSK2-related BR signaling response and regulates grain size by interacting with OsOFP19.

Short grain 5 controls grain length in rice by regulating cell expansion

Grain shape is a crucial determinant of grain weight and quality and plays a vital role in rice breeding. Although many grain shape-related genes have been reported, the regulatory relationship between them has not been well characterized in rice. In this study, we report the isolation of a short-grain-length mutant called sg5 from the heavy-panicle-type hybrid rice elite restorer line 'ShuhuiR498' (R498) after ethyl methanesulfonate (EMS) treatment. MutMap cloning revealed that SG5 encodes a Myb-like transcription factor. A missense mutation in the first exon of SG5 was found to cause an amino acid change from leucine to proline at position 197 in the mutant SG5 protein. Gene knockout and genetic complementation experiments confirmed that the point mutation in SG5 was responsible for the sg5 mutant phenotype. SG5 is mainly expressed in young panicles and hulls. In addition, the SG5 protein is found in the nucleus and does not affect subcellular localization. Histochemical observation and gene expression analysis indicated that SG5 regulates spikelet hull development by mediating cell expansion. Moreover, the expression levels of BG1, GS2, and DEP1 were reduced in sg5 plants, and dual-luciferase (LUC) assays showed that SG5 can bind to the BG1 gene promoter. The effect of pyramiding sg5 and GS3 suggests that sg5 and GS3 regulate grain length independently. The results of our study show that the missense mutation in sg5 is essential for the molecular function of SG5 and SG5 is involved in regulating cell expansion and expression of grain-shape-related genes to regulate grain length. This work provides new data to help study and understand the molecular function of SG5.

Natural variation of GhSI7 increases seed index in cotton

qSI07.1, a major QTL for seed index in cotton, was fine-mapped to a 17.45-kb region, and the candidate gene GhSI7 was verified in transgenic plants. Improving production to meet human needs is a vital objective in cotton breeding. The yield-related trait seed index is a complex quantitative trait, but few candidate genes for seed index have been characterized. Here, a major QTL for seed index qSI07.1 was fine-mapped to a 17.45-kb region by linkage analysis and substitutional mapping. Only GhSI7, encoding the transcriptional regulator STERILE APETALA, was contained in the candidate region. Association test and genetic analysis indicated that an 845-bp-deletion in its intron was responsible for the seed index variation. Origin analysis revealed that this variation was unique in Gossypium hirsutum and originated from race accessions. Overexpression of GhSI7 (haplotype 2) significantly increased the seed index and organ size in cotton plants. Our findings provided a diagnostic marker for breeding and selecting cotton varieties with high seed index, and laid a foundation for further studies to understand the molecular mechanism of cotton seed morphogenesis.

Rice co-expression network analysis identifies gene modules associated with agronomic traits

Identifying trait-associated genes is critical for rice (Oryza sativa) improvement, which usually relies on map-based cloning, quantitative trait locus analysis, or genome-wide association studies. Here we show that trait-associated genes tend to form modules within rice gene co-expression networks, a feature that can be exploited to discover additional trait-associated genes using reverse genetics. We constructed a rice gene co-expression network based on the graphical Gaussian model using 8,456 RNA-seq transcriptomes, which assembled into 1,286 gene co-expression modules functioning in diverse pathways. A number of the modules were enriched with genes associated with agronomic traits, such as grain size, grain number, tiller number, grain quality, leaf angle, stem strength, and anthocyanin content, and these modules are considered to be trait-associated gene modules. These trait-associated gene modules can be used to dissect the genetic basis of rice agronomic traits and to facilitate the identification of trait genes. As an example, we identified a candidate gene, OCTOPUS-LIKE 1 (OsOPL1), a homolog of the Arabidopsis (Arabidopsis thaliana) OCTOPUS gene, from a grain size module and verified it as a regulator of grain size via functional studies. Thus, our network represents a valuable resource for studying trait-associated genes in rice.

Florigen repression complexes involving rice CENTRORADIALIS2 regulate grain size

Grain size is one of the crucial factors determining grain yield. However, the genetic and molecular mechanisms of florigen repression complexes (FRCs) underlying grain size in rice (Oryza sativa L.) have not been reported. Here, we report that the rice CENTRORADIALIS (CEN) family member OsCEN2 (also known as Rice TFL1/CEN homolog, RCN1), a phosphatidylethanolamine-binding protein (PEBP) family protein, negatively controls grain size in rice. Overexpression of OsCEN2 led to small grains, and knockout of OsCEN2 resulted in large, heavy grains. OsCEN2 influenced grain size by restricting cell expansion in the spikelet hull and seed filling. In in vivo and in vitro experiments, OsCEN2 physically interacted with a G-box factor 14-3-3 homolog, GF14f, which negatively regulates grain size. Bimolecular fluorescence complementation and yeast two-hybrid assays revealed that GF14f directly interacts with the basic leucine zipper (bZIP) transcription factor, OsFD2. Plants overexpressing OsFD2 produced smaller and lighter grains than wild-type plants. We found that OsFD2 also influences grain size by controlling cell expansion and division in the spikelet hull. Our results reveal the molecular mechanisms of the OsCEN2-GF14f-OsFD2 regulatory module in controlling grain size. Additionally, our study provides insight into the functions of the FRC in rice and suggests a strategy for improving seed size and weight.

Control of Grain Shape and Size in Rice by Two Functional Alleles of OsPUB3 in Varied Genetic Background

Grain shape and size are key determinants of grain appearance quality and yield in rice. In our previous study, a grain shape QTL, qGS1-35.2, was fine-mapped using near-isogenic lines (NILs) derived from a cross between Zhenshan 97 (ZS97) and Milyang 46 (MY46). One annotated gene, OsPUB3, was found to be the most likely candidate gene. Here, knockout and overexpression experiments were performed to investigate the effects of OsPUB3 on grain shape and size. Four traits were tested, including grain length, grain width, grain weight, and the ratio of grain length to width. Knockout of OsPUB3 in NIL^ZS97, NIL^MY46, and another rice cultivar carrying the OsPUB3^MY46 allele all caused decreases in grain width and weight and increases in the ratio of grain length to width. Results also showed that the magnitude of the mutational effects varied depending on the target allele and the genetic background. Moreover, it was found that NIL^ZS97 and NIL^MY46 carried different functional alleles of OsPUB3, causing differences in grain shape rather than grain weight. In the overexpression experiment, significant differences between transgenic-positive and transgenic-negative plants were detected in all four traits. These results indicate that OsPUB3 regulates grain shape and size through a complex mechanism and is a good target for deciphering the regulatory network of grain shape. This gene could be used to improve grain appearance quality through molecular breeding as well.

Low grain weight, a new allele of BRITTLE CULM12, affects grain size through regulating GW7 expression in rice

Grain weight is a major determinant in rice yield, which is tightly associated with grain size. However, the underlying molecular mechanisms that control this trait remain unclear. Here, we report a rice (Oryza sativa) mutant, low grain weight (lgw), which shows that reduced grain length is caused by decreased cell elongation and proliferation. Map-based cloning revealed that all mutant phenotypes resulted from a nine-base pair (bp) deletion in LGW, which encodes the kinesin-like protein BRITTLE CULM12 (BC12). Protein sequence alignment analysis revealed that the mutation site was located at the nuclear localization signal (NLS) of LGW/BC12, resulting in the lgw protein not being located in the nucleus. LGW is preferentially expressed in both culms and roots, as well as in the early developing panicles. Overexpression of LGW increased the grain length, indicating that LGW is a positive regulator for regulating grain length. In addition, LGW/BC12 is directly bound to the promoter of GW7 and activates its expression. Elevating the GW7 expression levels in lgw plants rescued the small grain size phenotype. We conclude that LGW regulates grain development by directly binding to the GW7 promoter and activating its expression. Our findings revealed that LGW plays an important role in regulating grain size, and manipulation of this gene provides a new strategy for regulating grain weight in rice.

Discovery and Validation of Grain Shape Loci in U.S. Rice Germplasm Through Haplotype Characterization

Rice grain shape is a major determinant of rice market value and the end-use. We mapped quantitative trait loci (QTL) for grain shape traits in a bi-parental recombinant inbred line population (Trenasse/Jupiter) and discovered two major grain length QTLs—qGL3.1 and qGL7.1. Previously, a major grain shape gene GS3 was reported in the qGL3.1 region and grain length gene GL7 was reported to be encompassing qGL7.1 locus. The re-sequencing SNP data on the International Rice Research Institute (IRRI) 3K Rice Genome Project (RGP) panel were obtained from the IRRI SNP-Seek database for both genes and haplotype diversity was characterized for each gene in this diverse panel. United States rice germplasm was not well represented in the IRRI 3K RGP database. Therefore, a minimum SNP set was identified for each gene that could differentiate all the characterized haplotypes. These haplotypes in the 3K RGP panel were screened across 323 elite U.S. genotypes using the minimum SNP set. The screening of haplotypes and phenotype association confirmed the role of GS3 under qGL3.1. However, screening of the GL7 haplotypes in the U.S. germplasm panel showed that GL7 did not play a role in qGL7.1, and in addition, GL7.1 did not segregate in the Trenasse/Jupiter RIL population. This concluded that qGL7.1 is a novel QTL discovered on chr7 for grain shape in the Trenasse/Jupiter RIL population. A high-throughput KASP-based SNP marker for each locus (GS3 and qGL7.1) was identified and validated in elite U.S. rice germplasm to be used in an applied rice breeding program.

Genome-wide analysis of the JAZ subfamily of transcription factors and functional verification of BnC08.JAZ1-1 in Brassica napus

BACKGROUND

JAZ subfamily plays crucial roles in growth and development, stress, and hormone responses in various plant species. Despite its importance, the structural and functional analyses of the JAZ subfamily in Brassica napus are still limited.

RESULTS

Comparing to the existence of 12 JAZ genes (AtJAZ1-AtJAZ12) in Arabidopsis, there are 28, 31, and 56 JAZ orthologues in the reference genome of B. rapa, B. oleracea, and B. napus, respectively, in accordance with the proven triplication events during the evolution of Brassicaceae. The phylogenetic analysis showed that 127 JAZ proteins from A. thaliana, B. rapa, B. oleracea, and B. napus could fall into five groups. The structure analysis of all 127 JAZs showed that these proteins have the common motifs of TIFY and Jas, indicating their conservation in Brassicaceae species. In addition, the cis-element analysis showed that the main motif types are related to phytohormones, biotic and abiotic stresses. The qRT-PCR of the representative 11 JAZ genes in B. napus demonstrated that different groups of BnJAZ individuals have distinct patterns of expression under normal conditions or treatments with distinctive abiotic stresses and phytohormones. Especially, the expression of BnJAZ52 (BnC08.JAZ1-1) was significantly repressed by abscisic acid (ABA), gibberellin (GA), indoleacetic acid (IAA), polyethylene glycol (PEG), and NaCl treatments, while induced by methyl jasmonate (MeJA), cold and waterlogging. Expression pattern analysis showed that BnC08.JAZ1-1 was mainly expressed in the vascular bundle and young flower including petal, pistil, stamen, and developing ovule, but not in the stem, leaf, and mature silique and seed. Subcellular localization showed that the protein was localized in the nucleus, in line with its orthologues in Arabidopsis. Overexpression of BnC08.JAZ1-1 in Arabidopsis resulted in enhanced seed weight, likely through regulating the expression of the downstream response genes involved in the ubiquitin-proteasome pathway and phospholipid metabolism pathway.

CONCLUSIONS

The systematic identification, phylogenetic, syntenic, and expression analyses of BnJAZs subfamily improve our understanding of their roles in responses to stress and phytohormone in B. napus. In addition, the preliminary functional validation of BnC08.JAZ1-1 in Arabidopsis demonstrated that this subfamily might also play a role in regulating seed weight.

OsBSK2, a putative brassinosteroid-signalling kinase, positively controls grain size in rice

Grain size is an important trait that directly affects grain yield in rice; however, the genetic and molecular mechanisms regulating grain size remain unclear. In this study, we identified a mutant, grain length and grain weight 10 (glw10), which exhibited significantly reduced grain length and grain weight. Histological analysis demonstrated that GLW10 affects cell expansion, which regulates grain size. MutMap-based gene mapping and transgenic experiments demonstrated that GLW10 encodes a putative brassinosteroid (BR) signalling kinase, OsBSK2. OsBSK2 is a plasma membrane protein, and an N-myristoylation site is needed for both membrane localization and function. OsBSK2 directly interacts with the BR receptor kinase OsBRI1; however, genetic experiments have demonstrated that OsBSK2 may regulate grain size independent of the BR signalling pathway. OsBSK2 can form a homodimer or heterodimer with OsBSK3 and OsBSK4, and silencing OsBSK2, OsBSK3, and OsBSK4 reduce grain size. This indicates that OsBSKs seem to function as homodimers or heterodimers to positively regulate grain size in rice. OsBSK2/3/4 are all highly expressed in young panicles and spikelet hulls, suggesting that they control grain size. In summary, our results provide novel insights into the function of BSKs in rice, and identify novel targets for improving grain size during crop breeding.

cis-Regulatory variation affecting gene expression contributes to the improvement of maize kernel size

cis-Regulatory variations contribute to trait evolution and adaptation during crop domestication and improvement. As the most important harvested organ in maize (Zea mays L.), kernel size has undergone intensive selection for size. However, the associations between maize kernel size and cis-regulatory variations remain unclear. We chose two independent association populations to dissect the genetic architecture of maize kernel size together with transcriptomic and genotypic data. The resulting phenotypes reflected a strong influence of population structure on kernel size. Compared with genome-wide association studies (GWASs), which accounted for population structure and relatedness, GWAS based on a naïve or simple linear model revealed additional associated single-nucleotide polymorphisms significantly involved in the conserved pathways controlling seed size in plants. Regulation analyses through expression quantitative trait locus mapping revealed that cis-regulatory variations likely control kernel size by fine-tuning the expression of proximal genes, among which ZmKL1 (GRMZM2G098305) was transgenically validated. We also proved that the pyramiding of the favorable cis-regulatory variations has contributed to the improvement of maize kernel size. Collectively, our results demonstrate that cis-regulatory variations, together with their regulatory genes, provide excellent targets for future maize improvement.

Emerging roles of the ubiquitin-proteasome pathway in enhancing crop yield by optimizing seed agronomic traits

Ubiquitin-proteasome pathway has the potential to modulate crop productivity by influencing agronomic traits. Being sessile, the plant often uses the ubiquitin-proteasome pathway to maintain the stability of different regulatory proteins to survive in an ever-changing environment. The ubiquitin system influences plant reproduction, growth, development, responses to the environment, and processes that control critical agronomic traits. E3 ligases are the major players in this pathway, and they are responsible for recognizing and tagging the targets/substrates. Plants have a variety of E3 ubiquitin ligases, whose functions have been studied extensively, ranging from plant growth to defense strategies. Here we summarize three agronomic traits influenced by ubiquitination: seed size and weight, seed germination, and accessory plant agronomic traits particularly panicle architecture, tillering in rice, and tassels branch number in maize. This review article highlights some recent progress on how the ubiquitin system influences the stability/modification of proteins that determine seed agronomic properties like size, weight, germination and filling, and ultimately agricultural productivity and quality. Further research into the molecular basis of the aforementioned processes might lead to the identification of genes that could be modified or selected for crop development. Likewise, we also propose advances and future perspectives in this regard.

Fine mapping and characterization of a major QTL for grain weight on wheat chromosome arm 5DL

We fine mapped QTL QTKW.caas-5DL for thousand kernel weight in wheat, predicted candidate genes and developed a breeding-applicable marker. Thousand kernel weight (TKW) is an important yield component trait in wheat, and identification of the underlying genetic loci is helpful for yield improvement. We previously identified a stable quantitative trait locus (QTL) QTKW.caas-5DL for TKW in a Doumai/Shi4185 recombinant inbred line (RIL) population. Here we performed fine mapping of QTKW.caas-5DL using secondary populations derived from 15 heterozygous recombinants and delimited the QTL to an approximate 3.9 Mb physical interval from 409.9 to 413.8 Mb according to the Chinese Spring (CS) reference genome. Analysis of genomic synteny showed that annotated genes in the physical interval had high collinearity among CS and eight other wheat genomes. Seven genes with sequence variation and/or differential expression between parents were predicted as candidates for QTKW.caas-5DL based on whole-genome resequencing and transcriptome assays. A kompetitive allele-specific PCR (KASP) marker for QTKW.caas-5DL was developed, and genotyping confirmed a significant association with TKW but not with other yield component traits in a panel of elite wheat cultivars. The superior allele of QTKW.caas-5DL was frequent in a panel of cultivars, suggesting that it had undergone positive selection. These findings not only lay a foundation for map-based cloning of QTKW.caas-5DL but also provide an efficient tool for marker-assisted selection.

Natural allelic variation of GmST05 controlling seed size and quality in soybean

Seed size is one of the most important agronomic traits determining the yield of crops. Cloning the key genes controlling seed size and pyramiding their elite alleles will facilitate yield improvement. To date, few genes controlling seed size have been identified in soybean, a major crop that provides half of the plant oil and one quarter of the plant protein globally. Here, through a genome-wide association study of over 1800 soybean accessions, we determined that natural allelic variation at GmST05 (Seed Thickness 05) predominantly controlled seed thickness and size in soybean germplasm. Further analyses suggested that the two major haplotypes of GmST05 differed significantly at the transcriptional level. Transgenic experiments demonstrated that GmST05 positively regulated seed size and influenced oil and protein contents, possibly by regulating the transcription of GmSWEET10a. Population genetic diversity analysis suggested that allelic variations of GmST05 were selected during geographical differentiation but have not been fixed. In summary, natural variation in GmST05 determines transcription levels and influences seed size and quality in soybean, making it an important gene resource for soybean molecular breeding.

FIND-IT: Accelerated trait development for a green evolution

Improved agricultural and industrial production organisms are required to meet the future global food demands and minimize the effects of climate change. A new resource for crop and microbe improvement, designated FIND-IT (Fast Identification of Nucleotide variants by droplet DigITal PCR), provides ultrafast identification and isolation of predetermined, targeted genetic variants in a screening cycle of less than 10 days. Using large-scale sample pooling in combination with droplet digital PCR (ddPCR) greatly increases the size of low-mutation density and screenable variant libraries and the probability of identifying the variant of interest. The method is validated by screening variant libraries totaling 500,000 barley (Hordeum vulgare) individuals and isolating more than 125 targeted barley gene knockout lines and miRNA or promoter variants enabling functional gene analysis. FIND-IT variants are directly applicable to elite breeding pipelines and minimize time-consuming technical steps to accelerate the evolution of germplasm.

Identification of QTLs for rice grain size and weight by high-throughput SNP markers in the IR64 x Sadri population

Grain appearance is one of the most important attributes of rice. It is determined by grain size, shape, and weight, which in turn influences the rice yield and market value. In this study, QTLs for grain length, grain width, grain length/width ratio, and grain weight were mapped using the high-throughput indica/indica SNP platforms. The population of the mega indica variety IR64 and the high-quality aromatic variety Sadri from Iran was phenotyped. Based on this phenotypic data, plants of 94 F2:3 families including both parents were selected. A linkage map analysis of 210 SNP markers identified 14 QTLs controlling the grain length, grain width, length/width ratio, and 1,000 grain weight. Among these 14, one important region containing the QTLs for all the four studies’ traits was mapped on chromosome 8. It was derived from Sadri for the decreased length/width ratio and increased grain weight. This study demonstrated the speed and efficiency in using multiplex SNP genotyping for QTL analysis. Moreover, this study identified four novel QTLs (qGL8, qTGW8, qLWR8, and qGW8) sharing the same position on chromosome 8 which were linked with grain quality characteristics between one indica and one aromatic variety. It will enable more precise marker-assisted selection for grain weight, shape, and size. Further in-depth studies are required to dissect this region of interest and identify the related gene(s).

GLW7.1, a Strong Functional Allele of Ghd7, Enhances Grain Size in Rice

Grain size is a key determinant of both grain weight and grain quality. Here, we report the map-based cloning of a novel quantitative trait locus (QTL), GLW7.1 (Grain Length, Width and Weight 7.1), which encodes the CCT motif family protein, GHD7. The QTL is located in a 53 kb deletion fragment in the cultivar Jin23B, compared with the cultivar CR071. Scanning electron microscopy analysis and expression analysis revealed that GLW7.1 promotes the transcription of several cell division and expansion genes, further resulting in a larger cell size and increased cell number, and finally enhancing the grain size as well as grain weight. GLW7.1 could also increase endogenous GA content by up-regulating the expression of GA biosynthesis genes. Yeast two-hybrid assays and split firefly luciferase complementation assays revealed the interactions of GHD7 with seven grain-size-related proteins and the rice DELLA protein SLR1. Haplotype analysis and transcription activation assay revealed the effect of six amino acid substitutions on GHD7 activation activity. Additionally, the NIL with GLW7.1 showed reduced chalkiness and improved cooking and eating quality. These findings provide a new insight into the role of Ghd7 and confirm the great potential of the GLW7.1 allele in simultaneously improving grain yield and quality.