Association Studies Identify Natural Variation at PHYC Linked to Flowering Time and Morphological Variation in Pearl Millet

The identification of genes selected during and after plant domestication is an important research topic to enhance knowledge on adaptative evolution. Adaptation to different climates was a key factor in the spread of domesticated crops. We conducted a study to identify genes responsible for these adaptations in pearl millet and developed an association framework to identify genetic variations associated with the phenotype in this species. A set of 90 inbred lines genotyped using microsatellite loci and AFLP markers was used. The population structure was assessed using two different Bayesian approaches that allow inbreeding or not. Association studies were performed using a linear mixed model considering both the population structure and familial relationships between inbred lines. We assessed the ability of the method to limit the number of false positive associations on the basis of the two different Bayesian methods, the number of populations considered and different morphological traits while also assessing the power of the methodology to detect given additive effects. Finally, we applied this methodology to a set of eight pearl millet genes homologous to cereal flowering pathway genes. We found significant associations between several polymorphisms of the pearl millet PHYC gene and flowering time, spike length, and stem diameter in the inbred line panel. To validate this association, we performed a second association analysis in a different set of pearl millet individuals from Niger. We confirmed a significant association between genetic variation in this gene and these characters.DOMESTICATION and dispersion of cultivated plants were associated with their adaptation to the agricultural environment. These adaptations led to genetic changes shared by all individuals of a cultivated species (domestication genes) or to variations between varieties within a cultivated species (genes controlling varietal differences). Domestication genes like tb1 (Doebley et al. 1997; Wang et al. 1999) in maize (Zea mays) were selected very early by human populations (Jaenicke-Després et al. 2003). After the first early selection, adaptation of the flowering phenotype to different climatic conditions was certainly a key innovation that enabled colonization of new environments. One of the most well-known examples was the adaptation of maize—a tropical plant—to northern climates. Maize cultivation spread late to northeastern America. By 1000 YBP, only maize was an established staple crop (Fritz 1995). A genetic variant of the Dwarf8 gene led to an earlier flowering phenotype (Thornsberry et al. 2001). This early allele was present at a high frequency in North America and was certainly selected after the domestication of maize under northern climatic conditions (Camus-Kulandaivelu et al. 2006).Pearl millet (Pennisetum glaucum [(L.) R. Br.]), one of the most important West African cereals, was most likely domesticated once in the Sahelian zone of West Africa (Oumar et al. 2008). By 3500 YBP, it was already being cultivated throughout Sahelian and tropical West African countries (D''Andrea et al. 2001; D''Andrea and Casey 2002). The adaptation of pearl millet in West Africa was also associated with an environmental gradient (Haussmann et al. 2006). Pearl millet varieties from tropical coastal West Africa flower very late (up to 160 days from planting to female flowering) as compared to varieties from Sahelian West Africa, which may have a flowering time as short as 45 days (Haussmann et al. 2006). The genetic factors underlying the differences between these varieties are still unknown.Association studies offer new opportunities for assessing the role of a particular gene on a phenotype. Contrary to QTL analysis, association studies have the challenging task of taking an unknown evolutionary history of studied individuals into account. For example, population structure is a common confounding effect in association studies (Pritchard et al. 2000a). Allele frequencies evolve between divergent structured populations via drift, mutation, and selection. Differences in allele frequencies may be correlated with any morphological traits that differentiate two populations. Then a statistical correlation between a gene and a trait is not necessarily associated with a “causative” relationship between the gene and the morphology, which can lead to a high number of false positives. The use of population structure to correct the number of false positives was a significant breakthrough in plant studies (Thornsberry et al. 2001). This approach was recently further refined by also using a matrix of kinship coefficients, which proves efficient when there is a complex structure and familial relationship between individuals (Yu et al. 2006; Kang et al. 2008; Stich et al. 2008). Complex structures and familial relationships are common in inbred cultivated crop material. In the current association study framework (Thornsberry et al. 2001; Yu et al. 2006; Casa et al. 2008; Kang et al. 2008; Stich et al. 2008), population structure was assessed using STRUCTURE software (Pritchard et al. 2000b). This tool is not implemented to deal with selfed inbred materials or inbred species (Pritchard et al. 2000b). Through new methodological developments, population structure analysis can now be performed using Bayesian methods in these particular cases (Gao et al. 2007). The extent to which the power of association studies will differ when dealing with inbred material or selfing species using either Bayesian method has yet to be evaluated.In this study, we developed an association framework for pearl millet to assess the role of flowering pathway genes. We assessed the ability of the method to control the number of false positives, while taking different methodological inferences of population structure that allow inbreeding or not into account. We also assessed the power of the association framework to detect given additive genetic effects. Finally, we applied this method to a set of eight flowering time gene homologs sequenced in pearl millet. We assessed sequence variation in light perception genes (PHYA, PHYB, PHYC, and CRY2) and downstream regulators of flowering (GI, Hd6, Hd1, and FLORICAULA). Variation was detected in the PHYC gene associated with variations in flowering time and morphological traits. This association was noted in two different data sets.     

Selection for earlier flowering crop associated with climatic variations in the Sahel

Climate changes will have an impact on food production and will require costly adaptive responses. Adapting to a changing environment will be particularly challenging in sub-Saharan Africa where climate change is expected to have a major impact. However, one important phenomenon that is often overlooked and is poorly documented is the ability of agro-systems to rapidly adapt to environmental variations. Such an adaptation could proceed by the adoption of new varieties or by the adaptation of varieties to a changing environment. In this study, we analyzed these two processes in one of the driest agro-ecosystems in Africa, the Sahel. We performed a detailed study in Niger where pearl millet is the main crop and covers 65% of the cultivated area. To assess how the agro-system is responding to recent recurrent drought, we analyzed samples of pearl millet landraces collected in the same villages in 1976 and 2003 throughout the entire cultivated area of Niger. We studied phenological and morphological differences in the 1976 and 2003 collections by comparing them over three cropping seasons in a common garden experiment. We found no major changes in the main cultivated varieties or in their genetic diversity. However, we observed a significant shift in adaptive traits. Compared to the 1976 samples, samples collected in 2003 displayed a shorter lifecycle, and a reduction in plant and spike size. We also found that an early flowering allele at the PHYC locus increased in frequency between 1976 and 2003. The increase exceeded the effect of drift and sampling, suggesting a direct effect of selection for earliness on this gene. We conclude that recurrent drought can lead to selection for earlier flowering in a major Sahelian crop. Surprisingly, these results suggest that diffusion of crop varieties is not the main driver of short term adaptation to climatic variation.     

Etiology and Factors Associated with Pneumonia in Children under 5 Years of Age in Mali: A Prospective Case-Control Study

Background

There are very limited data on children with pneumonia in Mali. The objective was to assess the etiology and factors associated with community-acquired pneumonia in hospitalized children <5 years of age in Mali.

Methods

A prospective hospital-based case-control study was implemented in the Pediatric department of Gabriel Touré University Hospital at Bamako, Mali, between July 2011-December 2012. Cases were children with radiologically-confirmed pneumonia; Controls were hospitalized children without respiratory features, matched for age and period. Respiratory specimens, were collected to identify 19 viruses and 5 bacteria. Whole blood was collected from cases only. Factors associated with pneumonia were assessed by multivariate logistic regression.

Results

Overall, 118 cases and 98 controls were analyzed; 44.1% were female, median age was 11 months. Among pneumonia cases, 30.5% were hypoxemic at admission, mortality was 4.2%. Pneumonia cases differed from the controls regarding clinical signs and symptoms but not in terms of past medical history. Multivariate analysis of nasal swab findings disclosed that S. pneumoniae (adjusted odds ratio [aOR] = 3.4, 95% confidence interval [95% CI]: 1.6–7.0), human metapneumovirus (aOR = 17.2, 95% CI: 2.0–151.4), respiratory syncytial virus [RSV] (aOR = 7.4, 95% CI: 2.3–23.3), and influenza A virus (aOR = 10.7, 95% CI: 1.0–112.2) were associated with pneumonia, independently of patient age, gender, period, and other pathogens. Distribution of S. pneumoniae and RSV differed by season with higher rates of S. pneumoniae in January-June and of RSV in July-September. Pneumococcal serotypes 1 and 5 were more frequent in pneumonia cases than in the controls (P = 0.009, and P = 0.04, respectively).

Conclusions

In this non-PCV population from Mali, pneumonia in children was mainly attributed to S. pneumoniae, RSV, human metapneumovirus, and influenza A virus. Increased pneumococcal conjugate vaccine coverage in children could significantly reduce the burden of pneumonia in sub-Saharan African countries.     

Enteric Bacterial Pathogens in Children with Diarrhea in Niger: Diversity and Antimicrobial Resistance

Background

Although rotavirus is the leading cause of severe diarrhea among children in sub-Saharan Africa, better knowledge of circulating enteric pathogenic bacteria and their antimicrobial resistance is crucial for prevention and treatment strategies.

Methodology/Principal Findings

As a part of rotavirus gastroenteritis surveillance in Maradi, Niger, we performed stool culture on a sub-population of children under 5 with moderate-to-severe diarrhea between April 2010 and March 2012. Campylobacter, Shigella and Salmonella were sought with conventional culture and biochemical methods. Shigella and Salmonella were serotyped by slide agglutination. Enteropathogenic Escherichia coli (EPEC) were screened by slide agglutination with EPEC O-typing antisera and confirmed by detection of virulence genes. Antimicrobial susceptibility was determined by disk diffusion. We enrolled 4020 children, including 230 with bloody diarrhea. At least one pathogenic bacterium was found in 28.0% of children with watery diarrhea and 42.2% with bloody diarrhea. Mixed infections were found in 10.3% of children. EPEC, Salmonella and Campylobacter spp. were similarly frequent in children with watery diarrhea (11.1%, 9.2% and 11.4% respectively) and Shigella spp. were the most frequent among children with bloody diarrhea (22.1%). The most frequent Shigella serogroup was S. flexneri (69/122, 56.5%). The most frequent Salmonella serotypes were Typhimurimum (71/355, 20.0%), Enteritidis (56/355, 15.8%) and Corvallis (46/355, 13.0%). The majority of putative EPEC isolates was confirmed to be EPEC (90/111, 81.1%). More than half of all Enterobacteriaceae were resistant to amoxicillin and co-trimoxazole. Around 13% (46/360) Salmonella exhibited an extended-spectrum beta-lactamase phenotype.

Conclusions

This study provides updated information on enteric bacteria diversity and antibiotic resistance in the Sahel region, where such data are scarce. Whether they are or not the causative agent of diarrhea, bacterial infections and their antibiotic resistance profiles should be closely monitored in countries like Niger where childhood malnutrition pre-disposes to severe and invasive infections.     

Preventing Acute Malnutrition among Young Children in Crises: A Prospective Intervention Study in Niger

Background

Finding the most appropriate strategy for the prevention of moderate acute malnutrition (MAM) and severe acute malnutrition (SAM) in young children is essential in countries like Niger with annual “hunger gaps.” Options for large-scale prevention include distribution of supplementary foods, such as fortified-blended foods or lipid-based nutrient supplements (LNSs) with or without household support (cash or food transfer). To date, there has been no direct controlled comparison between these strategies leading to debate concerning their effectiveness. We compared the effectiveness of seven preventive strategies—including distribution of nutritious supplementary foods, with or without additional household support (family food ration or cash transfer), and cash transfer only—on the incidence of SAM and MAM among children aged 6–23 months over a 5-month period, partly overlapping the hunger gap, in Maradi region, Niger. We hypothesized that distributions of supplementary foods would more effectively reduce the incidence of acute malnutrition than distributions of household support by cash transfer.

Methods and Findings

We conducted a prospective intervention study in 48 rural villages located within 15 km of a health center supported by Forum Santé Niger (FORSANI)/Médecins Sans Frontières in Madarounfa. Seven groups of villages (five to 11 villages) were allocated to different strategies of monthly distributions targeting households including at least one child measuring 60 cm–80 cm (at any time during the study period whatever their nutritional status): three groups received high-quantity LNS (HQ-LNS) or medium-quantity LNS (MQ-LNS) or Super Cereal Plus (SC+) with cash (€38/month [US$52/month]); one group received SC+ and family food ration; two groups received HQ-LNS or SC+ only; one group received cash only (€43/month [US$59/month]). Children 60 cm–80 cm of participating households were assessed at each monthly distribution from August to December 2011. Primary endpoints were SAM (weight-for-length Z-score [WLZ]<−3 and/or mid-upper arm circumference [MUAC]<11.5 cm and/or bipedal edema) and MAM (−3≤WLZ<−2 and/or 11.5≤MUAC<12.5 cm). A total of 5,395 children were included in the analysis (615 to 1,054 per group). Incidence of MAM was twice lower in the strategies receiving a food supplement combined with cash compared with the cash-only strategy (cash versus HQ-LNS/cash adjusted hazard ratio [HR] = 2.30, 95% CI 1.60–3.29; cash versus SC+/cash HR = 2.42, 95% CI 1.39–4.21; cash versus MQ-LNS/cash HR = 2.07, 95% CI 1.52–2.83) or with the supplementary food only groups (HQ-LNS versus HQ-LNS/cash HR = 1.84, 95% CI 1.35–2.51; SC+ versus SC+/cash HR = 2.53, 95% CI 1.47–4.35). In addition, the incidence of SAM was three times lower in the SC+/cash group compared with the SC+ only group (SC+ only versus SC+/cash HR = 3.13, 95% CI 1.65–5.94). However, non-quantified differences between groups, may limit the interpretation of the impact of the strategies.

Conclusions

Preventive distributions combining a supplementary food and cash transfer had a better preventive effect on MAM and SAM than strategies relying on cash transfer or supplementary food alone. As a result, distribution of nutritious supplementary foods to young children in conjunction with household support should remain a pillar of emergency nutritional interventions. Additional rigorous research is vital to evaluate the effectiveness of these and other nutritional interventions in diverse settings.

Trial registration

ClinicalTrials.gov {"type":"clinical-trial","attrs":{"text":"NCT01828814","term_id":"NCT01828814"}}NCT01828814 Please see later in the article for the Editors'' Summary     

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

Key message

Linkage analysis confirmed the association in the region of PHYC in pearl millet. The comparison of genes found in this region suggests that PHYC is the best candidate.

Abstract

Major efforts are currently underway to dissect the phenotype–genotype relationship in plants and animals using existing populations. This method exploits historical recombinations accumulated in these populations. However, linkage disequilibrium sometimes extends over a relatively long distance, particularly in genomic regions containing polymorphisms that have been targets for selection. In this case, many genes in the region could be statistically associated with the trait shaped by the selected polymorphism. Statistical analyses could help in identifying the best candidate genes into such a region where an association is found. In a previous study, we proposed that a fragment of the PHYTOCHROME C gene (PHYC) is associated with flowering time and morphological variations in pearl millet. In the present study, we first performed linkage analyses using three pearl millet F₂ families to confirm the presence of a QTL in the vicinity of PHYC. We then analyzed a wider genomic region of ~100 kb around PHYC to pinpoint the gene that best explains the association with the trait in this region. A panel of 90 pearl millet inbred lines was used to assess the association. We used a Markov chain Monte Carlo approach to compare 75 markers distributed along this 100-kb region. We found the best candidate markers on the PHYC gene. Signatures of selection in this region were assessed in an independent data set and pointed to the same gene. These results foster confidence in the likely role of PHYC in phenotypic variation and encourage the development of functional studies.