-
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.
-
The role of the Arabidopsis FUSCA3 transcription factor during inhibition of seed germination at high temperature
Background:
Imbibed seeds integrate environmental and endogenous signals to break dormancy and initiate growth under optimal conditions. Seed maturation plays an important role in determining the survival of germinating seeds, for example one of the roles of dormancy is to stagger germination to prevent mass growth under suboptimal conditions. The B3-domain transcription factor FUSCA3 (FUS3) is a master regulator of seed development and an important node in hormonal interaction networks in Arabidopsis thaliana. Its function has been mainly characterized during embryonic development, where FUS3 is highly expressed to promote seed maturation and dormancy by regulating ABA/GA levels.
Results:
In this study, we present evidence for a role of FUS3 in delaying seed germination at supraoptimal temperatures that would be lethal for the developing seedlings. During seed imbibition at supraoptimal temperature, the FUS3 promoter is reactivated and induces de novo synthesis of FUS3 mRNA, followed by FUS3 protein accumulation. Genetic analysis shows that FUS3 contributes to the delay of seed germination at high temperature. Unlike WT, seeds overexpressing FUS3 (ML1:FUS3-GFP) during imbibition are hypersensitive to high temperature and do not germinate, however, they can fully germinate after recovery at control temperature reaching 90% seedling survival. ML1:FUS3-GFP hypersensitivity to high temperature can be partly recovered in the presence of fluridone, an inhibitor of ABA biosynthesis, suggesting this hypersensitivity is due in part to higher ABA level in this mutant. Transcriptomic analysis shows that WT seeds imbibed at supraoptimal temperature activate seed-specific genes and ABA biosynthetic and signaling genes, while inhibiting genes that promote germination and growth, such as GA biosynthetic and signaling genes.
Conclusion:
In this study, we have uncovered a novel function for the master regulator of seed maturation, FUS3, in delaying germination at supraoptimal temperature. Physiologically, this is important since delaying germination has a protective role at high temperature. Transcriptomic analysis of seeds imbibed at supraoptimal temperature reveal that a complex program is in place, which involves not only the regulation of heat and dehydration response genes to adjust cellular functions, but also the activation of seed-specific programs and the inhibition of germination-promoting programs to delay germination.
-
Combining SNP discovery from next-generation sequencing data with bulked segregant analysis (BSA) to fine-map genes in polyploid wheat
Background:
Next generation sequencing (NGS) technologies are providing new ways to accelerate fine-mapping and gene isolation in many species. To date, the majority of these efforts have focused on diploid organisms with readily available whole genome sequence information. In this study, as a proof of concept, we tested the use of NGS for SNP discovery in tetraploid wheat lines differing for the previously cloned grain protein content (GPC) gene GPC-B1. Bulked segregant analysis (BSA) was used to define a subset of putative SNPs within the candidate gene region, which were then used to fine-map GPC-B1.
Results:
We used Illumina paired end technology to sequence mRNA (RNAseq) from near isogenic lines differing across a ~30-cM interval including the GPC-B1 locus. After discriminating for SNPs between the two homoeologous wheat genomes and additional quality filtering, we identified inter-varietal SNPs in wheat unigenes between the parental lines. The relative frequency of these SNPs was examined by RNAseq in two bulked samples made up of homozygous recombinant lines differing for their GPC phenotype. SNPs that were enriched at least 3-fold in the corresponding pool (6.5% of all SNPs) were further evaluated. Marker assays were designed for a subset of the enriched SNPs and mapped using DNA from individuals of each bulk. Thirty nine new SNP markers, corresponding to 67% of the validated SNPs, mapped across a 12.2-cM interval including GPC-B1. This translated to 1 SNP marker per 0.31 cM defining the GPC-B1 gene to within 13-18 genes in syntenic cereal genomes and to a 0.4 cM interval in wheat.
Conclusions:
This study exemplifies the use of RNAseq for SNP discovery in polyploid species and supports the use of BSA as an effective way to target SNPs to specific genetic intervals to fine-map genes in unsequenced genomes.
-
Characterization of resistance to pine wood nematode infection in Pinus thunbergii using suppression subtractive hybridization
Background:
Pine wilt disease is caused by the pine wood nematode, Bursaphelenchus xylophilus, which threatens pine forests and forest ecosystems worldwide and causes serious economic losses. In the 40 years since the pathogen was identified, the physiological changes occurring as the disease progresses have been characterized using anatomical and biochemical methods, and resistant trees have been selected via breeding programs. However, no studies have assessed the molecular genetics, e.g. transcriptional changes, associated with infection-induced physiological changes in resistant or susceptible trees.
Results:
We constructed seven subtractive suppression hybridization (SSH) cDNA libraries using time-course sampling of trees inoculated with pine wood nematode at 1, 3, or 7 days post-inoculation (dpi) in susceptible trees and at 1, 3, 7, or 14 dpi in resistant trees. A total of 3,299 sequences was obtained from these cDNA libraries, including from 138 to 315 non-redundant sequences in susceptible SSH libraries and from 351 to 435 in resistant SSH libraries. Using Gene Ontology hierarchy, those non-redundant sequences were classified into 15 subcategories of the biological process Gene Ontology category and 17 subcategories of the molecular function category. The transcriptional components revealed by the Gene Ontology classification clearly differed between resistant and susceptible libraries. Some transcripts were discriminative: expression of antimicrobial peptide and putative pathogenesis-related genes (e.g., PR-1b, 2, 3, 4, 5, 6) was much higher in susceptible trees than in resistant trees at every time point, whereas expression of PR-9, PR-10, and cell wall-related genes (e.g., for hydroxyproline-rich glycoprotein precursor and extensin) was higher in resistant trees than in susceptible trees at 7 and 14 dpi.
Conclusions:
Following inoculation with pine wood nematode, there were marked differences between resistant and susceptible trees in transcript diversity and the timing and level of transcripts expressed in common; in particular, expression of stress response and defense genes differed. This study provided new insight into the differences in the physiological changes between resistant and susceptible trees that have been observed in anatomical and biochemical studies.
-
QTL and candidate gene mapping for polyphenolic composition in apple fruit
Background:
The polyphenolic products of the phenylpropanoid pathway, including proanthocyanidins, anthocyanins and flavonols, possess antioxidant properties that may provide health benefits. To investigate the genetic architecture of control of their biosynthesis in apple fruit, various polyphenolic compounds were quantified in progeny from a 'Royal Gala' x 'Braeburn' apple population segregating for antioxidant content, using ultra high performance liquid chromatography of extracts derived from fruit cortex and skin.
Results:
Construction of genetic maps for 'Royal Gala' and 'Braeburn' enabled detection of 79 quantitative trait loci (QTL) for content of 17 fruit polyphenolic compounds. Seven QTL clusters were stable across two years of harvest and included QTLs for content of flavanols, flavonols, anthocyanins and hydroxycinnamic acids. Alignment of the parental genetic maps with the apple whole genome sequence in silico enabled screening for co-segregation with the QTLs of a range of candidate genes coding for enzymes in the polyphenolic biosynthetic pathway. This co-location was confirmed by genetic mapping of markers derived from the gene sequences. Leucoanthocyanidin reductase (LAR1) co-located with a QTL cluster for the fruit flavanols catechin, epicatechin, procyanidin dimer and five unknown procyanidin oligomers identified near the top of linkage group (LG) 16, while hydroxy cinnamate/quinate transferase (HCT/HQT) co-located with a QTL for chlorogenic acid concentration mapping near the bottom of LG 17.
Conclusion:
We conclude that LAR1 and HCT/HQT are likely to influence the concentration of these compounds in apple fruit and provide useful allele-specific markers for marker assisted selection of trees bearing fruit with healthy attributes.
-
Plasma membrane protein OsMCA1 is involved in regulation of hypo-osmotic shock-induced Ca2+ influx and modulates generation of reactive oxygen species in cultured rice cells
Background:
Mechanosensing and its downstream responses are speculated to involve sensory complexes containing Ca2+-permeable mechanosensitive channels. On recognizing osmotic signals, plant cells initiate activation of a widespread signal transduction network that induces second messengers and triggers inducible defense responses. Characteristic early signaling events include Ca2+ influx, protein phosphorylation and generation of reactive oxygen species (ROS). Pharmacological analyses show Ca2+ influx mediated by mechanosensitive Ca2+ channels to influence induction of osmotic signals, including ROS generation. However, molecular bases and regulatory mechanisms for early osmotic signaling events remain poorly elucidated.
Results:
We here identified and investigated OsMCA1, the sole rice homolog of putative Ca2+-permeable mechanosensitive channels in Arabidopsis (MCAs). OsMCA1 was specifically localized at the plasma membrane. A promoter-reporter assay suggested that OsMCA1 mRNA is widely expressed in seed embryos, proximal and apical regions of shoots, and mesophyll cells of leaves and roots in rice. Ca2+ uptake was enhanced in OsMCA1-overexpressing suspension-cultured cells, suggesting that OsMCA1 is involved in Ca2+ influx across the plasma membrane. Hypo-osmotic shock-induced ROS generation mediated by NADPH oxidases was also enhanced in OsMCA1-overexpressing cells. We also generated and characterized OsMCA1-RNAi transgenic plants and cultured cells; OsMCA1-suppressed plants showed retarded growth and shortened rachises, while OsMCA1-suppressed cells carrying Ca2+-sensitive photoprotein aequorin showed partially impaired changes in cytosolic free Ca2+ concentration ([Ca2+]cyt) induced by hypo-osmotic shock and trinitrophenol, an activator of mechanosensitive channels.
Conclusions:
We have identified a sole MCA ortholog in the rice genome and developed both overexpression and suppression lines. Analyses of cultured cells with altered levels of this putative Ca2+-permeable mechanosensitive channel indicate that OsMCA1 is involved in regulation of plasma membrane Ca2+ influx and ROS generation induced by hypo-osmotic stress in cultured rice cells. These findings shed light on our understanding of mechanical sensing pathways.
-
Development and characterization of BAC-end sequence derived SSRs, and their incorporation into a new higher density genetic map for cultivated peanut (Arachis hypogaea L.)
Background:
Cultivated peanut (Arachis hypogaea L.) is an important crop worldwide, valued for its edible oil and digestible protein. It has a very narrow genetic base that may well derive from a relatively recent single polyploidization event. Accordingly molecular markers have low levels of polymorphism and the number of polymorphic molecular markers available for cultivated peanut is still limiting.
Results:
Here, we report a large set of BAC-end sequences (BES), use them for developing SSR (BES-SSR) markers, and apply them in genetic linkage mapping. The majority of BESs had no detectable homology to known genes (49.5%) followed by sequences with similarity to known genes (44.3%), and miscellaneous sequences (6.2%) such as transposable element, retroelement, and organelle sequences. A total of 1,424 SSRs were identified from 36,435 BESs. Among these identified SSRs, dinucleotide (47.4%) and trinucleotide (37.1%) SSRs were predominant. The new set of 1,152 SSRs as well as about 4,000 published or unpublished SSRs were screened against two parents of a mapping population, generating 385 polymorphic loci. A genetic linkage map was constructed, consisting of 318 loci onto 21 linkage groups and covering a total of 1,674.4 cM, with an average distance of 5.3 cM between adjacent loci. Two markers related to resistance gene homologs (RGH) were mapped to two different groups, thus anchoring 1 RGH-BAC contig and 1 singleton.
Conclusions:
The SSRs mined from BESs will be of use in further molecular analysis of the peanut genome, providing a novel set of markers, genetically anchoring BAC clones, and incorporating gene sequences into a linkage map. This will aid in the identification of markers linked to genes of interest and map-based cloning.
-
Integration of gene-based markers in a pearl millet genetic map for identification of candidate genes underlying drought tolerance quantitative trait loci
Background:
Identification of genes underlying drought tolerance (DT) quantitative trait loci (QTLs) will facilitate understanding of molecular mechanisms of drought tolerance, and also will accelerate genetic improvement of pearl millet through marker-assisted selection. We report a map based on genes with assigned functional roles in plant adaptation to drought and other abiotic stresses and demonstrate its use in identifying candidate genes underlying a major DT-QTL.
Results:
Seventy five single nucleotide polymorphism (SNP) and conserved intron spanning primer (CISP) markers were developed from available expressed sequence tags (ESTs) using four genotypes, H 77/833-2, PRLT 2/89-33, ICMR 01029 and ICMR 01004, representing parents of two mapping populations. A total of 228 SNPs were obtained from 30.5 kb sequenced region resulting in a SNP frequency of 1/134 bp. The positions of major pearl millet linkage group (LG) 2 DT-QTLs (reported from crosses H 77/833-2 x PRLT 2/89-33 and 841B x 863B) were added to the present consensus function map which identified 18 genes, coding for PSI reaction center subunit III, PHYC, actin, alanine glyoxylate aminotransferase, uridylate kinase, acyl-CoA oxidase, dipeptidyl peptidase IV, MADS-box, serine/threonine protein kinase, ubiquitin conjugating enzyme, zinc finger C-x8-Cx5-Cx3-H type, Hd3, acetyl CoA carboxylase, chlorophyll a/b binding protein, photolyase, protein phosphatase1 regulatory subunit SDS22 and two hypothetical proteins, co-mapping in this DT-QTL interval. Many of these candidate genes were found to have significant association with QTLs of grain yield, flowering time and leaf rolling under drought stress conditions.
Conclusions:
We have exploited available pearl millet EST sequences to generate a mapped resource of seventy five new gene-based markers for pearl millet and demonstrated its use in identifying candidate genes underlying a major DT-QTL in this species. The reported gene-based markers represent an important resource for identification of candidate genes for other mapped abiotic stress QTLs in pearl millet. They also provide a resource for initiating association studies using candidate genes and also for comparing the structure and function of distantly related plant genomes such as other Poaceae members.
-
Defects in leaf carbohydrate metabolism compromise acclimation to high light and lead to a high chlorophyll fluorescence phenotype in Arabidopsis thaliana
Background:
We have studied the impact of carbohydrate-starvation on the acclimation response to high light using Arabidopsis thaliana double mutants strongly impaired in the day- and night path of photoassimilate export from the chloroplast. A complete knock-out mutant of the triose phosphate/phosphate translocator (TPT; tpt-2 mutant) was crossed to mutants defective in (i) starch biosynthesis (adg1-1, pgm1 and pgi1-1; knock-outs of ADP-glucose pyrophosphorylase, plastidial phosphoglucomutase and phosphoglucose isomerase) or (ii) starch mobilization (sex1-3, knock-out of glucan water dikinase) as well as in (iii) maltose export from the chloroplast (mex1-2).
Results:
All double mutants were viable and indistinguishable from the wild type when grown under low light conditions, but - except for sex1-3/tpt-2 - developed a high chlorophyll fluorescence (HCF) phenotype and growth retardation when grown in high light. Immunoblots of thylakoid proteins, Blue-Native gel electrophoresis and chlorophyll fluorescence emission analyses at 77 Kelvin with the adg1-1/tpt-2 double mutant revealed that HCF was linked to a specific decrease in plastome-encoded core proteins of both photosystems (with the exception of the PSII component cytochrome b559), whereas nuclear-encoded antennae (LHCs) accumulated normally, but were predominantly not attached to their photosystems. Uncoupled antennae are the major cause for HCF of dark-adapted plants. Feeding of sucrose or glucose to high light-grown adg1-1/tpt-2 plants rescued the HCF- and growth phenotypes. Elevated sugar levels induce the expression of the glucose-6-phosphate/phosphate translocator2 (GPT2), which in principle could compensate for the deficiency in the TPT. A triple mutant with an additional defect in GPT2 (adg1-1/tpt-2/gpt2-1) exhibited an identical rescue of the HCF- and growth phenotype in response to sugar feeding as the adg1-1/tpt-2 double mutant, indicating that this rescue is independent from the sugar-triggered induction of GPT2.
Conclusions:
We propose that cytosolic carbohydrate availability modulates acclimation to high light in A. thaliana. It is conceivable that the strong relationship between the chloroplast and nucleus with respect to a co-ordinated expression of photosynthesis genes is modified in carbohydrate-starved plants. Hence carbohydrates may be considered as a novel component involved in chloroplast-to-nucleus retrograde signaling, an aspect that will be addressed in future studies.
-
A genomics approach to understanding the role of auxin in apple (Malus x domestica) fruit size control
Background:
Auxin is an important phytohormone for fleshy fruit development, having been shown to be involved in the initial signal for fertilisation, fruit size through the control of cell division and cell expansion, and ripening related events. There is considerable knowledge of auxin-related genes, mostly from work in model species. With the apple genome now available, it is possible to carry out genomics studies on auxin-related genes to identify genes that may play roles in specific stages of apple fruit development.
Results:
High amounts of auxin in the seed compared with the fruit cortex were observed in 'Royal Gala' apples, with amounts increasing through fruit development. Injection of exogenous auxin into developing apples at the start of cell expansion caused an increase in cell size. An expression analysis screen of auxin-related genes involved in auxin reception, homeostasis, and transcriptional regulation showed complex patterns of expression in each class of gene. Two mapping populations were phenotyped for fruit size over multiple seasons, and multiple quantitative trait loci (QTLs) were observed. One QTL mapped to a region containing an Auxin Response Factor (ARF106). This gene is expressed during cell division and cell expansion stages, consistent with a potential role in the control of fruit size.
Conclusions:
The application of exogenous auxin to apples increased cell expansion, suggesting that endogenous auxin concentrations are at least one of the limiting factors controlling fruit size. The expression analysis of ARF106 linked to a strong QTL for fruit weight suggests that the auxin signal regulating fruit size could partially be modulated through the function of this gene. One class of gene (GH3) removes free auxin by conjugation to amino acids. The lower expression of these GH3 genes during rapid fruit expansion is consistent with the apple maximising auxin concentrations at this point.
|
