Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase.

BACKGROUND: Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1. PURPOSE: Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system. RESULTS: Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification. CONCLUSION: Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

Promoter insertion leads to polyembryony in mango - a case of convergent evolution with citrus.

Sexual reproduction in plants is the main pathway for creating new genetic combinations in modern agriculture. In heterozygous plants, after the identification of a plant with desired traits, vegetative propagation (cloning) is the primary path to create genetically uniform plants. Another natural plant mechanism that creates genetically uniform plants (clones) is apomixis. In fruit crops like citrus and mango, sporophytic apomixis results in polyembryony, where seeds contain multiple embryos, one of which is sexually originated and the others are vegetative clones of the parent mother tree. Utilizing the mango genome and genetic analysis of a diverse germplasm collection, we identified MiRWP as the gene that causes polyembryony in mango. There is a strong correlation between a specific insertion in the gene's promoter region and altered expression in flowers and developing fruitlets, inducing multiple embryos. The MiRWP gene is an ortholog of CitRWP that causes polyembryony in citrus. Based on the data, we speculate that promoter insertion events, which occurred independently in citrus and mango, induced nucellar embryogenesis. The results suggest convergent evolution of polyembryony in the two species. Further work is required to demonstrate the utility of these genes (mango and citrus) in other biological systems as a tool for the clonal production of other crops.

Genotype-by-environment interactions for starch, mineral, and agronomic traits in pearl millet hybrids evaluated across five locations in West Africa.

Introduction: Pearl millet is a staple cereal grown in the harshest environments of arid and semi-arid regions of Asia and sub-Saharan Africa. It is the primary source of calories for millions of people in these regions because it has better adaptation to harsh environmental conditions and better nutritional traits than many other cereals. By screening the pearl millet inbred germplasm association panel (PMiGAP), we earlier reported the best genotypes with the highest concentration of slowly digestible and resistant starch in their grains. Methods: In the current study, we tested these 20 top-performing pearl millet hybrids, identified based on starch data, in a randomised block design with three replications at five locations in West Africa, viz. Sadore and Konni (Niger), Bambey (Senegal), Kano (Nigeria), and Bawku (Ghana). Phenotypic variability was assessed for agronomic traits and mineral traits (Fe and Zn). Results and discussion: Analysis of variance demonstrated significant genotypic, environmental, and GEI effects among five testing environments for agronomic traits (days to 50% flowering, panicle length, and grain yield), starch traits (rapidly digestible starch, slowly digestible starch, resistant starch, and total starch), and mineral trait (iron and zinc). Starch traits, such as rapidly digestible starch (RDS) and slowly digestible starch (SDS), showed nonsignificant genotypic and environmental interactions but high heritability, indicating the lower environmental influence on these traits in the genotype × testing environments. Genotype stability and mean performance across all the traits were estimated by calculating the multi-trait stability index (MTSI), which showed that genotypes G3 (ICMX207070), G8 (ICMX207160), and G13 (ICMX207184) were the best performing and most stable among the five test environments.

Genotype-by-Environment Interaction Analysis of Metabolites in Pearl Millet Genotypes with High Concentrations of Slowly Digestible and Resistant Starch in Their Grains.

Genotype × environment interactions (GEIs) should play an important role in the selection of suitable germplasm in breeding programmes. We here assessed GEI effects on pearl millet (Pennisetum glaucum L.) genotypes, selected to possess a high concentration of slowly digestible starch (SDS) and resistant starch (RS) in their grains. Entries were grown in a randomized complete block design with three replications at locations in Bawku-Ghana, Sadore-Niger, Bamako-Mali, Konni-Nigeria, and Gampella-Burkina Faso across West Africa. Harvested grains from these locations were metabolomically profiled using flow injection ionization-high-resolution mass spectrometry (FIE-HRMS). A total of 3144 mass features (m/z) (1560 negative ion mode and 1584 positive ion mode) were detected, of which, 475 m/z were linked to metabolites be involved in starch, antioxidant and lipid biosynthesis, and vitamin metabolism. Combined ANOVA revealed that the GEI was significantly evident for 54 health-benefiting metabolites, many associated with sugar, especially galactose, metabolism. Additive main effects and multiplicative interaction (AMMI) analysis examined genotype variation and GEI effects, which, when combined with principal component analysis (PCA), found that m/z 171.14864 (positive ionisation, propenyl heptanoate) accounted for 89% of the GEI variation along PC1. The AMMI-based stability parameter (ASTAB), modified AMMI stability value (MASV), and modified AMMI stability index (MASI) were then applied to identify stable and high-performing genotypes for all the health-benefiting metabolites. Similarly, the best-linear-unbiased-prediction (BLUP)-based stability estimation was also performed using the harmonic mean of genotypic values (HMGV), relative performance of genotypic values (RPGV), and harmonic mean of relative performance of genotypic values (HMRPGV), to identify genotype rankings across multiple environments. The multi-trait stability index (MTSI) was calculated and found that the genotypes G1 (ICMH-177111) and G24 (ICMX-207137) were the most stable and were the best mean performers across 52 health-benefiting metabolic traits. These findings demonstrate the potential of G × E assessments on the delivery of health-benefiting metabolite-rich grains in future varieties and hybrids of pearl millet.

Metabolite Diversity and Metabolic Genome-Wide Marker Association Studies (Mgwas) for Health Benefiting Nutritional Traits in Pearl Millet Grains.

As efforts are made to increase food security, millets are gaining increasing importance due to their excellent nutritional credentials. Among the millets, pearl millet is the predominant species possessing several health benefiting nutritional traits in its grain that are helpful in mitigating chronic illnesses such as type-2 diabetes and obesity. In this paper, we conducted metabolomic fingerprinting of 197 pearl millet inbred lines drawn randomly from within the world collection of pearl millet germplasm and report the extent of genetic variation for health benefitting metabolites in these genotypes. Metabolites were extracted from seeds and assessed using flow infusion high-resolution mass spectrometry (FIE-HRMS). Metabolite features (m/z), whose levels significantly differed among the germplasm inbred lines, were identified by ANOVA corrected for FDR and subjected to functional pathway analysis. A number of health-benefiting metabolites linked to dietary starch, antioxidants, vitamins, and lipid metabolism-related compounds were identified. Metabolic genome-wide association analysis (mGWAS) performed using the 396 m/z as phenotypic traits and the 76 K SNP as genotypic variants identified a total of 897 SNPs associated with health benefiting nutritional metabolite at the -log p-value ≤ 4.0. From these associations, 738 probable candidate genes were predicted to have an important role in starch, antioxidants, vitamins, and lipid metabolism. The mGWAS analysis focused on genes involved in starch branching (α-amylase, ß-amylase), vitamin-K reductase, UDP-glucuronosyl, and UDP-glucosyl transferase (UGTs), L-ascorbate oxidase, and isoflavone 2'-monooxygenase genes, which are known to be linked to increases in human health benefiting metabolites. We demonstrate how metabolomic, genomic, and statistical approaches can be utilized to pinpoint genetic variations and their functions linked to key nutritional properties in pearl millet, which in turn can be bred into millets and other cereals crops using plant breeding methods.

Identifying Anti-Oxidant Biosynthesis Genes in Pearl Millet [Pennisetum glaucum (L.) R. Br.] Using Genome-Wide Association Analysis.

Pearl millet [Pennisetum glaucum (L.) R Br.] is an important staple food crop in the semi-arid tropics of Asia and Africa. It is a cereal grain that has the prospect to be used as a substitute for wheat flour for celiac patients. It is an important antioxidant food resource present with a wide range of phenolic compounds that are good sources of natural antioxidants. The present study aimed to identify the total antioxidant content of pearl millet flour and apply it to evaluate the antioxidant activity of its 222 genotypes drawn randomly from the pearl millet inbred germplasm association panel (PMiGAP), a world diversity panel of this crop. The total phenolic content (TPC) significantly correlated with DPPH (1,1-diphenyl-2-picrylhydrazyl) radical scavenging activity (% inhibition), which ranged from 2.32 to 112.45% and ferric-reducing antioxidant power (FRAP) activity ranging from 21.68 to 179.66 (mg ascorbic acid eq./100 g). Genome-wide association studies (GWAS) were conducted using 222 diverse accessions and 67 K SNPs distributed across all the seven pearl millet chromosomes. Approximately, 218 SNPs were found to be strongly associated with DPPH and FRAP activity at high confidence [-log (p) > 3.0-7.4]. Furthermore, flanking regions of significantly associated SNPs were explored for candidate gene harvesting. This identified 18 candidate genes related to antioxidant pathway genes (flavanone 7-O-beta-glycosyltransferase, GDSL esterase/lipase, glutathione S-transferase) residing within or near the association signal that can be selected for further functional characterization. Patterns of genetic variability and the associated genes reported in this study are useful findings, which would need further validation before their utilization in molecular breeding for high antioxidant-containing pearl millet cultivars.

Epigenetics and Epigenomics of Plants.

The genetic material DNA in association with histone proteins forms the complex structure called chromatin, which is prone to undergo modification through certain epigenetic mechanisms including cytosine DNA methylation, histone modifications, and small RNA-mediated methylation. Alterations in chromatin structure lead to inaccessibility of genomic DNA to various regulatory proteins such as transcription factors, which eventually modulates gene expression. Advancements in high-throughput sequencing technologies have provided the opportunity to study the epigenetic mechanisms at genome-wide levels. Epigenomic studies using high-throughput technologies will widen the understanding of mechanisms as well as functions of regulatory pathways in plant genomes, which will further help in manipulating these pathways using genetic and biochemical approaches. This technology could be a potential research tool for displaying the systematic associations of genetic and epigenetic variations, especially in terms of cytosine methylation onto the genomic region in a specific cell or tissue. A comprehensive study of plant populations to correlate genotype to epigenotype and to phenotype, and also the study of methyl quantitative trait loci (QTL) or epiGWAS, is possible by using high-throughput sequencing methods, which will further accelerate molecular breeding programs for crop improvement. Graphical Abstract.

Salinity induced differential methylation patterns in contrasting cultivars of foxtail millet (Setaria italica L.).

KEY MESSAGE: Genome-wide methylation analysis of foxtail millet cultivars contrastingly differing in salinity tolerance revealed DNA demethylation events occurring in tolerant cultivar under salinity stress, eventually modulating the expression of stress-responsive genes. Reduced productivity and significant yield loss are the adverse effects of environmental conditions on physiological and biochemical pathways in crop plants. In this context, understanding the epigenetic machinery underlying the tolerance traits in a naturally stress tolerant crop is imperative. Foxtail millet (Setaria italica) is known for its better tolerance to abiotic stresses compared to other cereal crops. In the present study, methylation-sensitive amplified polymorphism (MSAP) technique was used to quantify the salt-induced methylation changes in two foxtail millet cultivars contrastingly differing in their tolerance levels to salt stress. The study highlighted that the DNA methylation level was significantly reduced in tolerant cultivar compared to sensitive cultivar. A total of 86 polymorphic MSAP fragments were identified, sequenced and functionally annotated. These fragments showed sequence similarity to several genes including ABC transporter, WRKY transcription factor, serine threonine-protein phosphatase, disease resistance, oxidoreductases, cell wall-related enzymes and retrotransposon and transposase like proteins, suggesting salt stress-induced methylation in these genes. Among these, four genes were chosen for expression profiling which showed differential expression pattern between both cultivars of foxtail millet. Altogether, the study infers that salinity stress induces genome-wide DNA demethylation, which in turn, modulates expression of corresponding genes.

Comprehensive analysis of SET domain gene family in foxtail millet identifies the putative role of SiSET14 in abiotic stress tolerance.

SET domain-containing genes catalyse histone lysine methylation, which alters chromatin structure and regulates the transcription of genes that are involved in various developmental and physiological processes. The present study identified 53 SET domain-containing genes in C4 panicoid model, foxtail millet (Setaria italica) and the genes were physically mapped onto nine chromosomes. Phylogenetic and structural analyses classified SiSET proteins into five classes (I-V). RNA-seq derived expression profiling showed that SiSET genes were differentially expressed in four tissues namely, leaf, root, stem and spica. Expression analyses using qRT-PCR was performed for 21 SiSET genes under different abiotic stress and hormonal treatments, which showed differential expression of these genes during late phase of stress and hormonal treatments. Significant upregulation of SiSET gene was observed during cold stress, which has been confirmed by over-expressing a candidate gene, SiSET14 in yeast. Interestingly, hypermethylation was observed in gene body of highly differentially expressed genes, whereas methylation event was completely absent in their transcription start sites. This suggested the occurrence of demethylation events during various abiotic stresses, which enhance the gene expression. Altogether, the present study would serve as a base for further functional characterization of SiSET genes towards understanding their molecular roles in conferring stress tolerance.

Genome-wide SNP identification and characterization in two soybean cultivars with contrasting Mungbean Yellow Mosaic India Virus disease resistance traits.

Mungbean yellow mosaic India virus (MYMIV) is a bipartite Geminivirus, which causes severe yield loss in soybean (Glycine max). Considering this, the present study was conducted to develop large-scale genome-wide single nucleotide polymorphism (SNP) markers and identify potential markers linked with known disease resistance loci for their effective use in genomics-assisted breeding to impart durable MYMIV tolerance. The whole-genome re-sequencing of MYMIV resistant cultivar 'UPSM-534' and susceptible Indian cultivar 'JS-335' was performed to identify high-quality SNPs and InDels (insertion and deletions). Approximately 234 and 255 million of 100-bp paired-end reads were generated from UPSM-534 and JS-335, respectively, which provided ~98% coverage of reference soybean genome. A total of 3083987 SNPs (1559556 in UPSM-534 and 1524431 in JS-335) and 562858 InDels (281958 in UPSM-534 and 280900 in JS-335) were identified. Of these, 1514 SNPs were found to be present in 564 candidate disease resistance genes. Among these, 829 non-synonymous and 671 synonymous SNPs were detected in 266 and 286 defence-related genes, respectively. Noteworthy, a non-synonymous SNP (in chromosome 18, named 18-1861613) at the 149th base-pair of LEUCINE-RICH REPEAT RECEPTOR-LIKE PROTEIN KINASE gene responsible for a G/C transversion [proline (CCC) to alanine(GCC)] was identified and validated in a set of 12 soybean cultivars. Taken together, the present study generated a large-scale genomic resource such as, SNPs and InDels at a genome-wide scale that will facilitate the dissection of various complex traits through construction of high-density linkage maps and fine mapping. In the present scenario, these markers can be effectively used to design high-density SNP arrays for their large-scale validation and high-throughput genotyping in diverse natural and mapping populations, which could accelerate genomics-assisted MYMIV disease resistance breeding in soybean.

Genome-wide development of transposable elements-based markers in foxtail millet and construction of an integrated database.

Transposable elements (TEs) are major components of plant genome and are reported to play significant roles in functional genome diversity and phenotypic variations. Several TEs are highly polymorphic for insert location in the genome and this facilitates development of TE-based markers for various genotyping purposes. Considering this, a genome-wide analysis was performed in the model plant foxtail millet. A total of 30,706 TEs were identified and classified as DNA transposons (24,386), full-length Copia type (1,038), partial or solo Copia type (10,118), full-length Gypsy type (1,570), partial or solo Gypsy type (23,293) and Long- and Short-Interspersed Nuclear Elements (3,659 and 53, respectively). Further, 20,278 TE-based markers were developed, namely Retrotransposon-Based Insertion Polymorphisms (4,801, ∼24%), Inter-Retrotransposon Amplified Polymorphisms (3,239, ∼16%), Repeat Junction Markers (4,451, ∼22%), Repeat Junction-Junction Markers (329, ∼2%), Insertion-Site-Based Polymorphisms (7,401, ∼36%) and Retrotransposon-Microsatellite Amplified Polymorphisms (57, 0.2%). A total of 134 Repeat Junction Markers were screened in 96 accessions of Setaria italica and 3 wild Setaria accessions of which 30 showed polymorphism. Moreover, an open access database for these developed resources was constructed (Foxtail millet Transposable Elements-based Marker Database; http://59.163.192.83/ltrdb/index.html). Taken together, this study would serve as a valuable resource for large-scale genotyping applications in foxtail millet and related grass species.

Identification and molecular characterization of MYB Transcription Factor Superfamily in C4 model plant foxtail millet (Setaria italica L.).

MYB proteins represent one of the largest transcription factor families in plants, playing important roles in diverse developmental and stress-responsive processes. Considering its significance, several genome-wide analyses have been conducted in almost all land plants except foxtail millet. Foxtail millet (Setaria italica L.) is a model crop for investigating systems biology of millets and bioenergy grasses. Further, the crop is also known for its potential abiotic stress-tolerance. In this context, a comprehensive genome-wide survey was conducted and 209 MYB protein-encoding genes were identified in foxtail millet. All 209 S. italica MYB (SiMYB) genes were physically mapped onto nine chromosomes of foxtail millet. Gene duplication study showed that segmental- and tandem-duplication have occurred in genome resulting in expansion of this gene family. The protein domain investigation classified SiMYB proteins into three classes according to number of MYB repeats present. The phylogenetic analysis categorized SiMYBs into ten groups (I-X). SiMYB-based comparative mapping revealed a maximum orthology between foxtail millet and sorghum, followed by maize, rice and Brachypodium. Heat map analysis showed tissue-specific expression pattern of predominant SiMYB genes. Expression profiling of candidate MYB genes against abiotic stresses and hormone treatments using qRT-PCR revealed specific and/or overlapping expression patterns of SiMYBs. Taken together, the present study provides a foundation for evolutionary and functional characterization of MYB TFs in foxtail millet to dissect their functions in response to environmental stimuli.