Bulk segregant analysis coupled with transcriptomics and metabolomics revealed key regulators of bacterial leaf blight resistance in rice.

BACKGROUND: Bacterial leaf blight (BLB) is a highly destructive disease, causing significant yield losses in rice (Oryza sativa). Genetic variation is contemplated as the most effective measure for inducing resistance in plants. The mutant line T1247 derived from R3550 (BLB susceptible) was highly resistant to BLB. Therefore, by utilizing this valuable source, we employed bulk segregant analysis (BSA) and transcriptome profiling to identify the genetic basis of BLB resistance in T1247. RESULTS: The differential subtraction method in BSA identified a quantitative trait locus (QTL) on chromosome 11 spanning a 27-27.45 Mb region with 33 genes and 4 differentially expressed genes (DEGs). Four DEGs (P < 0.01) with three putative candidate genes, OsR498G1120557200, OsR498G1120555700, and OsR498G1120563600,0.01 in the QTL region were identified with specific regulation as a response to BLB inoculation. Moreover, transcriptome profiling identified 37 resistance analogs genes displaying differential regulation. CONCLUSIONS: Our study provides a substantial addition to the available information regarding QTLs associated with BLB, and further functional verification of identified candidate genes can broaden the scope of understanding the BLB resistance mechanism in rice.

Transcriptomic and methylomic analyses provide insights into the molecular mechanism and prediction of heterosis in rice.

Heterosis has been widely used in multiple crops. However, the molecular mechanism and prediction of heterosis remains elusive. We generated five F1 hybrids [four showing better-parent heterosis (BPH) and one showing mid-parent heterosis], and performed the transcriptomic and methylomic analyses to identify the candidate genes for BPH and explore the molecular mechanism of heterosis and the potential predictors for heterosis. Transcriptomic results showed that most of the differentially expressed genes shared in the four better-parent hybrids were significantly enriched into the terms of molecular function, and the additive and dominant effects played crucial roles for BPH. DNA methylation level, especially in CG context, significantly and positively correlated with grain yield per plant. The ratios of differentially methylated regions in CG context in exons to transcription start sites between the parents exhibited significantly negative correlation with the heterosis levels of their hybrids, as was further confirmed in 24 pairwise comparisons of other rice lines, implying that this ratio could be a feasible predictor for heterosis level, and this ratio of less than 5 between parents in early growth stages might be a critical index for judging that their F1 hybrids would show BPH. Additionally, we identified some important genes showing differential expression and methylation, such as OsDCL2, Pi5, DTH2, DTH8, Hd1 and GLW7 in the four better-parent hybrids as the candidate genes for BPH. Our findings helped shed more light on the molecular mechanism and heterosis prediction.

Chromosomal-level genome assembly of the high-quality Xian/Indica rice (Oryza sativa L.) Xiangyaxiangzhan.

The indica rice variety XYXZ carries elite traits including appearance and eating quality. Here, we report the de novo assembly of XYXZ using Illumine paired-end whole-genome shotgun sequencing and Nanopore sequencing. We annotated 39,722 protein-coding genes in the 395.04 Mb assembly. In comparison to other cultivars, XYXZ showed a larger gene size including the transcripts and introns, and more exons per gene. And hundreds of ultra-long genes were also detected. A total of 4362 complete LTRs were annotated, and among them, many were located next to or in protein-coding genes including several genes related to rice quality. We observed the different distributions of LTRs in these genes among XYXZ, Nipponbare, and R498, implying these LTRs might potentially affect expressions of the proximal genes and rice quality. Overall, This chromosome-length genome assembly of XYXZ provides a valuable resource for gene discovery, genetic variation and evolution, and the breeding of high-quality rice.

The grain yield regulator NOG1 plays a dual role in latitudinal adaptation and cold tolerance during rice domestication.

Rice originated in tropical and subtropical regions and is distributed worldwide. Low temperature is one of the most critical abiotic stresses affecting grain yield and geographical distribution of rice. It is vital to elucidate the molecular mechanism of chilling tolerance in rice for ensuring cereals production. Previously we isolated the domestication-related gene NOG1 which affects rice grain number and yield. In this study, we specified that rice varieties harboring high-yielding NOG1 allele are more distributed in low-latitude regions. Additionally, we observed NOG1 influences the chilling tolerance of rice. Through genome-wide transcriptional analysis after cold treatment at 10°C, there were 717 differentially expressed genes (DEGs) in nog1 near-isogenic lines compared with the control Guichao 2, including 432 up-regulated DEGs and 284 down-regulated DEGs. Gene ontology annotations and KEGG enrichment analysis of DEGs showed that various biological processes and signaling pathways were related to cold stress, such as lipid metabolism and genetic information processing. These results provide new insights into the mechanism of chilling tolerance in rice and the molecular basis of environmental adaptation during rice domestication.

Identification of Heterotic Groups and Patterns Based on Genotypic and Phenotypic Characteristics Among Rice Accessions of Diverse Origins.

Identification of the right parental combinations to maximize heterosis is the major goal of hybrid breeding, which could be achieved through identification of heterotic groups. The main objective of this study was to identify promising heterotic groups for future rice breeding programs. A collection of 359 rice genotypes of diverse origins of China and abroad, composed of inbreds, maintainers, restorers, and temperature-sensitive genic male sterile (TGMS) lines were genotyped using 10K SNP chips. The SNP data set was subjected to genomic analyses for estimation of genetic divergence and diversity. Significant variations were observed in the germplasm with the identification of six different genetic groups. These lines were assigned to the genetic groups independent of their origin. Taking an account of commercially used heterotic groups present in each cluster, three cytoplasmic male sterile (CMS) lines and 14 inbred and restorer lines with moderate to high genetic distances selected from five heterotic patterns were crossed and obtained 42 F1 hybrids. A total of 14 hybrids were found with significant maximum mid- and better-parent heterosis, namely, TaifengA × Guang122, TaifengA × Wushansimiao, and TaifengA × Minghui63 for earliness; Guang8A × Huazhan for dwarf stature; and Guang8A × Huanghuzhan-1, TaifengA × Yuexiangzhan, Guang8A × Minhui3301, TianfengA × Guang122, Guang8A × Yahui2115, TianfengA × Huanghuazhan, TianfengA × Minghui63, TianfengA × Minhui3301, TaifengA × Gui99, and Guang8A × Yuenongsimiao for yield and yield-related traits. Mid-parent and better-parent heterotic F1 hybrids were in positive correlation with the genetic distances as that manifested by commercially used heterotic groups, encouraging the use of genotypic data for identification of heterotic groups. Our study provides an informative strategy for the development of early maturing, lodging resistant and high-yielding commercial hybrids and cultivars in future heterosis breeding programs.

Rice Brittle Culm19 Encoding Cellulose Synthase Subunit CESA4 Causes Dominant Brittle Phenotype But has No Distinct Influence on Growth and Grain Yield.

BACKGROUND: Mechanical strength is a crucial agronomic trait in rice (Oryza sativa), and brittle mutants are thought suitable materials to investigate the mechanism of cell wall formation. So far, almost all brittle mutants are recessive, and most of them are defected in multiple morphologies and/or grain yield, limiting their application in hybrid breeding and in rice straw recycling. RESULTS: We identified a semi-dominant brittle mutant Brittle culm19 (Bc19) isolated from the japonica variety Nipponbare through chemical mutagenesis. The mutant showed the same apparent morphologies and grain yield to the wild type plant except for its weak mechanical strength. Its development of secondary cell wall in sclerenchyma cells was affected, along with reduced contents of cellulose, hemicellulose, lignin and sugars in culms and leaves. Positional cloning suggested that the Bc19 gene was allelic to OsCESA4, encoding one of the cellulose synthase A (CESA) catalytic subunits. In this mutant, a C-to-T substitution occurred in the coding sequence of BC19, causing the P507S missense mutation in its encoded product, which was located in the second cytoplasmic region of the OsCESA4 protein. Furthermore, introducing mutant gene Bc19 into the wild-type plant resulted in brittle plants, confirming that the P507S point mutation in OsCESA4 protein was responsible for the semi-dominant brittle phenotype of Bc19 mutant. Reverse correlation was revealed between cellulose contents and expression levels of mutant gene Bc19 among the homozygous mutant, the hybrid F1 plant, and the Bc19 overexpression transgenic plants, implying that gene Bc19 might affect cellulose synthesis in a dosage-dependent manner. CONCLUSIONS: Bc19, a semi-dominant brittle mutant allele of gene OsCESA4, was identified using map-based cloning approach. The mutated protein of Bc19 possessing the P507S missense mutation behaved in a dosage-dependent semi-dominant manner. Unique brittle effect on phenotype and semi-dominant genetic quality of gene Bc19 indicated its potential application in grain-straw dual-purpose hybrid rice breeding.

Genomic and transcriptomic analysis reveal molecular basis of salinity tolerance in a novel strong salt-tolerant rice landrace Changmaogu.

BACKGROUND: Salt stress is an important factor that limits rice yield. We identified a novel, strongly salt tolerant rice landrace called Changmaogu (CMG) collected from a coastal beach of Zhanjiang, Guangdong Province, China. The salt tolerance of CMG was much better than that of the international recognized salt tolerant rice cultivar Pokkali in the germination and seedling stages. RESULTS: To understand the molecular basis of salt tolerance in CMG, we performed BSA-seq for two extreme bulks derived from the cross between CMG and a cultivar sensitive to salt, Zhefu802. Transcriptomic sequencing was conducted for CMG at the germination and young seedling stages. Six candidate regions for salt tolerance were mapped on Chromosome 1 by BSA-seq using the extreme populations. Based on the polymorphisms identified between both parents, we detected 32 genes containing nonsynonymous coding single nucleotide polymorphisms (SNPs) and frameshift mutations in the open reading frame (ORF) regions. With transcriptomic sequencing, we detected a large number of differentially expressed genes (DEGs) at the germination and seedling stages under salt stress. KEGG analysis indicated two of 69 DEGs shared at the germination and seedling stages were significantly enriched in the pathway of carotenoid biosynthesis. Of the 169 overlapping DEGs among three sample points at the seedling stage, 13 and six DEGs were clustered into the pathways of ABA signal transduction and carotenoid biosynthesis, respectively. Of the 32 genes carrying sequence variation, only OsPP2C8 (Os01g0656200) was differentially expressed in the young seedling stage under salt stress and also showed sequence polymorphism in the ORFs between CMG and Zhefu802. CONCLUSION: OsPP2C8 was identified as the target candidate gene for salinity tolerance in the seedling stage. This provides an important genetic resource for the breeding of novel salt tolerant rice cultivars.

Comparison of Translocation and Transformation from Soil to Rice and Metabolism in Rats for Four Arsenic Species.

Arsenic (As) is ubiquitously present in the environment. The toxicity of As is related to its forms. This study was designed to compare the translocation and transformation of four As species from soil to rice, and metabolism in rats for four arsenic species. A set of 26550 data was obtained from pot experiments of rice plants grown in soil fortified with four As species, and 4050 data were obtained from rat experiments in which 81 rats were administered with the four As species. The total As in grain from the methyl arsenate fortified soil was 6.1, 4.9, and 5.2 times that from As(III), As(V), and dimethyl arsenate fortified soil, respectively. The total As in husk was 1.2-7.8 times greater than that in grain. After oral administration of each As species to rats, 83-96% was accumulatively excreted via feces and urine, while 0.1-16% was detected in blood. The translocation, transformation, and metabolism of different forms of arsenic vary greatly.

QTL mapping using an ultra-high-density SNP map reveals a major locus for grain yield in an elite rice restorer R998.

To dissect the genetic basis of yield formation in restorer line of hybrid rice, we conducted QTL analysis for 6 yield traits including panicles per plant (PPP), grains per panicle (GPP), grain yield per plant (GY), thousand-grain weight (TGW), above-ground biomass (AGB), and harvest index (HI) using SNP markers in a recombinant inbred lines (RILs) population derived from a cross between a tropical japonica inbred Francis and an elite indica restorer Guanghui 998 (R998). A total of 26 QTLs were detected using a high density genetic map consisting of 3016 bin markers. Nineteen out of the 26 QTL alleles from R998 had a beneficial effect on yield traits. Most of the QTLs were co-located with previously reported rice QTLs. qAGB6 and qHI9, controlling AGB and HI respectively, were detected as novel QTLs. Four QTLs for GY were repeatedly detected across two years, with all the beneficial alleles from R998. Notably, qGY8 explained over 20% of the yield variance in both years. Moreover, qGY8 together with qTGW8 and qHI8 formed a QTL cluster. Markers tightly linked with qGY8 were developed. Cloning of qGY8 will facilitate its further exploitation in high-yield breeding.

Transcriptomic Analysis Reveals New Insights into High-Temperature-Dependent Glume-Unclosing in an Elite Rice Male Sterile Line.

Glume-unclosing after anthesis is a widespread phenomenon in hybrid rice and also a maternal hereditary trait. The character of Glume-unclosing in rice male sterile lines also seriously influences germination rate and the commercial quality of hybrid rice seeds. We validated that the type of glume-unclosing after anthesis in the elite rice thermo-sensitive genic male sterile (TGMS) line RGD-7S was caused by high temperature. Transcriptomic sequencing of rice panicles was performed to explore the change of transcript profiles under four conditions: pre- and post-anthesis under high temperature (HRGD0 and HRGD1), and pre- and post-anthesis under low temperature (LRGD0 and LRGD1). We identified a total of 14,540 differentially expressed genes (DEGs) including some heat shock factors (HSFs) across the four samples. We found that more genes were up-regulated than down-regulated in the sample pair HRGD1vsHRGD0. These up-regulated genes were significantly enriched in the three biological processes of carbohydrate metabolism, response to water and cell wall macromolecular metabolism. Simultaneously, we also found that the HSF gene OsHsfB1 was specially up-regulated in HRGD1vsHRGD0. However, the down-regulated DEGs in LRGD1vsLRGD0 were remarkably clustered in the biological process of carbohydrate metabolism. This suggests that carbohydrate metabolism may play a key role in regulation of glume-unclosing under high temperature in RGD-7S. We also analyzed the expression pattern of genes enriched in carbohydrate metabolism and several HSF genes under different conditions and provide new insights into the cause of rice glume-unclosing.

Genome-wide DNA polymorphism in the indica rice varieties RGD-7S and Taifeng B as revealed by whole genome re-sequencing.

Next-generation sequencing technologies provide opportunities to further understand genetic variation, even within closely related cultivars. We performed whole genome resequencing of two elite indica rice varieties, RGD-7S and Taifeng B, whose F1 progeny showed hybrid weakness and hybrid vigor when grown in the early- and late-cropping seasons, respectively. Approximately 150 million 100-bp pair-end reads were generated, which covered ∼86% of the rice (Oryza sativa L. japonica 'Nipponbare') reference genome. A total of 2,758,740 polymorphic sites including 2,408,845 SNPs and 349,895 InDels were detected in RGD-7S and Taifeng B, respectively. Applying stringent parameters, we identified 961,791 SNPs and 46,640 InDels between RGD-7S and Taifeng B (RGD-7S/Taifeng B). The density of DNA polymorphisms was 256.8 SNPs and 12.5 InDels per 100 kb for RGD-7S/Taifeng B. Copy number variations (CNVs) were also investigated. In RGD-7S, 1989 of 2727 CNVs were overlapped in 218 genes, and 1231 of 2010 CNVs were annotated in 175 genes in Taifeng B. In addition, we verified a subset of InDels in the interval of hybrid weakness genes, Hw3 and Hw4, and obtained some polymorphic InDel markers, which will provide a sound foundation for cloning hybrid weakness genes. Analysis of genomic variations will also contribute to understanding the genetic basis of hybrid weakness and heterosis.

Expression of MyHC genes, composition of muscle fiber type and their association with intramuscular fat, tenderness in skeletal muscle of Simmental hybrids.

In adult bovine skeletal muscle, it expressed four isoforms of Myosin heavy chain (MyHC) gene, MyHC-I, MyHC-IIa, MyHC-IIb, and MyHC-IIx that are translated into different structural protein myofibrils, and then further form different types of muscle fiber. In the studies, our objective is to reveal the expression patterns of MyHC genes in longissimus dorsi (Ld), semitendinosus (Se) and soleus (Sol) of Simmental hybrids cattle, and their association with intramuscular fat (IMF) content and meat shearing force (MSF). The muscle tissue of Ld, Se and Sol were collected from 6, 12 and 36-month old Simmental hybrids respectively, then the expression levels of MyHCs were examined by real-time PCR, at the same time, IMF, MSF and muscle type were measured with chemical assessment, shearing force measurer and immunostaining respectively. Our results showed that t Ld, Se, and Sol expressed MyHC-I, MyHC-IIa and MyHC-IIx isoforms but not MyHC-IIb, furthermore MyHC-I, MyHC-IIa and MyHC-IIx had different expression patterns in different skeletal muscle. The expression of MyHC-I in Se and Sol, MyHC-IIa in Ld, Se, and Sol, and MyHC-IIx in Sol was decreased with increasing age. The highest expression of MyHC-I in Ld, and MyHC-IIx in Ld and Se was observed in 12-month-old animals. The percentage of type-IIa fiber approximately occupied 70-80 % among various muscle fiber of Ld, Se and Sol. The percentage of different type fiber was not related to IMF content and MSF, but the expression levels of MyHC-I and MyHC-IIa were negatively related to IMF content (r = -0.724, and -0.681, respectively) and MSF (r = -0.672, and -0.641, respectively). The expression level of MyHC-IIx was also negatively related to MSF (r = -0.655). In conclusion, MyHC gene might be considered as a negative factor in genetic selection of IMF content and MSF.

Genetic and cytological analysis of a novel type of low temperature-dependent intrasubspecific hybrid weakness in rice.

Hybrid weakness (HW) is an important postzygotic isolation which occurs in both intra- and inter-specific crosses. In this study, we described a novel low temperature-dependent intrasubspecific hybrid weakness in the F1 plants derived from the cross between two indica rice varieties Taifeng A and V1134. HW plants showed growth retardation, reduced panicle number and pale green leaves with chlorotic spots. Cytological assay showed that there were reduced cell numbers, larger intercellular spaces, thicker cell walls, and abnormal development of chloroplast and mitochondria in the mature leaves from HW F1 plants in comparison with that from both of the parental lines. Genetic analysis revealed that HW was controlled by two complementary dominant genes Hw3 from V1134 and Hw4 from Taifeng A. Hw3 was mapped in a 136 kb interval between the markers Indel1118 and Indel1117 on chromosome 11, and Hw4 was mapped in the region of about 15 cM between RM182 and RM505 on chromosome 7, respectively. RT-PCR analysis revealed that only LOC_Os11g44310, encoding a putative calmodulin-binding protein (OsCaMBP), differentially expressed among Taifeng A, V1134 and their HW F1. No recombinant was detected using the markers designed based on the sequence of LOC_Os11g44310 in the BC1F2 (Taifeng A//Taifeng A/V1134) population. Hence, LOC_Os11g44310 was probably the candidate gene of Hw3. Gene amplification suggested that LOC_Os11g44310 was present in V1134 and absent in Taifeng A. BLAST search revealed that LOC_Os11g44310 had one copy in the japonica genomic sequence of Nipponbare, and no homologous sequence in the indica reference sequence of 9311. Our results indicate that Hw3 is a novel gene for inducing hybrid weakness in rice.

[Genetic analysis and mapping of a rice white midrib mutant Oswm2].

A rice (Oryza sativa L.) plant with white midrib (Oswm2) was selected from the T-DNA inserted mutant pool of rice. The basal midrib of the leaves in the mutant, except for the flag leaves, is white. The blade close to the white midrib is etiolated, and alterations of the agronomic traits, such as plant height, panicle length occur in the mutant. Genetic analysis indicated that this mutant trait was controlled by a recessive nuclear gene. To map the OsWM2 gene, an F2 population was constructed by crossing mutant Oswm2 with 02428. The OsWM2 locus was preliminarily located between the SSR molecular markers RM21478 and RM418 on chromosome 7 with the distances of 8.7 and 15.9 cM, respectively.