Association mapping, patterns of linkage disequilibrium and selection in the vicinity of the PHYTOCHROME C gene in pearl millet

Genome-wide SNPs and candidate genes underlying the genetic variations for protein and amino acids in pearl millet (Pennisetum glaucum) germplasm

MAIN CONCLUSION

A total of 544 significant marker-trait associations and 286 candidate genes associated with total protein and 18 amino acids were identified. Thirty-three candidate genes were found near the strong marker trait associations (- log10P ≥ 5.5). Pearl millet (Pennisetum glaucum) is largely grown as a subsistence crop in South Asia and sub-Saharan Africa. It serves as a major source of daily protein intake in these regions. Despite its importance, no systematic effort has been made to study the genetic variations of protein and amino acid content in pearl millet germplasm. The present study was undertaken to dissect the global genetic variations of total protein and 18 essential and non-essential amino acids in pearl millet, using a set of 435 K Single Nucleotide Polymorphisms (SNPs) and 161 genotypes of the Pearl Millet Inbred Germplasm Association Panel (PMiGAP). A total of 544 significant marker-trait associations (at P < 0.0001; - log10P ≥ 4) were detected and 23 strong marker-trait associations were identified using Bonferroni's correction method. Forty-eight pleiotropic loci were found in the genome for the studied traits. In total, 286 candidate genes associated with total protein and 18 amino acids were identified. Thirty-three candidate genes were found near strongly associated SNPs. The associated markers and the candidate genes provide an insight into the genetic architecture of the traits studied and are going to be useful in breeding improved pearl millet varieties in the future. Availabilities of improved pearl millet varieties possessing higher protein and amino acid compositions will help combat the rising malnutrition problem via diet.

Genome-wide association analysis identifies natural allelic variants associated with panicle architecture variation in African rice, Oryza glaberrima Steud

African rice (Oryza glaberrima Steud), a short-day cereal crop closely related to Asian rice (Oryza sativa L.), has been cultivated in Sub-Saharan Africa for ∼ 3,000 years. Although less cultivated globally, it is a valuable genetic resource in creating high-yielding cultivars that are better adapted to diverse biotic and abiotic stresses. While inflorescence architecture, a key trait for rice grain yield improvement, has been extensively studied in Asian rice, the morphological and genetic determinants of this complex trait are less understood in African rice. In this study, using a previously developed association panel of 162 O. glaberrima accessions and new SNP variants characterized through mapping to a new version of the O. glaberrima reference genome, we conducted a genome-wide association study of four major morphological panicle traits. We have found a total of 41 stable genomic regions that are significantly associated with these traits, of which 13 co-localized with previously identified QTLs in O. sativa populations and 28 were unique for this association panel. Additionally, we found a genomic region of interest on chromosome 3 that was associated with the number of spikelets and primary and secondary branches. Within this region was localized the O. sativa ortholog of the PHYTOCHROME B gene (Oglab_006903/OgPHYB). Haplotype analysis revealed the occurrence of natural sequence variants at the OgPHYB locus associated with panicle architecture variation through modulation of the flowering time phenotype, whereas no equivalent alleles were found in O. sativa. The identification in this study of genomic regions specific to O. glaberrima indicates panicle-related intra-specific genetic variation in this species, increasing our understanding of the underlying molecular processes governing panicle architecture. Identified candidate genes and major haplotypes may facilitate the breeding of new African rice cultivars with preferred panicle traits.

Development of a new AgriSeq 4K mid-density SNP genotyping panel and its utility in pearl millet breeding

Pearl millet is a crucial nutrient-rich staple food in Asia and Africa and adapted to the climate of semi-arid topics. Since the genomic resources in pearl millet are very limited, we have developed a brand-new mid-density 4K SNP panel and demonstrated its utility in genetic studies. A set of 4K SNPs were mined from 925 whole-genome sequences through a comprehensive in-silico pipeline. Three hundred and seventy-three genetically diverse pearl millet inbreds were genotyped using the newly-developed 4K SNPs through the AgriSeq Targeted Genotyping by Sequencing technology. The 4K SNPs were uniformly distributed across the pearl millet genome and showed considerable polymorphism information content (0.23), genetic diversity (0.29), expected heterozygosity (0.29), and observed heterozygosity (0.03). The SNP panel successfully differentiated the accessions into two major groups, namely B and R lines, through genetic diversity, PCA, and structure models as per their pedigree. The linkage disequilibrium (LD) analysis showed Chr3 had higher LD regions while Chr1 and Chr2 had more low LD regions. The genetic divergence between the B- and R-line populations was 13%, and within the sub-population variability was 87%. In this experiment, we have mined 4K SNPs and optimized the genotyping protocol through AgriSeq technology for routine use, which is cost-effective, fast, and highly reproducible. The newly developed 4K mid-density SNP panel will be useful in genomics and molecular breeding experiments such as assessing the genetic diversity, trait mapping, backcross breeding, and genomic selection in pearl millet.

Genomic footprints of selection in early-and late-flowering pearl millet landraces

Pearl millet is among the top three-cereal production in one of the most climate vulnerable regions, sub-Saharan Africa. Its Sahelian origin makes it adapted to grow in poor sandy soils under low soil water regimes. Pearl millet is thus considered today as one of the most interesting crops to face the global warming. Flowering time, a trait highly correlated with latitude, is one of the key traits that could be modulated to face future global changes. West African pearl millet landraces, can be grouped into early- (EF) and late-flowering (LF) varieties, each flowering group playing a specific role in the functioning and resilience of Sahelian smallholders. The aim of this study was thus to detect genes linked to flowering but also linked to relevant traits within each flowering group. We thus investigated genomic and phenotypic diversity in 109 pearl millet landrace accessions, i.e., 66 early-flowering and 43 late-flowering, grown in the groundnut basin, the first area of rainfed agriculture in Senegal dominated by dry cereals (millet, maize, and sorghum) and legumes (groundnuts, cowpeas). We were able to confirm the role of PhyC gene in pearl millet flowering and identify several other genes that appear to be as much as important, such as FSR12 and HAC1. HAC1 and two other genes appear to be part of QTLs previously identified and deserve further investigation. At the same time, we were able to highlight a several genes and variants that could contribute to the improvement of pearl millet yield, especially since their impact was demonstrated across flowering cycles.

Genetic insights in pearl millet breeding in the genomic era: challenges and prospects

Pearl millet, a vital staple food and an important cereal, is emerging as crop having various end-uses as feed, food as well as fodder. Advancement in high-throughput sequencing technology has boosted up pearl millet genomic research in past few years. The available draft genome of pearl millet providing an insight into the advancement of several breeding lines. Comparative and functional genomics have untangled several loci and genes regulating adaptive and agronomic traits in pearl millet. Additionally, the knowledge achieved has far away from being applicable in real breeding practices. We believe that the best path ahead is to adopt genome-based approaches for tailored designing of pearl millet as multi-functional crop with outstanding agronomic traits for various end uses. Presently review highlight several novel concepts and techniques in crop breeding, and summarize the recent advances in pearl millet genomic research, peculiarly genome-wide association dissections of several novel alleles and genes for agronomically important traits.

Pearl Millet: A Climate-Resilient Nutricereal for Mitigating Hidden Hunger and Provide Nutritional Security

Pearl millet [Pennisetum glaucum (L.) R. Br.] is the sixth most important cereal crop after rice, wheat, maize, barley and sorghum. It is widely grown on 30 million ha in the arid and semi-arid tropical regions of Asia and Africa, accounting for almost half of the global millet production. Climate change affects crop production by directly influencing biophysical factors such as plant and animal growth along with the various areas associated with food processing and distribution. Assessment of the effects of global climate changes on agriculture can be helpful to anticipate and adapt farming to maximize the agricultural production more effectively. Pearl millet being a climate-resilient crop is important to minimize the adverse effects of climate change and has the potential to increase income and food security of farming communities in arid regions. Pearl millet has a deep root system and can survive in a wide range of ecological conditions under water scarcity. It has high photosynthetic efficiency with an excellent productivity and growth in low nutrient soil conditions and is less reliant on chemical fertilizers. These attributes have made it a crop of choice for cultivation in arid and semi-arid regions of the world; however, fewer efforts have been made to study the climate-resilient features of pearl millet in comparison to the other major cereals. Several hybrids and varieties of pearl millet were developed during the past 50 years in India by both the public and private sectors. Pearl millet is also nutritionally superior and rich in micronutrients such as iron and zinc and can mitigate malnutrition and hidden hunger. Inclusion of minimum standards for micronutrients-grain iron and zinc content in the cultivar release policy-is the first of its kind step taken in pearl millet anywhere in the world, which can lead toward enhanced food and nutritional security. The availability of high-quality whole-genome sequencing and re-sequencing information of several lines may aid genomic dissection of stress tolerance and provide a good opportunity to further exploit the nutritional and climate-resilient attributes of pearl millet. Hence, more efforts should be put into its genetic enhancement and improvement in inheritance to exploit it in a better way. Thus, pearl millet is the next-generation crop holding the potential of nutritional richness and the climate resilience and efforts must be targeted to develop nutritionally dense hybrids/varieties tolerant to drought using different omics approaches.

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.

Genome-Wide Association Studies and Genomic Selection in Pearl Millet: Advances and Prospects

Pearl millet is a climate-resilient, drought-tolerant crop capable of growing in marginal environments of arid and semi-arid regions globally. Pearl millet is a staple food for more than 90 million people living in poverty and can address the triple burden of malnutrition substantially. It remained a neglected crop until the turn of the 21st century, and much emphasis has been placed since then on the development of various genetic and genomic resources for whole-genome scan studies, such as the genome-wide association studies (GWAS) and genomic selection (GS). This was facilitated by the advent of sequencing-based genotyping, such as genotyping-by-sequencing (GBS), RAD-sequencing, and whole-genome re-sequencing (WGRS) in pearl millet. To carry out GWAS and GS, a world association mapping panel called the Pearl Millet inbred Germplasm Association Panel (PMiGAP) was developed at ICRISAT in partnership with Aberystwyth University. This panel consisted of germplasm lines, landraces, and breeding lines from 27 countries and was re-sequenced using the WGRS approach. It has a repository of circa 29 million genome-wide SNPs. PMiGAP has been used to map traits related to drought tolerance, grain Fe and Zn content, nitrogen use efficiency, components of endosperm starch, grain yield, etc. Genomic selection in pearl millet was jump-started recently by WGRS, RAD, and tGBS (tunable genotyping-by-sequencing) approaches for the PMiGAP and hybrid parental lines. Using multi-environment phenotyping of various training populations, initial attempts have been made to develop genomic selection models. This mini review discusses advances and prospects in GWAS and GS for pearl millet.

Single nucleotide polymorphisms in a regulatory site of VRN-A1 first intron are associated with differences in vernalization requirement in winter wheat

Winter wheats require a long exposure to cold temperatures (vernalization) to accelerate flowering. However, varieties differ in the length of the period of cold required to saturate the vernalization response. Here we show that single nucleotide polymorphisms (SNP) at the binding site of the GRP2 protein in the VRN-A1 first intron (henceforth, RIP3) are associated with significant differences in heading time after a partial vernalization treatment. The ancestral winter VRN-A1 allele in ‘Triple Dirk C’ has one SNP in the RIP3 region (1_SNP) relative to the canonical RIP3 sequence, whereas the derived ‘Jagger’ allele has three SNPs (3_SNPs). Both varieties have a single VRN-A1 copy encoding identical proteins. In an F₂ population generated from a cross between these two varieties, plants with the 3_SNPs haplotype headed significantly earlier (P < 0.001) than those with the 1_SNP haplotype, both in the absence of vernalization (17 days difference) and after 3-weeks of vernalization (11 days difference). Plants with the 3_SNPs haplotype showed higher VRN-A1 transcript levels than those with the 1_SNP haplotype. The 3_SNPs haplotype was also associated with early heading in a panel of 127 winter wheat varieties grown in three separate controlled-environment experiments under partial vernalization (36 to 54 days, P < 0.001) and one experiment under field conditions (21 d, P < 0.0001). The RIP3 polymorphisms can be used by wheat breeders to develop winter wheat varieties adapted to regions with different duration or intensity of the cold season.

Myosin XI is associated with fitness and adaptation to aridity in wild pearl millet

Phenotypic changes in plants can be observed along many environmental gradients and are determined by both environmental and genetic factors. The identification of alleles associated with phenotypic variations is a rapidly developing area of research. We studied the genetic basis of phenotypic variations in 11 populations of wild pearl millet (Pennisetum glaucum) on two North-South aridity gradients, one in Niger and one in Mali. Most of the 11 phenotypic traits assessed in a common garden experiment varied between the populations studied. Moreover, the size of the inflorescence, the number of flowers and aboveground dry mass co-varied positively with a decrease in rainfall. To decipher the genetic basis of these phenotypes, we used an association mapping strategy with a mixed model. We found two SNPs on the same myosin XI contig significantly associated with variations in the average number of flowers. Both the allele frequency of the two SNPs and the average number of flowers co-varied with the rainfall gradient on the two gradients. Interestingly, this gene was also a target of selection during domestication. The Myosin XI gene is thus a good candidate for fitness-related adaptation in wild populations.

Deciphering Genomic Regions for High Grain Iron and Zinc Content Using Association Mapping in Pearl Millet

Micronutrient malnutrition, especially deficiency of two mineral elements, iron [Fe] and zinc [Zn] in the developing world needs urgent attention. Pearl millet is one of the best crops with many nutritional properties and is accessible to the poor. We report findings of the first attempt to mine favorable alleles for grain iron and zinc content through association mapping in pearl millet. An association mapping panel of 130 diverse lines was evaluated at Delhi, Jodhpur and Dharwad, representing all the three pearl millet growing agro-climatic zones of India, during 2014 and 2015. Wide range of variation was observed for grain iron (32.3-111.9 ppm) and zinc (26.6-73.7 ppm) content. Genotyping with 114 representative polymorphic SSRs revealed 0.35 mean gene diversity. STRUCTURE analysis revealed presence of three sub-populations which was further supported by Neighbor-Joining method of clustering and principal coordinate analysis (PCoA). Marker-trait associations (MTAs) were analyzed with 267 markers (250 SSRs and 17 genic markers) in both general linear model (GLM) and mixed linear model (MLM), however, MTAs resulting from MLM were considered for more robustness of the associations. After appropriate Bonferroni correction, Xpsmp 2261 (13.34% R²-value), Xipes 0180 (R²-value of 11.40%) and Xipes 0096 (R²-value of 11.38%) were consistently associated with grain iron and zinc content for all the three locations. Favorable alleles and promising lines were identified for across and specific environments. PPMI 1102 had highest number (7) of favorable alleles, followed by four each for PPMFeZMP 199 and PPMI 708 for across the environment performance for both grain Fe and Zn content, while PPMI 1104 had alleles specific to Dharwad for grain Fe and Zn content. When compared with the reference genome Tift 23D₂B₁-P1-P5, Xpsmp 2261 amplicon was identified in intergenic region on pseudomolecule 5, while the other marker, Xipes 0810 was observed to be overlapping with aspartic proteinase (Asp) gene on pseudomolecule 3. Thus, this study can help in breeding new lines with enhanced micronutrient content using marker-assisted selection (MAS) in pearl millet leading to improved well-being especially for women and children.

Spatial and Temporal Variation in Selection of Genes Associated with Pearl Millet Varietal Quantitative Traits In situ

Ongoing global climate changes imply new challenges for agriculture. Whether plants and crops can adapt to such rapid changes is still a widely debated question. We previously showed adaptation in the form of earlier flowering in pearl millet at the scale of a whole country over three decades. However, this analysis did not deal with variability of year to year selection. To understand and possibly manage plant and crop adaptation, we need more knowledge of how selection acts in situ. Is selection gradual, abrupt, and does it vary in space and over time? In the present study, we tracked the evolution of allele frequency in two genes associated with pearl millet phenotypic variation in situ. We sampled 17 populations of cultivated pearl millet over a period of 2 years. We tracked changes in allele frequencies in these populations by genotyping more than seven thousand individuals. We demonstrate that several allele frequencies changes are compatible with selection, by correcting allele frequency changes associated with genetic drift. We found marked variation in allele frequencies from year to year, suggesting a variable selection effect in space and over time. We estimated the strength of selection associated with variations in allele frequency. Our results suggest that the polymorphism maintained at the genes we studied is partially explained by the spatial and temporal variability of selection. In response to environmental changes, traditional pearl millet varieties could rapidly adapt thanks to this available functional variability.

Inference of domestication history and differentiation between early- and late-flowering varieties in pearl millet

Pearl millet (Pennisetum glaucum) is a staple crop in Sahelian Africa. Farmers usually grow varieties with different cycle lengths and complementary functions in Sahelian agrosystems. Both the level of genetic differentiation of these varieties and the domestication history of pearl millet have been poorly studied. We investigated the neutral genetic diversity and population genetic structure of early- and late-flowering domesticated and wild pearl millet populations using 18 microsatellite loci and 8 nucleotide sequences. Strikingly, early- and late-flowering domesticated varieties were not differentiated over their whole distribution area, despite a clear difference in their isolation-by-distance pattern. Conversely, our data brought evidence for two well-differentiated genetic pools in wild pearl millet, allowing us to test scenarios with different numbers and origins of domestication using approximate Bayesian computation (ABC). The ABC analysis showed the likely existence of asymmetric migration between wild and domesticated populations. The model choice procedure indicated that a single domestication from the eastern wild populations was the more likely scenario to explain the polymorphism patterns observed in cultivated pearl millet.

From phenotypes to causal sequences: using genome wide association studies to dissect the sequence basis for variation of plant development

Tremendous natural variation of growth and development exists within species. Uncovering the molecular mechanisms that tune growth and development promises to shed light on a broad set of biological issues including genotype to phenotype relations, regulatory mechanisms of biological processes and evolutionary questions. Recent progress in sequencing and data processing capabilities has enabled Genome Wide Association Studies (GWASs) to identify DNA sequence polymorphisms that underlie the variation of biological traits. In the last years, GWASs have proven powerful in revealing the complex genetic bases of many phenotypes in various plant species. Here we highlight successful recent GWASs that uncovered mechanistic and sequence bases of trait variation related to plant growth and development and discuss important considerations for conducting successful GWASs.

PHYTOCHROME C is an essential light receptor for photoperiodic flowering in the temperate grass, Brachypodium distachyon

We show that in the temperate grass, Brachypodium distachyon, PHYTOCHROME C (PHYC), is necessary for photoperiodic flowering. In loss-of-function phyC mutants, flowering is extremely delayed in inductive photoperiods. PHYC was identified as the causative locus by utilizing a mapping by sequencing pipeline (Cloudmap) optimized for identification of induced mutations in Brachypodium. In phyC mutants the expression of Brachypodium homologs of key flowering time genes in the photoperiod pathway such as GIGANTEA (GI), PHOTOPERIOD 1 (PPD1/PRR37), CONSTANS (CO), and florigen/FT are greatly attenuated. PHYC also controls the day-length dependence of leaf size as the effect of day length on leaf size is abolished in phyC mutants. The control of genes upstream of florigen production by PHYC was likely to have been a key feature of the evolution of a long-day flowering response in temperate pooid grasses.

Phytochrome C plays a major role in the acceleration of wheat flowering under long-day photoperiod

Phytochromes are dimeric proteins that function as red and far-red light sensors influencing nearly every phase of the plant life cycle. Of the three major phytochrome families found in flowering plants, phytochrome C (PHYC) is the least understood. In Arabidopsis and rice, PHYC is unstable and functionally inactive unless it heterodimerizes with another phytochrome. However, when expressed in an Arabidopsis phy-null mutant, wheat PHYC forms signaling active homodimers that translocate into the nucleus in red light to mediate photomorphogenic responses. Tetraploid wheat plants homozygous for loss-of-function mutations in all PHYC copies (phyC(AB)) flower on average 108 d later than wild-type plants under long days but only 19 d later under short days, indicating a strong interaction between PHYC and photoperiod. This interaction is further supported by the drastic down-regulation in the phyC(AB) mutant of the central photoperiod gene photoperiod 1 (PPD1) and its downstream target flowering locus T1, which are required for the promotion of flowering under long days. These results implicate light-dependent, PHYC-mediated activation of PPD1 expression in the acceleration of wheat flowering under inductive long days. Plants homozygous for the phyC(AB) mutations also show altered profiles of circadian clock and clock-output genes, which may also contribute to the observed differences in heading time. Our results highlight important differences in the photoperiod pathways of the temperate grasses with those of well-studied model plant species.