An eQTL analysis of partial resistance to Puccinia hordei in barley

Molecular mechanisms underpinning quantitative resistance to Phytophthora sojae in Glycine max using a systems genomics approach

Expression of quantitative disease resistance in many host-pathogen systems is controlled by genes at multiple loci, each contributing a small effect to the overall response. We used a systems genomics approach to study the molecular underpinnings of quantitative disease resistance in the soybean-Phytophthora sojae pathosystem, incorporating expression quantitative trait loci (eQTL) mapping and gene co-expression network analysis to identify the genes putatively regulating transcriptional changes in response to inoculation. These findings were compared to previously mapped phenotypic (phQTL) to identify the molecular mechanisms contributing to the expression of this resistance. A subset of 93 recombinant inbred lines (RILs) from a Conrad × Sloan population were inoculated with P. sojae isolate 1.S.1.1 using the tray-test method; RNA was extracted, sequenced, and the normalized read counts were genetically mapped from tissue collected at the inoculation site 24 h after inoculation from both mock and inoculated samples. In total, more than 100,000 eQTLs were mapped. There was a switch from predominantly cis-eQTLs in the mock treatment to an almost entirely nonoverlapping set of predominantly trans-eQTLs in the inoculated treatment, where greater than 100-fold more eQTLs were mapped relative to mock, indicating vast transcriptional reprogramming due to P. sojae infection occurred. The eQTLs were organized into 36 hotspots, with the four largest hotspots from the inoculated treatment corresponding to more than 70% of the eQTLs, each enriched for genes within plant-pathogen interaction pathways. Genetic regulation of trans-eQTLs in response to the pathogen was predicted to occur through transcription factors and signaling molecules involved in plant-pathogen interactions, plant hormone signal transduction, and MAPK pathways. Network analysis identified three co-expression modules that were correlated with susceptibility to P. sojae and associated with three eQTL hotspots. Among the eQTLs co-localized with phQTLs, two cis-eQTLs with putative functions in the regulation of root architecture or jasmonic acid, as well as the putative master regulators of an eQTL hotspot nearby a phQTL, represent candidates potentially underpinning the molecular control of these phQTLs for resistance.

Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis-Botrytis pathosystem

Bidirectional flow of information shapes the outcome of the host-pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens' genetic variation and both organisms' transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host's genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.

Pathogen Genetic Control of Transcriptome Variation in the Arabidopsis thaliana - Botrytis cinerea Pathosystem

In plant-pathogen relations, disease symptoms arise from the interaction of the host and pathogen genomes. Host-pathogen functional gene interactions are well described, whereas little is known about how the pathogen genetic variation modulates both organisms' transcriptomes. To model and generate hypotheses on a generalist pathogen control of gene expression regulation, we used the Arabidopsis thaliana-Botrytis cinerea pathosystem and the genetic diversity of a collection of 96 B. cinerea isolates. We performed expression-based genome-wide association (eGWA) for each of 23,947 measurable transcripts in Arabidopsis (host), and 9267 measurable transcripts in B. cinerea (pathogen). Unlike other eGWA studies, we detected a relative absence of locally acting expression quantitative trait loci (cis-eQTL), partly caused by structural variants and allelic heterogeneity hindering their identification. This study identified several distantly acting trans-eQTL linked to eQTL hotspots dispersed across Botrytis genome that altered only Botrytis transcripts, only Arabidopsis transcripts, or transcripts from both species. Gene membership in the trans-eQTL hotspots suggests links between gene expression regulation and both known and novel virulence mechanisms in this pathosystem. Genes annotated to these hotspots provide potential targets for blocking manipulation of the host response by this ubiquitous generalist necrotrophic pathogen.

Transcriptomic and presence/absence variation in the barley genome assessed from multi-tissue mRNA sequencing and their power to predict phenotypic traits

Background

Barley is the world’s fourth most cultivated cereal and is an important crop model for genetic studies. One layer of genomic information that remains poorly explored in barley is presence/absence variation (PAV), which has been suggested to contribute to phenotypic variation of agronomic importance in various crops.

Results

An mRNA sequencing approach was used to study genomic PAV and transcriptomic variation in 23 spring barley inbreds. 1502 new genes identified here were physically absent from the Morex reference sequence, and 11,523 previously unannotated genes were not expressed in Morex. The procedure applied to detect expression PAV revealed that more than 50% of all genes of our data set are not expressed in all inbreds. Interestingly, expression PAV were not in strong linkage disequilibrium with neighboring sequence variants (SV), and therefore provided an additional layer of genetic information. Optimal combinations of expression PAV, SV, and gene abundance data could enhance the prediction accuracy of predicting three different agronomic traits.

Conclusions

Our results highlight the advantage of mRNA sequencing for genomic prediction over other technologies, as it allows extracting multiple layers of genomic data from a single sequencing experiment. Finally, we propose low coverage mRNA sequencing based characterization of breeding material harvested as seedlings in petri dishes as a powerful and cost efficient approach to replace current single nucleotide polymorphism (SNP) based characterizations.

Development of Genome-Wide SNP Markers for Barley via Reference- Based RNA-Seq Analysis

Marker-assisted selection of crop plants requires DNA markers that can distinguish between the closely related strains often used in breeding. The availability of reference genome sequence facilitates the generation of markers, by elucidating the genomic positions of new markers as well as of their neighboring sequences. In 2017, a high quality genome sequence was released for the six-row barley (Hordeum vulgare) cultivar Morex. Here, we developed a de novo RNA-Seq-based genotyping procedure for barley strains used in Japanese breeding programs. Using RNA samples from the seedling shoot, seedling root, and immature flower spike, we mapped next-generation sequencing reads onto the transcribed regions, which correspond to ∼590 Mb of the whole ∼4.8-Gbp reference genome sequence. Using 150 samples from 108 strains, we detected 181,567 SNPs and 45,135 indels located in the 28,939 transcribed regions distributed throughout the Morex genome. We evaluated the quality of this polymorphism detection approach by analyzing 387 RNA-Seq-derived SNPs using amplicon sequencing. More than 85% of the RNA-Seq SNPs were validated using the highly redundant reads from the amplicon sequencing, although half of the indels and multiple-allele loci showed different polymorphisms between the platforms. These results demonstrated that our RNA-Seq-based de novo polymorphism detection system generates genome-wide markers, even in the closely related barley genotypes used in breeding programs.

OsSYP121 Accumulates at Fungal Penetration Sites and Mediates Host Resistance to Rice Blast

Magnaporthe oryzae is a fungal pathogen that causes rice (Oryza sativa) blast. SNAREs (soluble N-ethylmaleimide-sensitive factor attachment protein receptors) are key components in vesicle trafficking in eukaryotic cells and are known to contribute to fungal pathogen resistance. Syntaxin of Plants121 (SYP121), a Qa-SNARE, has been reported to function in nonhost resistance in Arabidopsis (Arabidopsis thaliana). However, the functions of SYP121 in host resistance to rice blast are largely unknown. Here, we report that the rice SYP121 protein, OsSYP121, accumulates at fungal penetration sites and mediates host resistance to rice blast. OsSYP121 is plasma membrane localized and its expression was obviously induced by the rice blast in both the blast-resistant rice landrace Heikezijing and the blast-susceptible landrace Suyunuo (Su). Overexpression of OsSYP121 in Su resulted in enhanced resistance to blast. Knockdown of OsSYP121 expression in Su resulted in a more susceptible phenotype. However, knockdown of OsSYP121 expression in the resistant landrace Heikezijing resulted in susceptibility to the blast fungus. The P_OsSYP121 ::GFP-OsSYP121 accumulated at rice blast penetration sites in transgenic rice, as observed by confocal microscopy. Yeast two-hybrid results showed that OsSYP121 can interact with OsSNAP32 (Synaptosome-associated protein of 32 kD) and Vesicle-associated membrane protein714/724. The interaction between OsSYP121 and OsSNAP32 may contribute to host resistance to rice blast. Our study reveals that OsSYP121 plays an important role in rice blast resistance as it is a key component in vesicle trafficking.

Triticale Improvement for Forage and Cover Crop Uses in the Southern Great Plains of the United States

Triticale (×Triticosecale Wittmack) is a man-made species developed by crossing wheat (Triticum spp.) and rye (Secale cereale L.). It incorporates favorable alleles from both progenitor species (wheat and rye), enabling adaptation to environments that are less favorable for wheat yet providing better biomass yield and forage quality. Triticale has huge potential for both grain and forage production, though research to improve the crop for better adaptation and grain quality is lagging behind that of other small grains. It is also gaining popularity as a cover crop to improve soil health and reduce nutrient leaching. Because of its genetic and flower structure, triticale is suitable for both line and hybrid breeding methods. Advances in the areas of molecular biology and the wealth of genomic resources from both wheat and rye can be exploited for triticale improvement. Gene mapping and genomic selection will facilitate triticale breeding by increasing selection precision and reducing time and cost. The objectives of this review are to summarize current triticale production status, breeding, and genetics research achievements and to highlight gaps for future research.

Characterization of Botrytis cinerea isolates collected on pepper in Southern Turkey by using molecular markers, fungicide resistance genes and virulence assay

Botrytis cinerea is a polyphagous fungal pathogen causing gray mold disease. Moreover, it is one of the most destructive infections of small fruit crops such as pepper (Capsicum annnum L.). C. sativum is a species belonging to the Solanaceae family and Turkey is one of the main producers in the World. In the present work, aiming to obtain information useful for pest management, fifty B. cinerea isolates collected from Turkey and a reference isolate (B05.10) were characterized using molecular markers and fungicide resistance genes. Morphological and molecular (ITS1-ITS4) identification of B. cinerea isolates, the degree of virulence and mating types were determined. Since one or several allelic mutations in the histidine kinase (Bos1) and β-tubulin genes generally confer the resistance to fungicides, the sequences of these target genes were investigated in the selected isolates, which allowed the identification of two different haplotypes. Mating types were also determined by PCR assays using primer specific for MAT1-1 alpha gene (MAT1-1-1) and MAT1-2 HMG (MAT1-2-1) of B. cinerea. Twenty-two out of 50 isolates (44%) were MAT1-2, while 38% were MAT1-1. Interestingly, out of whole studied samples, 9 isolates (18%) were heterokaryotic or mixed colonies. In addition, cluster and population structure analyses identified five main groups and two genetic pools, respectively, underlining a good level of variability in the analysed panel. The results highlighted the presence of remarkable genetic diversity in B. cinerea isolates collected in a crucial economical area for pepper cultivation in Turkey and the data will be beneficial in view of future gray mold disease management.

Verticillium dahliae Disease Resistance and the Regulatory Pathway for Maturity and Tuberization in Potato

Kleb. is a pathogenic fungus causing wilting, chlorosis, and early dying in potato ( L.). Genetic mapping of resistance to was done using a diploid population of potato. The major quantitative trait locus (QTL) for resistance was found on chromosome 5. The gene, controlling earliness of maturity and tuberization, was mapped within the interval. Another QTL on chromosome 9 co-localized with the wilt resistance gene marker. Epistasis analysis indicated that the loci on chromosomes 5 and 9 had a highly significant interaction, and that functioned downstream of The alleles were sequenced and found to encode StCDF1.1 and StCDF1.3. Interaction between the resistance allele and the was demonstrated, but not for Genome-wide expression QTL (eQTL) analysis was performed and genes with eQTL at the and loci were both found to have similar functions involving the chloroplast, including photosynthesis, which declines in both maturity and wilt. Among the gene ontology (GO) terms that were specific to genes with eQTL at the , but not the locus, were those associated with fungal defense. These results suggest that controls fungal defense and reduces early dying in wilt through affecting genetic pathway controlling tuberization timing.

Genetical genomics of quality related traits in potato tubers using proteomics

BACKGROUND

Recent advances in ~omics technologies such as transcriptomics, metabolomics and proteomics along with genotypic profiling have permitted the genetic dissection of complex traits such as quality traits in non-model species. To get more insight into the genetic factors underlying variation in quality traits related to carbohydrate and starch metabolism and cold sweetening, we determined the protein content and composition in potato tubers using 2D-gel electrophoresis in a diploid potato mapping population. Upon analyzing we made sure that the proteins from the patatin family were excluded to ensure a better representation of the other proteins.

RESULTS

We subsequently performed pQTL analyses for all other proteins with a sufficient representation in the population and established a relationship between proteins and 26 potato tuber quality traits (e.g. flesh colour, enzymatic discoloration) by co-localization on the genetic map and a direct correlation study of protein abundances and phenotypic traits. Over 1643 unique protein spots were detected in total over the two harvests. We were able to map pQTLs for over 300 different protein spots some of which co-localized with traits such as starch content and cold sweetening. pQTLs were observed on every chromosome although not evenly distributed over the chromosomes. The largest number of pQTLs was found for chromosome 8 and the lowest for chromosome number 10. For some 20 protein spots multiple QTLs were observed.

CONCLUSIONS

From this analysis, hotspot areas for protein QTLs were identified on chromosomes three, five, eight and nine. The hotspot on chromosome 3 coincided with a QTL previously identified for total protein content and had more than 23 pQTLs in the region from 70 to 80 cM. Some of the co-localizing protein spots associated with some of the most interesting tuber quality traits were identified, albeit far less than we had anticipated at the onset of the experiments.

RNA-Seq analysis of resistant and susceptible sub-tropical maize lines reveals a role for kauralexins in resistance to grey leaf spot disease, caused by Cercospora zeina

BACKGROUND

Cercospora zeina is a foliar pathogen responsible for maize grey leaf spot in southern Africa that negatively impacts maize production. Plants use a variety of chemical and structural mechanisms to defend themselves against invading pathogens such as C. zeina, including the production of secondary metabolites with antimicrobial properties. In maize, a variety of biotic and abiotic stressors induce the accumulation of the terpenoid phytoalexins, zealexins and kauralexins.

RESULTS

C. zeina-susceptible line displayed pervasive rectangular grey leaf spot lesions, running parallel with the leaf veins in contrast to C. zeina-resistant line that had restricted disease symptoms. Analysis of the transcriptome of both lines indicated that genes involved in primary and secondary metabolism were up-regualted, and although different pathways were prioritized in each line, production of terpenoid compounds were common to both. Targeted phytoalexin analysis revealed that C. zeina-inoculated leaves accumulated zealexins and kauralexins. The resistant line shows a propensity toward accumulation of the kauralexin B series metabolites in response to infection, which contrasts with the susceptible line that preferentially accumulates the kauralexin A series. Kauralexin accumulation was correlated to expression of the kauralexin biosynthetic gene, ZmAn2 and a candidate biosynthetic gene, ZmKSL2. We report the expression of a putative copalyl diphosphate synthase gene that is induced by C. zeina in the resistant line exclusively.

DISCUSSION

This study shows that zealexins and kauralexins, and expression of their biosynthetic genes, are induced by C. zeina in both resistant and susceptible germplasm adapted to the southern African climate. The data presented here indicates that different forms of kauralexins accumulate in the resistant and susceptible maize lines in response to C. zeina, with the accumulation of kauralexin B compounds in a resistant maize line and kauralexin A compounds accumulating in the susceptible line.

Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum

Fusarium head blight is a devastating disease of small grain cereals such as bread wheat (Triticum aestivum). The pathogen switches from a biotrophic to a nectrotrophic lifestyle in course of disease development forcing its host to adapt its defence strategies. Using a genetical genomics approach, we illustrate genome-wide reconfigurations of genetic control over transcript abundances between two decisive time points after inoculation with the causative pathogen Fusarium graminearum. Whole transcriptome measurements have been recorded for 163 lines of a wheat doubled haploid population segregating for several resistance genes yielding 15 552 at 30 h and 15 888 eQTL at 50 h after inoculation. The genetic map saturated with transcript abundance-derived markers identified of a novel QTL on chromosome 6A, besides the previously reported QTL Fhb1 and Qfhs.ifa-5A. We find a highly different distribution of eQTL between time points with about 40% of eQTL being unique for the respective assessed time points. But also for more than 20% of genes governed by eQTL at either time point, genetic control changes in time. These changes are reflected in the dynamic compositions of three major regulatory hotspots on chromosomes 2B, 4A and 5A. In particular, control of defence-related biological mechanisms concentrated in the hotspot at 4A shift to hotspot 2B as the disease progresses. Hotspots do not colocalize with phenotypic QTL, and within their intervals no higher than expected number of eQTL was detected. Thus, resistance conferred by either QTL is mediated by few or single genes.

Interchromosomal Transfer of Immune Regulation During Infection of Barley with the Powdery Mildew Pathogen

Powdery mildew pathogens colonize over 9500 plant species, causing critical yield loss. The Ascomycete fungus, Blumeria graminis f. sp. hordei (Bgh), causes powdery mildew disease in barley (Hordeum vulgare L.). Successful infection begins with penetration of host epidermal cells, culminating in haustorial feeding structures, facilitating delivery of fungal effectors to the plant and exchange of nutrients from host to pathogen. We used expression Quantitative Trait Locus (eQTL) analysis to dissect the temporal control of immunity-associated gene expression in a doubled haploid barley population challenged with Bgh. Two highly significant regions possessing trans eQTL were identified near the telomeric ends of chromosomes (Chr) 2HL and 1HS. Within these regions reside diverse resistance loci derived from barley landrace H. laevigatum (MlLa) and H. vulgare cv. Algerian (Mla1), which associate with the altered expression of 961 and 3296 genes during fungal penetration of the host and haustorial development, respectively. Regulatory control of transcript levels for 299 of the 961 genes is reprioritized from MlLa on 2HL to Mla1 on 1HS as infection progresses, with 292 of the 299 alternating the allele responsible for higher expression, including Adaptin Protein-2 subunit μ AP2M and Vesicle Associated Membrane Protein VAMP72 subfamily members VAMP721/722. AP2M mediates effector-triggered immunity (ETI) via endocytosis of plasma membrane receptor components. VAMP721/722 and SNAP33 form a Soluble N-ethylmaleimide-sensitive factor Attachment Protein REceptor (SNARE) complex with SYP121 (PEN1), which is engaged in pathogen associated molecular pattern (PAMP)-triggered immunity via exocytosis. We postulate that genes regulated by alternate chromosomal positions are repurposed as part of a conserved immune complex to respond to different pathogen attack scenarios.

Quantitative Resistance: More Than Just Perception of a Pathogen

Molecular plant pathology has focused on studying large-effect qualitative resistance loci that predominantly function in detecting pathogens and/or transmitting signals resulting from pathogen detection. By contrast, less is known about quantitative resistance loci, particularly the molecular mechanisms controlling variation in quantitative resistance. Recent studies have provided insight into these mechanisms, showing that genetic variation at hundreds of causal genes may underpin quantitative resistance. Loci controlling quantitative resistance contain some of the same causal genes that mediate qualitative resistance, but the predominant mechanisms of quantitative resistance extend beyond pathogen recognition. Indeed, most causal genes for quantitative resistance encode specific defense-related outputs such as strengthening of the cell wall or defense compound biosynthesis. Extending previous work on qualitative resistance to focus on the mechanisms of quantitative resistance, such as the link between perception of microbe-associated molecular patterns and growth, has shown that the mechanisms underlying these defense outputs are also highly polygenic. Studies that include genetic variation in the pathogen have begun to highlight a potential need to rethink how the field considers broad-spectrum resistance and how it is affected by genetic variation within pathogen species and between pathogen species. These studies are broadening our understanding of quantitative resistance and highlighting the potentially vast scale of the genetic basis of quantitative resistance.

High-resolution mapping of genes involved in plant stage-specific partial resistance of barley to leaf rust

Partial resistance quantitative trait loci (QTLs) Rphq11 and rphq16 against Puccinia hordei isolate 1.2.1 were previously mapped in seedlings of the mapping populations Steptoe/Morex and Oregon Wolfe Barleys, respectively. In this study, QTL mapping was performed at adult plant stage for the two mapping populations challenged with the same rust isolate. The results suggest that Rphq11 and rphq16 are effective only at seedling stage, and not at adult plant stage. The cloning of several genes responsible for partial resistance of barley to P. hordei will allow elucidation of the molecular basis of this type of plant defence. A map-based cloning approach requires to fine-map the QTL in a narrow genetic window. In this study, Rphq11 and rphq16 were fine-mapped using an approach aiming at speeding up the development of plant material and simplifying its evaluation. The plant materials for fine-mapping were identified from early plant materials developed to produce QTL-NILs. The material was first selected to carry the targeted QTL in heterozygous condition and susceptibility alleles at other resistance QTLs in homozygous condition. This strategy took four to five generations to obtain fixed QTL recombinants (i.e., homozygous resistant at the Rphq11 or rphq16 QTL alleles, homozygous susceptible at the non-targeted QTL alleles). In less than 2 years, Rphq11 was fine-mapped into a 0.2-cM genetic interval and a 1.4-cM genetic interval for rphq16. The strongest candidate gene for Rphq11 is a phospholipid hydroperoxide glutathione peroxidase. Thus far, no candidate gene was identified for rphq16.

Systems genetics reveals a transcriptional network associated with susceptibility in the maize-grey leaf spot pathosystem

We used a systems genetics approach to elucidate the molecular mechanisms of the responses of maize to grey leaf spot (GLS) disease caused by Cercospora zeina, a threat to maize production globally. Expression analysis of earleaf samples in a subtropical maize recombinant inbred line population (CML444 × SC Malawi) subjected in the field to C. zeina infection allowed detection of 20 206 expression quantitative trait loci (eQTLs). Four trans-eQTL hotspots coincided with GLS disease QTLs mapped in the same field experiment. Co-expression network analysis identified three expression modules correlated with GLS disease scores. The module (GY-s) most highly correlated with susceptibility (r = 0.71; 179 genes) was enriched for the glyoxylate pathway, lipid metabolism, diterpenoid biosynthesis and responses to pathogen molecules such as chitin. The GY-s module was enriched for genes with trans-eQTLs in hotspots on chromosomes 9 and 10, which also coincided with phenotypic QTLs for susceptibility to GLS. This transcriptional network has significant overlap with the GLS susceptibility response of maize line B73, and may reflect pathogen manipulation for nutrient acquisition and/or unsuccessful defence responses, such as kauralexin production by the diterpenoid biosynthesis pathway. The co-expression module that correlated best with resistance (TQ-r; 1498 genes) was enriched for genes with trans-eQTLs in hotspots coinciding with GLS resistance QTLs on chromosome 9. Jasmonate responses were implicated in resistance to GLS through co-expression of COI1 and enrichment of genes with the Gene Ontology term 'cullin-RING ubiquitin ligase complex' in the TQ-r module. Consistent with this, JAZ repressor expression was highly correlated with the severity of GLS disease in the GY-s susceptibility network.

eQTLs Regulating Transcript Variations Associated with Rapid Internode Elongation in Deepwater Rice

To avoid low oxygen, oxygen deficiency or oxygen deprivation, deepwater rice cultivated in flood planes can develop elongated internodes in response to submergence. Knowledge of the gene regulatory networks underlying rapid internode elongation is important for an understanding of the evolution and adaptation of major crops in response to flooding. To elucidate the genetic and molecular basis controlling their deepwater response we used microarrays and performed expression quantitative trait loci (eQTL) and phenotypic QTL (phQTL) analyses of internode samples of 85 recombinant inbred line (RIL) populations of non-deepwater (Taichung 65)- and deepwater rice (Bhadua). After evaluating the phenotypic response of the RILs exposed to submergence, confirming the genotypes of the populations, and generating 188 genetic markers, we identified 10,047 significant eQTLs comprised of 2,902 cis-eQTLs and 7,145 trans-eQTLs and three significant eQTL hotspots on chromosomes 1, 4, and 12 that affect the expression of many genes. The hotspots on chromosomes 1 and 4 located at different position from phQTLs detected in this study and other previous studies. We then regarded the eQTL hotspots as key regulatory points to infer causal regulatory networks of deepwater response including rapid internode elongation. Our results suggest that the downstream regulation of the eQTL hotspots on chromosomes 1 and 4 is independent, and that the target genes are partially regulated by SNORKEL1 and SNORKEL2 genes (SK1/2), key ethylene response factors. Subsequent bioinformatic analyses, including gene ontology-based annotation and functional enrichment analysis and promoter enrichment analysis, contribute to enhance our understanding of SK1/2-dependent and independent pathways. One remarkable observation is that the functional categories related to photosynthesis and light signaling are significantly over-represented in the candidate target genes of SK1/2. The combined results of these investigations together with genetical genomics approaches using structured populations with a deepwater response are also discussed in the context of current molecular models concerning the rapid internode elongation in deepwater rice. This study provides new insights into the underlying genetic architecture of gene expression regulating the response to flooding in deepwater rice and will be an important community resource for analyses on the genetic basis of deepwater responses.

Genome-Editing Technologies for Enhancing Plant Disease Resistance

One of the greatest challenges for agricultural science in the 21st century is to improve yield stability through the progressive development of superior cultivars. The increasing numbers of infectious plant diseases that are caused by plant-pathogens make it ever more necessary to develop new strategies for plant disease resistance breeding. Targeted genome engineering allows the introduction of precise modifications directly into a commercial variety, offering a viable alternative to traditional breeding methods. Genome editing is a powerful tool for modifying crucial players in the plant immunity system. In this work, we propose and discuss genome-editing strategies and targets for improving resistance to phytopathogens. First of all, we present the opportunities to rewrite the effector-target sequence for avoiding effector-target molecular interaction and also to modify effector-target promoters for increasing the expression of target genes involved in the resistance process. In addition, we describe potential approaches for obtaining synthetic R-genes through genome-editing technologies (GETs). Finally, we illustrate a genome editing flowchart to modify the pathogen recognition sites and engineer an R-gene that mounts resistance to some phylogenetically divergent pathogens. GETs potentially mark the beginning of a new era, in which synthetic biology affords a basis for obtaining a reinforced plant defense system. Nowadays it is conceivable that by modulating the function of the major plant immunity players, we will be able to improve crop performance for a sustainable agriculture.

eQTL Regulating Transcript Levels Associated with Diverse Biological Processes in Tomato

Variation in gene expression, in addition to sequence polymorphisms, is known to influence developmental, physiological, and metabolic traits in plants. Genetic mapping populations have facilitated identification of expression quantitative trait loci (eQTL), the genetic determinants of variation in gene expression patterns. We used an introgression population developed from the wild desert-adapted Solanum pennellii and domesticated tomato (Solanum lycopersicum) to identify the genetic basis of transcript level variation. We established the effect of each introgression on the transcriptome and identified approximately 7,200 eQTL regulating the steady-state transcript levels of 5,300 genes. Barnes-Hut t-distributed stochastic neighbor embedding clustering identified 42 modules revealing novel associations between transcript level patterns and biological processes. The results showed a complex genetic architecture of global transcript abundance pattern in tomato. Several genetic hot spots regulating a large number of transcript level patterns relating to diverse biological processes such as plant defense and photosynthesis were identified. Important eQTL regulating transcript level patterns were related to leaf number and complexity as well as hypocotyl length. Genes associated with leaf development showed an inverse correlation with photosynthetic gene expression, but eQTL regulating genes associated with leaf development and photosynthesis were dispersed across the genome. This comprehensive eQTL analysis details the influence of these loci on plant phenotypes and will be a valuable community resource for investigations on the genetic effects of eQTL on phenotypic traits in tomato.

The Quantitative Basis of the Arabidopsis Innate Immune System to Endemic Pathogens Depends on Pathogen Genetics

The most established model of the eukaryotic innate immune system is derived from examples of large effect monogenic quantitative resistance to pathogens. However, many host-pathogen interactions involve many genes of small to medium effect and exhibit quantitative resistance. We used the Arabidopsis-Botrytis pathosystem to explore the quantitative genetic architecture underlying host innate immune system in a population of Arabidopsis thaliana. By infecting a diverse panel of Arabidopsis accessions with four phenotypically and genotypically distinct isolates of the fungal necrotroph B. cinerea, we identified a total of 2,982 genes associated with quantitative resistance using lesion area and 3,354 genes associated with camalexin production as measures of the interaction. Most genes were associated with resistance to a specific Botrytis isolate, which demonstrates the influence of pathogen genetic variation in analyzing host quantitative resistance. While known resistance genes, such as receptor-like kinases (RLKs) and nucleotide-binding site leucine-rich repeat proteins (NLRs), were found to be enriched among associated genes, they only account for a small fraction of the total genes associated with quantitative resistance. Using publically available co-expression data, we condensed the quantitative resistance associated genes into co-expressed gene networks. GO analysis of these networks implicated several biological processes commonly connected to disease resistance, including defense hormone signaling and ROS production, as well as novel processes, such as leaf development. Validation of single gene T-DNA knockouts in a Col-0 background demonstrate a high success rate (60%) when accounting for differences in environmental and Botrytis genetic variation. This study shows that the genetic architecture underlying host innate immune system is extremely complex and is likely able to sense and respond to differential virulence among pathogen genotypes.

Exploiting Differential Gene Expression and Epistasis to Discover Candidate Genes for Drought-Associated QTLs in Arabidopsis thaliana

Soil water availability represents one of the most important selective agents for plants in nature and the single greatest abiotic determinant of agricultural productivity, yet the genetic bases of drought acclimation responses remain poorly understood. Here, we developed a systems-genetic approach to characterize quantitative trait loci (QTLs), physiological traits and genes that affect responses to soil moisture deficit in the TSUxKAS mapping population of Arabidopsis thaliana. To determine the effects of candidate genes underlying QTLs, we analyzed gene expression as a covariate within the QTL model in an effort to mechanistically link markers, RNA expression, and the phenotype. This strategy produced ranked lists of candidate genes for several drought-associated traits, including water use efficiency, growth, abscisic acid concentration (ABA), and proline concentration. As a proof of concept, we recovered known causal loci for several QTLs. For other traits, including ABA, we identified novel loci not previously associated with drought. Furthermore, we documented natural variation at two key steps in proline metabolism and demonstrated that the mitochondrial genome differentially affects genomic QTLs to influence proline accumulation. These findings demonstrate that linking genome, transcriptome, and phenotype data holds great promise to extend the utility of genetic mapping, even when QTL effects are modest or complex.

Quantitative resistance to biotrophic filamentous plant pathogens: concepts, misconceptions, and mechanisms

Quantitative resistance (QR) refers to a resistance that is phenotypically incomplete and is based on the joined effect of several genes, each contributing quantitatively to the level of plant defense. Often, QR remains durably effective, which is the primary driver behind the interest in it. The various terms that are used to refer to QR, such as field resistance, adult plant resistance, and basal resistance, reflect the many properties attributed to it. In this article, we discuss aspects connected to those attributions, in particular the hypothesis that much of the QR to biotrophic filamentous pathogens is basal resistance, i.e., poor suppression of PAMP-triggered defense by effectors. We discuss what role effectors play in suppressing defense or improving access to nutrients. Based on the functions of the few plant proteins identified as involved in QR, vesicle trafficking and protein/metabolite transportation are likely to be common physiological processes relevant to QR.

Expression quantitative trait loci infer the regulation of isoflavone accumulation in soybean (Glycine max L. Merr.) seed

BACKGROUND

Mapping expression quantitative trait loci (eQTL) of targeted genes represents a powerful and widely adopted approach to identify putative regulatory variants. Linking regulation differences to specific genes might assist in the identification of networks and interactions. The objective of this study is to identify eQTL underlying expression of four gene families encoding isoflavone synthetic enzymes involved in the phenylpropanoid pathway, which are phenylalanine ammonia-lyase (PAL; EC 4.3.1.5), chalcone synthase (CHS; EC 2.3.1.74), 2-hydroxyisoflavanone synthase (IFS; EC1.14.13.136) and flavanone 3-hydroxylase (F3H; EC 1.14.11.9). A population of 130 recombinant inbred lines (F5:11), derived from a cross between soybean cultivar 'Zhongdou 27' (high isoflavone) and 'Jiunong 20' (low isoflavone), and a total of 194 simple sequence repeat (SSR) markers were used in this study. Overlapped loci of eQTLs and phenotypic QTLs (pQTLs) were analyzed to identify the potential candidate genes underlying the accumulation of isoflavone in soybean seed.

RESULTS

Thirty three eQTLs (thirteen cis-eQTLs and twenty trans-eQTLs) underlying the transcript abundance of the four gene families were identified on fifteen chromosomes. The eQTLs between Satt278-Sat_134, Sat_134-Sct_010 and Satt149-Sat_234 underlie the expression of both IFS and CHS genes. Five eQTL intervals were overlapped with pQTLs. A total of eleven candidate genes within the overlapped eQTL and pQTL were identified.

CONCLUSIONS

These results will be useful for the development of marker-assisted selection to breed soybean cultivars with high or low isoflavone contents and for map-based cloning of new isoflavone related genes.

A first step toward the development of a barley NAM population and its utilization to detect QTLs conferring leaf rust seedling resistance

We suggest multi-parental nested association mapping as a valuable innovation in barley genetics, which increases the power to map quantitative trait loci and assists in extending genetic diversity of the elite barley gene pool. Plant genetic resources are a key asset to further improve crop species. The nested association mapping (NAM) approach was introduced to identify favorable genes in multi-parental populations. Here, we report toward the development of the first explorative barley NAM population and demonstrate its usefulness in a study on mapping quantitative trait loci (QTLs) for leaf rust resistance. The NAM population HEB-5 was developed from crossing and backcrossing five exotic barley donors with the elite barley cultivar 'Barke,' resulting in 295 NAM lines in generation BC1S1. HEB-5 was genetically characterized with 1,536 barley SNPs. Across HEB-5 and within the NAM families, no deviation from the expected genotype and allele frequencies was detected. Genetic similarity between 'Barke' and the NAM families ranged from 78.6 to 83.1 %, confirming the backcrossing step during population development. To explore its usefulness, a screen for leaf rust (Puccinia hordei) seedling resistance was conducted. Resistance QTLs were mapped to six barley chromosomes, applying a mixed model genome-wide association study. In total, four leaf rust QTLs were detected across HEB-5 and four QTLs within family HEB-F23. Favorable exotic QTL alleles reduced leaf rust symptoms on two chromosomes by 33.3 and 36.2 %, respectively. The located QTLs may represent new resistance loci or correspond to new alleles of known resistance genes. We conclude that the exploratory population HEB-5 can be applied to mapping and utilizing exotic QTL alleles of agronomic importance. The NAM concept will foster the evaluation of the genetic diversity, which is present in our primary barley gene pool.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Selection of candidate genes for grape proanthocyanidin pathway by an integrative approach

Proanthocyanidins (PA) play a major role in plant protection against biotic and abiotic stresses. Moreover these molecules are known to be beneficial for human health and are responsible for astringency of foods and beverages such as wine and thus have a great impact on the final quality of the product. Genes playing a role in the PA pathway are only partially known. The amount of available transcriptomic and genetic data to select candidate genes without a priori knowledge from orthologous function increases every day. However, the methods used so far generate so many candidate genes that it is impossible to validate all of them. In this study, we used an integrative strategy based on different screening methods to select a reduced list of candidate genes. We have crossed results from different screening methods including QTL mapping and three transcriptomic studies to select 20 candidate genes, located in QTL intervals and fulfilling at least two transcriptomic screenings. This list includes three glucosyltransferases, already suspected to have a role in the PA biosynthetic pathway. Among the 17 remaining genes, we selected three genes to perform further analysis by association genetic studies. For each of these genes, we found a polymorphism linked to PA variation. The three genes (VvMybC2-L1, VvGAT-like and VvCob-like), not previously known to play a role in PA synthesis, are promising candidates for further molecular physiology studies.

Transgenic barley: a prospective tool for biotechnology and agriculture

Barley (Hordeum vulgare L.) is one of the founder crops of agriculture, and today it is the fourth most important cereal grain worldwide. Barley is used as malt in brewing and distilling industry, as an additive for animal feed, and as a component of various food and bread for human consumption. Progress in stable genetic transformation of barley ensures a potential for improvement of its agronomic performance or use of barley in various biotechnological and industrial applications. Recently, barley grain has been successfully used in molecular farming as a promising bioreactor adapted for production of human therapeutic proteins or animal vaccines. In addition to development of reliable transformation technologies, an extensive amount of various barley genetic resources and tools such as sequence data, microarrays, genetic maps, and databases has been generated. Current status on barley transformation technologies including gene transfer techniques, targets, and progeny stabilization, recent trials for improvement of agricultural traits and performance of barley, especially in relation to increased biotic and abiotic stress tolerance, and potential use of barley grain as a protein production platform have been reviewed in this study. Overall, barley represents a promising tool for both agricultural and biotechnological transgenic approaches, and is considered an ancient but rediscovered crop as a model industrial platform for molecular farming.

Silicon era of carbon-based life: application of genomics and bioinformatics in crop stress research

Abiotic and biotic stresses lead to massive reprogramming of different life processes and are the major limiting factors hampering crop productivity. Omics-based research platforms allow for a holistic and comprehensive survey on crop stress responses and hence may bring forth better crop improvement strategies. Since high-throughput approaches generate considerable amounts of data, bioinformatics tools will play an essential role in storing, retrieving, sharing, processing, and analyzing them. Genomic and functional genomic studies in crops still lag far behind similar studies in humans and other animals. In this review, we summarize some useful genomics and bioinformatics resources available to crop scientists. In addition, we also discuss the major challenges and advancements in the "-omics" studies, with an emphasis on their possible impacts on crop stress research and crop improvement.

Prospects for advancing defense to cereal rusts through genetical genomics

Rusts are one of the most severe threats to cereal crops because new pathogen races emerge regularly, resulting in infestations that lead to large yield losses. In 1999, a new race of stem rust, Puccinia graminis f. sp. tritici (Pgt TTKSK or Ug99), was discovered in Uganda. Most of the wheat and barley cultivars grown currently worldwide are susceptible to this new race. Pgt TTKSK has already spread northward into Iran and will likely spread eastward throughout the Indian subcontinent in the near future. This scenario is not unique to stem rust; new races of leaf rust (Puccinia triticina) and stripe rust (Puccinia striiformis) have also emerged recently. One strategy for countering the persistent adaptability of these pathogens is to stack complete- and partial-resistance genes, which requires significant breeding efforts in order to reduce deleterious effects of linkage drag. These varied resistance combinations are typically more difficult for the pathogen to defeat, since they would be predicted to apply lower selection pressure. Genetical genomics or expression Quantitative Trait Locus (eQTL) analysis enables the identification of regulatory loci that control the expression of many to hundreds of genes. Integrated deployment of these technologies coupled with efficient phenotyping offers significant potential to elucidate the regulatory nodes in genetic networks that orchestrate host defense responses. The focus of this review will be to present advances in genetical genomic experimental designs and analysis, particularly as they apply to the prospects for discovering partial disease resistance alleles in cereals.

A major quantitative trait locus for cold-responsive gene expression is linked to frost-resistance gene Fr-A2 in common wheat

Low temperature induces expression of Cor (cold-responsive)/Lea (late embryogenesis-abundant) gene family members through C-repeat binding factor (CBF) transcription factors in common wheat. However, the relationship between the genetic loci controlling cold-responsive gene expression and freezing tolerance is unclear. In expression quantitative trait locus (eQTL) analysis, accumulated transcripts of Cor/Lea and CBF genes were quantified in recombinant inbred lines derived from a cross between two common wheat cultivars with different levels of freezing tolerance. Four eQTLs controlling five cold-responsive genes were found, and the major eQTL with the greatest effect was located on the long arm of chromosome 5A. At least the 1D and 5A eQTLs played important roles in development of freezing tolerance in common wheat. The chromosomal location of the 5A eQTL, controlling four cold-responsive genes, coincided with a region homoeologous to a frost-tolerance locus (Fr-A (m) 2) reported as a CBF cluster region in einkorn wheat. The 5A eQTL plays a significant role through Cor/Lea gene expression in cold acclimation of wheat. In addition, our results suggest that one or more CBF copies at the Fr-2 region positively regulate other copies, which might amplify the positive effects of the CBF cluster on downstream Cor/Lea gene activation.

Proteomic characterization of the Rph15 barley resistance gene-mediated defence responses to leaf rust

Background

Leaf rust, caused by the biotrophic fungal pathogen Puccinia hordei, is one of the most important foliar disease of barley (Hordeum vulgare) and represents a serious threat in many production regions of the world. The leaf rust resistance gene Rph15 is of outstanding interest for resistance breeding because it confers resistance to over 350 Puccinia hordei isolates collected from around the world. Molecular and biochemical mechanisms responsible for the Rph15 effectiveness are currently not investigated. The aim of the present work was to study the Rph15-based defence responses using a proteomic approach.

Results

Protein pattern changes in response to the leaf rust pathogen infection were investigated in two barley near isogenic lines (NILs), Bowman (leaf rust susceptible) and Bowman-Rph15 (leaf rust resistant), differing for the introgression of the leaf rust resistance gene Rph15. Two infection time points, 24 hours and four days post inoculation (dpi), were analysed. No statistically significant differences were identified at the early time point, while at 4 dpi eighteen protein spots were significantly up or down regulated with a fold-change equal or higher than two in response to pathogen infection. Almost all the pathogen-responsive proteins were identified in the Bowman-Rph15 resistant NIL. Protein spots were characterized by LC-MS/MS analysis and found to be involved in photosynthesis and energy metabolism, carbohydrate metabolism, protein degradation and defence. Proteomic data were complemented by transcriptional analysis of the respective genes. The identified proteins can be related to modulation of the photosynthetic apparatus components, re-direction of the metabolism to sustain defence responses and deployment of defence proteins.

Conclusions

The identification of leaf rust infection-modulated defence responses restricted to the resistant NIL support the hypothesis that basal defence responses of Bowman, but not the Rph15 resistance gene-based ones, are suppressed or delayed by pathogen effectors to levels below the detection power of the adopted proteomic approach. Additionally, Rph15-mediated resistance processes identified mainly resides on a modulation of primary metabolism, affecting photosyntesis and carbohydrate pool.

Identification of positive and negative regulators of disease resistance to rice blast fungus using constitutive gene expression patterns

Elevated constitutive expression of components of the defence arsenal is associated with quantitative resistance to the rice blast fungus, a phenomenon called preformed defence. While the role of many disease regulators in inducible defence systems has been extensively studied, little attention has been paid so far to genes that regulate preformed defence. In this study, we show by microarray analysis across rice diversity that the preformed defence phenomenon impacts on a large number of defence-related genes without apparently affecting other biological processes. Using a guilt-by-association strategy, we identified two positive regulators that promote constitutive expression of known defence markers and partial resistance to rice blast. The HSF23 gene encodes for a putative member of the heat shock transcription factor family, while CaMBP encodes for a putative Calmodulin-binding protein. Both HSF23 and CaMBP strongly affect preformed defence and also plant growth. Additionally, we identified the OB-fold gene as a negative regulator of blast resistance, which could be involved in RNA stabilization. The OB-fold mutants do not suffer from obvious developmental defects. Taken together, our results prove that our strategy of combining analysis of gene expression diversity with guilt-by-association is a powerful way to identify disease resistance regulators in rice.

TaWIR1 contributes to post-penetration resistance to Magnaporthe oryzae, but not Blumeria graminis f. sp. tritici, in wheat

Members of the Wheat-Induced Resistance 1 (TaWIR1) gene family are highly induced in response to a wide range of pathogens. Homologues have been identified in barley, but not in Brachypodium, whereas, in rice, only distant WIR1 candidates are known. Phylogenetic analysis placed TaWIR1a and TaWIR1b within a distinct clade of wheat transcripts, whereas TaWIR1c clustered with HvWIR1 genes. Transcripts of all three TaWIR1 genes were strongly induced by a wheat-adapted isolate of Magnaporthe oryzae. Virus-induced gene silencing of the TaWIR1 gene family had no effect on the initial penetration of epidermal cells by M. oryzae. However, following the establishment of an infection site, the fungus was able to grow more extensively within the leaf tissue, relative to control leaves, indicating a role for the TaWIR1 gene family in the cell-to-cell movement of M. oryzae. In contrast, the silencing of TaWIR1 transcripts had no effect on epidermal cell penetration by a wheat-adapted isolate of Blumeria graminis, or on the subsequent growth of hyphae. Differential transcription of TaWIR1 genes was also seen in epidermal peels, relative to the remaining leaf tissue, following inoculation with M. oryzae.

Origin of leaf rust adult plant resistance gene Rph20 in barley

Rph20 is the only reported, simply inherited gene conferring moderate to high levels of adult plant resistance (APR) to leaf rust ( Puccinia hordei Otth) in barley ( Hordeum vulgare L.). Key parental genotypes were examined to determine the origin of Rph20 in two-rowed barley. The Dutch cultivar ‘Vada’ (released in the 1950s) and parents, ‘Hordeum laevigatum’ and ‘Gull’ (‘Gold’), along with the related cultivar ‘Emir’ (a derivative of ‘Delta’), were assessed for APR to P. hordei in a disease screening nursery. The marker bPb-0837-PCR, co-located with Rph20 on the short arm of chromosome 5H (5HS), was used to screen genotypes for the resistance allele, Rph20.ai. Results from phenotypic assessment and DNA analysis confirmed that Rph20 originated from the landrace ‘H. laevigatum’ (i.e., Hordeum vulgare subsp. vulgare). Tracing back this gene through the pedigrees of two-rowed barley cultivars, indicated that Rph20 has contributed APR to P. hordei for more than 60 years. Although there have been no reports of an Rph20-virulent pathotype, the search for alternative sources of APR should continue to avoid widespread reliance upon a single resistance factor.

Lessons from eQTL mapping studies: non-coding regions and their role behind natural phenotypic variation in plants

Even if considerable progress has been achieved towards the understanding of natural variation in plant systems, the contribution of transcript abundance variation to phenotypic diversity remains unappreciated. Over the last decade, efforts to characterise the genome-wide expression variation in natural accessions, structured populations and hybrids have improved our knowledge of the contribution of non-coding polymorphisms to gene expression regulation. Moreover, new studies are helping to unravel the role of expression polymorphisms and their orchestrated performance. Recent advances involving classical linkage analysis, GWAS and improved eQTL mapping strategies will provide a greater resolution to determine the genetic variants shaping the broad diversity in plant systems.

cDNA-AFLP-based genetical genomics in cotton fibers

Genetical genomics, or genetic analysis applied to gene expression data, has not been widely used in plants. We used quantitative cDNA-AFLP to monitor the variation in the expression level of cotton fiber transcripts among a population of inter-specific Gossypium hirsutum × G. barbadense recombinant inbred lines (RILs). Two key fiber developmental stages, elongation (10 days post anthesis, dpa), and secondary cell wall thickening (22 dpa), were studied. Normalized intensity ratios of 3,263 and 1,201 transcript-derived fragments (TDFs) segregating over 88 RILs were analyzed for quantitative trait loci (QTL) mapping for the 10 and 22 dpa fibers, respectively. Two-thirds of all TDFs mapped between 1 and 6 eQTLs (LOD > 3.5). Chromosome 21 had a higher density of eQTLs than other chromosomes in both data sets and, within chromosomes, hotspots of presumably trans-acting eQTLs were identified. The eQTL hotspots were compared to the location of phenotypic QTLs for fiber characteristics among the RILs, and several cases of co-localization were detected. Quantitative RT-PCR for 15 sequenced TDFs showed that 3 TDFs had at least one eQTL at a similar location to those identified by cDNA-AFLP, while 3 other TDFs mapped an eQTL at a similar location but with opposite additive effect. In conclusion, cDNA-AFLP proved to be a cost-effective and highly transferable platform for genome-wide and population-wide gene expression profiling. Because TDFs are anonymous, further validation and interpretation (in silico analysis, qPCR gene profiling) of the eQTL and eQTL hotspots will be facilitated by the increasing availability of cDNA and genomic sequence resources in cotton.

Genetical genomics for evolutionary studies

Genetical genomics combines acquired high-throughput genomic data with genetic analysis. In this chapter, we discuss the application of genetical genomics for evolutionary studies, where new high-throughput molecular technologies are combined with mapping quantitative trait loci (QTL) on the genome in segregating populations.The recent explosion of high-throughput data--measuring thousands of proteins and metabolites, deep sequencing, chromatin, and methyl-DNA immunoprecipitation--allows the study of the genetic variation underlying quantitative phenotypes, together termed xQTL. At the same time, mining information is not getting easier. To deal with the sheer amount of information, powerful statistical tools are needed to analyze multidimensional relationships. In the context of evolutionary computational biology, a well-designed experiment may help dissect a complex evolutionary trait using proven statistical methods for associating phenotypical variation with genomic locations.Evolutionary expression QTL (eQTL) studies of the last years focus on gene expression adaptations, mapping the gene expression landscape, and, tentatively, eQTL networks. Here, we discuss the possibility of introducing an evolutionary prior, in the form of gene families displaying evidence of positive selection, and using that in the context of an eQTL experiment for elucidating host-pathogen protein-protein interactions. Through the example of an experimental design, we discuss the choice of xQTL platform, analysis methods, and scope of results. The resulting eQTL can be matched, resulting in putative interacting genes and their regulators. In addition, a prior may help distinguish QTL causality from reactivity, or independence of traits, by creating QTL networks.

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast- and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.

Quantitative and qualitative stem rust resistance factors in barley are associated with transcriptional suppression of defense regulons

Stem rust (Puccinia graminis f. sp. tritici; Pgt) is a devastating fungal disease of wheat and barley. Pgt race TTKSK (isolate Ug99) is a serious threat to these Triticeae grain crops because resistance is rare. In barley, the complex Rpg-TTKSK locus on chromosome 5H is presently the only known source of qualitative resistance to this aggressive Pgt race. Segregation for resistance observed on seedlings of the Q21861 × SM89010 (QSM) doubled-haploid (DH) population was found to be predominantly qualitative, with little of the remaining variance explained by loci other than Rpg-TTKSK. In contrast, analysis of adult QSM DH plants infected by field inoculum of Pgt race TTKSK in Njoro, Kenya, revealed several additional quantitative trait loci that contribute to resistance. To molecularly characterize these loci, Barley1 GeneChips were used to measure the expression of 22,792 genes in the QSM population after inoculation with Pgt race TTKSK or mock-inoculation. Comparison of expression Quantitative Trait Loci (eQTL) between treatments revealed an inoculation-dependent expression polymorphism implicating Actin depolymerizing factor3 (within the Rpg-TTKSK locus) as a candidate susceptibility gene. In parallel, we identified a chromosome 2H trans-eQTL hotspot that co-segregates with an enhancer of Rpg-TTKSK-mediated, adult plant resistance discovered through the Njoro field trials. Our genome-wide eQTL studies demonstrate that transcript accumulation of 25% of barley genes is altered following challenge by Pgt race TTKSK, but that few of these genes are regulated by the qualitative Rpg-TTKSK on chromosome 5H. It is instead the chromosome 2H trans-eQTL hotspot that orchestrates the largest inoculation-specific responses, where enhanced resistance is associated with transcriptional suppression of hundreds of genes scattered throughout the genome. Hence, the present study associates the early suppression of genes expressed in this host-pathogen interaction with enhancement of R-gene mediated resistance.

Combining genetical genomics and bulked segregant analysis-based differential expression: an approach to gene localization

Positional gene isolation in unsequenced species generally requires either a reference genome sequence or an inference of gene content based on conservation of synteny with a genomic model. In the large unsequenced genomes of the Triticeae cereals the latter, i.e. conservation of synteny with the rice and Brachypodium genomes, provides a powerful proxy for establishing local gene content and order. However, efficient exploitation of conservation of synteny requires 'homology bridges' between the model genome and the target region that contains a gene of interest. As effective homology bridges are generally the sequences of genetically mapped genes, increasing the density of these genes around a target locus is an important step in the process. We used bulked segregant analysis (BSA) of transcript abundance data to identify genes located in a specific region of the barley genome. The approach is valuable because only a relatively small proportion of barley genes are currently placed on a genetic map. We analyzed eQTL datasets from the reference Steptoe × Morex doubled haploid population and showed a strong association between differential gene expression and cis-regulation, with 83% of differentially expressed genes co-locating with their eQTL. We then performed BSA by assembling allele-specific pools based on the genotypes of individuals at the partial resistance QTL Rphq11. BSA identified a total of 411 genes as differentially expressed, including HvPHGPx, a gene previously identified as a promising candidate for Rphq11. The genetic location of 276 of these genes could be determined from both eQTL datasets and conservation of synteny, and 254 (92%) of these were located on the target chromosome. We conclude that the identification of differential expression by BSA constitutes a novel method to identify genes located in specific regions of interest. The datasets obtained from such studies provide a robust set of candidate genes for the analysis and serve as valuable resources for targeted marker development and comparative mapping with other grass species.

Transcriptional responses of Italian ryegrass during interaction with Xanthomonas translucens pv. graminis reveal novel candidate genes for bacterial wilt resistance

Xanthomonas translucens pv. graminis (Xtg) causes bacterial wilt, a severe disease of forage grasses such as Italian ryegrass (Lolium multiflorum Lam.). In order to gain a more detailed understanding of the genetic control of resistance mechanisms and to provide prerequisites for marker assisted selection, the partial transcriptomes of two Italian ryegrass genotypes, one resistant and one susceptible to bacterial wilt were compared at four time points after Xtg infection. A cDNA microarray developed from a perennial ryegrass (Lolium perenne) expressed sequence tag set consisting of 9,990 unique genes was used for transcriptome analysis in Italian ryegrass. An average of 4,487 (45%) of the perennial ryegrass sequences spotted on the cDNA microarray were detected by cross-hybridisation to Italian ryegrass. Transcriptome analyses of the resistant versus the susceptible genotype revealed substantial gene expression differences (>1,200) indicating that great gene expression differences between different Italian ryegrass genotypes exist which potentially contribute to the observed phenotypic divergence in Xtg resistance between the two genotypes. In the resistant genotype, several genes differentially expressed after Xtg inoculation were identified which revealed similarities to transcriptional changes triggered by pathogen-associated molecular patterns in other plant-pathogen interactions. These genes represent candidate genes of particular interest for the development of tools for marker assisted resistance breeding.

Early response mechanisms of perennial ryegrass (Lolium perenne) to phosphorus deficiency

BACKGROUND AND AIMS

Improving phosphorus (P) nutrient efficiency in Lolium perenne (perennial ryegrass) is likely to result in considerable economic and ecological benefits. To date, research into the molecular and biochemical response of perennial ryegrass to P deficiency has been limited, particularly in relation to the early response mechanisms. This study aimed to identify molecular mechanisms activated in response to the initial stages of P deficiency.

METHODS

A barley microarray was successfully used to study gene expression in perennial ryegrass and this was complemented with gas chromatography-mass spectrometry metabolic profiling to obtain an overview of the plant response to early stages of P deficiency.

KEY RESULTS

After 24 h of P deficiency, internal phosphate concentrations were reduced and significant alterations were detected in the metabolome and transcriptome of two perennial ryegrass genotypes. Results indicated a replacement of phospholipids with sulfolipids and the utilization of glycolytic bypasses in response to P deficiency in perennial ryegrass.

CONCLUSIONS

The transcriptome and metabolome of perennial ryegrass undergo changes in response to reductions in P supply after 24 h.

Pathogenesis-associated transcriptional patterns in Triticeae

The Triticeae tribe of the plant Poaceae family contains some of the most important cereal crop plants for nutrition of humans and livestock such as wheat and barley. Despite the agronomical relevance of plant immunity, knowledge on mechanisms of disease or resistance in Triticeae is limited. It is hardly understood what actually stops a microbial invader when restricted by the plant and in how far a susceptible host plant contributes to pathogenesis. Transcriptional reprogramming of the host plant may be involved in both immunity and disease. This paper gives an overview about recent analyses of global pathogenesis-related transcriptional patterns in response of Triticeae to biotrophic or non-biotrophic fungal pathogens and their toxins. It highlights enriched biological functions in association with successful plant defence or disease as well as experiments that successfully translated gene expression data into analysis of gene functions.

Specific patterns of gene space organisation revealed in wheat by using the combination of barley and wheat genomic resources

Background

Because of its size, allohexaploid nature and high repeat content, the wheat genome has always been perceived as too complex for efficient molecular studies. We recently constructed the first physical map of a wheat chromosome (3B). However gene mapping is still laborious in wheat because of high redundancy between the three homoeologous genomes. In contrast, in the closely related diploid species, barley, numerous gene-based markers have been developed. This study aims at combining the unique genomic resources developed in wheat and barley to decipher the organisation of gene space on wheat chromosome 3B.

Results

Three dimensional pools of the minimal tiling path of wheat chromosome 3B physical map were hybridised to a barley Agilent 15K expression microarray. This led to the fine mapping of 738 barley orthologous genes on wheat chromosome 3B. In addition, comparative analyses revealed that 68% of the genes identified were syntenic between the wheat chromosome 3B and barley chromosome 3 H and 59% between wheat chromosome 3B and rice chromosome 1, together with some wheat-specific rearrangements. Finally, it indicated an increasing gradient of gene density from the centromere to the telomeres positively correlated with the number of genes clustered in islands on wheat chromosome 3B.

Conclusion

Our study shows that novel structural genomics resources now available in wheat and barley can be combined efficiently to overcome specific problems of genetic anchoring of physical contigs in wheat and to perform high-resolution comparative analyses with rice for deciphering the organisation of the wheat gene space.

Differential gene expression in nearly isogenic lines with QTL for partial resistance to Puccinia hordei in barley

BACKGROUND

The barley-Puccinia hordei (barley leaf rust) pathosystem is a model for investigating partial disease resistance in crop plants and genetic mapping of phenotypic resistance has identified several quantitative trait loci (QTL) for partial resistance. Reciprocal QTL-specific near-isogenic lines (QTL-NILs) have been developed that combine two QTL, Rphq2 and Rphq3, the largest effects detected in a recombinant-inbred-line (RIL) population derived from a cross between the super-susceptible line L94 and partially-resistant line Vada. The molecular mechanism underpinning partial resistance in these QTL-NILs is unknown.

RESULTS

An Agilent custom microarray consisting of 15,000 probes derived from barley consensus EST sequences was used to investigate genome-wide and QTL-specific differential expression of genes 18 hours post-inoculation (hpi) with Puccinia hordei. A total of 1,410 genes were identified as being significantly differentially expressed across the genome, of which 55 were accounted for by the genetic differences defined by QTL-NILs at Rphq2 and Rphq3. These genes were predominantly located at the QTL regions and are, therefore, positional candidates. One gene, encoding the transcriptional repressor Ethylene-Responsive Element Binding Factor 4 (HvERF4) was located outside the QTL at 71 cM on chromosome 1H, within a previously detected eQTL hotspot for defence response. The results indicate that Rphq2 or Rphq3 contains a trans-eQTL that modulates expression of HvERF4. We speculate that HvERF4 functions as an intermediate that conveys the response signal from a gene(s) contained within Rphq2 or Rphq3 to a host of down-stream defense responsive genes. Our results also reveal that barley lines with extreme or intermediate partial resistance phenotypes exhibit a profound similarity in their spectrum of Ph-responsive genes and that hormone-related signalling pathways are actively involved in response to Puccinia hordei.

CONCLUSIONS

Differential gene expression between QTL-NILs identifies genes predominantly located within the target region(s) providing both transcriptional and positional candidate genes for the QTL. Genetically mapping the differentially expressed genes relative to the QTL has the potential to discover trans-eQTL mediated regulatory relays initiated from genes within the QTL regions.

Data-driven assessment of eQTL mapping methods

Background

The analysis of expression quantitative trait loci (eQTL) is a potentially powerful way to detect transcriptional regulatory relationships at the genomic scale. However, eQTL data sets often go underexploited because legacy QTL methods are used to map the relationship between the expression trait and genotype. Often these methods are inappropriate for complex traits such as gene expression, particularly in the case of epistasis.

Results

Here we compare legacy QTL mapping methods with several modern multi-locus methods and evaluate their ability to produce eQTL that agree with independent external data in a systematic way. We found that the modern multi-locus methods (Random Forests, sparse partial least squares, lasso, and elastic net) clearly outperformed the legacy QTL methods (Haley-Knott regression and composite interval mapping) in terms of biological relevance of the mapped eQTL. In particular, we found that our new approach, based on Random Forests, showed superior performance among the multi-locus methods.

Conclusions

Benchmarks based on the recapitulation of experimental findings provide valuable insight when selecting the appropriate eQTL mapping method. Our battery of tests suggests that Random Forests map eQTL that are more likely to be validated by independent data, when compared to competing multi-locus and legacy eQTL mapping methods.