Genome-wide association analysis of seedling root development in maize (Zea mays L.)

Background

Plants rely on the root system for anchorage to the ground and the acquisition and absorption of nutrients critical to sustaining productivity. A genome wide association analysis enables one to analyze allelic diversity of complex traits and identify superior alleles. 384 inbred lines from the Ames panel were genotyped with 681,257 single nucleotide polymorphism markers using Genotyping-by-Sequencing technology and 22 seedling root architecture traits were phenotyped.

Results

Utilizing both a general linear model and mixed linear model, a GWAS study was conducted identifying 268 marker trait associations (p ≤ 5.3×10^-7). Analysis of significant SNP markers for multiple traits showed that several were located within gene models with some SNP markers localized within regions of previously identified root quantitative trait loci. Gene model GRMZM2G153722 located on chromosome 4 contained nine significant markers. This predicted gene is expressed in roots and shoots.

Conclusion

This study identifies putatively associated SNP markers associated with root traits at the seedling stage. Some SNPs were located within or near (<1 kb) gene models. These gene models identify possible candidate genes involved in root development at the seedling stage. These and respective linked or functional markers could be targets for breeders for marker assisted selection of seedling root traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1226-9) contains supplementary material, which is available to authorized users.     

Genetic‐based dissection of arsenic accumulation in maize using a genome‐wide association analysis method

Understanding the mechanism of arsenic (As) accumulation in plants is important in reducing As's toxicity to plants and its potential risks to human health. Here, we performed a genome‐wide association study to dissect the genetic basis of the As contents of different maize tissues in Xixian, which was irrigated with As‐rich surface water, and Changge using an association population consisting of 230 representative maize inbred lines. Phenotypic data revealed a wide normal distribution and high repeatability for the As contents in maize tissues. The As concentrations in maize tissues followed the same trend in the two locations: kernels < axes < stems < bracts < leaves. In total, 15, 16 and 15 non‐redundant quantitative trait loci (QTL s) associated with As concentrations were identified (P  ≤ 2.04 × 10^?6) in five tissues from Xixian, Changge, and the combination of the locations, respectively, explaining 9.70%–24.65% of the phenotypic variation for each QTL , on average. Additionally, four QTL s [involving 15 single nucleotide polymorphisms (SNP s)] were detected in the single and the combined locations, indicating that these loci/SNP s might be stable across different environments. The candidate genes associated with these four loci were predicted. In addition, four non‐redundant QTL s (6 SNP s), including a QTL that was detected in multiple locations according to the genome‐wide association study, were found to co‐localize with four previously reported QTL intervals. These results are valuable to understand the genetic architecture of As mechanism in maize and facilitate the genetic improvement of varieties without As toxicity.     

Distribution of growth and proton efflux in gravireactive roots of maize (Zea mays L.)

Roots of Zea mays were maintained in a vertical orhorizontal position and the local elongation rate and H⁺ fluxes were measured using Sephadex beads containing a pH indicator. When the roots were kept horizontally, the growth of the lower side was strongly inhibited and that of the upper side slightly stimulated as compared with vertical roots. The H⁺ extrusion, which was greatest in the elongation zone, was strongly inhibited on the lower side and slightly stimulated on the upper side as compared with vertical roots.     

Inducing salt tolerance in maize (Zea mays L.) through seed priming with chloride salts: Growth and ion transport at early growth stages

The study was conducted to determine whether salt tolerance could be induced in maize at germination stage by soaking of seeds for 8 h in distilled water or in 200 meq·L⁻¹ of NaCl, KCl, CaCl₂·2H₂O. Both primed and un-primed seeds were subjected for 14 days to 0, 100 or 200 mol·m⁻³ NaCl under controlled conditions. Although all priming agents were effective in alleviating adverse effects of salt stress on maize at germination stage, CaCl₂·2H₂O proved to be more effective since the seeds primed with this salt had significantly higher final germination, rate of germination and fresh and dry weights of plumules and radicles than those treated with other salts or distilled water. Concentration of Na⁺, K⁺ and Ca²⁺ increased significantly in all parts of germinating seeds of maize seeds primed with NaCl, KCl, or CaCl₂·2H₂O, respectively. In addition, seeds primed with CaCl₂·2H₂O were the highest in Cl⁻ accumulation in all parts of the germinating seeds, followed by seeds treated with NaCl and KCl. Most of the Ca²⁺ was retained in seeds and mesocotyl, because of which, transport of this ion to plumules and radicles was low.     

Novel temporal, fine-scale and growth variation phenotypes in roots of adult-stage maize (Zea mays L.) in response to low nitrogen stress

There is interest in discovering root traits associated with acclimation to nutrient stress. Large root systems, such as in adult maize, have proven difficult to be phenotyped comprehensively and over time, causing target traits to be missed. These challenges were overcome here using aeroponics, a system where roots grow in the air misted with a nutrient solution. Applying an agriculturally relevant degree of low nitrogen (LN) stress, 30-day-old plants responded by increasing lengths of individual crown roots (CRs) by 63%, compensated by a 40% decline in CR number. LN increased the CR elongation rate rather than lengthening the duration of CR growth. Only younger CR were significantly responsive to LN stress, a novel finding. LN shifted the root system architectural balance, increasing the lateral root (LR)-to-CR ratio, adding ～70 m to LR length. LN caused a dramatic increase in second-order LR density, not previously reported in adult maize. Despite the near-uniform aeroponics environment, LN induced increased variation in the relative lengths of opposing LR pairs. Large-scale analysis of root hairs (RHs) showed that LN decreased RH length and density. Time-course experiments suggested the RH responses may be indirect consequences of decreased biomass/demand under LN. These results identify novel root traits for genetic dissection.     

Genetic association with high‐resolution climate data reveals selection footprints in the genomes of barley landraces across the Iberian Peninsula

Landraces are local populations of crop plants adapted to a particular environment. Extant landraces are surviving genetic archives, keeping signatures of the selection processes experienced by them until settling in their current niches. This study intends to establish relationships between genetic diversity of barley (Hordeum vulgare L.) landraces collected in Spain and the climate of their collection sites. A high‐resolution climatic data set (5 × 5 km spatial, 1‐day temporal grid) was computed from over 2,000 temperature and 7,000 precipitation stations across peninsular Spain. This data set, spanning the period 1981–2010, was used to derive agroclimatic variables meaningful for cereal production at the collection sites of 135 barley landraces. Variables summarize temperature, precipitation, evapotranspiration, potential vernalization and frost probability at different times of the year and time scales (season and month). SNP genotyping of the landraces was carried out combining Illumina Infinium assays and genotyping‐by‐sequencing, yielding 9,920 biallelic markers (7,479 with position on the barley reference genome). The association of these SNPs with agroclimatic variables was analysed at two levels of genetic diversity, with and without taking into account population structure. The whole data sets and analysis pipelines are documented and available at https://eead-csic-compbio.github.io/barley-agroclimatic-association . We found differential adaptation of the germplasm groups identified to be dominated by reactions to cold temperature and late‐season frost occurrence, as well as to water availability. Several significant associations pointing at specific adaptations to agroclimatic features related to temperature and water availability were observed, and candidate genes underlying some of the main regions are proposed.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

Genome‐wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield‐limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome‐wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific‐locus amplified fragment sequencing (SLAF‐seq) approach. Trait–SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10^–5 (P = 1/85 630, –log₁₀P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P‐value (–log₁₀P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372‐kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR‐NBS‐LRR domains that may be involved in disease resistance response. Real‐time quantitative PCR and virus‐induced gene silencing (VIGS) analysis showed that CG02 was specific to up‐regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.     

Methane yield and species diversity dynamics of perennial wild plant mixtures established alone,under cover crop maize (Zea mays L.), and after spring barley (Hordeum vulgare L.)

The cultivation of perennial wild plant mixtures (WPMs) in biogas cropping systems dominated by maize (Zea mays L.) restores numerous ecosystem functions and improves both spatial and temporal agrobiodiversity. In addition, the colorful appearance of WPM can help enhance landscape beauty. However, their methane yield per hectare (MYH) varies greatly and amounts to only about 50% that of maize. This study aimed at decreasing MYH variability and increasing accumulated MYH of WPM by optimizing the establishment method. A field trial was established in southwest Germany in 2014, and is still running. It tested the effects of three WPM establishment procedures (E1: alone [without maize, in May], E2: undersown in cover crop maize [in May], E3: WPM sown after whole‐crop harvest of spring barley [Hordeum vulgare L.] in June) on both MYH and species diversity of two WPMs [S1, S2]). Mono‐cropped maize and cup plant (Silphium perfoliatum L.) were used as reference crops. Of the WPM treatments tested, S2E2 achieved the highest (19,296 , 60.5% of maize) and S1E1 the lowest accumulated MYH (8,156 , 25.6% of maize) in the years 2014–2018. Cup plant yielded slightly higher than S2E2 (19,968 , 62.6% of maize). In 2014, the WPM sown under maize did not significantly affect the cover crop performance. From 2015 onward, E1 and E2 had comparable average annual MYH and average annual number of WPM species. With a similar accumulated MYH but significantly higher number of species (3.5–10.2), WPM S2E2 outperformed cup plant. Overall, the long‐term MYH performance of WPM cultivation for biogas production can be significantly improved by undersowing with maize as cover crop. This improved establishment method could help facilitate the implementation of WPM cultivation for biogas production and thus reduce the trade‐off between bioenergy and biodiversity.     

Changes in root-shoot ratio and ion uptake of maize (Zea mays L.) from soil as influenced by a plant growth regulator

A pot experiment with maize cv. Limac was conducted to investigate the influence of BAS 110.. W, a plant growth regulator (PGR), on root and shoot development and nutrient uptake. The PGR was applied via the soil with 0, 5, 10, 20, and 40 mg a.i. per pot. Shoot dry matter production was reduced to a higher degree than root length, resulting in a higher root-shoot ratio (RSR) of the treated plants. Shoots of treated plants contained higher concentrations of N, P, Ca, Mg, and unchanged K concentrations. The alterations in concentration could be explained by the changes in RSR induced by the plant growth retardant. The effect was strongest with P (+40%) which was limited by soil supply. N, Ca, and Mgconcentrations were positively influenced (+20%), there was no increase for Kvs RSR.     

Effect of plant growth promoting rhizobacteria containing ACC-deaminase on maize (Zea mays L.) growth under axenic conditions and on nodulation in mung bean (Vigna radiata L.)

AIMS: This study was conducted to test the hypothesis that the bacterial strains possessing 1-aminocyclopropane-1-carboxylic acid (ACC)-deaminase activity may also promote growth of inoculated plants and could increase nodulation in legumes upon co-inoculation with rhizobia. METHODS AND RESULTS: Several rhizobacteria were isolated from maize rhizosphere through enrichment on ACC as a sole N source. Purified isolates were screened for growth promotion in maize under axenic conditions and for in vitro ACC-deaminase activity. A significant positive correlation was observed between in vitro ACC-deaminase activity of bacterial cells and root elongation. None of the isolates produced auxins. Bradyrhizobium japonicum produced less amount of auxins but did not carry ACC-deaminase activity. Results of pot experiment revealed that co-inoculation with Bradyrhizobium and plant growth promoting rhizobacteria (PGPR) isolates enhanced the nodulation in mung bean compared with inoculation with Bradyrhizobium alone. CONCLUSIONS: It is highly expected that inoculation with rhizobacteria containing ACC-deaminase hydrolysed endogenous ACC into ammonia and alpha-ketobutyrate instead of ethylene. Consequently, root and shoot growth as well as nodulation were promoted. SIGNIFICANCE AND IMPACT OF THE STUDY: The ACC-deaminase trait could be employed as an efficient tool to screen effective PGPR, which could be successfully used as biofertilizers to increase the growth of inoculated plants as well as nodulation in legumes.     

High‐throughput two‐dimensional root system phenotyping platform facilitates genetic analysis of root growth and development

High‐throughput phenotyping of root systems requires a combination of specialized techniques and adaptable plant growth, root imaging and software tools. A custom phenotyping platform was designed to capture images of whole root systems, and novel software tools were developed to process and analyse these images. The platform and its components are adaptable to a wide range root phenotyping studies using diverse growth systems (hydroponics, paper pouches, gel and soil) involving several plant species, including, but not limited to, rice, maize, sorghum, tomato and Arabidopsis. The RootReader2D software tool is free and publicly available and was designed with both user‐guided and automated features that increase flexibility and enhance efficiency when measuring root growth traits from specific roots or entire root systems during large‐scale phenotyping studies. To demonstrate the unique capabilities and high‐throughput capacity of this phenotyping platform for studying root systems, genome‐wide association studies on rice (Oryza sativa) and maize (Zea mays) root growth were performed and root traits related to aluminium (Al) tolerance were analysed on the parents of the maize nested association mapping (NAM) population.     

Use of computed tomography as a non‐invasive method for diagnosing cephenemyiosis in roe deer (Capreolus capreolus)

This study was conducted to assess the reliability of computed tomography (CT) for diagnosing bot fly infestations by Cephenemyia stimulator (Clark) (Diptera: Oestridae) in roe deer (Capreolus capreolus L.) (Artiodactyla: Cervidae). For this purpose, the heads of 30 animals were analysed, firstly by CT and then by necropsy, which was used as the reference standard method. The prevalence values obtained by both methods were identical; the prevalence of infestation was 40.0% overall, and was higher in males (45.5%) than in females (25.0%). These results highlight the usefulness of CT as an alternative or non‐invasive method for diagnosing cephenemyiosis in live‐captured roe deer and in hunting trophies or museum collections that cannot be destroyed or damaged.     

Roothairless5, which functions in maize (Zea mays L.) root hair initiation and elongation encodes a monocot‐specific NADPH oxidase

Root hairs are instrumental for nutrient uptake in monocot cereals. The maize (Zea mays L.) roothairless5 (rth5) mutant displays defects in root hair initiation and elongation manifested by a reduced density and length of root hairs. Map‐based cloning revealed that the rth5 gene encodes a monocot‐specific NADPH oxidase. RNA‐Seq, in situ hybridization and qRT‐PCR experiments demonstrated that the rth5 gene displays preferential expression in root hairs but also accumulates to low levels in other tissues. Immunolocalization detected RTH5 proteins in the epidermis of the elongation and differentiation zone of primary roots. Because superoxide and hydrogen peroxide levels are reduced in the tips of growing rth5 mutant root hairs as compared with wild‐type, and Reactive oxygen species (ROS) is known to be involved in tip growth, we hypothesize that the RTH5 protein is responsible for establishing the high levels of ROS in the tips of growing root hairs required for elongation. Consistent with this hypothesis, a comparative RNA‐Seq analysis of 6‐day‐old rth5 versus wild‐type primary roots revealed significant over‐representation of only two gene ontology (GO) classes related to the biological functions (i.e. oxidation/reduction and carbohydrate metabolism) among 893 differentially expressed genes (FDR <5%). Within these two classes the subgroups ‘response to oxidative stress’ and ‘cellulose biosynthesis’ were most prominently represented.     

Genome‐wide SNP analysis reveals a genetic basis for sea‐age variation in a wild population of Atlantic salmon (Salmo salar)

Delaying sexual maturation can lead to larger body size and higher reproductive success, but carries an increased risk of death before reproducing. Classical life history theory predicts that trade‐offs between reproductive success and survival should lead to the evolution of an optimal strategy in a given population. However, variation in mating strategies generally persists, and in general, there remains a poor understanding of genetic and physiological mechanisms underlying this variation. One extreme case of this is in the Atlantic salmon (Salmo salar), which can show variation in the age at which they return from their marine migration to spawn (i.e. their ‘sea age’). This results in large size differences between strategies, with direct implications for individual fitness. Here, we used an Illumina Infinium SNP array to identify regions of the genome associated with variation in sea age in a large population of Atlantic salmon in Northern Europe, implementing individual‐based genome‐wide association studies (GWAS) and population‐based F_ST outlier analyses. We identified several regions of the genome which vary in association with phenotype and/or selection between sea ages, with nearby genes having functions related to muscle development, metabolism, immune response and mate choice. In addition, we found that individuals of different sea ages belong to different, yet sympatric populations in this system, indicating that reproductive isolation may be driven by divergence between stable strategies. Overall, this study demonstrates how genome‐wide methodologies can be integrated with samples collected from wild, structured populations to understand their ecology and evolution in a natural context.     

Effect of sub-optimal root zone temperatures at varied nutrient supply and shoot meristem temperature on growth and nutrient concentrations in maize seedlings (Zea mays L.)

Maize seedlings were grown for 10 to 20 days in either nutrient solution or in soils with or without fertilizer supply. Air temperature was kept uniform for all treatments, while root zone temperature (RZT) was varied between 12 and 24°C. In some treatments the basal part of the shoot (with apical shoot meristem and zone of leaf elongation) was lifted up to separate the indirect effects of root zone temperature on shoot growth from the direct effects of temperature on the shoot meristem.Shoot and root growth were decreased by low RZT to a similar extent irrespective of the growth medium (i.e. nutrient solution, fertilized or unfertilized soil). In all culture media Ca concentration was similar or even higher in plants grown at 12 as compared to 24°. At lower RZT concentrations of N, P and K in the shoot dry matter decreased in unfertilized soil, whereas in nutrient solution and fertilized soil only the K concentration decreased.When direct temperature effects on the shoot meristem were reduced by lifting the basal part of the shoot above the temperature-controlled root zone, shoot growth at low RZT was significantly increased in nutrient solution and fertilized soil, but not in unfertilized soil. In fertilized soil and nutrient solution at low RZT the uptake of K increased to a similar extent as plant growth, and thus shoot K concentration was not reduced by increasing shoot growth rates. In contrast, uptake of N and P was not increased, resulting in significantly decreased shoot concentrations.It is concluded that shoot growth at suboptimal RZT was limited both by a direct temperature effect on shoot activity and by a reduced nutrient supply through the roots. Nutrient concentrations in the shoot tissue at low RZT were not only influenced by availability in the substrate and dilution by growth, but also by the internal demand for growth.     

Genome‐wide identification,expression analysis of auxin‐responsive GH3 family genes in maize (Zea mays L.) under abiotic stresses

Auxin is involved in different aspects of plant growth and development by regulating the expression of auxinresponsive family genes. As one of the three major auxinresponsive families, GH3(Gretchen Hagen3) genes participate in auxin homeostasis by catalyzing auxin conjugation and bounding free indole-3-acetic acid(IAA) to amino acids.However, how GH3 genes function in responses to abiotic stresses and various hormones in maize is largely unknown.Here, the latest updated maize(Zea mays L.) reference genome sequence was used to characterize and analyze the Zm GH3 family genes from maize. The results showed that 13 Zm GH3 genes were mapped on fi ve maize chromosomes(total10 chromosomes). Highly diversi fi ed gene structures and tissue-speci fi c expression patterns suggested the possibility of function diversi fi cation for these genes in response to environmental stresses and hormone stimuli. The expression patterns of Zm GH3 genes are responsive to several abiotic stresses(salt, drought and cadmium) and major stress-relatedhormones(abscisic acid, salicylic acid and jasmonic acid)Various environmental factors suppress auxin free IAA contents in maize roots suggesting that these abiotic stresses and hormones might alter GH3-mediated auxin levels. The respon siveness of Zm GH3 genes to a wide range of abiotic stresses and stress-related hormones suggested that Zm GH3 s are involved in maize tolerance to environmental stresses.