DNA polymorphisms and haplotype patterns of transcription factors involved in barley endosperm development are associated with key agronomic traits

Genetic loci mediating circadian clock output plasticity and crop productivity under barley domestication

Circadian clock rhythms are shown to be intertwined with crop adaptation. To realize the adaptive value of changes in these rhythms under crop domestication and improvement, there is a need to compare the genetics of clock and yield traits. We compared circadian clock rhythmicity based on Chl leaf fluorescence and transcriptomics among wild ancestors, landraces, and breeding lines of barley under optimal and high temperatures. We conducted a genome scan to identify pleiotropic loci regulating the clock and field phenotypes. We also compared the allelic diversity in wild and cultivated barley to test for selective sweeps. We found significant loss of thermal plasticity in circadian rhythms under domestication. However, transcriptome analysis indicated that this loss was only for output genes and that temperature compensation in the core clock machinery was maintained. Drivers of the circadian clock (DOC) loci were identified via genome-wide association study. Notably, these loci also modified growth and reproductive outputs in the field. Diversity analysis indicated selective sweep in these pleiotropic DOC loci. These results indicate a selection against thermal clock plasticity under barley domestication and improvement and highlight the importance of identifying genes underlying for understanding the biochemical basis of crop adaptation to changing environments.

Deep sampling and pooled amplicon sequencing reveals hidden genic variation in heterogeneous rye accessions

BACKGROUND

Loss of genetic variation negatively impacts breeding efforts and food security. Genebanks house over 7 million accessions representing vast allelic diversity that is a resource for sustainable breeding. Discovery of DNA variations is an important step in the efficient use of these resources. While technologies have improved and costs dropped, it remains impractical to consider resequencing millions of accessions. Candidate genes are known for most agronomic traits, providing a list of high priority targets. Heterogeneity in seed stocks means that multiple samples from an accession need to be evaluated to recover available alleles. To address this we developed a pooled amplicon sequencing approach and applied it to the out-crossing cereal rye (Secale cereale L.).

RESULTS

Using the amplicon sequencing approach 95 rye accessions of different improvement status and worldwide origin, each represented by a pooled sample comprising DNA of 96 individual plants, were evaluated for sequence variation in six candidate genes with significant functions on biotic and abiotic stress resistance, and seed quality. Seventy-four predicted deleterious variants were identified using multiple algorithms. Rare variants were recovered including those found only in a low percentage of seed.

CONCLUSIONS

We conclude that this approach provides a rapid and flexible method for evaluating stock heterogeneity, probing allele diversity, and recovering previously hidden variation. A large extent of within-population heterogeneity revealed in the study provides an important point for consideration during rye germplasm conservation and utilization efforts.

Development of Decreased-Gluten Wheat Enabled by Determination of the Genetic Basis of lys3a Barley

Celiac disease is the most common food-induced enteropathy in humans, with a prevalence of approximately 1% worldwide. It is induced by digestion-resistant, proline- and glutamine-rich seed storage proteins, collectively referred to as gluten, found in wheat (Triticum aestivum). Related prolamins are present in barley (Hordeum vulgare) and rye (Secale cereale). The incidence of both celiac disease and a related condition called nonceliac gluten sensitivity is increasing. This has prompted efforts to identify methods of lowering gluten in wheat, one of the most important cereal crops. Here, we used bulked segregant RNA sequencing and map-based cloning to identify the genetic lesion underlying a recessive, low-prolamin mutation (lys3a) in diploid barley. We confirmed the mutant identity by complementing the lys3a mutant with a transgenic copy of the wild-type barley gene and then used targeting-induced local lesions in genomes to identify induced single-nucleotide polymorphisms in the three homeologs of the corresponding wheat gene. Combining inactivating mutations in the three subgenomes of hexaploid bread wheat in a single wheat line lowered gliadin and low-molecular-weight glutenin accumulation by 50% to 60% and increased free and protein-bound lysine by 33%.

Construction of a network describing asparagine metabolism in plants and its application to the identification of genes affecting asparagine metabolism in wheat under drought and nutritional stress

A detailed network describing asparagine metabolism in plants was constructed using published data from Arabidopsis (Arabidopsis thaliana) maize (Zea mays), wheat (Triticum aestivum), pea (Pisum sativum), soybean (Glycine max), lupin (Lupus albus), and other species, including animals. Asparagine synthesis and degradation is a major part of amino acid and nitrogen metabolism in plants. The complexity of its metabolism, including limiting and regulatory factors, was represented in a logical sequence in a pathway diagram built using yED graph editor software. The network was used with a Unique Network Identification Pipeline in the analysis of data from 18 publicly available transcriptomic data studies. This identified links between genes involved in asparagine metabolism in wheat roots under drought stress, wheat leaves under drought stress, and wheat leaves under conditions of sulfur and nitrogen deficiency. The network represents a powerful aid for interpreting the interactions not only between the genes in the pathway but also among enzymes, metabolites and smaller molecules. It provides a concise, clear understanding of the complexity of asparagine metabolism that could aid the interpretation of data relating to wider amino acid metabolism and other metabolic processes.

A Catalog of Regulatory Sequences for Trait Gene for the Genome Editing of Wheat

Wheat has been cultivated for 10000 years and ever since the origin of hexaploid wheat it has been exempt from natural selection. Instead, it was under the constant selective pressure of human agriculture from harvest to sowing during every year, producing a vast array of varieties. Wheat has been adopted globally, accumulating variation for genes involved in yield traits, environmental adaptation and resistance. However, one small but important part of the wheat genome has hardly changed: the regulatory regions of both the x- and y-type high molecular weight glutenin subunit (HMW-GS) genes, which are alone responsible for approximately 12% of the grain protein content. The phylogeny of the HMW-GS regulatory regions of the Triticeae demonstrates that a genetic bottleneck may have led to its decreased diversity during domestication and the subsequent cultivation. It has also highlighted the fact that the wild relatives of wheat may offer an unexploited genetic resource for the regulatory region of these genes. Significant research efforts have been made in the public sector and by international agencies, using wild crosses to exploit the available genetic variation, and as a result synthetic hexaploids are now being utilized by a number of breeding companies. However, a newly emerging tool of genome editing provides significantly improved efficiency in exploiting the natural variation in HMW-GS genes and incorporating this into elite cultivars and breeding lines. Recent advancement in the understanding of the regulation of these genes underlines the needs for an overview of the regulatory elements for genome editing purposes.

A Genome-wide Combinatorial Strategy Dissects Complex Genetic Architecture of Seed Coat Color in Chickpea

The study identified 9045 high-quality SNPs employing both genome-wide GBS- and candidate gene-based SNP genotyping assays in 172, including 93 cultivated (desi and kabuli) and 79 wild chickpea accessions. The GWAS in a structured population of 93 sequenced accessions detected 15 major genomic loci exhibiting significant association with seed coat color. Five seed color-associated major genomic loci underlying robust QTLs mapped on a high-density intra-specific genetic linkage map were validated by QTL mapping. The integration of association and QTL mapping with gene haplotype-specific LD mapping and transcript profiling identified novel allelic variants (non-synonymous SNPs) and haplotypes in a MATE secondary transporter gene regulating light/yellow brown and beige seed coat color differentiation in chickpea. The down-regulation and decreased transcript expression of beige seed coat color-associated MATE gene haplotype was correlated with reduced proanthocyanidins accumulation in the mature seed coats of beige than light/yellow brown seed colored desi and kabuli accessions for their coloration/pigmentation. This seed color-regulating MATE gene revealed strong purifying selection pressure primarily in LB/YB seed colored desi and wild Cicer reticulatum accessions compared with the BE seed colored kabuli accessions. The functionally relevant molecular tags identified have potential to decipher the complex transcriptional regulatory gene function of seed coat coloration and for understanding the selective sweep-based seed color trait evolutionary pattern in cultivated and wild accessions during chickpea domestication. The genome-wide integrated approach employed will expedite marker-assisted genetic enhancement for developing cultivars with desirable seed coat color types in chickpea.

Genome-wide high-throughput SNP discovery and genotyping for understanding natural (functional) allelic diversity and domestication patterns in wild chickpea

We identified 82489 high-quality genome-wide SNPs from 93 wild and cultivated Cicer accessions through integrated reference genome- and de novo-based GBS assays. High intra- and inter-specific polymorphic potential (66-85%) and broader natural allelic diversity (6-64%) detected by genome-wide SNPs among accessions signify their efficacy for monitoring introgression and transferring target trait-regulating genomic (gene) regions/allelic variants from wild to cultivated Cicer gene pools for genetic improvement. The population-specific assignment of wild Cicer accessions pertaining to the primary gene pool are more influenced by geographical origin/phenotypic characteristics than species/gene-pools of origination. The functional significance of allelic variants (non-synonymous and regulatory SNPs) scanned from transcription factors and stress-responsive genes in differentiating wild accessions (with potential known sources of yield-contributing and stress tolerance traits) from cultivated desi and kabuli accessions, fine-mapping/map-based cloning of QTLs and determination of LD patterns across wild and cultivated gene-pools are suitably elucidated. The correlation between phenotypic (agromorphological traits) and molecular diversity-based admixed domestication patterns within six structured populations of wild and cultivated accessions via genome-wide SNPs was apparent. This suggests utility of whole genome SNPs as a potential resource for identifying naturally selected trait-regulating genomic targets/functional allelic variants adaptive to diverse agroclimatic regions for genetic enhancement of cultivated gene-pools.

Genome-wide association mapping for kernel and malting quality traits using historical European barley records

Malting quality is an important trait in breeding barley (Hordeum vulgare L.). It requires elaborate, expensive phenotyping, which involves micro-malting experiments. Although there is abundant historical information available for different cultivars in different years and trials, that historical information is not often used in genetic analyses. This study aimed to exploit historical records to assist in identifying genomic regions that affect malting and kernel quality traits in barley. This genome-wide association study utilized information on grain yield and 18 quality traits accumulated over 25 years on 174 European spring and winter barley cultivars combined with diversity array technology markers. Marker-trait associations were tested with a mixed linear model. This model took into account the genetic relatedness between cultivars based on principal components scores obtained from marker information. We detected 140 marker-trait associations. Some of these associations confirmed previously known quantitative trait loci for malting quality (on chromosomes 1H, 2H, and 5H). Other associations were reported for the first time in this study. The genetic correlations between traits are discussed in relation to the chromosomal regions associated with the different traits. This approach is expected to be particularly useful when designing strategies for multiple trait improvements.

Natural allelic diversity, genetic structure and linkage disequilibrium pattern in wild chickpea

Characterization of natural allelic diversity and understanding the genetic structure and linkage disequilibrium (LD) pattern in wild germplasm accessions by large-scale genotyping of informative microsatellite and single nucleotide polymorphism (SNP) markers is requisite to facilitate chickpea genetic improvement. Large-scale validation and high-throughput genotyping of genome-wide physically mapped 478 genic and genomic microsatellite markers and 380 transcription factor gene-derived SNP markers using gel-based assay, fluorescent dye-labelled automated fragment analyser and matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) mass array have been performed. Outcome revealed their high genotyping success rate (97.5%) and existence of a high level of natural allelic diversity among 94 wild and cultivated Cicer accessions. High intra- and inter-specific polymorphic potential and wider molecular diversity (11-94%) along with a broader genetic base (13-78%) specifically in the functional genic regions of wild accessions was assayed by mapped markers. It suggested their utility in monitoring introgression and transferring target trait-specific genomic (gene) regions from wild to cultivated gene pool for the genetic enhancement. Distinct species/gene pool-wise differentiation, admixed domestication pattern, and differential genome-wide recombination and LD estimates/decay observed in a six structured population of wild and cultivated accessions using mapped markers further signifies their usefulness in chickpea genetics, genomics and breeding.

The BdGAMYB protein from Brachypodium distachyon interacts with BdDOF24 and regulates transcription of the BdCathB gene upon seed germination

BdDOF24 interacting with BdGAMYB regulates the BdCathB gene upon germination. During barley seed germination, hydrolytic enzymes (α-amylases, proteases, etc.) synthesized in the aleurone layer in response to gibberellins (GA), catalyse the mobilization of storage reserves accumulated in the endosperm during seed maturation. In Brachypodium distachyon, the BdCathB gene that encodes a Cathepsin B-like thiol-protease, orthologous to the wheat Al21 and barley HvCathB, is highly induced in germinating seeds and its expression is regulated by transcription factors (TFs) encoded by genes BdGamyb and BdDof24, orthologous to the barley HvGamyb and BPBF-HvDof24, respectively. Transcripts of both TF genes increase during germination and treatments with abscisic acid (ABA) or paclobutrazol (PAC, an inhibitor of GA biosynthesis) decrease mRNA expression of BdGamyb but do not affect that of BdDof24. Besides, proteins BdDOF24 and BdGAMYB interact in yeast-2 hybrid systems and in plant nuclei, and in transient expression assays in aleurone layers BdDOF24 is a transcriptional repressor and BdGAMYB is an activator of the BdCathB promoter, as occurs with the putative orthologous in barley BPBF-HvDOF24 and HvGAMYB. However, when both TFs are co-bombarded, BdDOF24 enhances the activation driven by BdGAMYB while BPBF-HvDOF24 strongly decreases the HvGAMYB-mediated activation of the BdCathB promoter. The different results obtained when BdDOF24 and BPBF-HvDOF24 interact with BdGAMYB and HvGAMYB are discussed.

Genome-wide association studies of agronomic and quality traits in a set of German winter barley (Hordeum vulgare L.) cultivars using Diversity Arrays Technology (DArT)

A set of about 100 winter barley (Hordeum vulgare L.) cultivars, comprising diverse and economically important German barley elite germplasm released during the last six decades, was previously genotypically characterized by single nucleotide polymorphism (SNP) markers using the Illumina GoldenGate BeadArray Technology to detect associations with phenotypic data estimated in three-year field trials at 12 locations. In order to identify further associations and to obtain information on whether the marker type influences the outcome of association genetics studies, the set of winter barley cultivars was re-analyzed using Diversity Arrays Technology (DArT) markers. As with the analysis of the SNPs, only polymorphic markers present at an allele frequency >5% were included to detect associations in a mixed linear model (MLM) approach using the TASSEL software (P ≤ 0.001). The population structure and kinship matrix were estimated on 72 simple sequence repeats (SSRs) covering the whole barley genome. The respective average linkage disequilibrium (LD) analyzed with DArT markers was estimated at 5.73 cM. A total of 52 markers gave significant associations with at least one of the traits estimated which, therefore, may be suitable for marker-assisted breeding. In addition, by comparing the results to those generated using the Illumina GoldenGate BeadArray Technology, it turned out that a different number of associations for respective traits is detected, depending on the marker system. However, as only a few of the respective DArT and Illumina markers are present in a common map, no comprehensive comparison of the detected associations was feasible, but some were probably detected in the same chromosomal regions. Because of the identification of additional marker-trait associations, it may be recommended to use both marker techniques in genome-wide association studies.

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana. Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r (2) = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r (2) = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana. Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r (2) = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r (2) = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.

Development of SNP markers for genes of the phenylpropanoid pathway and their association to kernel and malting traits in barley

BACKGROUND

Flavonoids are an important class of secondary compounds in angiosperms. Next to certain biological functions in plants, they play a role in the brewing process and have an effect on taste, color and aroma of beer. The aim of this study was to reveal the haplotype diversity of candidate genes involved in the phenylpropanoid biosynthesis pathway in cultivated barley varieties (Hordeum vulgare L.) and to determine associations to kernel and malting quality parameters.

RESULTS

Five genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS), flavanone 3-hydroxylase (F3H) and dihydroflavonol reductase (DFR) of the phenylpropanoid biosynthesis pathway were partially resequenced in 16 diverse barley reference genotypes. Their localization in the barley genome, their genetic structure, and their genetic variation e.g. single nucleotide polymorphism (SNP) and Insertion/Deletion (InDel) patterns were revealed. In total, 130 SNPs and seven InDels were detected. Of these, 21 polymorphisms were converted into high-throughput pyrosequencing markers. The resulting SNP and haplotype patterns were used to calculate associations with kernel and malting quality parameters.

CONCLUSIONS

SNP patterns were found to be highly variable for the investigated genes. The developed high-throughput markers are applicable for assessing the genetic variability and for the determination of haplotype patterns in a set of barley accessions. The candidate genes PAL, C4H and F3H were shown to be associated to several malting properties like glassiness (PAL), viscosity (C4H) or to final attenuation (F3H).

Nitrogen-metabolism related genes in barley - haplotype diversity, linkage mapping and associations with malting and kernel quality parameters

Background

Several studies report about intra-specific trait variation of nitrogen-metabolism related traits, such as N(itrogen)-use efficiency, protein content, N-storage and remobilization in barley and related grass species. The goal of this study was to assess the intra-specific genetic diversity present in primary N-metabolism genes of barley and to investigate the associations of the detected haplotype diversity with malting and kernel quality related traits.

Results

Partial sequences of five genes related to N-metabolism in barley (Hordeum vulgare L.) were obtained, i.e. nitrate reductase 1, glutamine synthetase 2, ferredoxin-dependent glutamate synthase, aspartate aminotransferase and asparaginase. Two to five haplotypes in each gene were discovered in a set of 190 various varieties. The development of 33 SNP markers allowed the genotyping of all these barley varieties consisting of spring and winter types. Furthermore, these markers could be mapped in several doubled haploid populations. Cluster analysis based on haplotypes revealed a more uniform pattern of the spring barleys as compared to the winter barleys. Based on linear model approaches associations to several malting and kernel quality traits including soluble N and protein were identified.

Conclusions

A study was conducted to investigate the presence of sequence variation of several genes related to the primary N-metabolism in barley. The detected diversity could be related to particular phenotypic traits. Specific differences between spring and winter barleys most likely reflect different breeding aims. The developed markers can be used as tool for further genetic studies and marker-assisted selection during breeding of barley.

Association study of wheat grain protein composition reveals that gliadin and glutenin composition are trans-regulated by different chromosome regions

Wheat grain storage protein (GSP) content and composition are the main determinants of the end-use value of bread wheat (Triticum aestivum L.) grain. The accumulation of glutenins and gliadins, the two main classes of GSP in wheat, is believed to be mainly controlled at the transcriptional level through a network of transcription factors. This regulation network could lead to stable cross-environment allometric scaling relationships between the quantity of GSP classes/subunits and the total quantity of nitrogen per grain. This work conducted a genetic mapping study of GSP content and composition and allometric scaling parameters of grain N allocation using a bread wheat worldwide core collection grown in three environments. The core collection was genotyped with 873 markers for genome-wide association and 167 single nucleotide polymorphism markers in 51 candidate genes for candidate association. The candidate genes included 35 transcription factors (TFs) expressed in grain. This work identified 74 loci associated with 38 variables, of which 19 were candidate genes or were tightly linked with candidate genes. Besides structural GSP genes, several loci putatively trans-regulating GSP accumulation were identified. Seven candidate TFs, including four wheat orthologues of barley TFs that control hordein gene expression, were associated or in strong linkage disequilibrium with markers associated with the composition or quantity of glutenin or gliadin, or allometric grain N allocation parameters, confirming the importance of the transcriptional control of GSP accumulation. Genome-wide association results suggest that the genes regulating glutenin and gliadin compositions are mostly distinct from each other and operate differently.

Adaptation of maize to temperate climates: mid-density genome-wide association genetics and diversity patterns reveal key genomic regions, with a major contribution of the Vgt2 (ZCN8) locus

The migration of maize from tropical to temperate climates was accompanied by a dramatic evolution in flowering time. To gain insight into the genetic architecture of this adaptive trait, we conducted a 50K SNP-based genome-wide association and diversity investigation on a panel of tropical and temperate American and European representatives. Eighteen genomic regions were associated with flowering time. The number of early alleles cumulated along these regions was highly correlated with flowering time. Polymorphism in the vicinity of the ZCN8 gene, which is the closest maize homologue to Arabidopsis major flowering time (FT) gene, had the strongest effect. This polymorphism is in the vicinity of the causal factor of Vgt2 QTL. Diversity was lower, whereas differentiation and LD were higher for associated loci compared to the rest of the genome, which is consistent with selection acting on flowering time during maize migration. Selection tests also revealed supplementary loci that were highly differentiated among groups and not associated with flowering time in our panel, whereas they were in other linkage-based studies. This suggests that allele fixation led to a lack of statistical power when structure and relatedness were taken into account in a linear mixed model. Complementary designs and analysis methods are necessary to unravel the architecture of complex traits. Based on linkage disequilibrium (LD) estimates corrected for population structure, we concluded that the number of SNPs genotyped should be at least doubled to capture all QTLs contributing to the genetic architecture of polygenic traits in this panel. These results show that maize flowering time is controlled by numerous QTLs of small additive effect and that strong polygenic selection occurred under cool climatic conditions. They should contribute to more efficient genomic predictions of flowering time and facilitate the dissemination of diverse maize genetic resources under a wide range of environments.

Molecular footprints of domestication and improvement in soybean revealed by whole genome re-sequencing

BACKGROUND

Artificial selection played an important role in the origin of modern Glycine max cultivars from the wild soybean Glycine soja. To elucidate the consequences of artificial selection accompanying the domestication and modern improvement of soybean, 25 new and 30 published whole-genome re-sequencing accessions, which represent wild, domesticated landrace, and Chinese elite soybean populations were analyzed.

RESULTS

A total of 5,102,244 single nucleotide polymorphisms (SNPs) and 707,969 insertion/deletions were identified. Among the SNPs detected, 25.5% were not described previously. We found that artificial selection during domestication led to more pronounced reduction in the genetic diversity of soybean than the switch from landraces to elite cultivars. Only a small proportion (2.99%) of the whole genomic regions appear to be affected by artificial selection for preferred agricultural traits. The selection regions were not distributed randomly or uniformly throughout the genome. Instead, clusters of selection hotspots in certain genomic regions were observed. Moreover, a set of candidate genes (4.38% of the total annotated genes) significantly affected by selection underlying soybean domestication and genetic improvement were identified.

CONCLUSIONS

Given the uniqueness of the soybean germplasm sequenced, this study drew a clear picture of human-mediated evolution of the soybean genomes. The genomic resources and information provided by this study would also facilitate the discovery of genes/loci underlying agronomically important traits.

microRNA profiling of root tissues and root forming explant cultures in Medicago truncatula

Plant root architecture is regulated by the initiation and modulation of cell division in regions containing pluripotent stem cells known as meristems. In roots, meristems are formed early in embryogenesis, in the case of the root apical meristem (RAM), and during organogenesis at the site of lateral root or, in legumes, nodule formation. Root meristems can also be generated in vitro from leaf explants cultures supplemented with auxin. microRNAs (miRNAs) have emerged as regulators of many key biological functions in plants including root development. To identify key miRNAs involved in root meristem formation in Medicago truncatula, we used deep sequencing to compare miRNA populations. Comparisons were made between: (1) the root tip (RT), containing the RAM and the elongation zone (EZ) tissue and (2) root forming callus (RFC) and non-root forming callus (NRFC). We identified 83 previously reported miRNAs, 24 new to M. truncatula, in 44 families. For the first time in M. truncatula, members of conserved miRNA families miR165, miR181 and miR397 were found. Bioinformatic analysis identified 38 potential novel miRNAs. Selected miRNAs and targets were validated using Taqman miRNA assays and 5' RACE. Many miRNAs were differentially expressed between tissues, particularly RFC and NRFC. Target prediction revealed a number of miRNAs to target genes previously shown to be differentially expressed between RT and EZ or RFC and NRFC and important in root development. Additionally, we predict the miRNA/target relationships for miR397 and miR160 to be conserved in M. truncatula. Amongst the predictions, were AUXIN RESPONSE FACTOR 10, targeted by miR160 and a LACCASE-like gene, targeted by miR397, both are miRNA/target pairings conserved in other species.

Arabidopsis thaliana bZIP44: a transcription factor affecting seed germination and expression of the mannanase-encoding gene AtMAN7

Endo-β-mannanases (MAN; EC. 3.2.1.78) catalyze the cleavage of β1→4 bonds in mannan polymers and have been associated with the process of weakening the tissues surrounding the embryo during seed germination. In germinating Arabidopsis thaliana seeds, the most highly expressed MAN gene is AtMAN7 and its transcripts are restricted to the micropylar endosperm and to the radicle tip just before radicle emergence. Mutants with a T-DNA insertion in AtMAN7 have a slower germination than the wild type. To gain insight into the transcriptional regulation of the AtMAN7 gene, a bioinformatic search for conserved non-coding cis-elements (phylogenetic shadowing) within the Brassicaceae MAN7 gene promoters has been done, and these conserved motifs have been used as bait to look for their interacting transcription factors (TFs), using as a prey an arrayed yeast library from A. thaliana. The basic-leucine zipper TF AtbZIP44, but not the closely related AtbZIP11, has thus been identified and its transcriptional activation upon AtMAN7 has been validated at the molecular level. In the knock-out lines of AtbZIP44, not only is the expression of the AtMAN7 gene drastically reduced, but these mutants have a significantly slower germination than the wild type, being affected in the two phases of the germination process, both in the rupture of the seed coat and in the breakage of the micropylar endosperm cell walls. In the over-expression lines the opposite phenotype is observed.

Comparative analysis of kabuli chickpea transcriptome with desi and wild chickpea provides a rich resource for development of functional markers

Chickpea (Cicer arietinum L.) is an important crop legume plant with high nutritional value. The transcriptomes of desi and wild chickpea have already been sequenced. In this study, we sequenced the transcriptome of kabuli chickpea, C. arietinum (genotype ICCV2), having higher commercial value, using GS-FLX Roche 454 and Illumina technologies. The assemblies of both Roche 454 and Illumina datasets were optimized using various assembly programs and parameters. The final optimized hybrid assembly generated 43,389 transcripts with an average length of 1065 bp and N50 length of 1653 bp representing 46.2 Mb of kabuli chickpea transcriptome. We identified a total of 5409 simple sequence repeats (SSRs) in these transcript sequences. Among these, at least 130 and 493 SSRs were polymorphic with desi (ICC4958) and wild (PI489777) chickpea, respectively. In addition, a total of 1986 and 37,954 single nucleotide polymorphisms (SNPs) were predicted in kabuli/desi and kabuli/wild genotypes, respectively. The SNP frequency was 0.043 SNP per kb for kabuli/desi and 0.821 SNP per kb for kabuli/wild, reflecting very low genetic diversity in chickpea. Further, SSRs and SNPs present in tissue-specific and transcription factor encoding transcripts have been identified. The experimental validation of a selected set of polymorphic SSRs and SNPs exhibited high intra-specific polymorphism potential between desi and kabuli chickpea, suggesting their utility in large-scale genotyping applications. The kabuli chickpea gene index assembled, and SSRs and SNPs identified in this study will serve as useful genomic resource for genetic improvement of chickpea.

Phenotypic landscapes: phenological patterns in wild and cultivated barley

Phenotypic variation in natural populations is the outcome of the joint effects of environmentally induced adaptations and neutral processes on the genetic architecture of quantitative traits. In this study, we examined the role of adaptation in shaping wild barley phenotypic variation along different environmental gradients. Detailed phenotyping of 164 wild barley (Hordeum spontaneum) accessions from Israel (of the Barley1K collection) and 18 cultivated barley (H. vulgare) varieties was conducted in common garden field trials. Cluster analysis based on phenotypic data indicated that wild barley in this region can be differentiated into three ecotypes in accordance with their ecogeographical distribution: north, coast and desert. Population differentiation (Q(ST) ) for each trait was estimated using a hierarchical Bayesian model and compared to neutral differentiation (F(ST) ) based on 42 microsatellite markers. This analysis indicated that the three clusters diverged in morphological but not in reproductive characteristics. To address the issue of phenotypic variation along environmental gradients, climatic and soil gradients were compared with each of the measured traits given the geographical distance between sampling sites using a partial Mantel test. Flowering time and plant growth were found to be differentially correlated with climatic and soil characteristic gradients, respectively. The H. vulgare varieties were superior to the H. spontaneum accessions in yield components, yet resembled the Mediterranean types in vegetative characteristics and flowering time, which may indicate the geographical origin of domesticated barley.

The family of DOF transcription factors in Brachypodium distachyon: phylogenetic comparison with rice and barley DOFs and expression profiling

BACKGROUND

Transcription factors (TFs) are proteins that have played a central role both in evolution and in domestication, and are major regulators of development in living organisms. Plant genome sequences reveal that approximately 7% of all genes encode putative TFs. The DOF (DNA binding with One Finger) TF family has been associated with vital processes exclusive to higher plants and to their close ancestors (algae, mosses and ferns). These are seed maturation and germination, light-mediated regulation, phytohormone and plant responses to biotic and abiotic stresses, etc. In Hordeum vulgare and Oryza sativa, 26 and 30 different Dof genes, respectively, have been annotated. Brachypodium distachyon has been the first Pooideae grass to be sequenced and, due to its genomic, morphological and physiological characteristics, has emerged as the model system for temperate cereals, such as wheat and barley.

RESULTS

Through searches in the B. distachyon genome, 27 Dof genes have been identified and a phylogenetic comparison with the Oryza sativa and the Hordeum vulgare DOFs has been performed. To explore the evolutionary relationship among these DOF proteins, a combined phylogenetic tree has been constructed with the Brachypodium DOFs and those from rice and barley. This phylogenetic analysis has classified the DOF proteins into four Major Cluster of Orthologous Groups (MCOGs). Using RT-qPCR analysis the expression profiles of the annotated BdDof genes across four organs (leaves, roots, spikes and seeds) has been investigated. These results have led to a classification of the BdDof genes into two groups, according to their expression levels. The genes highly or preferentially expressed in seeds have been subjected to a more detailed expression analysis (maturation, dry stage and germination).

CONCLUSIONS

Comparison of the expression profiles of the Brachypodium Dof genes with the published functions of closely related DOF sequences from the cereal species considered here, deduced from the phylogenetic analysis, indicates that although the expression profile has been conserved in many of the putative orthologs, in some cases duplication followed by subsequent divergence may have occurred (neo-functionalization).

Genome-wide association mapping of flowering time and northern corn leaf blight (Setosphaeria turcica) resistance in a vast commercial maize germplasm set

BACKGROUND

Setosphaeria turcica is a fungal pathogen that causes northern corn leaf blight (NCLB) which is a serious foliar disease in maize. In order to unravel the genetic architecture of the resistance against this disease, a vast association mapping panel comprising 1487 European maize inbred lines was used to (i) identify chromosomal regions affecting flowering time (FT) and northern corn leaf blight (NCLB) resistance, (ii) examine the epistatic interactions of the identified chromosomal regions with the genetic background on an individual molecular marker basis, and (iii) dissect the correlation between NCLB resistance and FT.

RESULTS

The single marker analyses performed for 8 244 single nucleotide polymorphism (SNP) markers revealed seven, four, and four SNP markers significantly (α=0.05, amplicon wise Bonferroni correction) associated with FT, NCLB, and NCLB resistance corrected for FT, respectively. These markers explained individually between 0.36 and 14.29% of the genetic variance of the corresponding trait.

CONCLUSIONS

The very well interpretable pattern of SNP associations observed for FT suggested that data from applied plant breeding programs can be used to dissect polygenic traits. This in turn indicates that the associations identified for NCLB resistance might be successfully used in marker-assisted selection programs. Furthermore, the associated genes are also of interest for further research concerning the mechanism of resistance to NCLB and plant diseases in general, because some of the associated genes have not been mentioned in this context so far.

Genome-wide association studies for agronomical traits in a world wide spring barley collection

BACKGROUND

Genome-wide association studies (GWAS) based on linkage disequilibrium (LD) provide a promising tool for the detection and fine mapping of quantitative trait loci (QTL) underlying complex agronomic traits. In this study we explored the genetic basis of variation for the traits heading date, plant height, thousand grain weight, starch content and crude protein content in a diverse collection of 224 spring barleys of worldwide origin. The whole panel was genotyped with a customized oligonucleotide pool assay containing 1536 SNPs using Illumina's GoldenGate technology resulting in 957 successful SNPs covering all chromosomes. The morphological trait "row type" (two-rowed spike vs. six-rowed spike) was used to confirm the high level of selectivity and sensitivity of the approach. This study describes the detection of QTL for the above mentioned agronomic traits by GWAS.

RESULTS

Population structure in the panel was investigated by various methods and six subgroups that are mainly based on their spike morphology and region of origin. We explored the patterns of linkage disequilibrium (LD) among the whole panel for all seven barley chromosomes. Average LD was observed to decay below a critical level (r2-value 0.2) within a map distance of 5-10 cM. Phenotypic variation within the panel was reasonably large for all the traits. The heritabilities calculated for each trait over multi-environment experiments ranged between 0.90-0.95. Different statistical models were tested to control spurious LD caused by population structure and to calculate the P-value of marker-trait associations. Using a mixed linear model with kinship for controlling spurious LD effects, we found a total of 171 significant marker trait associations, which delineate into 107 QTL regions. Across all traits these can be grouped into 57 novel QTL and 50 QTL that are congruent with previously mapped QTL positions.

CONCLUSIONS

Our results demonstrate that the described diverse barley panel can be efficiently used for GWAS of various quantitative traits, provided that population structure is appropriately taken into account. The observed significant marker trait associations provide a refined insight into the genetic architecture of important agronomic traits in barley. However, individual QTL account only for a small portion of phenotypic variation, which may be due to insufficient marker coverage and/or the elimination of rare alleles prior to analysis. The fact that the combined SNP effects fall short of explaining the complete phenotypic variance may support the hypothesis that the expression of a quantitative trait is caused by a large number of very small effects that escape detection. Notwithstanding these limitations, the integration of GWAS with biparental linkage mapping and an ever increasing body of genomic sequence information will facilitate the systematic isolation of agronomically important genes and subsequent analysis of their allelic diversity.

NGS technologies for analyzing germplasm diversity in genebanks

More than 70 years after the first ex situ genebanks have been established, major efforts in this field are still concerned with issues related to further completion of individual collections and securing of their storage. Attempts regarding valorization of ex situ collections for plant breeders have been hampered by the limited availability of phenotypic and genotypic information. With the advent of molecular marker technologies first efforts were made to fingerprint genebank accessions, albeit on a very small scale and mostly based on inadequate DNA marker systems. Advances in DNA sequencing technology and the development of high-throughput systems for multiparallel interrogation of thousands of single nucleotide polymorphisms (SNPs) now provide a suite of technological platforms facilitating the analysis of several hundred of Gigabases per day using state-of-the-art sequencing technology or, at the same time, of thousands of SNPs. The present review summarizes recent developments regarding the deployment of these technologies for the analysis of plant genetic resources, in order to identify patterns of genetic diversity, map quantitative traits and mine novel alleles from the vast amount of genetic resources maintained in genebanks around the world. It also refers to the various shortcomings and bottlenecks that need to be overcome to leverage the full potential of high-throughput DNA analysis for the targeted utilization of plant genetic resources.

Haplotyping, linkage mapping and expression analysis of barley genes regulated by terminal drought stress influencing seed quality

BACKGROUND

The increasingly narrow genetic background characteristic of modern crop germplasm presents a challenge for the breeding of cultivars that require adaptation to the anticipated change in climate. Thus, high priority research aims at the identification of relevant allelic variation present both in the crop itself as well as in its progenitors. This study is based on the characterization of genetic variation in barley, with a view to enhancing its response to terminal drought stress.

RESULTS

The expression patterns of drought regulated genes were monitored during plant ontogeny, mapped and the location of these genes was incorporated into a comprehensive barley SNP linkage map. Haplotypes within a set of 17 starch biosynthesis/degradation genes were defined, and a particularly high level of haplotype variation was uncovered in the genes encoding sucrose synthase (types I and II) and starch synthase. The ability of a panel of 50 barley accessions to maintain grain starch content under terminal drought conditions was explored.

CONCLUSION

The linkage/expression map is an informative resource in the context of characterizing the response of barley to drought stress. The high level of haplotype variation among starch biosynthesis/degradation genes in the progenitors of cultivated barley shows that domestication and breeding have greatly eroded their allelic diversity in current elite cultivars. Prospective association analysis based on core drought-regulated genes may simplify the process of identifying favourable alleles, and help to understand the genetic basis of the response to terminal drought.

Allele mining in barley genetic resources reveals genes of race-non-specific powdery mildew resistance

Race-non-specific, or quantitative, pathogen resistance is of high importance to plant breeders due to its expected durability. However, it is usually controlled by multiple quantitative trait loci (QTL) and therefore difficult to handle in practice. Knowing the genes that underlie race-non-specific resistance (NR) would allow its exploitation in a more targeted manner. Here, we performed an association-genetic study in a customized worldwide collection of spring barley accessions for candidate genes of race-NR to the powdery mildew fungus Blumeria graminis f. sp. hordei (Bgh) and combined data with results from QTL mapping as well as functional-genomics approaches. This led to the identification of 11 associated genes with converging evidence for an important role in race-NR in the presence of the Mlo gene for basal susceptibility. Outstanding in this respect was the gene encoding the transcription factor WRKY2. The results suggest that unlocking plant genetic resources and integrating functional-genomic with genetic approaches can accelerate the discovery of genes underlying race-NR in barley and other crop plants.