Epigenetic modifications and miRNAs determine the transition of somatic cells into somatic embryos

KEY MESSAGE

This review discusses the epigenetic changes during somatic embryo (SE) development, highlights the genes and miRNAs involved in the transition of somatic cells into SEs as a result of epigenetic changes, and draws insights on biotechnological opportunities to study SE development. Somatic embryogenesis from somatic cells occurs in a series of steps. The transition of somatic cells into somatic embryos (SEs) is the most critical step under genetic and epigenetic regulations. Major regulatory genes such as SERK, WUS, BBM, FUS3/FUSA3, AGL15, and PKL, control SE steps and development by turning on and off other regulatory genes. Gene transcription profiles of somatic cells during SE development is the result of epigenetic changes, such as DNA and histone protein modifications, that control and decide the fate of SE formation. Depending on the type of somatic cells and the treatment with plant growth regulators, epigenetic changes take place dynamically. Either hypermethylation or hypomethylation of SE-related genes promotes the transition of somatic cells. For example, the reduced levels of DNA methylation of SERK and WUS promotes SE initiation. Histone modifications also promote SE induction by regulating SE-related genes in somatic cells. In addition, miRNAs contribute to the various stages of SE by regulating the expression of auxin signaling pathway genes (TIR1, AFB2, ARF6, and ARF8), transcription factors (CUC1 and CUC2), and growth-regulating factors (GRFs) involved in SE formation. These epigenetic and miRNA functions are unique and have the potential to regenerate bipolar structures from somatic cells when a pluripotent state is induced. However, an integrated overview of the key regulators involved in SE development and downstream processes is lacking. Therefore, this review discusses epigenetic modifications involved in SE development, SE-related genes and miRNAs associated with epigenetics, and common cis-regulatory elements in the promoters of SE-related genes. Finally, we highlight future biotechnological opportunities to alter epigenetic pathways using the genome editing tool and to study the transition mechanism of somatic cells.

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

Genome-wide association study reveals novel genomic regions governing agronomic and grain quality traits and superior allelic combinations for Basmati rice improvement

BACKGROUND

Basmati is a speciality segment in the rice genepool characterised by explicit grain quality. For the want of suitable populations, genome-wide association study (GWAS) in Basmati rice has not been attempted.

MATERIALS

To address this gap, we have performed a GWAS on a panel of 172 elite Basmati multiparent population comprising of potential restorers and maintainers. Phenotypic data was generated for various agronomic and grain quality traits across seven different environments during two consecutive crop seasons. Based on the observed phenotypic variation, three agronomic traits namely, days to fifty per cent flowering, plant height and panicle length, and three grain quality traits namely, kernel length before cooking, length breadth ratio and kernel length after cooking were subjected to GWAS. Genotyped with 80K SNP array, the population was subjected to principal component analysis to stratify the underlying substructure and subjected to the association analysis using Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) model.

RESULTS

We identified 32 unique MTAs including 11 robust MTAs for the agronomic traits and 25 unique MTAs including two robust MTAs for the grain quality traits. Six out of 13 robust MTAs were novel. By genome annotation, six candidate genes associated with the robust MTAs were identified. Further analysis of the allelic combinations of the robust MTAs enabled the identification of superior allelic combinations in the population. This information was utilized in selecting 77 elite Basmati rice genotypes from the panel.

CONCLUSION

This is the first ever GWAS study in Basmati rice which could generate valuable information usable for further breeding through marker assisted selection, including enhancing of heterosis.

From farm to plate: Spatio-temporal characterization revealed compositional changes and reduced retention of γ-oryzanol upon processing in rice

Background

Antioxidants detain the development and proliferation of various non-communicable diseases (NCDs). γ-oryzanol, a group of steryl ferulates and caffeates, is a major antioxidant present in rice grain with proven health benefits. The present study evaluated the distribution and dynamics of γ-oryzanol and its components in spatial and temporal scales and also delineated the effect of processing and cooking on its retention.

Methods

Six rice varieties (four Basmati and two non-Basmati) belonging to indica group were analyzed at spatial scale in four different tissues (leaf blades, leaf sheaths, peduncle and spikelets) and temporal scale at three developmental stages (booting, milky and dough). Additionally, the matured grains were fractioned into husk, embryo, bran, and endosperm to assess differential accumulation in these tissues. Further, milling and cooking of the samples was done to assess the retention upon processing. After extraction of γ-oryzanol by solvent extraction method, individual components were identified by UPLC-QToF-ESI-MS and quantified by RP-HPLC.

Results

The non-seed tissues were significantly different from the seed tissues for composition and quantitative variation of γ-oryzanol. Cycloartenyl caffeate was predominant in all the non-seed tissues during the three developmental stages while it showed significant reduction during the growth progression toward maturity and was totally absent in the matured grains. In contrary, the 24-methylenecycloartanyl ferulate, campesteryl ferulate and β-sitosteryl ferulate showed significant increment toward the growth progression to maturity. Milling caused significant reduction, retaining only an average of 58.77% γ-oryzanol. Cooking of brown rice in excess water showed relatively lower average retention (43.31%) to samples cooked in minimal water (54.42%). Cooked milled rice showed least mean retention of 21.66%.

Conclusion

The results demonstrate prominent compositional variation of γ-oryzanol during different growth stages. For the first time, the study demonstrated that ferulate esters of γ-oryzanol were predominant in the seed tissues while caffeate esters were dominant in non-seed tissues. Basmati cultivars show differential expression of γ-oryzanol and its components compared to non-Basmati cultivars. Cooking in excess water causes maximum degradation of γ-oryzanol. Post-harvest losses due to milling and cooking indicate the necessity of biofortification for γ-oryzanol content in rice grain.

Molecular mapping of QTLs for grain dimension traits in Basmati rice

Basmati rice is known for its extra-long slender grains, exceptional kernel dimensions after cooking, high volume expansion, and strong aroma. Developing high yielding Basmati rice varieties with good cooking quality is a gigantic task. Therefore, identifying the genomic regions governing the grain and cooked kernel dimension traits is of utmost importance for its use in marker-assisted breeding. Although several QTLs governing grain dimension traits have been reported, limited attempts have been made to map QTLs for grain and cooked kernel dimension traits of Basmati rice. In the current study, a population of recombinant inbred lines (RIL) was generated from a cross of Sonasal and Pusa Basmati 1121 (PB1121). In the RIL population, there was a significant positive correlation among the length (RRL: rough rice length, MRL: milled rice length, CKL: cooked kernel length) and breadth (RRB: rough rice breadth, MRB: milled rice breadth and CKB: cooked kernel breadth) of the related traits, while there was significant negative correlation between them. QTL mapping has led to the identification of four major genomic regions governing MRL and CKL. Two QTLs co-localize with the earlier reported major gene GS3 and a QTL qGRL7.1, while the remaining two QTLs viz., qCKL3.2 (qMRL3.2) and qCKL4.1 (qMRL4.1) were novel. The QTL qCKL3.2 has been bracketed to a genomic region of 0.78 Mb between the markers RM15247 and RM15281. Annotation of this region identified 18 gene models, of which the genes predicted to encode pentatricopeptides and brassinosteroid insensitive 1-associated receptor kinase 1 precursor may be the putative candidate genes. Furthermore, we identified a novel QTL qKER2.1 governing kernel elongation ratio (KER) in Basmati rice.

The plant epitranscriptome: revisiting pseudouridine and 2'-O-methyl RNA modifications

There is growing evidence that post-transcriptional RNA modifications are highly dynamic and can be used to improve crop production. Although more than 172 unique types of RNA modifications have been identified throughout the kingdom of life, we are yet to leverage upon the understanding to optimize RNA modifications in crops to improve productivity. The contributions of internal mRNA modifications such as N6-methyladenosine (m⁶ A) and 5-methylcytosine (m⁵ C) methylations to embryonic development, root development, leaf morphogenesis, flowering, fruit ripening and stress response are sufficiently known, but the roles of the two most abundant RNA modifications, pseudouridine (Ψ) and 2'-O-methylation (Nm), in the cell remain unclear due to insufficient advances in high-throughput technologies in plant development. Therefore, in this review, we discuss the latest methods and insights gained in mapping internal Ψ and Nm and their unique properties in plants and other organisms. In addition, we discuss the limitations that remain in high-throughput technologies for qualitative and quantitative mapping of these RNA modifications and highlight future challenges in regulating the plant epitranscriptome.

Genome-Wide Association Mapping Reveals Novel Putative Gene Candidates Governing Reproductive Stage Heat Stress Tolerance in Rice

Temperature rise predicted for the future will severely affect rice productivity because the crop is highly sensitive to heat stress at the reproductive stage. Breeding tolerant varieties is an economically viable option to combat heat stress, for which the knowledge of target genomic regions associated with the reproductive stage heat stress tolerance (RSHT) is essential. A set of 192 rice genotypes of diverse origins were evaluated under natural field conditions through staggered sowings for RSHT using two surrogate traits, spikelet fertility and grain yield, which showed significant reduction under heat stress. These genotypes were genotyped using a 50 k SNP array, and the association analysis identified 10 quantitative trait nucleotides (QTNs) for grain yield, of which one QTN (qHTGY8.1) was consistent across the different models used. Only two out of 10 MTAs coincided with the previously reported QTLs, making the remaing eight novel. A total of 22 QTNs were observed for spikelet fertility, among which qHTSF5.1 was consistently found across three models. Of the QTNs identified, seven coincided with previous reports, while the remaining QTNs were new. The genes near the QTNs were found associated with the protein–protein interaction, protein ubiquitination, stress signal transduction, and so forth, qualifying them to be putative for RSHT. An in silico expression analysis revealed the predominant expression of genes identified for spikelet fertility in reproductive organs. Further validation of the biological relevance of QTNs in conferring heat stress tolerance will enable their utilization in improving the reproductive stage heat stress tolerance in rice.

Genome-wide identification and expression characterization of the DoG gene family of moso bamboo (Phyllostachys edulis)

BACKGROUND

The DoG (Delay of Germination1) family plays a key regulatory role in seed dormancy and germination. However, to date, there is no complete genomic overview of the DoG gene family of any economically valuable crop, including moso bamboo (Phyllostachys edulis), and no studies have been conducted to characterize its expression profile. To identify the DoG gene members of moso bamboo (PeDoG) and to investigate their family structural features and tissue expression profile characteristics, a study was conducted. Based on the whole genome and differential transcriptome data, in this investigation, we have scrutinized the physicochemical properties, gene structure, cis-acting elements, phylogenetic relationships, conserved structural (CS) domains, CS motifs and expression patterns of the PeDoG1 family of moso bamboo.

RESULTS

The DoG family genes of moso bamboo were found distributed across 16 chromosomal scaffolds with 24 members. All members were found to carry DoG1 structural domains, while 23 members additionally possessed basic leucine zipper (bZIP) structural domains. We could divide the PeDoG genes into three subfamilies based on phylogenetic relationships. Covariance analysis revealed that tandem duplication was the main driver of amplification of the PeDoG genes. The upstream promoter of these genes containing several cis-acting elements indicates a plausible role in abiotic stress and hormone induction. Gene expression pattern according to transcriptome data revealed participation of the PeDoG genes in tissue and organ development. Analysis using Short Time-series Expression Miner (STEM) tool revealed that the PeDoG gene family is also associated with rapid early shoot growth. Gene ontology (GO) and KEGG analyses showed a dual role of the PeDoG genes. We found that PeDoGs has a possible role as bZIP transcription factors by regulating Polar like1 (PL1) gene expression, and thereby playing a disease response role in moso bamboo. Quantitative gene expression of the PeDoG genes revealed that they were abundantly expressed in roots and leaves, and could be induced in response to gibberellin (GA).

CONCLUSION

In this study, we found that the PeDoG genes are involved in a wide range of activities such as growth and development, stress response and transcription. This forms the first report of PeDoG genes and their potential roles in moso bamboo.

Transcriptome analysis of a near-isogenic line and its recurrent parent reveals the role of Pup1 QTL in phosphorus deficiency tolerance of rice at tillering stage

Phosphorus (P) is essential for cellular processes like respiration, photosynthesis, biosynthesis of membrane phospholipids, etc. To cope with P deficiency stress, plants adopt reprograming of the expression of genes involved in different metabolic/signaling pathways for survival, growth, and development. Plants use transcriptional, post-transcriptional, and/or post-translational machinery to achieve P homeostasis. Several transcription factors (TFs), miRNAs, and P transporters play important roles in P deficiency tolerance; however, the underlying mechanisms responsible for P deficiency tolerance remain poorly understood. Studies on P starvation/deficiency responses in plants at early (seedling) stage of growth have been reported but only a few of them focused on molecular responses of the plant at advanced (tillering or reproductive) stage of growth. To decipher the strategies adopted by rice at tillering stage under P deficiency stress, a pair of contrasting genotypes [Pusa-44 (a high-yielding, P deficiency sensitive cultivar) and its near-isogenic line (NIL-23, P deficiency tolerant) for Pup1 QTL] was used for morphophysiological, biochemical, and molecular analyses. Comparative analyses of shoot and root tissues from 45-day-old plants grown hydroponically under P sufficient (16 ppm) or P deficient (4 ppm) medium confirmed some of the known morphophysiological responses. Moreover, RNA-seq analysis revealed the important roles of phosphate transporters, TFs, auxin-responsive proteins, modulation in the cell wall, fatty acid metabolism, and chromatin architecture/epigenetic modifications in providing P deficiency tolerance to NIL-23, which were brought in due to the introgression of the Pup1 QTL in Pusa-44. This study provides insights into the molecular functions of Pup1 for P deficiency tolerance, which might be utilized to improve P-use efficiency of rice for better productivity in P deficient soils. KEY MESSAGE: Introgression of Pup1 QTL in high-yielding rice cultivar modulates mainly phosphate transporters, TFs, auxin-responsive proteins, cell wall structure, fatty acid metabolism, and chromatin architecture/epigenetic modifications at tillering stage of growth under phosphorus deficiency stress.

Redox status of the plant cell determines epigenetic modifications under abiotic stress conditions and during developmental processes

BACKGROUND

The oxidation-reduction (redox) status of the cell influences or regulates transcription factors and enzymes involved in epigenetic changes, such as DNA methylation, histone protein modifications, and chromatin structure and remodeling. These changes are crucial regulators of chromatin architecture, leading to differential gene expression in eukaryotes. But the cell's redox homeostasis is difficult to sustain since the production of reactive oxygen species (ROS) and reactive nitrogen species (RNS) is not equal in plants at different developmental stages and under abiotic stress conditions. Exceeding optimum ROS and RNS levels leads to oxidative stress and thus alters the redox status of the cell. Consequently, this alteration modulates intracellular epigenetic modifications that either mitigate or mediate the plant growth and stress response.

AIM OF REVIEW

Recent studies suggest that the altered redox status of the cell reform the cellular functions and epigenetic changes. Recent high-throughput techniques have also greatly advanced redox-mediated gene expression discovery, but the integrated view of the redox status, and its associations with epigenetic changes and subsequent gene expression in plants are still scarce. In this review, we accordingly focus on how the redox status of the cell affects epigenetic modifications in plants under abiotic stress conditions and during developmental processes. This is a first comprehensive review on the redox status of the cell covering the redox components and signaling, redox status alters the post-translational modification of proteins, intracellular epigenetic modifications, redox interplay during DNA methylation, redox regulation of histone acetylation and methylation, redox regulation of miRNA biogenesis, redox regulation of chromatin structure and remodeling and conclusion, future perspectives and biotechnological opportunities for the future development of the plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

The interaction of redox mediators such as ROS, RNS and antioxidants regulates redox homeostasis and redox-mediated epigenetic changes. We discuss how redox mediators modulate epigenetic changes and show the opportunities for smart use of the redox status of the cell in plant development and abiotic stress adaptation. However, how a redox mediator triggers epigenetic modification without activating other redox mediators remains yet unknown.

Identification of potential MTAs and candidate genes for juice quality- and yield-related traits in Saccharum clones: a genome-wide association and comparative genomic study

Sugarcane is an economically important commercial crop which provides raw material for the production of sugar, jaggery, bioethanol, biomass and other by-products. Sugarcane breeding till today heavily relies on conventional breeding approaches which is time consuming, laborious and costly. Integration of marker-assisted selection (MAS) in sugarcane genetic improvement programs for difficult to select traits like sucrose content, resistance to pests and diseases and tolerance to abiotic stresses will accelerate varietal development. In the present study, association mapping approach was used to identify QTLs and genes associated with sucrose and other important yield-contributing traits. A mapping panel of 110 diverse sugarcane genotypes and 148 microsatellite primers were used for structured association mapping study. An optimal subpopulation number (ΔK) of 5 was identified by structure analysis. GWAS analysis using TASSEL identified a total of 110 MTAs which were localized into 27 QTLs by GLM and MLM (Q + K, PC + K) approaches. Among the 24 QTLs sequenced, 12 were able to identify potential candidate genes, viz., starch branching enzyme, starch synthase 4, sugar transporters and G3P-DH related to carbohydrate metabolism and hormone pathway-related genes ethylene insensitive 3-like 1, reversion to ethylene sensitive1-like, and auxin response factor associated to juice quality- and yield-related traits. Six markers, NKS 5_185, SCB 270_144, SCB 370_256, NKS 46_176 and UGSM 648_245, associated with juice quality traits and marker SMC31CUQ_304 associated with NMC were validated and identified as significantly associated to the traits by one-way ANOVA analysis. In conclusion, 24 potential QTLs identified in the present study could be used in sugarcane breeding programs after further validation in larger population. The candidate genes from carbohydrate and hormone response pathway presented in this study could be manipulated with genome editing approaches to further improve sugarcane crop.

A novel LOX3-null allele (lox3-b) originated in the aromatic Basmati rice cultivars imparts storage stability to rice bran

Rapid deterioration of rice bran due to the LOX3 enzyme catalysed oxidation of PUFA is the major bottleneck for its utilization in various downstream applications. In the present study, we have identified a set of nine novel LOX3-null rice accessions carrying a deletion of C residue in the exon2 causing a frameshift mutation resulting in a truncated non-functional LOX3 protein. Our study, further manifested the predominance of C deletion based LOX3-null allele, named lox3-b, in the aromatic rice germplasm particularly in the Indian Basmati rice group. The LOX3-null genotypes exhibited significantly reduced rancidity, after six months of storage. They also showed significantly lower percentage reduction of linoleic acid (LA), higher γ-oryzanol content and lower hexanal content. A functional dCAPS marker designed based on the deletion polymorphism clearly differentiated LOX3 and lox3-b alleles, and has the potential application in marker assisted rice breeding programmes to develop cultivars with better bran storability.

Long terminal repeats (LTR) and transcription factors regulate PHRE1 and PHRE2 activity in Moso bamboo under heat stress

BACKGROUND

LTR retrotransposons play a significant role in plant growth, genome evolution, and environmental stress response, but their regulatory response to heat stress remains unclear. We have investigated the activities of two LTR retrotransposons, PHRE1 and PHRE2, of moso bamboo (Phyllostachys edulis) in response to heat stress.

RESULTS

The differential overexpression of PHRE1 and PHRE2 with or without CaMV35s promoter showed enhanced expression under heat stress in transgenic plants. The transcriptional activity studies showed an increase in transposition activity and copy number among moso bamboo wild type and Arabidopsis transgenic plants under heat stress. Comparison of promoter activity in transgenic plants indicated that 5'LTR promoter activity was higher than CaMV35s promoter. Additionally, yeast one-hybrid (Y1H) system and in planta biomolecular fluorescence complementation (BiFC) assay revealed interactions of heat-dependent transcription factors (TFs) with 5'LTR sequence and direct interactions of TFs with pol and gag.

CONCLUSIONS

Our results conclude that the 5'LTR acts as a promoter and could regulate the LTR retrotransposons in moso bamboo under heat stress.

Characterization of contrasting rice (Oryza sativa L.) genotypes reveals the Pi-efficient schema for phosphate starvation tolerance

BACKGROUND

Phosphorus (P), being one of the essential components of nucleic acids, cell membranes and enzymes, indispensable for diverse cellular processes like photosynthesis/carbohydrate metabolism, energy production, redox homeostasis and signaling. Crop yield is severely affected due to Phosphate (Pi) deficiency; and to cope with Pi-deficiency, plants have evolved several strategies. Some rice genotypes are compatible with low Pi availability, whereas others are sensitive to Pi deficiency. However, the underlying molecular mechanism for low Pi tolerance remains largely unexplored.

RESULT

Several studies were carried out to understand Pi-deficiency responses in rice at seedling stage, but few of them targeted molecular aspects/responses of Pi-starvation at the advanced stage of growth. To delineate the molecular mechanisms for low Pi tolerance, a pair of contrasting rice (Oryza sativa L.) genotypes [viz. Pusa-44 (Pi-deficiency sensitive) and its near isogenic line (NIL-23, Pi-deficiency tolerant) harboring Phosphorus uptake 1 (Pup1) QTL from an aus landrace Kasalath] were used. Comparative morphological, physiological, and biochemical analyses confirmed some of the well-known findings. Transcriptome analysis of shoot and root tissues from 45-day-old rice plants grown hydroponically under P-sufficient (16 ppm Pi) or P-starved (0 ppm Pi) medium revealed that Pi-starvation stress causes global transcriptional reprogramming affecting several transcription factors, signaling pathways and other regulatory genes. We could identify several significantly up-regulated genes in roots of NIL-23 under Pi-starvation which might be responsible for the Pi starvation tolerance. Pathway enrichment analysis indicated significant role of certain phosphatases, transporters, transcription factors, carbohydrate metabolism, hormone-signaling, and epigenetic processes in improving P-starvation stress tolerance in NIL-23.

CONCLUSION

We report the important candidate mechanisms for Pi acquisition/solubilization, recycling, remobilization/transport, sensing/signalling, genetic/epigenetic regulation, and cell wall structural changes to be responsible for P-starvation tolerance in NIL-23. The study provides some of the novel information useful for improving phosphorus-use efficiency in rice cultivars.

Genome-Wide Association Study Reveals Marker-Trait Associations for Early Vegetative Stage Salinity Tolerance in Rice

Rice germplasm is a rich resource for discovering genes associated with salt tolerance. In the current study, a set of 96 accessions were evaluated for seedling stage salinity tolerance and its component traits. Significant phenotypic variation was observed among the genotypes for all the measured traits and eleven accessions with high level of salt tolerance at seedling stage were identified. The germplasm set comprised of three sub-populations and genome-wide association study (GWAS) identified a total of 23 marker-trait associations (MTAs) for traits studied. These MTAs were located on rice chromosomes 1, 2, 5, 6, 7, 9, and 12 and explained the trait phenotypic variances ranging from 13.98 to 29.88 %. Twenty-one MTAs identified in this study were located either in or near the previously reported quantitative trait loci (QTLs), while two MTAs namely, qSDW2.1 and qSNC5 were novel. A total of 18 and 13 putative annotated candidate genes were identified in a genomic region spanning ~200 kb around the MTAs qSDW2.1 and qSNC5, respectively. Some of the important genes underlying the novel MTAs were OsFBA1,OsFBL7, and mTERF which are known to be associated with salinity tolerance in crops. These MTAs pave way for combining salinity tolerance with high yield in rice genotypes through molecular breeding.

Multi-omics analysis of cellular pathways involved in different rapid growth stages of moso bamboo

Moso bamboo (Phyllostachys edulis (Carriere) J. Houzeau) is a rapidly growing grass of industrial and ecological importance. However, the molecular mechanisms of its remarkable growth are not well understood. In this study, we investigated the early-stage growth of moso bamboo shoots and defined three different growth stages based on histological and biochemical analyses, namely, starting of cell division (SD), rapid division (RD) and rapid elongation (RE). Further analyses on potentially relevant cellular pathways in these growth stages using multi-omics approaches such as transcriptomics and proteomics revealed the involvement of multiple cellular pathways, including DNA replication, repair and ribosome biogenesis. A total of 8045 differentially expressed genes (DEGs) and 1053 differentially expressed proteins (DEPs) were identified in our analyses. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses of detected DEGs identified several key biological pathways such as phytohormone metabolism, signal transduction, cell wall development and carbohydrate metabolism. The comparative analysis of proteins displayed that a total of 213 DEPs corresponded with DEGs and 3 significant expression profiles that could be promoting the fast growth of bamboo internodes. Moreover, protein-protein interaction network prediction analysis is suggestive of the involvement of five major proteins of signal transduction, DNA synthesis and RNA transcription, and may act as key elements responsible for the rapid shoot growth. Our work exploits multi-omics and bioinformatic approaches to unfurl the complexity of molecular networks involved in the rapid growth of moso bamboo and opens up questions related to the interactions between the functions played by individual molecular pathway.

Marker Assisted Development and Characterization of Herbicide Tolerant Near Isogenic Lines of a Mega Basmati Rice Variety, "Pusa Basmati 1121"

BACKGROUND

Direct-seeded rice (DSR) is a potential technology for sustainable rice farming as it saves water and labor. However, higher incidence of weed under DSR limits productivity. Therefore, there is a need to develop herbicide tolerant (HT) rice varieties.

RESULTS

We used marker assisted backcross breeding (MABB) to transfer a mutant allele of Acetohydroxy acid synthase (AHAS) gene, which confers tolerance to imidazolinone group of herbicides from the donor parent (DP) "Robin" into the genetic background of an elite popular Basmati rice variety, Pusa Basmati 1121 (PB 1121). Foreground selection was done using the AHAS gene linked Simple Sequence Repeat (SSR) marker RM6844 and background selection was performed using 112 genome-wide SSR markers polymorphic between PB 1121 and Robin. Phenotypic selection for agronomic, Basmati grain and cooking quality traits in each generation was carried out to improve the recovery of recurrent parent phenome (RPP). Finally, a set of 12 BC₄F₄ near isogenic lines (NILs), with recurrent parent genome (RPG) recovery ranging from 98.66 to 99.55% were developed and evaluated. PB 1121-HT NILs namely 1979-14-7-33-99-10, 1979-14-7-33-99-15 and 1979-14-7-33-99-66 were found superior to PB 1121 in yield with comparable grain and cooking quality traits and herbicide tolerance similar to Robin.

CONCLUSION

Overall, the present study reports successful development of HT NILs in the genetic background of popular Basmati rice variety, PB 1121 by introgression of mutated AHAS allele. This is the first report on the development of HT Basmati rice. Superior NILs are being evaluated in the national Basmati trials, the release of which is likely to provide a viable option for the adoption of DSR technology in Basmati rice cultivation.

Discovery of a Novel Induced Polymorphism in SD1 Gene Governing Semi-Dwarfism in Rice and Development of a Functional Marker for Marker-Assisted Selection

The semi-dwarfing allele, sd1-d, has been widely utilized in developing high-yielding rice cultivars across the world. Originally identified from the rice cultivar Dee-Geo-Woo-Gen (DGWG), sd1-d, derived from a spontaneous mutation, has a 383-bp deletion in the SD1 gene. To date, as many as seven alleles of the SD1 gene have been identified and used in rice improvement, either with a functional single-nucleotide polymorphism (SNP), with insertion–deletions (InDels), or both. Here, we report discovery of a novel SNP in the SD1 gene from the rice genotype, Pusa 1652. Genetic analysis revealed that the inheritance of the semi-dwarfism in Pusa 1652 is monogenic and recessive, but it did not carry the sd1-d allele. However, response to exogenous gibberellic acid (GA₃) application and the subsequent bulked segregant and linkage analyses confirmed that the SD1 gene is involved in the plant height reduction in Pusa 1652. Sequencing of the SD1 gene from Pusa 1652 revealed a novel transition in exon 3 (T/A) causing a nonsense mutation at the 300th codon. The stop codon leads to premature termination, resulting in a truncated protein of OsGA20ox2 obstructing the GA₃ biosynthesis pathway. This novel recessive allele, named sd1-bm, is derived from Bindli Mutant 34 (BM34), a γ-ray induced mutant of a short-grain aromatic landrace, Bindli. BM34 is the parent of an aromatic semi-dwarf cultivar, Pusa 1176, from which Pusa 1652 is derived. The semi-dwarfing allele, sd1-bm, was further validated by developing a derived cleaved amplified polymorphic sequence (dCAPS) marker, AKS-sd1. This allele provides an alternative to the most widely used sd1-d in rice improvement programs and the functional dCAPS marker will facilitate marker-assisted introgression of the semi-dwarf trait into tall genotypes.

Marker aided introgression of 'Saltol', a major QTL for seedling stage salinity tolerance into an elite Basmati rice variety 'Pusa Basmati 1509'

Marker assisted backcross breeding was used to transfer Saltol, a major QTL for seedling stage salinity tolerance from the donor FL478 to Pusa Basmati 1509 (PB 1509), a high yielding and early maturing Basmati rice variety. Foreground selection was carried out using three markers namely, AP3206f, RM3412b and RM10793, linked to Saltol. In addition, 105 genome-wide SSR markers polymorphic between FL478 and PB 1509 were used in background selection. Among the BC₃F₄ near isogenic lines (NILs) developed, recurrent parent genome recovery ranged from 96.67 to 98.57%. Multi-season evaluation identified some of the NILs showing significantly higher yield with grain and cooking quality comparable to PB 1509. All the NILs exhibited tolerance to salinity with significantly higher relative water content, membrane stability index and proline content as compared to PB 1509. The root and shoot concentration of Na⁺, K⁺ and Na⁺/K⁺ in NILs was at par with FL478 under stress conditions. The gene OsHKT1;5 located in the Saltol region showed higher expression levels under stress indicating its role in conferring salinity tolerance. Salt tolerant NILs of PB 1509 will be useful in stabilizing production in salt affected areas.

Genome-Wide Association Study Reveals Novel Marker-Trait Associations (MTAs) Governing the Localization of Fe and Zn in the Rice Grain

Micronutrient malnutrition due to Fe and Zn, affects around two billion people globally particularly in the developing countries. More than 90% of the Asian population is dependent on rice-based diets, which is low in these micronutrients. In the present study, a set of 192 Indian rice germplasm accessions, grown at two locations, were evaluated for Fe and Zn in brown rice (BR) and milled rice (MR). A significant variation was observed in the rice germplasm for these micronutrients. The grain Fe concentration was in the range of 6.2–23.1 ppm in BR and 0.8–12.3 ppm in MR, while grain Zn concentration was found to be in the range of 11.0–47.0 ppm and 8.2–40.8 ppm in the BR and MR, respectively. Grain Fe exhibited maximum loss upon milling with a mean retention of 24.9% in MR, while Zn showed a greater mean retention of 74.2% in MR. A genome-wide association study (GWAS) was carried out implementing the FarmCPU model to control the population structure and kinship, and resulted in the identification of 29 marker-trait associations (MTAs) with significant associations for traits viz. FeBR (6 MTAs), FeMR (7 MTAs), ZnBR (11 MTAs), and ZnMR (5 MTAs), which could explain the phenotypic variance from 2.1 to as high as 53.3%. The MTAs governing the correlated traits showed co-localization, signifying the possibility of their simultaneous improvement. The robust MTAs identified in the study could be valuable resource for enhancing Fe and Zn concentration in the rice grain and addressing the problem of Fe and Zn malnutrition among rice consumers.

Nuclear export signal (NES) of transposases affects the transposition activity of mariner-like elements Ppmar1 and Ppmar2 of moso bamboo

Ppmar1 and Ppmar2 are two active mariner-like elements (MLEs) cloned from moso bamboo (Phyllostachys edulis (Carrière) J. Houz) genome possessing transposases that harbour nuclear export signal (NES) domain, but not any nuclear localization signal (NLS) domain. To understand the functions of NES in transposon activity, we have conducted two experiments, fluorescence and excision frequency assays in the yeast system. For this, by site-directed mutagenesis, three NES mutants were developed from each of the MLE. In the fluorescence assay, the mutants, NES-1, 2 and 3 along with the wild types (NES-0) were fused with fluorescent proteins, enhanced yellow fluorescent protein (EYFP) and enhanced cyan fluorescent protein (ECFP) were co-transformed into yeast system. To differentiate protein localisation under the NES influence, ECFP alone was fused to wild and mutant NES domains either on N- or C-terminal and not to EYFP. Fluorescence assay revealed that blue fluorescence of ECFP was more intense than the red fluorescence of the EYFP in the yeast cell matrix. Further, ECFP had a wider localisation in the cellular matrix, but EYFP was largely located in the nucleus. The NES-1 domain was related to the comparatively high spread of ECFP, while NES-2 and NES-3 indicated a low spread, implying that NES activity on nuclear export increased when the NES is made leucine-rich, while the signalling activity was reduced when the leucine content was lowered in the NES domain. In the transposon excision assay, the mutant and wild type NES of both the Ppmar elements were integrated into an Ade2 vector, and within the Ade2 gene. Co-transformation of the vector together with non-autonomous Ppmar transposons and NES-lacking transposases was used to assess the differential excision frequencies of the mutants NES domains. In both the MLEs, NES-1 had the highest excision suppression, which was less than half of the excision frequency of the wild type. NES-2 and NES-3 elements showed, up to three times increase in transposon excision than the wild types. The results suggested that NES is an important regulator of nuclear export of transposase in Ppmar elements and the mutation of the NES domains can either increase or decrease the export signalling. We speculate that in moso bamboo, NESs regulates the transposition activity of MLEs to maintain the genome integrity.

WA-CMS-based iso-cytoplasmic restorers derived from commercial rice hybrids reveal distinct population structure and genetic divergence towards restorer diversification

One hundred diverse iso-cytoplasmic restorer (ICR) lines carrying WA cytoplasm indicated significant but moderate variability for agro-morphological traits as well as for the microsatellite-based allele patterns. There were two major groups of ICRs based on agro-morphological clustering. Simple sequence repeat (SSR) markers identified allelic variants with an average of 2.48 alleles per locus and the gene diversity (GD) ranged from 0.02 to 0.62 at different loci. ICR lines showed a genetic structure involving two sub-populations, POP1 and POP2. Both the subpopulations had the presence of admixture lines. Nearest ancestry-based grouping of ICRs by neighbour-joining (NJ) method showed near similar grouping as that of sub-population division. The POP2 was the largest group but with fewer admixed lines. POP1 was more distinct than POP2. Since the hybrid parents of the ICRs had limited diversity on maternal lineage, paternal lineage was concluded as the major contributor to the observed divergence and population differentiation. ICRs developed from certain hybrids were more genetically distinct than other hybrids. Even with the moderate variability, ICRs could be considered as a potential source of fertility restoration in hybrid development because of their distinct population structure and the full complement of restorer genes they contained. ICR lines with high per se performance can be utilized in hybrid rice development by estimating their combining ability.

Affinities of Terminal Inverted Repeats to DNA Binding Domain of Transposase Affect the Transposition Activity of Bamboo Ppmar2 Mariner-Like Element

Mariner-like elements (MLE) are a super-family of DNA transposons widespread in animal and plant genomes. Based on their transposition characteristics, such as random insertions and high-frequency heterogeneous transpositions, several MLEs have been developed to be used as tools in gene tagging and gene therapy. Two active MLEs, Ppmar1 and Ppmar2, have previously been identified in moso bamboo (Phyllostachys edulis). Both of these have a preferential insertion affinity to AT-rich region and their insertion sites are close to random in the host genome. In Ppmar2 element, we studied the affinities of terminal inverted repeats (TIRs) to DNA binding domain (DBD) and their influence on the transposition activity. We could identify two putative boxes in the TIRs which play a significant role in defining the TIR's affinities to the DBD. Seven mutated TIRs were constructed, differing in affinities based on similarities with those of other plant MLEs. Gel mobility shift assays showed that the TIR mutants with mutation sites G669A-C671A had significantly higher affinities than the mutants with mutation sites C657T-A660T. The high-affinity TIRs indicated that their transposition frequency was 1.5-2.0 times higher than that of the wild type TIRs in yeast transposition assays. The MLE mutants with low-affinity TIRs had relatively lower transposition frequency from that of wild types. We conclude that TIR affinity to DBD significantly affects the transposition activity of Ppmar2. The mutant MLEs highly active TIRs constructed in this study can be used as a tool for bamboo genetic studies.

Marker-Assisted Introgression of Saltol QTL Enhances Seedling Stage Salt Tolerance in the Rice Variety "Pusa Basmati 1"

Marker-assisted selection is an unequivocal translational research tool for crop improvement in the genomics era. Pusa Basmati 1 (PB1) is an elite Indian Basmati rice cultivar sensitive to salinity. Here, we report enhanced seedling stage salt tolerance in improved PB1 genotypes developed through marker-assisted transfer of a major QTL, Saltol. A highly salt tolerant line, FL478, was used as the Saltol donor. Parental polymorphism survey using 456 microsatellite (SSR)/QTL-linked markers revealed 14.3% polymorphism between PB1 and FL478. Foreground selection was carried out using three Saltol-linked polymorphic SSR markers RM8094, RM493, and RM10793 and background selection by 62 genome-wide polymorphic SSR markers. In every backcross generation, foreground selection was restricted to the triple heterozygotes of foreground markers, which was followed by phenotypic and background selections. Twenty-four near isogenic lines (NILs), with recurrent parent genome recovery of 96.0-98.4%, were selected after two backcrosses followed by three selfing generations. NILs exhibited agronomic traits similar to those of PB1 and additional improvement in the seedling stage salt tolerance. They are being tested for per se performance under salt-affected locations for release as commercial varieties. These NILs appear promising for enhancing rice production in salinity-affected pockets of Basmati Geographical Indication (GI) areas of India.

Marker-aided Incorporation of Xa38, a Novel Bacterial Blight Resistance Gene, in PB1121 and Comparison of its Resistance Spectrum with xa13 + Xa21

Basmati rice is preferred internationally because of its appealing taste, mouth feel and aroma. Pusa Basmati 1121 (PB1121) is a widely grown variety known for its excellent grain and cooking quality in the international and domestic market. It contributes approximately USD 3 billion to India's forex earning annually by being the most traded variety. However, PB1121 is highly susceptible to bacterial blight (BB) disease. A novel BB resistance gene Xa38 was incorporated in PB1121 from donor parent PR114-Xa38 using a modified marker-assisted backcross breeding (MABB) scheme. Phenotypic selection prior to background selection was instrumental in identifying the novel recombinants with maximum recovery of recurrent parent phenome. The strategy was effective in delimiting the linkage drag to <0.5 mb upstream and <1.9 mb downstream of Xa38 with recurrent parent genome recovery upto 96.9% in the developed NILs. The NILs of PB1121 carrying Xa38 were compared with PB1121 NILs carrying xa13 + Xa21 (developed earlier in our lab) for their resistance to BB. Both NILs showed resistance against the Xoo races 1, 2, 3 and 6. Additionally, Xa38 also resisted Xoo race 5 to which xa13 + Xa21 was susceptible. The PB1121 NILs carrying Xa38 gene will provide effective control of BB in the Basmati growing region.

Improvement of Basmati rice varieties for resistance to blast and bacterial blight diseases using marker assisted backcross breeding

Marker assisted backcross breeding was employed to incorporate the blast resistance genes, Pi2 and Pi54 and bacterial blight (BB) resistance genes xa13 and Xa21 into the genetic background of Pusa Basmati 1121 (PB1121) and Pusa Basmati 6. Foreground selection for target gene(s) was followed by arduous phenotypic and background selection which fast-tracked the recovery of recurrent parent genome (RPG) to an extent of 95.8% in one of the near-isogenic lines (NILs) namely, Pusa 1728-23-33-31-56, which also showed high degree of resemblance to recurrent parent, PB6 in phenotype. The phenotypic selection prior to background selection provided an additional opportunity for identifying the novel recombinants viz., Pusa 1884-9-12-14 and Pusa 1884-3-9-175, superior to parental lines in terms of early maturity, higher yield and improved quality parameters. There was no significant difference between the RPG recovery estimated based on SSR or SNP markers, however, the panel of SNPs markers was considered as the better choice for background selection as it provided better genome coverage and included SNPs in the genic regions. Multi-location evaluation of NILs depicted their stable and high mean performance in comparison to the respective recurrent parents. The Pi2+Pi54 carrying NILs were effective in combating a pan-India panel of Magnaporthe oryzae isolates with high level of field resistance in northern, eastern and southern parts of India. Alongside, the PB1121-NILs and PB6-NILs carrying BB resistance genes xa13+Xa21 were resistant against Xanthomonas oryzae pv. oryzae races of north-western, southern and eastern parts of the country. Three of NILs developed in this study, have been promoted to final stage of testing during the ​Kharif 2015 in the Indian National Basmati Trial.

Development and evaluation of near-isogenic lines for major blast resistance gene(s) in Basmati rice

A set of NILs carrying major blast resistance genes in a Basmati rice variety has been developed. Also, the efficacy of pyramids over monogenic NILs against rice blast pathogen Magnaporthe oryzae has been demonstrated. Productivity and quality of Basmati rice is severely affected by rice blast disease. Major genes and QTLs conferring resistance to blast have been reported only in non-Basmati rice germplasm. Here, we report incorporation of seven blast resistance genes from the donor lines DHMASQ164-2a (Pi54, Pi1, Pita), IRBLz5-CA (Pi2), IRBLb-B (Pib), IRBL5-M (Pi5) and IRBL9-W (Pi9) into the genetic background of an elite Basmati rice variety Pusa Basmati 1 (PB1). A total of 36 near-isogenic lines (NILs) comprising of 14 monogenic, 16 two-gene pyramids and six three-gene pyramids were developed through marker-assisted backcross breeding (MABB). Foreground, recombinant and background selection was used to identify the plants with target gene(s), minimize the linkage drag and increase the recurrent parent genome (RPG) recovery (93.5-98.6 %), respectively, in the NILs. Comparative analysis performed using 50,051 SNPs and 500 SSR markers revealed that the SNPs provided better insight into the RPG recovery. Most of the monogenic NILs showed comparable performance in yield and quality, concomitantly, Pusa1637-18-7-6-20 (Pi9), was significantly superior in yield and stable across four different environments as compared to recurrent parent (RP) PB1. Further, among the pyramids, Pusa1930-12-6 (Pi2+Pi5) showed significantly higher yield and Pusa1633-7-8-53-6-8 (Pi54+Pi1+Pita) was superior in cooking quality as compared to RP PB1. The NILs carrying gene Pi9 were found to be the most effective against the concoction of virulent races predominant in the hotspot locations for blast disease. Conversely, when analyzed under artificial inoculation, three-gene pyramids expressed enhanced resistance as compared to the two-gene and monogenic NILs.

Molecular evaluation of genetic diversity and association studies in rice (Oryza sativa L.)

In the present study, we tested rice genotypes that included un(der)exploited landraces of Tamil Nadu along with indica and japonica test cultivars to ascertain their genetic diversity structure. Highly polymorphic microsatellite markers were used for generating marker segregation data. A novel measure, allele discrimination index, was used to determine subpopulation differentiation power of each marker. Phenotypic data were collected for yield and component traits. Pattern of molecular differentiation separated indica and japonica genotypes; indica genotypes had two subpopulations within. Landraces were found to have indica genome, but formed a separate subgroup with low linkage disequilibrium. The landraces further separated into distinct group in both hierarchical clustering analysis using neighbour-joining method as well as in the model based population structure analysis. Japonica and the remaining indica cultivars formed two other distinct groups. Linkage disequilibrium observed in the whole population was considerably reduced in subpopulations. Low linkage disequilibrium of landforms suggests their narrow adaptation in local geographical niche. Many population specific alleles could be identified particularly for japonica cultivars and landraces. Association analysis revealed nine marker-trait associations with three agronomic traits, of which 67% were previously reported. Although the testing landraces together with known cultivars had permitted genomewide association mapping, the experiment offers scope to study more landraces collected from the entire geographical region for drawing more reliable information.

QTL and QTL x environment effects on agronomic and nitrogen acquisition traits in rice

Agricultural environments deteriorate due to excess nitrogen application. Breeding for low nitrogen responsive genotypes can reduce soil nitrogen input. Rice genotypes respond variably to soil available nitrogen. The present study attempted quantification of genotype x nitrogen level interaction and mapping of quantitative trait loci (QTLs) associated with nitrogen use efficiency (NUE) and other associated agronomic traits. Twelve parameters were observed across a set of 82 double haploid (DH) lines derived from IR64/Azucena. Three nitrogen regimes namely, native (0 kg/ha; no nitrogen applied), optimum (100 kg/ha) and high (200 kg/ha) replicated thrice were the environments. The parents and DH lines were significantly varying for all traits under different nitrogen regimes. All traits except plant height recorded significant genotype x environment interaction. Individual plant yield was positively correlated with nitrogen use efficiency and nitrogen uptake. Sixteen QTLs were detected by composite interval mapping. Eleven QTLs showed significant QTL x environment interactions. On chromosome 3, seven QTLs were detected associated with nitrogen use, plant yield and associated traits. A QTL region between markers RZ678, RZ574 and RZ284 was associated with nitrogen use and yield. This chromosomal region was enriched with expressed gene sequences of known key nitrogen assimilation genes.