Metabolic GWAS‐based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome‐wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass‐spectrometry‐based analytical systems and genome sequencing technologies. Metabolic genome‐wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS‐based multi‐dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS‐based multi‐dimensional analysis and emerging insights into the diversity of metabolism.     

Genetic dissection of drought and heat‐responsive agronomic traits in wheat

High yield and wide adaptation are principal targets of wheat breeding but are hindered by limited knowledge on genetic basis of agronomic traits and abiotic stress tolerances. In this study, 277 wheat accessions were phenotyped across 30 environments with non‐stress, drought‐stressed, heat‐stressed, and drought‐heat‐stressed treatments and were subjected to genome‐wide association study using 395 681 single nucleotide polymorphisms. We detected 295 associated loci including consistent loci for agronomic traits across different treatments and eurytopic loci for multiple abiotic stress tolerances. A total of 22 loci overlapped with quantitative trait loci identified by biparental quantitative trait loci mapping. Six loci were simultaneously associated with agronomic traits and abiotic stress tolerance, four of which fell within selective sweep regions. Selection in Chinese wheat has increased the frequency of superior marker alleles controlling yield‐related traits in the four loci during past decades, which conversely diminished favourable genetic variation controlling abiotic stress tolerance in the same loci; two promising candidate paralogous genes colocalized with such loci, thereby providing potential targets for studying the molecular mechanism of stress tolerance–productivity trade‐off. These results uncovering promising alleles controlling agronomic traits and/or multiple abiotic stress tolerances, providing insights into heritable covariation between yield and abiotic stress tolerance, will accelerate future efforts for wheat improvement.     

A genome‐wide association study of 23 agronomic traits in Chinese wheat landraces

Uncovering the genetic basis of agronomic traits in wheat landraces is important for ensuring global food security via the development of improved varieties. Here, 723 wheat landraces from 10 Chinese agro‐ecological zones were evaluated for 23 agronomic traits in six environments. All accessions could be clustered into five subgroups based on phenotypic data via discriminant function analysis, which was highly consistent with genotypic classification. A genome‐wide association study was conducted for these traits using 52 303 DArT‐seq markers to identify marker‐trait associations and candidate genes. Using both the general linear model and the mixed linear model, 149 significant markers were identified for 21 agronomic traits based on best linear unbiased prediction values. Considering the linkage disequilibrium decay distance in this study, significant markers within 10 cM were combined as a quantitative trait locus (QTL), with a total of 29 QTL identified for 15 traits. Of these, five QTL for heading date, flag leaf width, peduncle length, and thousand kernel weight had been reported previously. Twenty‐five candidate genes associated with significant markers were identified. These included the known vernalization genes VRN‐B1 and vrn‐B3 and the photoperiod response genes Ppd and PRR. Overall, this study should be helpful in elucidating the underlying genetic mechanisms of complex agronomic traits and performing marker‐assisted selection in wheat.     

Combining genome‐wide association study and FST‐based approaches to identify targets of Borrelia‐mediated selection in natural rodent hosts

Recent advances in high‐throughput sequencing technologies provide opportunities to gain novel insights into the genetic basis of phenotypic trait variation. Yet to date, progress in our understanding of genotype–phenotype associations in nonmodel organisms in general and natural vertebrate populations in particular has been hampered by small sample sizes typically available for wildlife populations and a resulting lack of statistical power, as well as a limited ability to control for false‐positive signals. Here we propose to combine a genome‐wide association study (GWAS) and F_ST‐based approach with population‐level replication to partly overcome these limitations. We present a case study in which we used this approach in combination with genotyping‐by‐sequencing (GBS) single nucleotide polymorphism (SNP) data to identify genomic regions associated with Borrelia afzelii resistance or susceptibility in the natural rodent host of this Lyme disease‐causing spirochete, the bank vole (Myodes glareolus). Using this combined approach we identified four consensus SNPs located in exonic regions of the genes Slc26a4, Tns3, Wscd1 and Espnl, which were significantly associated with the voles’ Borrelia infectious status within and across populations. Functional links between host responses to bacterial infections and most of these genes have previously been demonstrated in other rodent systems, making them promising new candidates for the study of evolutionary host responses to Borrelia emergence. Our approach is applicable to other systems and may facilitate the identification of genetic variants underlying disease resistance or susceptibility, as well as other ecologically relevant traits, in wildlife populations.     

Low‐gluten,nontransgenic wheat engineered with CRISPR/Cas9

Coeliac disease is an autoimmune disorder triggered in genetically predisposed individuals by the ingestion of gluten proteins from wheat, barley and rye. The α‐gliadin gene family of wheat contains four highly stimulatory peptides, of which the 33‐mer is the main immunodominant peptide in patients with coeliac. We designed two sgRNAs to target a conserved region adjacent to the coding sequence for the 33‐mer in the α‐gliadin genes. Twenty‐one mutant lines were generated, all showing strong reduction in α‐gliadins. Up to 35 different genes were mutated in one of the lines of the 45 different genes identified in the wild type, while immunoreactivity was reduced by 85%. Transgene‐free lines were identified, and no off‐target mutations have been detected in any of the potential targets. The low‐gluten, transgene‐free wheat lines described here could be used to produce low‐gluten foodstuff and serve as source material to introgress this trait into elite wheat varieties.     

iTRAQ‐coupled 2‐D LC‐MS/MS analysis of protein profile associated with HBV‐modulated DNA methylation

The development of hepatocellular carcinoma (HCC) is believed to be associated with multiple risk factors, including the infection of hepatitis B virus (HBV). Based on the analysis of individual genes, evidence has indicated the association between HCC and HBV and has also been expanded to epigenetic regulation, with an involvement of HBV in the DNA methylation of the promoter of cellular target genes leading to changes in their expression. Proteomic study has been widely used to map a comprehensive protein profile, which in turn could provide a better understanding of underlying mechanisms of disease onset. In the present study, we performed a proteomic profiling by using iTRAQ‐coupled 2‐D LC/MS‐MS analysis to identify cellular genes down‐regulated in HBV‐producing HepG2.2.15 cells compared with HepG2 cells. A total of 15 proteins including S100A6 and Annexin A2 were identified by our approach. The significance of these cellular proteins as target of HBV‐mediated epigenetic regulation was supported by our validation assays, including their reactivation in cells treated with 5‐aza‐2′‐deoxycytidine (a DNA methyltransferase inhibitor) by real‐time RT‐PCR and Western blot analysis, as well as the DNA methylation status analysis by bisulfite genome sequencing. Our approach provides a comprehensive analysis of cellular target proteins to HBV‐mediated epigenetic regulation and further analysis should facilitate a better understanding of its involvement in HCC development.     

Metabolomics analysis and metabolite‐agronomic trait associations using kernels of wheat (Triticum aestivum) recombinant inbred lines

Plants produce numerous metabolites that are important for their development and growth. However, the genetic architecture of the wheat metabolome has not been well studied. Here, utilizing a high‐density genetic map, we conducted a comprehensive metabolome study via widely targeted LC‐MS/MS to analyze the wheat kernel metabolism. We further combined agronomic traits and dissected the genetic relationship between metabolites and agronomic traits. In total, 1260 metabolic features were detected. Using linkage analysis, 1005 metabolic quantitative trait loci (mQTLs) were found distributed unevenly across the genome. Twenty‐four candidate genes were found to modulate the levels of different metabolites, of which two were functionally annotated by in vitro analysis to be involved in the synthesis and modification of flavonoids. Combining the correlation analysis of metabolite‐agronomic traits with the co‐localization of methylation quantitative trait locus (mQTL) and phenotypic QTL (pQTL), genetic relationships between the metabolites and agronomic traits were uncovered. For example, a candidate was identified using correlation and co‐localization analysis that may manage auxin accumulation, thereby affecting number of grains per spike (NGPS). Furthermore, metabolomics data were used to predict the performance of wheat agronomic traits, with metabolites being found that provide strong predictive power for NGPS and plant height. This study used metabolomics and association analysis to better understand the genetic basis of the wheat metabolism which will ultimately assist in wheat breeding.     

Testing a ‘genes‐to‐ecosystems’ approach to understanding aquatic–terrestrial linkages

A ‘genes‐to‐ecosystems’ approach has been proposed as a novel avenue for integrating the consequences of intraspecific genetic variation with the underlying genetic architecture of a species to shed light on the relationships among hierarchies of ecological organization (genes → individuals → communities → ecosystems). However, attempts to identify genes with major effect on the structure of communities and/or ecosystem processes have been limited and a comprehensive test of this approach has yet to emerge. Here, we present an interdisciplinary field study that integrated a common garden containing different genotypes of a dominant, riparian tree, Populus trichocarpa, and aquatic mesocosms to determine how intraspecific variation in leaf litter alters both terrestrial and aquatic communities and ecosystem functioning. Moreover, we incorporate data from extensive trait screening and genome‐wide association studies estimating the heritability and genes associated with litter characteristics. We found that tree genotypes varied considerably in the quality and production of leaf litter, which contributed to variation in phytoplankton abundances, as well as nutrient dynamics and light availability in aquatic mesocosms. These ‘after‐life’ effects of litter from different genotypes were comparable to the responses of terrestrial communities associated with the living foliage. We found that multiple litter traits corresponding with aquatic community and ecosystem responses differed in their heritability. Moreover, the underlying genetic architecture of these traits was complex, and many genes contributed only a small proportion to phenotypic variation. Our results provide further evidence that genetic variation is a key component of aquatic–terrestrial linkages, but challenge the ability to predict community or ecosystem responses based on the actions of one or a few genes.     

High‐density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool

In wheat, a lack of genetic diversity between breeding lines has been recognized as a significant block to future yield increases. Species belonging to bread wheat's secondary and tertiary gene pools harbour a much greater level of genetic variability, and are an important source of genes to broaden its genetic base. Introgression of novel genes from progenitors and related species has been widely employed to improve the agronomic characteristics of hexaploid wheat, but this approach has been hampered by a lack of markers that can be used to track introduced chromosome segments. Here, we describe the identification of a large number of single nucleotide polymorphisms that can be used to genotype hexaploid wheat and to identify and track introgressions from a variety of sources. We have validated these markers using an ultra‐high‐density Axiom^® genotyping array to characterize a range of diploid, tetraploid and hexaploid wheat accessions and wheat relatives. To facilitate the use of these, both the markers and the associated sequence and genotype information have been made available through an interactive web site.     

The genetics underlying metabolic signatures in a brown rice diversity panel and their vital role in human nutrition

Brown rice (Oryza sativa) possesses various nutritionally dense bioactive phytochemicals exhibiting a wide range of antioxidant, anti-cancer, and anti-diabetic properties known to promote various human health benefits. However, despite the wide claims made about the importance of brown rice for human nutrition the underlying metabolic diversity has not been systematically explored. Non-targeted metabolite profiling of developing and mature seeds of a diverse genetic panel of 320 rice cultivars allowed quantification of 117 metabolites. The metabolite genome-wide association study (mGWAS) detected genetic variants influencing diverse metabolic targets in developing and mature seeds. We further interlinked genetic variants on chromosome 7 (6.06–6.43 Mb region) with complex epistatic genetic interactions impacting multi-dimensional nutritional targets, including complex carbohydrate starch quality, the glycemic index, antioxidant catechin, and rice grain color. Through this nutrigenomics approach rare gene bank accessions possessing genetic variants in bHLH and IPT5 genes were identified through haplotype enrichment. These variants were associated with a low glycemic index, higher catechin levels, elevated total flavonoid contents, and heightened antioxidant activity in the whole grain with elevated anti-cancer properties being confirmed in cancer cell lines. This multi-disciplinary nutrigenomics approach thus allowed us to discover the genetic basis of human health-conferring diversity in the metabolome of brown rice.     

Discovery and development of exome‐based,co‐dominant single nucleotide polymorphism markers in hexaploid wheat (Triticum aestivum L.)

Globally, wheat is the most widely grown crop and one of the three most important crops for human and livestock feed. However, the complex nature of the wheat genome has, until recently, resulted in a lack of single nucleotide polymorphism (SNP)‐based molecular markers of practical use to wheat breeders. Recently, large numbers of SNP‐based wheat markers have been made available via the use of next‐generation sequencing combined with a variety of genotyping platforms. However, many of these markers and platforms have difficulty distinguishing between heterozygote and homozygote individuals and are therefore of limited use to wheat breeders carrying out commercial‐scale breeding programmes. To identify exome‐based co‐dominant SNP‐based assays, which are capable of distinguishing between heterozygotes and homozygotes, we have used targeted re‐sequencing of the wheat exome to generate large amounts of genomic sequences from eight varieties. Using a bioinformatics approach, these sequences have been used to identify 95 266 putative single nucleotide polymorphisms, of which 10 251 were classified as being putatively co‐dominant. Validation of a subset of these putative co‐dominant markers confirmed that 96% were true polymorphisms and 65% were co‐dominant SNP assays. The new co‐dominant markers described here are capable of genotypic classification of a segregating locus in polyploid wheat and can be used on a variety of genotyping platforms; as such, they represent a powerful tool for wheat breeders. These markers and related information have been made publically available on an interactive web‐based database to facilitate their use on genotyping programmes worldwide.     

Roles of mouse sperm‐associated alpha‐L‐fucosidases in fertilization

Sperm‐associated α‐L ‐fucosidases have been implicated in fertilization in many species. Previously, we documented the existence of α‐L ‐fucosidase in mouse cauda epididymal contents, and showed that sperm‐associated α‐L ‐fucosidase is cryptically stored within the acrosome and reappears within the sperm equatorial segment after the acrosome reaction. The enrichment of sperm membrane‐associated α‐L ‐fucosidase within the equatorial segment of acrosome‐reacted cells implicates its roles during fertilization. Here, we document the absence of α‐L ‐fucosidase in mouse oocytes and early embryos, and define roles of sperm associated α‐L ‐fucosidase in fertilization using specific inhibitors and competitors. Mouse sperm were pretreated with deoxyfuconojirimycin (DFJ, an inhibitor of α‐L ‐fucosidase) or with anti‐fucosidase antibody; alternatively, mouse oocytes were pretreated with purified human liver α‐L ‐fucosidase. Five‐millimolar DFJ did not inhibit sperm–zona pellucida (ZP) binding, membrane binding, or fusion and penetration, but anti‐fucosidase antibody and purified human liver α‐L ‐fucosidase significantly decreased the frequency of these events. To evaluate sperm‐associated α‐L ‐fucosidase enzyme activity in post‐fusion events, DFJ‐pretreated sperm were microinjected into oocytes, and 2‐pronuclear (2‐PN) embryos were treated with 5 mM DFJ with no significant effects, suggesting that α‐L ‐fucosidase enzyme activity does not play a role in post‐fusion events and/or early embryo development in mice. The recognition and binding of mouse sperm to the ZP and oolemma involves the glycoprotein structure of α‐L ‐fucosidase, but not its catalytic action. These observations suggest that deficits in fucosidase protein and/or the presence of anti‐fucosidase antibody may be responsible for some types of infertility. Mol. Reprod. Dev. 80: 273–285, 2013. © 2013 Wiley Periodicals, Inc.     

Molecular Characterization and Expression Analysis of Triticum aestivum Squamosa‐Promoter Binding Protein‐Box Genes Involved in Ear Development

Wheat (Triticum aestivum L.) is one of the most important crops in the world. Squamosa‐promoter binding protein (SBP)‐box genes play a critical role in regulating flower and fruit development. In this study, 10 novel SBP‐box genes (TaSPL genes) were isolated from wheat ((Triticum aestivum L.) cultivar Yanzhan 4110). Phylogenetic analysis classified the TaSPL genes into five groups (G1–G5). The motif combinations and expression patterns of the TaSPL genes varied among the five groups with each having own distinctive characteristics: TaSPL20/21 in G1 and TaSPL17 in G2 mainly expressed in the shoot apical meristem and the young ear, and their expression levels responded to development of the ear; TaSPL6/15 belonging to G3 were upregulated and TaSPL1/23 in G4 were downregulated during grain development; the gene in G5 (TaSPL3) expressed constitutively. Thus, the consistency of the phylogenetic analysis, motif compositions, and expression patterns of the TaSPL genes revealed specific gene structures and functions. On the other hand, the diverse gene structures and different expression patterns suggested that wheat SBP‐box genes have a wide range of functions. The results also suggest a potential role for wheat SBP‐box genes in ear development. This study provides a significant beginning of functional analysis of SBP‐box genes in wheat.     

ADP‐glucose pyrophosphorylase genes,associated with kernel weight,underwent selection during wheat domestication and breeding

ADP‐glucose pyrophosphorylase, comprising two small subunits and two large subunits, is considered a key enzyme in the endosperm starch synthesis pathway in wheat (Triticum aestivum L.). Two genes, TaAGP‐S1‐7A and TaAGP‐L‐1B, were investigated in this study. Haplotypes of these genes were associated with thousand kernel weight (TKW) in different populations. Mean TKWs of favoured haplotypes were significantly higher than those of nonfavoured ones. Two molecular markers developed to distinguish these haplotypes could be used in molecular breeding. Frequencies of favoured haplotypes were dramatically increased in cultivars released in China after the 1940s. These favoured haplotypes were also positively selected in six major wheat production regions globally. Selection of AGP‐S1 and AGP‐L‐1B in wheat mainly occurred during and after hexaploidization. Strong additive effects of the favoured haplotypes of with other genes for starch synthesis were also detected in different populations.     

Genome‐wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome‐wide association study (GWAS) of the six quality‐related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment‐significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1‐B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1‐B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker‐assisted selection to improve the process of breeding quality wheat in China.