Fine mapping of the panicle length QTL qPL5 in rice.

Panicle length is a crucial trait tightly associated with spikelets per panicle and grain yield in rice. To dissect the genetic basis of panicle length, a population of 161 recombinant inbred lines (RILs) was developed from the cross between an aus variety Chuan 7 (C7) and a tropical Geng variety Haoboka (HBK). C7 has a panicle length of 30 cm, 7 cm longer than that of HBK, and the panicle length was normally distributed in the RIL population. A total of six quantitative trait loci (QTLs) for panicle length were identified, and single QTLs explained the phenotypic variance from 4.9 to 18.1%. Among them, three QTLs were mapped to the regions harbored sd1, DLT, and Ehd1, respectively. To validate the genetic effect of a minor QTL qPL5, a near-isogenic F2 (NIF2) population segregated at qPL5 was developed. Interestingly, panicle length displayed bimodal distribution, and heading date also exhibited significant variation in the NIF2 population. qPL5 accounted for 66.5% of the panicle length variance. The C7 allele at qPL5 increased panicle length by 2.4 cm and promoted heading date by 5 days. Finally, qPL5 was narrowed down to an 80-kb region flanked by markers M2197 and M2205 using a large NIF2 population of 7600 plants. LOC_Os05g37540, encoding a phytochrome signal protein whose homolog in Arabidopsis enlarges panicle length, is regarded as the candidate gene because a single-nucleotide mutation (C1099T) caused a premature stop codon in HBK. The characterization of qPL5 with enlarging panicle length but promoting heading date makes its great value in breeding early mature varieties without yield penalty in rice. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-024-01443-2.

A long transcript mutant of the rubisco activase gene RCA upregulated by the transcription factor Ghd2 enhances drought tolerance in rice.

The transcription factor Ghd2 increases rice yield potential under normal conditions and accelerates leaf senescence under drought stress. However, its mechanism on the regulation of leaf senescence under drought stress remains unclear. In the present study, to unveil the mechanism, one target of Ghd2, the Rubisco activase gene RCA, was identified through the combined analysis of Ghd2-CRISPR transcriptome data and Ghd2-overexpression microarray data. Ghd2 binds to the 'CACA' motif in the RCA promoter by its CCT domain and upregulates RCA expression. RCA has alternative transcripts, RCAS and RCAL, which are predominantly expressed under normal conditions and drought stress, respectively. Similar to Ghd2-overexpressing plants, RCAL-overexpressing plants were more sensitive to drought stress than the wild-type. However, the plants overexpressing RCAS showed a weak drought-sensitive phenotype. Moreover, RCAL knockdown and knockout plants did not show yield loss under normal conditions, but exhibited enhanced drought tolerance and delayed leaf senescence. The chlorophyll content, the free amino acid content and the expression of senescence-related genes in the RCAL mutant were lower than those in the wild-type plants under drought stress. In summary, Ghd2 induces leaf senescence by upregulating RCAL expression under drought stress, and the RCAL mutant has important values in breeding drought-tolerant varieties.

A Ghd7-centered regulatory network provides a mechanistic approximation to optimal heterosis in an elite rice hybrid.

Heterosis refers to the superior performance of hybrids over their parents, which is a general phenomenon occurring in diverse organisms. Many commercial hybrids produce high yield without delayed flowering, which we refer to as optimal heterosis and is desired in hybrid breeding. Here, we attempted to illustrate the genomic basis of optimal heterosis by reinvestigating the single-locus quantitative trait loci and digenic interactions of two traits, the number of spikelets per panicle (SP) and heading date (HD), using recombinant inbred lines and 'immortalized F2 s' derived from the elite rice (Oryza sativa) hybrid Shanyou 63. Our analysis revealed a regulatory network that may provide an approximation to the genetic constitution of the optimal heterosis observed in this hybrid. In this network, Ghd7 works as the core element, and three other genes, Ghd7.1, Hd1, and Hd3a/RFT1, also have major roles. The effects of positive dominance by Ghd7 and Ghd7.1 and negative dominance by Hd1 and Hd3a/RFT1 in the hybrid background contribute the major part to the high SP without delaying HD; numerous epistatic interactions, most of which involve Ghd7, also play important roles collectively. The results expand our understanding of the genic interaction networks underlying hybrid rice breeding programs, which may be very useful in future crop genetic improvement.

Structural Insight into DNA Recognition by CCT/NF-YB/YC Complexes in Plant Photoperiodic Flowering.

CONSTANS, CONSTANS-LIKE, and TIMING OF CAB EXPRESSION1 (CCT) domain-containing proteins are a large family unique to plants. They transcriptionally regulate photoperiodic flowering, circadian rhythms, vernalization, and other related processes. Through their CCT domains, CONSTANS and HEADING DATE1 (HD1) coordinate with the NUCLEAR FACTOR Y (NF-Y) B/C dimer to specifically target a conserved 'CCACA' motif within the promoters of their target genes. However, the mechanism underlying DNA recognition by the CCT domain remains unclear. Here we determined the crystal structures of the rice (Oryza sativa) NF-YB/YC dimer and the florigen gene Heading date 3a (Hd3a)-bound HD1CCT/NF-YB/YC trimer with resolutions of 2.0 Å and 2.55 Å, respectively. The CCT domain of HD1 displays an elongated structure containing two α-helices and two loops, tethering Hd3a to the NF-YB/YC dimer. Helix α2 and loop 2 are anchored into the minor groove of the 'CCACA' motif, which determines the specific base recognition. Our structures reveal the interaction mechanism among the CCT domain, NF-YB/YC dimer, and the target DNA. These results not only provide insight into the network between the CCT proteins and NF-Y subunits, but also offer potential approaches for improving productivity and global adaptability of crops by manipulating florigen expression.

A Pseudo-near isogenic F2 population strategy for rapid QTL cloning.

Near isogenic F2 (NIF2) population frequently developed by conventional backcross has dramatically contributed to QTL identification in plants. Developing such a NIF2 population is time-consuming. Thus, it is urgent to rapidly produce a NIF2 population for QTL cloning. Here, we proposed a rapid QTL cloning strategy by generating a Pseudo-near isogenic F2 population (Pseudo-NIF2), which segregates at the target QTL but is fixed at other QTLs for the target trait. Nineteen QTLs for GL, GW, and TGW were detected in the F2 population from the cross between Zhenshan 97 and Egy316. To verify the efficiency of Pseudo-NIF2 in QTL quick cloning, the novel moderate QTL qGL10.1 which explained 9.1% and 5.6% of grain length variation in F2 and F2:3 populations was taken as an example. An F2 plant (F2-120), which segregated at qGL10.1 but fixed at other 8 QTLs for grain length, was screened to generate a Pseudo-NIF2 population by selfing cross. In the Pseudo-NIF2 population, the segregation ratio of plants with long grains to short grains fits 3:1, indicating that one gene controlled the variation of grain length. Based on the Pseudo-NIF2 and its progeny, qGL10.1 was fine mapped to a 19.3-kb region, where a gene OsMADS56 was verified as the candidate by functional polymorphism between parental alleles. Pseudo-NIF2 strategy is a rapid way for QTL cloning, which saves 3 to 4 cropping seasons compared to the conventional way. Applying the method for cloning QTL with moderate or major effects is promising. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01408-x.

The amino acid residue E96 of Ghd8 is crucial for the formation of the flowering repression complex Ghd7-Ghd8-OsHAP5C in rice.

Ghd7 is an important gene involved in the photoperiod flowering pathway in rice. A Ghd7-involved transcriptional regulatory network has been established, but its translational regulatory pathway is poorly understood. The mutant suppressor of overexpression of Ghd7 (sog7) was identified from EMS-induced mutagenesis on the background of ZH11 overexpressing Ghd7. MutMap analysis revealed that SOG7 is allelic to Ghd8 and delayed flowering under long-day (LD) conditions. Biochemical assays showed that Ghd8 interacts with OsHAP5C and Ghd7 both in vivo and in vitro. Surprisingly, a point mutation E96K in the α2 helix of the Ghd8 histone fold domain (HFD) destroyed its ability to interact with Ghd7. The prediction of the structure shows that mutated amino acid is located in the interaction region of CCT/NF-YB/YC complexes, which alter the structure of α4 of Ghd8. This structural difference prevents the formation of complex NF-YB/YC. The triple complex of Ghd8-OsHAP5C-Ghd7 directly bound to the promotor of Hd3a and downregulated the expression of Ehd1, Hd3a and RFT1, and finally resulted in a delayed heading. These findings are helpful in deeply understanding the Ghd7-involved photoperiod flowering pathway and promote the elucidation of rice heading.

Five plants per RIL for phenotyping traits of high or moderate heritability ensure the power of QTL mapping in a rice MAGIC population.

Currently, the power of QTL mapping is mainly dependent on the quality of phenotypic data in a given population, regardless of the statistical method, as the quality of genotypic data is easily guaranteed in the laboratory. Increasing the sample size per line used for phenotyping is a good way to improve the quality of phenotypic data. However, accommodating a large-scale mapping population takes a large area of rice field, which frequently results in high costs and extra environmental noises. To acquire a reasonable small sample size without a penalty in mapping power, we conducted three experiments with a 4-way MAGIC population and measured phenotypes of 5, 10, and 20 plants per RIL. Three traits including heading date, plant height, and tillers per plant were focused. With SNP- and bin-based QTL mapping, 3 major and 3 minor QTLs for heading date with high heritability and 2 major QTLs for plant height with moderate heritability were commonly detected across the three experiments, but no QTL for tillers per plant with low heritability were commonly identified. In addition, bin-based QTL mapping was more powerful than SNP-based mapping and able to rank the genetic effects of parental alleles. Thus, 5 plants per RIL for phenotyping ensure the power of QTL mapping for traits of high or moderate heritability, and bin-based QTL mapping is recommended for multiparent populations.

Natural variation in the HAN1 gene confers chilling tolerance in rice and allowed adaptation to a temperate climate.

Rice (Oryza sativa L.) is a chilling-sensitive staple crop that originated in subtropical regions of Asia. Introduction of the chilling tolerance trait enables the expansion of rice cultivation to temperate regions. Here we report the cloning and characterization of HAN1, a quantitative trait locus (QTL) that confers chilling tolerance on temperate japonica rice. HAN1 encodes an oxidase that catalyzes the conversion of biologically active jasmonoyl-L-isoleucine (JA-Ile) to the inactive form 12-hydroxy-JA-Ile (12OH-JA-Ile) and fine-tunes the JA-mediated chilling response. Natural variants in HAN1 diverged between indica and japonica rice during domestication. A specific allele from temperate japonica rice, which gained a putative MYB cis-element in the promoter of HAN1 during the divergence of the two japonica ecotypes, enhances the chilling tolerance of temperate japonica rice and allows it to adapt to a temperate climate. The results of this study extend our understanding of the northward expansion of rice cultivation and provide a target gene for the improvement of chilling tolerance in rice.

Assembly of yield heterosis of an elite rice hybrid is promising by manipulating dominant quantitative trait loci.

In the past, rice hybrids with strong heterosis have been obtained empirically, by developing and testing thousands of combinations. Here, we aimed to determine whether heterosis of an elite hybrid could be achieved by manipulating major quantitative trait loci. We used 202 chromosome segment substitution lines from the elite hybrid Shanyou 63 to evaluate single segment heterosis (SSH) of yield per plant and identify heterotic loci. All nine detected heterotic loci acted in a dominant fashion, and no SSH exhibited overdominance. Functional alleles of key yield-related genes Ghd7, Ghd7.1, Hd1, and GS3 were dispersed in both parents. No functional alleles of three investigated genes were expressed at higher levels in the hybrids than in the more desirable parents. A hybrid pyramiding eight heterotic loci in the female parent Zhenshan 97 background had a comparable yield to Shanyou 63 and much higher yield than Zhenshan 97. Five hybrids pyramiding eight or nine heterotic loci in the combined parental genome background showed similar yield performance to that of Shanyou 63. These results suggest that dominance underlying functional complementation is an important contributor to yield heterosis and that heterosis assembly might be successfully promised by manipulating several major dominant heterotic loci.

Wide Grain 7 increases grain width by enhancing H3K4me3 enrichment in the OsMADS1 promoter in rice (Oryza sativa L.).

Grain size is a major determinant of grain weight, a key component of grain yield of rice. Here, we identified the grain size gene WIDE GRAIN 7 (WG7) from a T-DNA insertion mutant. The grain size of WG7 knockout mutants and WG7 overexpression lines indicated that WG7 is a positive regulator of grain size. WG7 encodes a cysteine-tryptophan (CW) domain-containing transcriptional activator. EMSAs and ChIP-qPCR assay confirmed that WG7 directly bound to the promoter of OsMADS1, a grain size gene, and thereby significantly activated its expression. Point mutations showed that the cis-element CATTTC motif in the promoter was the binding site of WG7. Compared with the wild-type, deletion mutants of the cis-element motif exhibited lower expression of OsMADS1 and produced narrower grains, implicating the requirement of this motif for WG7 function. ChIP-qPCR assays showed that WG7 enhanced histone H3K4me3 enrichment in the promoter of OsMADS1. WG7 underwent directional selection due to the poor fertility of the non-functional mutant. These findings demonstrated that WG7 upregulated OsMADS1 expression by directly binding to its promoter, enhanced histone H3K4me3 enrichment in the promoter and ultimately increased grain width. This study will enrich the knowledge concerning the regulatory network of grain size formation in rice.

An ethyl methanesulfonate-induced neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice.

Spontaneous mutants are mainly obtained from tissue culture or natural occurrences in plants. The traditional strategy for identifying spontaneously mutated genes is to continuously backcross these mutants to another variety and develop a near-isogenic F2 population for map-based cloning or bulked segregant analysis. However, this strategy is time-consuming. Here, we have developed a new method to efficiently accelerate the identification process. The chemical mutagen ethyl methanesulfonate was first used to treat the wild type of the spontaneous mutants to induce thousands of neutral mutations. An induced individual without any statistically significant phenotypic changes which was compared with the wild type was chosen as the neutral mutant. The spontaneous mutant was then crossed with the neutral mutant to develop a pseudo-near-isogenic F2 population in which only the induced neutral mutations and the causal mutation were segregated in the genome. This population ensures that the variation of the mutated trait is controlled only by the spontaneously mutated gene. Finally, after sequencing the neutral mutant and the mutant-type DNA pool of the F2 population the spontaneous mutation will be identified quickly by bioinformatics analysis. Using this method, two spontaneously mutated genes were identified successfully. Therefore, the neutral mutant-bridging method efficiently identifies spontaneously mutated genes in rice, and its value in other plants is discussed.

OsMFT2 is involved in the regulation of ABA signaling-mediated seed germination through interacting with OsbZIP23/66/72 in rice.

Seed germination is a complex process involving various physical and biochemical cues, determined by exogenous and endogenous factors. Here, we identified a gene, OsMFT2, that negatively regulates seed germination in rice. OsMFT2 knock-out lines exhibited pre-harvest sprouting, whereas OsMFT2 overexpression lines showed delayed germination. RNA expression profiling showed that OsMFT2 was specifically expressed in seeds. Subcellular localization indicated that OsMFT2 was a nuclear protein. Exogenous abscisic acid (ABA) treatment of imbibed seeds and seedlings indicated that OsMFT2 altered ABA sensitivity during seed germination and post-germination growth. In vivo and in vitro assays showed that three bZIP transcription factors, OsbZIP23, OsbZIP66 and OsbZIP72, interacted with OsMFT2. OsbZIP23/66/72 bound to the promoter of Rab16A, a typical gene containing the ABA-responsive element, and OsMFT2 enhanced the binding to the Rab16A promoter. Moreover, several ABA-responsive genes were differentially expressed in the imbibed seeds of OsMFT2 transgenic lines and the wild type. The performance of the transgenic plants demonstrated that overexpressing OsbZIP23 rescued the pre-harvest sprouting phenotype and the decrease in ABA-signaling genes expression caused by OsMFT2 knock-out. All of these results demonstrate that OsMFT2 positively regulates ABA-responsive genes through interacting with OsbZIP23/66/72 and functions in seed germination.

Genetic analyses of lodging resistance and yield provide insights into post-Green-Revolution breeding in rice.

Lodging reduces grain yield in cereal crops. Understanding the genetic basis of lodging resistance (LR) benefits LR breeding. In the study, 524 accessions from a rice germplasm collection and 193 recombinant inbred lines were phenotyped for 17 LR-related traits. Height and culm strength (the magnitude of applied force necessary to break the culm) were two major factors affecting LR. We conducted genome-wide association study (GWAS) and identified 127 LR-associated loci. Significant phenotypic correlations between culm-strength traits and yield-related traits were observed. To reveal the genetic relationship between them, we conducted GWAS of culm-strength traits with adding yield-related trait as a covariate and detected 63 loci linking culm strength and yield. As a proof, a near-isogenic line for an association locus on chromosome 7 showed enhanced LR and yield. Strikingly, 58 additional loci were identified in the covariate-added GWAS. Several LR-associated loci had undergone divergent selection. Linkage analysis supported the GWAS results. We propose that introgression of alleles beneficial for both culm strength and panicle weight without negative effects on panicle number or pyramiding high-yielding alleles and lodging-resistant alleles without effects on yield can be employed for the post-Green-Revolution breeding.

The chromosome-scale reference genome of safflower (Carthamus tinctorius) provides insights into linoleic acid and flavonoid biosynthesis.

Safflower (Carthamus tinctorius L.), a member of the Asteraceae, is a popular crop due to its high linoleic acid (LA) and flavonoid (such as hydroxysafflor yellow A) contents. Here, we report the first high-quality genome assembly (contig N50 of 21.23 Mb) for the 12 pseudochromosomes of safflower using single-molecule real-time sequencing, Hi-C mapping technologies and a genetic linkage map. Phyloge nomic analysis showed that safflower diverged from artichoke (Cynara cardunculus) and sunflower (Helianthus annuus) approximately 30.7 and 60.5 million years ago, respectively. Comparative genomic analyses revealed that uniquely expanded gene families in safflower were enriched for those predicted to be involved in lipid metabolism and transport and abscisic acid signalling. Notably, the fatty acid desaturase 2 (FAD2) and chalcone synthase (CHS) families, which function in the LA and flavonoid biosynthesis pathways, respectively, were expanded via tandem duplications in safflower. CarFAD2-12 was specifically expressed in seeds and was vital for high-LA content in seeds, while tandemly duplicated CarFAD2 genes were up-regulated in ovaries compared to CarFAD2-12, which indicates regulatory divergence of FAD2 in seeds and ovaries. CarCHS1, CarCHS4 and tandem-duplicated CarCHS5˜CarCHS6, which were up-regulated compared to other CarCHS members at early stages, contribute to the accumulation of major flavonoids in flowers. In addition, our data reveal multiple alternative splicing events in gene families related to fatty acid and flavonoid biosynthesis. Together, these results provide a high-quality reference genome and evolutionary insights into the molecular basis of fatty acid and flavonoid biosynthesis in safflower.

OsHOX1 and OsHOX28 Redundantly Shape Rice Tiller Angle by Reducing HSFA2D Expression and Auxin Content.

Tiller angle largely determines plant architecture, which in turn substantially influences crop production by affecting planting density. A recent study revealed that HEAT STRESS TRANSCRIPTION FACTOR2D (HSFA2D) acts upstream of LAZY1 (LA1) to regulate tiller angle establishment in rice (Oryza sativa). However, the mechanisms underlying transcriptional regulation of HSFA2D remain unknown. In this study, two class II homeodomain-Leu zipper genes, OsHOX1 and OsHOX28, were identified as positive regulators of tiller angle by affecting shoot gravitropism. OsHOX1 and OsHOX28 showed strong transcriptional suppressive activity in rice protoplasts and formed intricate self- and mutual-transcriptional negative feedback loops. Moreover, OsHOX1 and OsHOX28 bound to the pseudopalindromic sequence CAAT(C/G)ATTG within the promoter of HSFA2D, thus suppressing its expression. In contrast to HSFA2D and LA1, OsHOX1 and OsHOX28 attenuated lateral auxin transport, thus repressing the expression of WUSCHEL-RELATED HOMEOBOX 6 (WOX6) and WOX11 in the lower side of the shoot base of plants subjected to gravistimulation. Genetic analysis further confirmed that OsHOX1 and OsHOX28 act upstream of HSFA2D Additionally, both OsHOX1 and OsHOX28 inhibit the expression of multiple OsYUCCA genes and decrease auxin biosynthesis. Taken together, these results demonstrated that OsHOX1 and OsHOX28 regulate the local distribution of auxin, and thus tiller angle establishment, through suppression of the HSFA2D-LA1 pathway and reduction of endogenous auxin content. Our finding increases the knowledge concerning fine tuning of tiller angles to optimize plant architecture in rice.

The transcriptional repressor OsPRR73 links circadian clock and photoperiod pathway to control heading date in rice.

The phase transition from vegetative to reproductive growth is triggered by internal and external signals that participate in circadian clock in plants. We identified a rice floral inhibitor OsPRR73 encoding a CONSTANS protein. Overexpression of OsPRR73 resulted in late heading under both long-day (LD) and short-day (SD) conditions. Knockout mutants led to early heading under LD conditions but no change under SD. OsPRR73 mRNA accumulated at noon and exhibited a robust oscillation under constant light (LL) and constant darkness (DD) conditions. OsPRR73 overexpression exerted negative feedback on endogenous OsPRR73 expression and altered diurnal expressions of key flowering genes and circadian clock genes. OsPRR73 bound to the promoters of the floral gene Ehd1 and the circadian gene OsLHY, and significantly suppressed their expression at dawn. In LL and DD, the oscillatory patterns of the circadian genes OsLHY, OsTOC1, OsGI and OsELF3 were varied in OsPRR73OX and osprr73 mutants. OsPRR73 expression was decreased in osphyb mutants, and overexpression of OsPRR73 complemented the early heading date phenotype of osphyb, indicating OsPRR73 works downstream of OsPhyB. Therefore, OsPRR73 is involved in a feedback loop of the rice clock and connects the photoperiod flowering pathway by binding to the Ehd1 promoter in rice.

Combinations of Ghd7, Ghd8, and Hd1 determine strong heterosis of commercial rice hybrids in diverse ecological regions.

Heterosis of grain yield is closely associated with heading date in crops. Gene combinations of the major heading date genes Ghd7, Ghd8, and Hd1 play important roles in enhancing grain yield and adaptation to ecological regions in rice. However, the predominant three-gene combinations for a specific ecological region remain unclear in both three-line and two-line hybrids. In this study, we sequenced these three genes of 50 cytoplasmic male sterile/maintainer lines, 31 photo-thermo-sensitive genic male sterile lines, and 109 restorer lines. Sequence analysis showed that hybrids carrying strong functional alleles of Ghd7 and Hd1 and non-functional Ghd8 are predominant in three-line hybrids and are recommended for rice production in the subtropics around 30°N/S. Hybrids carrying strong functional Ghd7 and Ghd8 and non-functional Hd1 are predominant in two-line hybrids and are recommended for low latitude areas around 23.5°N/S rich in photothermal resources. Hybrids carrying strong functional Ghd7 and Ghd8 and functional Hd1 were not identified in commercial hybrids in the middle and lower reaches of the Yangtze River, but they have high yield potential in tropical regions because they have the strongest photoperiod sensitivity. Based on these findings, two genic sterile lines, Xiangling 628S and C815S, whose hybrids often head very late, were diagnosed with these three genes, and Hd1 was targeted to be knocked out in Xiangling 628S and replaced with hd1 in C815S. The hybrids developed from both modified sterile lines in turn had appropriate heading dates and significantly improved grain yield. This study provides new insights for breeding design to develop hybrids for various regions.

Genome-Wide Association Studies Reveal the Genetic Basis of Ionomic Variation in Rice.

Rice (Oryza sativa) is an important dietary source of both essential micronutrients and toxic trace elements for humans. The genetic basis underlying the variations in the mineral composition, the ionome, in rice remains largely unknown. Here, we describe a comprehensive study of the genetic architecture of the variation in the rice ionome performed using genome-wide association studies (GWAS) of the concentrations of 17 mineral elements in rice grain from a diverse panel of 529 accessions, each genotyped at ∼6.4 million single nucleotide polymorphism loci. We identified 72 loci associated with natural ionomic variations, 32 that are common across locations and 40 that are common within a single location. We identified candidate genes for 42 loci and provide evidence for the causal nature of three genes, the sodium transporter gene Os-HKT1;5 for sodium, Os-MOLYBDATE TRANSPORTER1;1 for molybdenum, and Grain number, plant height, and heading date7 for nitrogen. Comparison of GWAS data from rice versus Arabidopsis (Arabidopsis thaliana) also identified well-known as well as new candidates with potential for further characterization. Our study provides crucial insights into the genetic basis of ionomic variations in rice and serves as an important foundation for further studies on the genetic and molecular mechanisms controlling the rice ionome.

Genome-wide association studies reveal that members of bHLH subfamily 16 share a conserved function in regulating flag leaf angle in rice (Oryza sativa).

As a major component of ideal plant architecture, leaf angle especially flag leaf angle (FLA) makes a large contribution to grain yield in rice. We utilized a worldwide germplasm collection to elucidate the genetic basis of FLA that would be helpful for molecular design breeding in rice. Genome-wide association studies (GWAS) identified a total of 40 and 32 QTLs for FLA in Wuhan and Hainan, respectively. Eight QTLs were commonly detected in both conditions. Of these, 2 and 3 QTLs were identified in the indica and japonica subpopulations, respectively. In addition, the candidates of 5 FLA QTLs were verified by haplotype-level association analysis. These results indicate diverse genetic bases for FLA between the indica and japonica subpopulations. Three candidates, OsbHLH153, OsbHLH173 and OsbHLH174, quickly responded to BR and IAA involved in plant architecture except for OsbHLH173, whose expression level was too low to be detected; their overexpression in plants increased rice leaf angle. Together with previous studies, it was concluded that all 6 members in bHLH subfamily 16 had the conserved function in regulating FLA in rice. A comparison with our previous GWAS for tiller angle (TA) showed only one QTL had pleiotropic effects on FLA and TA, which explained low similarity of the genetic basis between FLA and TA. An ideal plant architecture is expected to be efficiently developed by combining favorable alleles for FLA from indica with favorable alleles for TA from japonica by inter-subspecies hybridization.

Root-to-Shoot Long-Distance Mobile miRNAs Identified from Nicotiana Rootstocks.

Root-derived mobile signals play critical roles in coordinating a shoot's response to underground conditions. However, the identification of root-to-shoot long-distance mobile signals has been scant. In this study, we aimed to characterize root-to-shoot endogenous mobile miRNAs by using an Arabidopsis/Nicotiana interfamilial heterograft in which these two taxonomically distant species with clear genetic backgrounds had sufficient diversity in differentiating miRNA sources. Small RNA deep sequencing analysis revealed that 82 miRNAs from the Arabidopsis scion could travel through the graft union to reach the rootstock, whereas only a very small subset of miRNA (6 miRNAs) preferred the root-to-shoot movement. We demonstrated in an ex vivo RNA imaging experiment that the root-to-shoot mobile Nb-miR164, Nb-miR395 and Nb-miR397 were targeted to plasmodesmata using the bacteriophage coat protein MS2 system. Furthermore, the Nb-miR164 was shown to move from the roots to the shoots to induce phenotypic changes when its overexpressing line was used as rootstock, strongly supporting that root-derived Nb-miR164 was able to modify the scion trait via its long-distance movement.