- 1Agricultural Research Center-Hays, Kansas State University, Hays, KS, USA
- 2Institute of Biological, Environmental and Rural Sciences, Aberystwyth University, Aberystwyth, UK
Pearl millet [Penisetum glaucum (L) R. Br.] is a hardy cereal crop grown in the arid and semiarid tropics where other cereals are likely to fail to produce economic yields due to drought and heat stresses. Adaptive evolution, a form of natural selection shaped the crop to grow and yield satisfactorily with limited moisture supply or under periodic water deficits in the soil. Drought tolerance is a complex polygenic trait that various morphological and physiological responses are controlled by 100s of genes and significantly influenced by the environment. The development of genomic tools will have enormous potential to improve the efficiency and precision of conventional breeding. The apparent independent domestication events, highly outcrossing nature and traditional cultivation in stressful environments maintained tremendous amount of polymorphism in pearl millet. This high polymorphism of the crop has been revealed by genome mapping that in turn stimulated the mapping and tagging of genomic regions controlling important traits such as drought tolerance. Mapping of a major QTL for terminal drought tolerance in independent populations envisaged the prospect for the development of molecular breeding in pearl millet. To accelerate genetic gains for drought tolerance targeted novel approaches such as establishment of marker-trait associations, genomic selection tools, genome sequence and genotyping-by-sequencing are still limited. Development and application of high throughput genomic tools need to be intensified to improve the breeding efficiency of pearl millet to minimize the impact of climate change on its production.
Introduction
The unfolding climate change; a rise in global average temperature caused mostly by increasing concentrations of greenhouse gasses in the atmosphere is expected to increasingly distress crop production in many parts of the world. An estimated 41% of the land area of the world constitutes arid and semiarid zones characterized by a very stochastic environment (Safriel et al., 2005), suboptimal for crop production due to drought. A crop species or genotype that is resilient to low rainfall amount and erratic distribution and elevated temperature would have a vital importance for sustainable food supply to the ever increasing world population. Developing a cultivar for such suboptimal environment through a carefully crafted breeding strategy is inevitable to sustainable food security in inhospitable climates.
Pearl millet [Penisetum glaucum (L) R. Br.] is a hardy cereal crop grown in the adverse climatic conditions of the arid and semiarid tropics where other cereals are likely to fail to produce economic yields. It is the staple food and fodder crop of millions of rural families in the hottest and driest environments of Africa and Indian subcontinent. As its cultivation is predominantly practiced on marginal lands under unirrigated condition, environmental stresses mainly drought and heat are the two most important production constraints that it faces.
A better understanding of the effects of drought on pearl millet and mechanisms of plant tolerance are vital for breeding new drought tolerant cultivars. However, complex mechanisms are involved in drought tolerance of crops and progress in developing drought tolerant varieties of pearl millet and other orphan crops have been slow (Bidinger et al., 1987; Howarth and Yadav, 2002; Yadav et al., 2002). Finding key functional components or machineries of drought tolerance would have utmost importance in overcoming this longstanding problem in breeding for drought tolerance. Application of molecular tools is proven to be quite important in dissecting complex quantitative traits such as drought tolerance into contributing traits, each controlled by QTL that follow Mendelian inheritance. With the advancements in next-generation sequencing (NGS) and molecular profiling technologies that is increasingly becoming more affordable and efficient to record DNA sequence variations over 1000s of genomic loci, identification and tagging of genes underling a trait in the genome have become possible.
In this review, the progress of development and application of genomic tools in breeding for drought tolerance in pearl millet are discussed. Some important morphological and physiological traits that can help in indirect selection, tools for molecular breeding approaches, and different approaches attempted are reviewed. The gaps in relation to breeding for drought tolerance, future prospects with respect to the application of genomic tools and potential platforms that can be used to accelerate the improvement of pearl millet for drought tolerance are examined.
General Effects of Drought on Plant Growth and Development
Water is required by plants for transpiration, surface cooling and carbon assimilation and drought imposes an adverse effect on plant growth and development and the effect is manifold in nature. One of the key processes affected by water deficit in plants is photosynthesis due to decreased CO2 diffusion to the chloroplast (Pinheiro and Chaves, 2011) as stomata is closed to conserve available water in the plant system. Reduced photosynthesis can also be worsened by heat stress aggravated by water deficit (Crafts-Brandner and Salvucci, 2002). Decreased water in the plant system also damages chloroplast and affects photosynthetic outputs due to reduced light interception (Stone, 2001). Furthermore, the drought accompanying increased air temperature results in reduced reproduction due to heat damage on the pollen grain (Bita and Gerats, 2013), and increased sterility (Schoper et al., 1987; Jagadish et al., 2012).
Drought stress also prompt plants to accumulate metabolites such as proline that plays important role as protein precursors in plant metabolism and development. Proline accumulated in the plant cell serves as an excellent osmolyte, as a metal chelator, an anti-oxidative defense molecule, and a signaling molecule (Hayat et al., 2012). It is known to impart stress tolerance by maintaining cell turgor or osmotic balance, preventing electrolyte leakage by stabilizing membranes and proteins such as RUBISCO (Ribulose 1,5-bisphosphate carboxylase) (Hayat et al., 2012) and mitochondrial electron transport complex II (Hamilton and Heckathorn, 2001), and controlling reactive oxygen species (ROS) concentrations in plants (Paleg et al., 1984; Hare et al., 1998). Drought screening and molecular and biochemical evaluations conducted so far in pearl millet never considered the level of proline accumulation and its correlation with drought tolerance of genotypes. Some other metabolic changes associated with water deficit in the plant system include increased concentration of total soluble sugars, reducing sugars, sucrose, and amylase activity (Comino et al., 1997) which result in decreased starch grains. In general, reduction or inhibition of enzymatic activity, inhibition of cell expansion due to cell wall hardening to maintain turgor through lignification (Fan and Neumann, 2004), metabolic changes, ionic imbalance, disturbances in solute accumulation or a combination of two or more of these mechanisms are some of the effects of water shortfall in the plant system.
Pearl Millet As A Drought Tolerant Cereal
Drought is the major constraint to pearl millet as it is grown in the drier semiarid and arid regions. However, adaptive evolution and natural selection made pearl millet relatively the most drought and heat tolerant among other cereals. Traditional landraces from drier regions are good sources for breeding drought tolerance (Kusaka et al., 2005; Yadav, 2010). Though, the yield potential of many of the landraces is lower as compared to improved cultivars, genes for drought tolerance in the landraces can be utilized to combine with high yield potential (Yadav et al., 2011).
Nevertheless, drought tolerance is a complex polygenic trait and influenced significantly by environments. The various morphological and physiological responses to drought are controlled by 100s of genes (Hu and Xiong, 2014). Plant stress responses involves adjustment of metabolism for physiological and morphological adaptation occurring as a result of dynamic and complex cross-talk between different regulatory networks (Saito and Matsuda, 2010). Several technical limitations prohibit the study of plant responses to environmental stress that appear to be complex and pose difficulty on the classical plant-breeding program. Apparently, the elucidation of the components of the complex mechanisms underlying drought tolerance in crops will accelerate breeding for drought tolerant cultivars.
The mechanism how plants endure the effect of drought on their growth and development takes various forms. The ability of plants to grow and yield satisfactorily with limited moisture supply or under periodic soil water deficits is termed as ‘drought resistance.’ It is principally the constitutive plant traits of maintaining high plant water status (dehydration avoidance) that enables plants to survive in water limited environments (Blum, 2005). This mechanism has virtuous relevance in terminal drought tolerance of pearl millet in which tolerant genotypes fill their grain and offset any drastic effect on their yield. Since phenotyping dehydration avoidance is difficult, no much study is devoted to understand the genetic bases of drought resistance in pearl millet and other related crops.
The other mechanism is the ability of plant to sense environmental factors and trigger development shift to mature before water stress becomes a serious limiting factor, which is termed as drought escape. The initiation of flowering in higher plants is a complex trait regulated by an integrated gene network that responds to the dynamics of both autonomous and environmental cues (Andres and Coupland, 2012). Flowering time is among the most common adaptive transitions (Zuellig et al., 2014) where ecologically relevant adaptive shifts taking place in plant development. Few recent studies shade light on the genetic factors underlying a shift in vegetative to reproductive stage and cycle length in pearl millet. A study on MADS-box gene family identified a polymorphism at a gene named PgMADS11 was associated with flowering time variation in pearl millet (Mariac et al., 2011). The early shift from the vegetative to the reproductive stage warrant the survival of the species through successful completion of the life cycle before moisture deficit becomes worse, but with significant negative effect on the yield performance.
On the other hand, ability of a plant to withstand water deficit as measured by degree and duration of low plant water potential is termed as drought avoidance. This is characterized by water saving and water spending physiological and morphological mechanisms. Both root and shoot architecture are important in the drought avoidance mechanisms of plants. The water savers are plants which conserve water by closing of stomata when they sense water deficit in the soil. The water spenders have specialized root system to extract more water from deeper soil horizon or by adjusting their osmotic balance. Pearl millet root system is characterized by fast primary root development at early growth stage and quickly colonizing deeper soil horizon (Passot et al., 2016) that enables water acquisition from deeper root zone.
Plants that can recover from a dry period by producing new leaves or tillers from buds and able to survive drought by slowing total water loss are regarded as drought tolerant. In a highly drought prone environments pearl millet farmers specifically grow high tillering landraces (van Oosterom et al., 1996) that can grow secondary tiller when the main culm is affected by drought. These high tillering landraces usually have small panicle and small individual grain sizes (Aparna et al., 2014) to minimize the impairment of grain filling in drought conditions. As pearl millet grain yield is highly correlated with grain number (Bidinger and Raju, 2000), and an adaptive trait to the arid environment (DeWet et al., 1992), preference for high tillering landraces with small grain sizes that takes short period to fill the grain is a wise drought coping strategy helpful in combating extreme aridity.
Morpho-Physiological Drought Adaptive Traits
Through evolution, the sessile plants evolved morpho-physiological traits that help them adapt to the dry environment. The morpho-physiological traits of importance in drought tolerance include both shoot and root characteristics naturally modified for efficient water absorption from the soil and conservation in the system. The morpho-physiological traits of drought tolerance effects relevant for yield are stomatal conductance, photosynthetic capacity, timing of phenological phases, starch availability during ovary/embryo development, stem reserve mobilization, stay green, single plant leaf area, rooting depth and density, cuticular resistance and surface roughness, osmotic adjustment, membrane composition, antioxidative defense, and accumulation of stress-related proteins (Cattivelli et al., 2008). More open stomata allows greater conductance and lower leaf temperature, a proxy for higher potential photosynthetic and transpiration rates. However, genotypes with higher stomatal conductance are relatively sensitive to drought stress as they lose water rapidly than those close their stomata and reduce their photosynthetic CO2 assimilation.
Increased mobilization of stem reserved assimilates to compensate for reduced current leaf photosynthesis in drought stress (Slafer et al., 2005), reduced leaf area through leaf rolling (Walter and Schurr, 2005), and leaf surface roughness that reduce water loss and induce greater light reflectance (Kerstiens, 1996) are some of the traits that have vital importance for drought stress tolerance.
Stay-green is a characteristic of some genotypes to extend the duration of active photosynthesis by delaying their leaves senescence for a longer period of time than other genotypes. Stay-green phenotype is induced by a complex signaling network that involves mutations suppressing ethylene, abscisic acid, brassinosteroid, and strigolactone signal transduction, or those activating cytokinin signaling otherwise regulating leaf senescence (Kusaba et al., 2013). The semi-dwarf pearl millet inbred lines developed in USA for hybrid breeding have stay-green characteristics. However, a detailed study of the genes involved or the region of the genome where such genes are residing is not conducted. A plant that maintains the integrity of its chlorophyll despite the soil water status during reproductive stage can continue photosynthesis and can fill the grain and maintain a reasonable grain yield performance in drought stress.
Root anatomy, architecture, depth, and density are also important traits in drought tolerance. Detailed characterization of pearl millet root system may have a significant contribution in pearl millet breeding for drought tolerance. Small fine root diameters and root length density are some of the root traits known to contributing to productivity under drought condition (Comas et al., 2013). Such type of genotypes are potential terminal drought tolerant for they can exploit the water available deeper in the soil profile when the upper root zone is dried. Neither the genes underlying such traits nor the genomic regions hosting the QTL for such root characteristics were investigated in pearl millet. Further, any genetic variation in the germplasm for such traits and its correlation with terminal drought tolerance was not studied yet.
History of Genomic Tools Development and Application in Pearl Millet
The application of genomic tools has become an essential component of breeding programs and proven to be useful for identifying novel genes for traits of agronomic importance and stress tolerance. Marker assisted-selection (MAS) and gene introgression into desirable genetic backgrounds have a great potential to facilitate crops improvement by cutting down lengthy phenotypic evaluation and selection. Alternative molecular plant breeding strategies have been developed for different crops based on the advancement in population genotyping (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/220903/fpls-07-01724-HTML/image_m/fpls-07-01724-g001.jpg)
FIGURE 1. Conceptual flow of molecular breeding using advanced genomic tools that can be applicable to pearl millet (F2, second filial generation; BC, backcross population; BILs, backcross inbred lines; RILs, recombinant inbred lines; DH, doubled haploid; AM, association mapping population; NAM, nested association mapping population; MAGIC, multi-parent advanced generation intercrossing population; GEBV, genetic and breeding value).
The development of molecular markers for the improvement of pearl millet was conceived in 1991 by DFID-JIC-ICRISAT project (Gale et al., 2005). Consequently, the first genetic linkage map of pearl millet using restriction fragment-length polymorphisms (RFLPs) was constructed in Liu et al. (1994). Seven linkage groups corresponding to the haploid chromosome number of the species were formed.
The straight forward diploid genetics and high levels of polymorphism in pearl millet prompted several subsequent studies including mapping and tagging of genes controlling important traits in the pearl millet genome. The prominent work was the mapping of QTL to understand the genetic and physiological basis of terminal drought tolerance in pearl millet (Yadav et al., 2002). This has provided a more-targeted approach to improving the drought tolerance and yield of this crop in water-limited environments (Howarth and Yadav, 2002). A major QTL for terminal drought tolerance that explains more than 30% of grain yield variability in severe drought conditions was mapped on LG2 (linkage group two) using two independent mapping populations (Yadav et al., 2002; Bidinger et al., 2007).
Later, three major QTL for grain yield with low QTL × environment interactions were identified across a range of post-flowering moisture environments (Bidinger et al., 2007). The effect of this major QTL explaining more than up to 32% of the variations under drought was validated in two independent marker-assisted backcross programs (Serraj et al., 2005) in which 30% improvement in general combining ability for grain yield expected from this QTL under terminal drought stress was recovered in introgression lines (Yadav et al., 2011). Being predominantly climate resilient crops, millets may serve as valuable sources of novel genes for stress tolerance, which need to be identified and characterized. The close phylogenetic relationships between millets and other cereals could enable the introgression of novel alleles, genes or QTL identified in millets for better agronomic traits into other cereals toward ensuring food security under changing climate.
The major drought tolerance QTL identified by Yadav et al. (2002) has currently been found to be associated with reduced salt uptake (Sharma et al., 2014). A QTL for low transpiration rates that contributes to terminal drought tolerance through conserving soil moisture for grain filling by minimizing water lose during vegetative stage was co-mapped with the terminal drought tolerance QTL (Kholová et al., 2012). This co-mapping of the transpiration rate (a proxy trait to canopy conductance) and drought tolerance indicates that water conservation at a vegetative stage is important for the post-flowering stage performance of pearl millet. Reduced transpiration rate is an interplay of physiological traits such as reduced stomatal density and morphological traits of the canopy that plays important role in the rate of canopy conductance.
Development of PCR-based markers that can be used readily in diversity assessment and MAS programs have been conducted at different times. One way of developing SSRs was isolation from BACs through targeted use of 3′ end-anchored SSR primers and flanking sequences suppression PCR (Qi et al., 2001). This way of developing markers enables anchoring individual BACs to the genetic maps that in turn facilitate the construction of BAC contigs. EST database mining with the aim of developing PCR-based codominant markers found EST sequences containing di- and trinucleotide repeats at a density of 1.75 kb EST sequences (Senthilvel et al., 2008). Some of these EST-SSRs markers were found polymorphic among several set of parental lines used to develop mapping population for pearl millet genetic linkage mapping (Rajaram et al., 2013).
To integrate gene-based markers in the pearl millet genetic map, SNP markers developed from genes with known functions in plant adaptation to drought and associated abiotic stresses were used to identify candidate genes underlying major drought-tolerance QTL (Sehgal et al., 2012). Accordingly, genes coding for transferases, kinases, oxidases, and other proteins mostly associated with QTL of grain yield, flowering time and leaf rolling in drought stress conditions were co-mapped with the drought-tolerance QTL. This demonstrated the importance of functional markers in candidate gene identification and functional comparison of distantly related genomes in the study of drought tolerance and other traits. The candidate gene approach is preferred for QTL characterization than map-based cloning and insertional mutagenesis (Pflieger et al., 2001; Zhu and Zhao, 2007) as involvement of several genes with minor effect on the phenotype affect the efficiency of the later. In the identification of genes responsible for quantitative traits, candidate-gene approach using linkage maps based on functional gene markers has become essential (Sehgal et al., 2012). This approach of mapping functional gene markers in the QTL region was proven practical in putative candidate gene identification for drought tolerance in rice (Nguyen et al., 2004), barley (Tondelli et al., 2006), and durum wheat (Diab et al., 2008).
Development of high throughput markers has become quick and inexpensive through the development of NGS technology. Discovery of high throughput markers enabled genotyping of larger mapping population and construction of high and ultrahigh density maps in crops with no reference genome sequence available (Poland et al., 2012). As compared to other orphan crops, a sizable progress has been made for pearl millet in the development of molecular tools that are proven efficient in expediting the breeding of improved cultivars in the last two decades (Vadez et al., 2012). Development of molecular markers, EST libraries, BACs, and genome mapping were some of the areas in which progress has been observed. Nevertheless, to harness the wealth of genetic variability in the germplasm for improvement in yield, drought tolerance, and other desired traits much is remaining to be done (Vadez et al., 2012).
Second-Generation Platforms For Molecular Breeding
The technological advances in NGS (Davey et al., 2011) also brought about a shift in linkage mapping population from the bi-parental to the second-generation multi-parent populations (Bohra, 2013) such as nested association mapping population (NAM) (Yu et al., 2008), multi-parent advanced generation intercross (MAGIC) population (Kover et al., 2009) to exploit multi-parent allelic diversity to a maximum advantage in fine mapping of QTL for desirable traits using an integrated mapping population. Use of second-generation populations with cutting-edge computing resources rendered a new opportunities in utilizing the advantages of both the linkage and association mapping for crop improvement (Bevan and Uauy, 2013).
The nested-association mapping utilizes the advantages of both linkage analysis and association mapping and eliminates disadvantages of both by taking into consideration recent and historical recombination events and facilitate high resolution mapping using an integrated multi-parent population (Yu et al., 2008). It has been utilized to identify functional markers or polymorphic sites that provide information about the genetic architecture underlying complex quantitative traits variability in many plant species (Andersen and Lübberstedt, 2003).
The NAM population should be developed in such a manner that it incorporates diverse parents crossed to a common parent (Figure 2). The information about the development and utilization of NAM in pearl millet to decipher allelic variations underlying agronomic and yield traits as well as tolerance to stresses including drought is lacking. As pearl millet exhibits high diversity, novel alleles for important traits can be captured through the utilization of such multi-parent population. As the cultivated pearl millet is classified into four races (typhoides, nigritarium, globosum, and leonis) (Brunken et al., 1977), second-generation mapping populations such as NAM, MAGIC and germplasm population needs to follow a comprehensive approach to incorporate such morphological variations and ethnographic isolates.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/220903/fpls-07-01724-HTML/image_m/fpls-07-01724-g002.jpg)
FIGURE 2. Schematic diagram of pearl millet NAM population development. The genome of each founder is color coded to show that the genome of RIL is a mosaic of founders and common parent genome segments (adapted from Yu et al., 2008). High-density genotyping of founders permits the linkage disequilibrium information captured in these diverse founders to be exploited for high-resolution mapping (SSD, single seed descent).
A germplasm based population [pearl millet inbred germplasm association panel (PMiGAP)] has been developed in pearl millet to capture wider genetic variations of traits (Sehgal et al., 2015). A structure analysis on this systematically developed association panel from ICRISAT global collection revealed six subpopulations corresponding predominantly to pedigree or similar agronomic traits (Sehgal et al., 2015). This result evidenced the importance of including different morphological forms in a molecular study for genome-wide marker discovery, QTL detection, or estimate of genomic based breeding value. Both within and between races assessment would likely help to discover novel allelic variants and extensive development of genomic tools for marker-assisted breeding.
Multi-parent advanced generation intercross (MAGIC) population was proposed for genome-wide association mapping for the identification of major genes, QTL, and potentially novel loci associated with important traits (Kover et al., 2009; Bandillo et al., 2013). The analytical methodology involves reconstructing the genome of each line developed as a mosaic of the founders used to develop the population. It is expected to improve the genetic dissection of complex traits and precision of QTL mapping (Kover et al., 2009) for MAS and map-based cloning of candidate genes. The fixed multi-parent population in pearl millet will have advantages in assessing larger genetic variation in the germplasm than association analysis on the sample of naturally occurring haplotypes (Rakshit et al., 2012). Across species, a standard of 1000 fixed individuals developed from intercrossing of eight founders is advocated as the most appropriate for MAGIC (Figure 3). However, MAGIC has disadvantage mainly in resource utilization as it takes longer time and more resources to develop. The population is also likely to show extensive segregation for developmental traits, limiting their use in the analysis of complex traits.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/220903/fpls-07-01724-HTML/image_m/fpls-07-01724-g003.jpg)
FIGURE 3. Proposed crossing scheme for the development of multi-parent advanced generation inter-cross (MAGIC) population development for pearl millet. The scheme involves half-diallel mating of eight molecularly diverse inbred germplasm, intercrossing of the recombinants, and repeated selfing to develop multi-parent recombinant inbred lines.
The MAS uses molecular markers in linkage disequilibrium with major QTL that are tightly linked with the QTL and inherited together. The recent extension of marker-assisted selection is genomic selection (Meuwissen et al., 2001; Goddard and Hayes, 2007; Heffner et al., 2009) where genome-wide molecular markers are targeted in a large breeding population to improve quantitative traits. In genomic selection promising individuals are selected based on genomics estimated breeding values. It combines marker genotype data with phenotype from several environments and pedigree or kinship to estimate the breeding values. This method of selection based on genetic values predicted from molecular markers in the whole genome was anticipated to substantially speed up genetic gain than that achieved through selection based on major QTL (Meuwissen et al., 2001).
Genomic selection has been used for improvement of quantitative traits in several different crops with effective training population design. Being an orphan crop, pearl millet has not been evaluated yet for the selection strategy using breeding values. It was preferred to major QTL to accurately predict the genetic value of breeding lines using genome-wide high-throughput DNA markers information (Heffner et al., 2009). The major advantage of genomic selection lies in its potential to allow selection decisions during the off-season, thereby increasing genetic gain on an annual basis (Heffner et al., 2010). Development and utilization of genomic selection in pearl millet breeding can be applied in accurate prediction of the hybrid performance with ideal resource allocation. As has been surveyed in different crops (Zhao et al., 2015) it slightly improved prediction accuracies of hybrid performance when dominance effects were exploited.
NGS Based Genotyping Platforms
Genotyping-by-sequencing (GBS) is the most widely used low-cost whole-genome marker profiling platform for a robust simultaneous marker discovery and genotyping (Elshire et al., 2011). This inexpensive highly multiplex protocol has been used in many crops for constructing libraries for NGS of genomes. A large number of high throughput single nucleotide polymorphisms (SNPs) markers were discovered and used to study within-species diversity and population structure, construct genome maps, and conduct genome-wide association studies (GWAS) (Metzker, 2010). Few works that have been conducted on pearl millet using this marker discovery and genotyping platform discovered important findings that encourage the application of the technology. High SNP polymorphism was observed between pairs of parental lines in gene based DNA sequences in pearl millet (Sehgal et al., 2012) indicating high allelic variation between the genotypes. Further, these markers are distributed across all the seven linkage groups that reflect the likely phenomenon of the distribution of the markers in the genome. Fine mapping of the DT-QTL on LG2 detected several markers mapped to the region that are associated with traits having prime importance in drought tolerance mainly delayed leaf senescence and leaf rolling under drought stress (Sehgal et al., 2012). GBS enabled construction of genetic map with high marker density and more uniform distribution than previously developed maps (Moumouni et al., 2015). As constructing a genetic map of a species is a key step in identifying and positioning a gene or QTL underlying a trait of interest (Pedraza-Garcia et al., 2010), marker density and uniform distribution play important role in the process of mapping and marker-trait association.
A recent molecular diversity and population structure analysis of 500 pearl millet accessions in two diversity panels (252 of global accessions from Africa, Asia, and the Americas; and 248 landraces from multiple agro-ecological zones in Senegal) using GBS discovered 83,875 SNPs that revealed the prevalence of high genetic diversity in accessions from Senegal as compared to the world collections (Hu et al., 2015). Rapid LD decay and lack of confounding population structure along agro-ecological zones in Senegal is predicted to simplify future pearl millet association mapping studies.
A genotyping method that allows an unbiased investigation of all exons or the complete protein-coding portion of the genome is exome sequencing. As only a small fraction of the genome is coding for protein (1–2%) (Warr et al., 2015), exome sequencing minimizes the sequencing space and associated cost. Either with de novo transcriptome assemblies for the development of array-based (Bi et al., 2012) or the solution-based whole-exome sequencing would have a practical application in pearl millet molecular breeding. The array-based is currently applicable to the pearl millet as no genome sequence available yet. One of the advantages of exome-sequencing commonly called exome-capture is more targeted development of SNP markers associated with variants of a gene controlling a phenotype of interest. Dealing with a small portion of the genome by excluding the repetitive parts also circumvents the complexity of data analysis and allow to handling larger numbers of samples than with other genotyping methods covering the whole-genome sequencing. It also enables to directly map mutations responsible for phenotypes of interest through direct sequencing of the exons.
Pearl Millet Genome Sequence
The pearl millet genome sequence provides a genetic blueprint of the species and presumed to facilitate the development of genomic tools that would expedite the development of improved cultivars through the application of marker-assisted breeding. It also facilitate marker-trait association studies and identification of candidate genes for important quantitative traits such as yield and tolerance to stresses like drought. With the overall objective to sequence the pearl millet genome, ICRISAT spearheaded an informal International Pearl Millet Genome Sequencing Consortium (IPMGSC) that brought together several leading research organizations in genome sequencing, crop genomics, and pearl millet research1. The project of sequencing the pearl millet genome selected an important ancestral genotype of many seed parents of pearl millet hybrids in use for forage and grain in the USA as well as dual-purpose hybrids in India, Tift23D2B1.
The project successfully assembled large clean sequence data from whole genome shotgun sequencing (WGS) raw data (Varshney, 2014). It is expected that the first version of genome sequence will be published soon and apparently brings a paradigm shift to pearl millet molecular breeding by providing a foundation for accelerating gene mapping and characterization. The genome sequence needs to be complemented with transcriptome analysis to predict gene functions, their position in the genome, and expression patterns in different tissues, developmental stages and stress factors.
In pearl millet, there were several independent domestication events resulted in four different races (Brunken et al., 1977), enriching the genome sequence through resequencing of the populations used in the QTL mapping and in entirety of independently sequencing a genotype from other race thus becomes highly desirable. To conduct inclusive characterization of the genetic factors responsible for the seed morphological variations used as a basis for the classification, integrating the known races into the genome sequence would broaden the application of the resource.
Future Prospects
The progress toward developing genomic resources for the application of NGS in pearl millet is promising but much awaits to implement genomics-assisted breeding to expedite the breeding for drought tolerance. With the development of such genomic tools, identification as well as development of improved genotypes is destined to become easier and faster. High density maps construction for fine mapping and candidate gene identification, integration with the genome sequence, and QTL detection needs to be strengthened. A shift from bi-parental to multi-parental and germplasm based association mapping platforms with the application of genomic tools can accelerate discovery of allelic variants that fast pace pearl millet improvement for drought tolerance.
Cost effective high throughput genotyping platforms are readily applicable in a diploid pearl millet genome. However, the respectable genome size should be taken into consideration in the process of selection for a genotyping method. The newly popularized exome-capture may reduce the extent of DNA sequence data handled and also leads to the identification of functional markers. Further, the advanced marker-assisted selection strategy that captures small-effect quantitative trait loci and breeding value needs to be given due emphasis for hybrid breeding. A release of the pending genome sequence would have paramount importance in the development and application of the aforementioned genomic tools for mapping and QTL detection. Integration of genome sequences with the consensus linkage map would increase the identification of genes underlying the QTL detected for drought tolerance.
Author Contributions
DS conceived the idea and prepared the draft of the manuscript. RY edited the draft manuscript and contributed relevant information included in the review.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
References
Andersen, J. R., and Lübberstedt, T. (2003). Functional markers in plants. Trends Plant Sci. 8, 554–560. doi: 10.1016/j.tplants.2003.09.010
Andres, F., and Coupland, G. (2012). The genetic basis of flowering responses to seasonal cues. Nat. Rev. Genet. 13, 627–639. doi: 10.1038/nrg3291
Aparna, K., Hash, C. T., Yadav, R. S., and Vadez, V. (2014). Seed number and 100-seed weight of pearl millet (Pennisetum glaucum L.) respond differently to low soil moisture in genotypes contrasting for drought tolerance. J. Agron. Crop Sci. 200, 119–131. doi: 10.1111/jac.12052
Bandillo, N., Raghavan, C., Muyco, P. A., Sevilla, M. A. L., Lobina, I. T., Dilla-Ermita, C. J., et al. (2013). Multi-parent advanced generation inter-cross (MAGIC) populations in rice: progress and potential for genetics research and breeding. Rice (N. Y). 6:11. doi: 10.1186/1939-8433-6-11
Bevan, M. W., and Uauy, C. (2013). Genomics reveals new landscapes for crop improvement. Genome Biol. 14:206. doi: 10.1186/gb-2013-14-6-206
Bi, K., Vanderpool, D., Singhal, S., Linderoth, T., Moritz, C., and Good, J. M. (2012). Transcriptome-based exon capture enables highly cost-effective comparative genomic data collection at moderate evolutionary scales. BMC Genomics 13:403. doi: 10.1186/1471-2164-13-403
Bidinger, F. R., Mahalakshmi, V., and Rao, G. D. P. (1987). Assessment of drought resistance in pearl millet [Pennisetum americanum (L.) Leeke]. 1. Factors affecting yield under stress. Aust. J. Agric. Res. 38, 37–48.
Bidinger, F. R., Nepolean, T., Hash, C. T., Yadav, R. S., and Howarth, C. J. (2007). Quantitative trait loci for grain yield in pearl millet under variable postflowering moisture conditions. Crop Sci. 47, 969–980. doi: 10.2135/cropsci2006.07.0465
Bidinger, F. R., and Raju, D. S. (2000). Response to selection for increased individual grain mass in pearl millet. Crop Sci. 40, 68–71. doi: 10.2135/cropsci2000.40168x
Bita, C. E., and Gerats, T. (2013). Plant tolerance to high temperature in a changing environment: scientific fundamentals and production of heat stress-tolerant crops. Front. Plant Sci. 4:273. doi: 10.3389/fpls.2013.00273
Blum, A. (2005). Drought resistance, water-use efficiency, and yield potential – Are they compatible, dissonant, or mutually exclusive? Aust. J. Agric. Res. 56, 1159–1168. doi: 10.1071/AR05069
Bohra, A. (2013). Emerging paradigms in genomics-based crop improvement. Sci. World J. 2013:585467. doi: 10.1155/2013/585467
Brunken, J., de Wet, J. M. J., and Harlan, J. R. (1977). The morphology and domestication of pearl millet. Econ. Bot. 31, 163–174. doi: 10.1007/BF02866587
Cattivelli, L., Rizza, F., Badeck, F. W., Mazzucotelli, E., Mastrangelo, A. M., Francia, E., et al. (2008). Drought tolerance improvement in crop plants: an integrated view from breeding to genomics. Field Crop Res. 105, 1–14. doi: 10.1016/j.fcr.2007.07.004
Comas, L. H., Becker, S. R., Cruz, V. M. V., Byrne, P. F., and Dierig, D. A. (2013). Root traits contributing to plant productivity under drought. Front. Plant Sci. 4:442. doi: 10.3389/fpls.2013.00442
Comino, M. L., de Felipe, M. R., Fernández-Pascual, M., and Martin, L. (1997). “Effect of drought stress on carbohydrate metabolism in nodules of Lupinus angustifolius,” in Eukaryotism and Symbiosis, eds H. E. A. Schenk, R. G. Herrmann, K. W. Jeon, N. E. Müller, and W. Schwemmler (Berlin: Springer), 449–456.
Crafts-Brandner, S. J., and Salvucci, M. E. (2002). Sensitivity of photosynthesis in a C4 plant, maize, to heat stress. Plant Physiol. 129, 1773–1780. doi: 10.1104/pp.002170
Davey, J. W., Hohenlohe, P. A., Etter, P. D., Boone, J. Q., Catchen, J. M., and Blaxter, M. L. (2011). Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat. Rev. Genet. 12, 499–510. doi: 10.1038/nrg3012
DeWet, J. M., Bidinger, F. R., and Peacock, J. M. (1992). “Pearl millet (Pennisetum glaucum) a cereal of the Sahel,” in Desertified Grasslands: Their Biology and Management, ed. G. P. Chappman (London: Academic Press), 259–267.
Diab, A. A., Kantety, R. V., Ozturk, N. Z., Benscher, D., Nachit, M. M., and Sorrells, M. E. (2008). Drought – inducible genes and differentially expressed sequence tags associated with components of drought tolerance in durum wheat. Sci. Res. Essay 3, 9–26.
Elshire, R. J., Glaubitz, J. C., Sun, Q., Poland, J. A., Kawamoto, K., Buckler, E. S., et al. (2011). A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6:e19379. doi: 10.1371/journal.pone.0019379
Fan, L., and Neumann, P. M. (2004). The spatially variable inhibition by water deficit of maize root growth correlates with altered profiles of proton flux and cell wall pH. Growth (Lakeland) 135, 2291–2300. doi: 10.1104/pp.104.041426.1
Gale, M. D., Devos, K. M., Zhu, J. H., Allouis, S., Couchman, M. S., Liu, H., et al. (2005). New molecular marker technologies for pearl millet improvement. SAT eJournal 1, 1–7. doi: 10.1186/1471-2229-8-119
Goddard, M. E., and Hayes, B. J. (2007). Genomic selection. J. Anim. Breed. Genet. 124, 323–330. doi: 10.1111/j.1439-0388.2007.00702.x
Hamilton, E. W., and Heckathorn, S. A. (2001). Mitochondrial adaptations to NaCl. Complex I is protected by anti-oxidants and small heat shock proteins, whereas complex II is protected by proline and betaine. Plant Physiol. 126, 1266–1274. doi: 10.1104/pp.126.3.1266
Hare, P. D., Cress, W. A., and Van Staden, J. (1998). Dissecting the roles of osmolyte accumulation during stress. Plant Cell Environ. 21, 535–553. doi: 10.1046/j.1365-3040.1998.00309.x
Hayat, S., Hayat, Q., Alyemeni, M. N., Wani, A. S., Pichtel, J., and Ahmad, A. (2012). Role of proline under changing environments: a review. Plant Signal. Behav. 7, 1456–1466. doi: 10.4161/psb.21949
Heffner, E. L., Lorenz, A. J., Jannink, J. L., and Sorrells, M. E. (2010). Plant breeding with genomic selection: gain per unit time and cost. Crop Sci. 50, 1681–1690. doi: 10.2135/cropsci2009.11.0662
Heffner, E. L., Sorrells, M. E., and Jannink, J. (2009). Genomic selection for crop improvement. Crop Sci. 49, 1–12. doi: 10.2135/cropsci2008.08.0512
Howarth, C. J., and Yadav, R. S. (2002). Successful marker assisted selection for drought tolerance and disease resistance in pearl millet. IGER Innov. 6, 18–21.
Hu, H., and Xiong, L. (2014). Genetic engineering and breeding of drought-resistant crops. Annu. Rev. Plant Biol. 65, 715–741. doi: 10.1146/annurev-arplant-050213-040000
Hu, Z., Mbacké, B., Perumal, R., Guèye, M. C., Sy, O., Bouchet, S., et al. (2015). Population genomics of pearl millet (Pennisetum glaucum (L.) R. Br.): comparative analysis of global accessions and Senegalese landraces. BMC Genomics 16:1048. doi: 10.1186/s12864-015-2255-0
Jagadish, S. V. K., Septiningsih, E. M., Kohli, A., Thomson, M. J., Ye, C., Redona, E., et al. (2012). Genetic advances in adapting rice to a rapidly changing climate. J. Agron. Crop Sci. 198, 360–373. doi: 10.1111/j.1439-037X.2012.00525.x
Kerstiens, G. (1996). Cuticular water permeability and its physiological significance. J. Exp. Bot. 47, 1813–1832. doi: 10.1093/jxb/47.12.1813
Kholová, J., Nepolean, T., Tom Hash, C., Supriya, A., Rajaram, V., Senthilvel, S., et al. (2012). Water saving traits co-map with a major terminal drought tolerance quantitative trait locus in pearl millet [Pennisetum glaucum (L.) R. Br.]. Mol. Breed. 30, 1337–1353. doi: 10.1007/s11032-012-9720-0
Kover, P. X., Valdar, W., Trakalo, J., Scarcelli, N., Ehrenreich, I. M., Purugganan, M. D., et al. (2009). A multiparent advanced generation inter-cross to fine-map quantitative traits in Arabidopsis thaliana. PLoS Genet. 5:e1000551. doi: 10.1371/journal.pgen.1000551
Kusaba, M., Tanaka, A., and Tanaka, R. (2013). Stay-green plants: what do they tell us about the molecular mechanism of leaf senescence. Photosynth. Res. 117, 221–234. doi: 10.1007/s11120-013-9862-x
Kusaka, M., Lalusin, A. G., and Fujimura, T. (2005). The maintenance of growth and turgor in pearl millet (Pennisetum glaucum [L.] Leeke) cultivars with different root structures and osmo-regulation under drought stress. Plant Sci. 168, 1–14. doi: 10.1016/j.plantsci.2004.06.021
Liu, C. J., Witcombe, J. R., Pittaway, T. S., Nash, M., Hash, C. T., Busso, C. S., et al. (1994). An RFLP-based genetic map of pearl millet (Pennisetum glaucum). Theor. Appl. Genet. 89, 481–487. doi: 10.1007/BF00225384
Mariac, C., Jehin, L., Saïdou, A. A., Thuillet, A. C., Couderc, M., Sire, P., et al. (2011). Genetic basis of pearl millet adaptation along an environmental gradient investigated by a combination of genome scan and association mapping. Mol. Ecol. 20, 80–91. doi: 10.1111/j.1365-294X.2010.04893.x
Metzker, M. L. (2010). Sequencing technologies – the next generation. Nat. Rev. Genet. 11, 31–46. doi: 10.1038/nrg2626
Meuwissen, T. H. E., Hayes, B. J., and Goddard, M. E. (2001). Prediction of total genetic value using genome-wide dense marker maps. Genetics 157, 1819–1829.
Moumouni, K. H., Kountche, B. A., Jean, M., Hash, C. T., Vigouroux, Y., Haussmann, B. I. G., et al. (2015). Construction of a genetic map for pearl millet, Pennisetum glaucum (L.) R. Br., using a genotyping-by-sequencing (GBS) approach. Mol. Breed. 35, 1–10. doi: 10.1007/s11032-015-0212-x
Nguyen, T. T. T., Klueva, N., Chamareck, V., Aarti, A., Magpantay, G., Millena, A. C. M., et al. (2004). Saturation mapping of QTL regions and identification of putative candidate genes for drought tolerance in rice. Mol. Genet. Genomics 272, 35–46. doi: 10.1007/s00438-004-1025-5
Paleg, L. G., Stewart, G. R., and Bradbeer, J. W. (1984). Proline and glycine betaine influence protein solvation. Plant Physiol. 75, 974–978. doi: 10.1104/pp.75.4.974
Passot, S., Gnacko, F., Moukouanga, D., Lucas, M., Guyomarc’h, S., Ortega, B. M., et al. (2016). Characterization of pearl millet root architecture and anatomy reveals three types of lateral roots. Front. Plant Sci. 7:829. doi: 10.3389/fpls.2016.00829
Pedraza-Garcia, F., Specht, J. E., and Dweikat, I. (2010). A new PCR-based linkage map in pearl millet. Crop Sci. 50, 1754–1760. doi: 10.2135/cropsci2009.10.0560
Pflieger, S., Lefebvre, V., and Causse, M. (2001). The candidate gene approach in plant genetics: a review. Mol. Breed. 7, 275–291. doi: 10.1023/A:1011605013259
Pinheiro, C., and Chaves, M. M. (2011). Photosynthesis and drought: can we make metabolic connections from available data? J. Exp. Bot. 62, 869–882. doi: 10.1093/jxb/erq340
Poland, J. A., Brown, P. J., Sorrells, M. E., and Jannink, J. L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. PLoS ONE 7:e32253. doi: 10.1371/journal.pone.0032253
Qi, X., Lindup, S., Pittaway, T. S., Allouis, S., Gale, M. D., and Devos, K. M. (2001). Development of simple sequence repeat markers from bacterial artificial chromosomes without subcloning. Biotechniques 31, 355–362.
Rajaram, V., Nepolean, T., Senthilvel, S., Varshney, R. K., Vadez, V., Srivastava, R. K., et al. (2013). Pearl millet [Pennisetum glaucum (L.) R. Br.] consensus linkage map constructed using four RIL mapping populations and newly developed EST-SSRs. BMC Genomics 14:159. doi: 10.1186/1471-2164-14-159
Rakshit, S., Rakshit, A., and Patil, J. V. (2012). Multiparent intercross populations in analysis of quantitative traits. J. Genet. 91, 111–117. doi: 10.1007/s12041-012-0144-8
Safriel, U., Adeel, Z., Niemeijer, D., Puigdefabregas, J., White, R., Lal, R., et al. (2005). “Dryland systems,” in Millennium Ecosystem Assessment: Ecosystems and Human Well-Being: Current State and Trends: Findings of the Condition and Trends Working Group, Vol. 1, eds R. Hassan, R. Scholes, and N. Ash (Washington, DC: Island Press), 623–662.
Saito, K., and Matsuda, F. (2010). Metabolomics for functional genomics, systems biology, and biotechnology. Annu. Rev. Plant Biol. 61, 463–489. doi: 10.1146/annurev.arplant.043008.092035
Schoper, J. B., Lambert, R. J., Vasilas, B. L., and Westgate, M. E. (1987). Plant factors controlling seed set in maize: the influence of silk, pollen, and ear-leaf water status and tassel heat treatment at pollination. Plant Physiol. 83, 121–125. doi: 10.1104/pp.83.1.121
Sehgal, D., Rajaram, V., Armstead, I. P., Vadez, V., Yadav, Y. P., Hash, C. T., et al. (2012). Integration of gene-based markers in a pearl millet genetic map for identification of candidate genes underlying drought tolerance quantitative trait loci. BMC Plant Biol. 12:9. doi: 10.1186/1471-2229-12-9
Sehgal, D., Skot, L., Singh, R., Srivastava, R. K., Das, S. P., Taunk, J., et al. (2015). Exploring potential of pearl millet germplasm association panel for association mapping of drought tolerance traits. PLoS ONE 10:e0122165. doi: 10.1371/journal.pone.0122165
Senthilvel, S., Jayashree, B., Mahalakshmi, V., Kumar, P. S., Nakka, S., Nepolean, T., et al. (2008). Development and mapping of simple sequence repeat markers for pearl millet from data mining of expressed sequence tags. BMC Plant Biol. 8:119. doi: 10.1186/1471-2229-8-119
Serraj, R., Tom Hash, C., Rizvi, S. M. H., Sharma, A., Yadav, R. S., and Bidinger, F. R. (2005). Recent advances in marker-assisted selection for drought tolerance in pearl millet. Plant Prod. Sci. 8, 334–337. doi: 10.1626/pps.8.334
Sharma, P. C., Singh, D., Sehgal, D., Singh, G., Hash, C. T., and Yadav, R. S. (2014). Further evidence that a terminal drought tolerance QTL of pearl millet is associated with reduced salt uptake. Environ. Exp. Bot. 102, 48–57. doi: 10.1016/j.envexpbot.2014.01.013
Slafer, G. A., Araus, J. L., Royo, C., and García Del Moral, L. F. (2005). Promising eco-physiological traits for genetic improvement of cereal yields in Mediterranean environments. Ann. Appl. Biol. 146, 61–70. doi: 10.1111/j.1744-7348.2005.04048.x
Stone, P. (2001). “The effects of heat stress on cereal yield and quality,” in Crop Responses and Adaptations to Temperature Stress, ed. A. S. Basra (Binghamton, NY: Food Products Press), 243–291.
Tondelli, A., Francia, E., Barabaschi, D., Aprile, A., Skinner, J. S., Stockinger, E. J., et al. (2006). Mapping regulatory genes as candidates for cold and drought stress tolerance in barley. Theor. Appl. Genet. 112, 445–454. doi: 10.1007/s00122-005-0144-7
Vadez, V., Hash, T., Bidinger, F. R., and Kholova, J. (2012). Phenotyping pearl millet for adaptation to drought. Front. Physiol. 3:386. doi: 10.3389/fphys.2012.00386
van Oosterom, E. J., Whitaker, M. L., and Weltzien, R. E. (1996). “Integrating genotype by environment interaction analysis, characterization of drought patterns, and farmer preferences to identify adaptive plant traits for pearl millet,” in Plant Adaptation and Crop Improvement, eds M. Cooper and G. L. Hammer (Wallingford: CAB Inernational), 383–402.
Varshney, R. (2014). “Towards developing a reference genome sequence of pearl millet (Pennisetum glaucum L.),” in Proceedings of the Plant and Animal Genome Conference XXII Poster Abstract, San Diego, CA.
Walter, A., and Schurr, U. (2005). Dynamics of leaf and root growth: endogenous control versus environmental impact. Ann. Bot. 95, 891–900. doi: 10.1093/aob/mci103
Warr, A., Robert, C., Hume, D., Archibald, A., Deeb, N., and Watson, M. (2015). Exome sequencing: current and future perspectives. G3 (Bethesda) 5, 1543–1550. doi: 10.1534/g3.115.018564
Yadav, O. P. (2010). Drought response of pearl millet landrace-based populations and their crosses with elite composites. Field Crop Res. 118, 51–56. doi: 10.1016/j.fcr.2010.04.005
Yadav, R. S., Hash, C. T., Bidinger, F. R., Cavan, G. P., and Howarth, C. J. (2002). Quantitative trait loci associated with traits determining grain and stover yield in pearl millet under terminal drought-stress conditions. Theor. Appl. Genet. 104, 67–83. doi: 10.1007/s001220200008
Yadav, R. S., Sehgal, D., and Vadez, V. (2011). Using genetic mapping and genomics approaches in understanding and improving drought tolerance in pearl millet. J. Exp. Bot. 62, 397–408. doi: 10.1093/jxb/erq265
Yu, J., Holland, J. B., McMullen, M. D., and Buckler, E. S. (2008). Genetic design and statistical power of nested association mapping in maize. Genetics 178, 539–551. doi: 10.1534/genetics.107.074245
Zhao, Y., Mette, M. F., and Reif, J. C. (2015). Genomic selection in hybrid breeding. Plant Breed. 134, 1–10. doi: 10.1111/pbr.12231
Zhu, M., and Zhao, S. (2007). Candidate gene identification approach: progress and challenges. Int. J. Biol. Sci. 3, 420–427. doi: 10.7150/ijbs.3.420
Keywords: pearl millet, Pennisetum glaucum, genomic tools, drought tolerance, GBS, QTL
Citation: Serba DD and Yadav RS (2016) Genomic Tools in Pearl Millet Breeding for Drought Tolerance: Status and Prospects. Front. Plant Sci. 7:1724. doi: 10.3389/fpls.2016.01724
Received: 05 August 2016; Accepted: 02 November 2016;
Published: 22 November 2016.
Edited by:
Xiaowu Wang, Biotechnology Research Institute – Chinese Academy of Agricultural Sciences (CAAS), ChinaReviewed by:
Xiyin Wang, North China Institute of Science and Technology, ChinaErli Pang, Beijing Normal University, China
Xianmin Diao, Institute of Crop Sciences – Chinese Academy of Agricultural Sciences (CAAS), China
Copyright © 2016 Serba and Yadav. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Desalegn D. Serba, ddserba@ksu.edu