Sequence divergence between spelt and common wheat.

KEY MESSAGE: Sequence comparison between spelt and common wheat reveals that the former has huge potential in enriching the genetic variation of the latter. Genetic variation is the foundation of crop improvement. By comparing genome sequences of a Triticum spelta accession and one of its derived hexaploid lines with the sequences of the international reference genotype Chinese Spring, we detected variants more than tenfold higher than those present among common wheat (T. aestivum L) genotypes. Furthermore, different from the typical 'V-shaped' pattern of variant distribution often observed along wheat chromosomes, the sequence variation detected in this study was more evenly distributed along the 3B chromosome. This was also the case between T. spelta and the wild emmer genome. Genetic analysis showed that T. spelta and common wheat formed discrete groups. These results showed that, although it is believed that the spelt and common wheat are evolutionarily closely related and belong to the same species, a significant sequence divergence exists between them. Thus, the values of T. spelta in enriching the genetic variation of common wheat can be huge.

Transcriptome analysis of Brachypodium during fungal pathogen infection reveals both shared and distinct defense responses with wheat.

Fusarium crown rot (FCR) of wheat and barley, predominantly caused by the fungal pathogen Fusarium pseudograminearum, is a disease of economic significance. The quantitative nature of FCR resistance within cultivated wheat germplasm has significantly limited breeding efforts to enhanced FCR resistance in wheat. In this study, we characterized the molecular responses of Brachypodium distachyon (Brachypodium hereafter) to F. pseudograminearum infection using RNA-seq to determine whether Brachypodium can be exploited as a model system towards better understanding of F. pseudograminearum-wheat interaction. The transcriptional response to infection in Brachypodium was strikingly similar to that previously reported in wheat, both in shared expression patterns of wheat homologs of Brachypodium genes and functional overlap revealed through comparative gene ontology analysis in both species. Metabolites produced by various biosynthetic pathways induced in both wheat and Brachypodium were quantified, revealing a high degree of overlap between these two species in metabolic response to infection but also showed Brachypodium does not produce certain defence-related metabolites found in wheat. Functional analyses of candidate genes identified in this study will improve our understanding of resistance mechanisms and may lead to the development of new strategies to protect cereal crops from pathogen infection.

MEDIATOR18 and MEDIATOR20 confer susceptibility to Fusarium oxysporum in Arabidopsis thaliana.

The conserved protein complex known as Mediator conveys transcriptional signals by acting as an intermediary between transcription factors and RNA polymerase II. As a result, Mediator subunits play multiple roles in regulating developmental as well as abiotic and biotic stress pathways. In this report we identify the head domain subunits MEDIATOR18 and MEDIATOR20 as important susceptibility factors for Fusarium oxysporum infection in Arabidopsis thaliana. Mutants of MED18 and MED20 display down-regulation of genes associated with jasmonate signaling and biosynthesis while up-regulation of salicylic acid associated pathogenesis related genes and reactive oxygen producing and scavenging genes. We propose that MED18 and MED20 form a sub-domain within Mediator that controls the balance of salicylic acid and jasmonate associated defense pathways.

The Fusarium crown rot pathogen Fusarium pseudograminearum triggers a suite of transcriptional and metabolic changes in bread wheat (Triticum aestivum L.).

Background and Aims: Fusarium crown rot caused by the fungal pathogen Fusarium pseudograminearum is a disease of wheat and barley, bearing significant economic cost. Efforts to develop effective resistance to this disease have been hampered by the quantitative nature of resistance and a lack of understanding of the factors associated with resistance and susceptibility. Here, we aimed to dissect transcriptional responses triggered in wheat by F. pseudograminearum infection. Methods: We used an RNA-seq approach to analyse host responses during a compatible interaction and identified >2700 wheat genes differentially regulated after inoculation with F. pseudograminearum . The production of a few key metabolites and plant hormones in the host during the interaction was also analysed. Key Results: Analysis of gene ontology enrichment showed that a disproportionate number of genes involved in primary and secondary metabolism, signalling and transport were differentially expressed in infected seedlings. A number of genes encoding pathogen-responsive uridine-diphosphate glycosyltransferases (UGTs) potentially involved in detoxification of the Fusarium mycotoxin deoxynivalenol (DON) were differentially expressed. Using a F. pseudograminearum DON-non-producing mutant, DON was shown to play an important role in virulence during Fusarium crown rot. An over-representation of genes involved in the phenylalanine, tryptophan and tyrosine biosynthesis pathways was observed. This was confirmed through metabolite analyses that demonstrated tryptamine and serotonin levels are induced after F. pseudograminearum inoculation. Conclusions: Overall, the observed host response in bread wheat to F. pseudograminearum during early infection exhibited enrichment of processes related to pathogen perception, defence signalling, transport and metabolism and deployment of chemical and enzymatic defences. Additional functional analyses of candidate genes should reveal their roles in disease resistance or susceptibility. Better understanding of host responses contributing to resistance and/or susceptibility will aid the development of future disease improvement strategies against this important plant pathogen.

The defence-associated transcriptome of hexaploid wheat displays homoeolog expression and induction bias.

Bread wheat (Triticum aestivum L.) is an allopolyploid species containing three ancestral genomes. Therefore, three homoeologous copies exist for the majority of genes in the wheat genome. Whether different homoeologs are differentially expressed (homoeolog expression bias) in response to biotic and abiotic stresses is poorly understood. In this study, we applied a RNA-seq approach to analyse homoeolog-specific global gene expression patterns in wheat during infection by the fungal pathogen Fusarium pseudograminearum, which causes crown rot disease in cereals. To ensure specific detection of homoeologs, we first optimized read alignment methods and validated the results experimentally on genes with known patterns of subgenome-specific expression. Our global analysis identified widespread patterns of differential expression among homoeologs, indicating homoeolog expression bias underpins a large proportion of the wheat transcriptome. In particular, genes differentially expressed in response to Fusarium infection were found to be disproportionately contributed from B and D subgenomes. In addition, we found differences in the degree of responsiveness to pathogen infection among homoeologous genes with B and D homoeologs exhibiting stronger responses to pathogen infection than A genome copies. We call this latter phenomenon as 'homoeolog induction bias'. Understanding how homoeolog expression and induction biases operate may assist the improvement of biotic stress tolerance in wheat and other polyploid crop species.

Enzyme-driven metabolomic screening: a proof-of-principle method for discovery of plant defence compounds targeted by pathogens.

Plants produce a variety of secondary metabolites to defend themselves from pathogen attack, while pathogens have evolved to overcome plant defences by producing enzymes that degrade or modify these defence compounds. However, many compounds targeted by pathogen enzymes currently remain enigmatic. Identifying host compounds targeted by pathogen enzymes would enable us to understand the potential importance of such compounds in plant defence and modify them to make them insensitive to pathogen enzymes. Here, a proof of concept metabolomics-based method was developed to discover plant defence compounds modified by pathogens using two pathogen enzymes with known targets in wheat and tomato. Plant extracts treated with purified pathogen enzymes were subjected to LC-MS, and the relative abundance of metabolites before and after treatment were comparatively analysed. Using two enzymes from different pathogens the in planta targets could be found by combining relatively simple enzymology with the power of untargeted metabolomics. Key to the method is dataset simplification based on natural isotope occurrence and statistical filtering, which can be scripted. The method presented here will aid in our understanding of plant-pathogen interactions and may lead to the development of new plant protection strategies.

Characterization of a JAZ7 activation-tagged Arabidopsis mutant with increased susceptibility to the fungal pathogen Fusarium oxysporum.

In Arabidopsis, jasmonate (JA)-signaling plays a key role in mediating Fusarium oxysporum disease outcome. However, the roles of JASMONATE ZIM-domain (JAZ) proteins that repress JA-signaling have not been characterized in host resistance or susceptibility to this pathogen. Here, we found most JAZ genes are induced following F. oxysporum challenge, and screening T-DNA insertion lines in Arabidopsis JAZ family members identified a highly disease-susceptible JAZ7 mutant (jaz7-1D). This mutant exhibited constitutive JAZ7 expression and conferred increased JA-sensitivity, suggesting activation of JA-signaling. Unlike jaz7 loss-of-function alleles, jaz7-1D also had enhanced JA-responsive gene expression, altered development and increased susceptibility to the bacterial pathogen PstDC3000 that also disrupts host JA-responses. We also demonstrate that JAZ7 interacts with transcription factors functioning as activators (MYC3, MYC4) or repressors (JAM1) of JA-signaling and contains a functional EAR repressor motif mediating transcriptional repression via the co-repressor TOPLESS (TPL). We propose through direct TPL recruitment, in wild-type plants JAZ7 functions as a repressor within the JA-response network and that in jaz7-1D plants, misregulated ectopic JAZ7 expression hyper-activates JA-signaling in part by disturbing finely-tuned COI1-JAZ-TPL-TF complexes.

The Fdb3 transcription factor of the Fusarium Detoxification of Benzoxazolinone gene cluster is required for MBOA but not BOA degradation in Fusarium pseudograminearum.

A number of cereals produce the benzoxazolinone class of phytoalexins. Fusarium species pathogenic towards these hosts can typically degrade these compounds via an aminophenol intermediate, and the ability to do so is encoded by a group of genes found in the Fusarium Detoxification of Benzoxazolinone (FDB) cluster. A zinc finger transcription factor encoded by one of the FDB cluster genes (FDB3) has been proposed to regulate the expression of other genes in the cluster and hence is potentially involved in benzoxazolinone degradation. Herein we show that Fdb3 is essential for the ability of Fusarium pseudograminearum to efficiently detoxify the predominant wheat benzoxazolinone, 6-methoxy-benzoxazolin-2-one (MBOA), but not benzoxazoline-2-one (BOA). Furthermore, additional genes thought to be part of the FDB gene cluster, based upon transcriptional response to benzoxazolinones, are regulated by Fdb3. However, deletion mutants for these latter genes remain capable of benzoxazolinone degradation, suggesting that they are not essential for this process.

EffectorP: predicting fungal effector proteins from secretomes using machine learning.

Eukaryotic filamentous plant pathogens secrete effector proteins that modulate the host cell to facilitate infection. Computational effector candidate identification and subsequent functional characterization delivers valuable insights into plant-pathogen interactions. However, effector prediction in fungi has been challenging due to a lack of unifying sequence features such as conserved N-terminal sequence motifs. Fungal effectors are commonly predicted from secretomes based on criteria such as small size and cysteine-rich, which suffers from poor accuracy. We present EffectorP which pioneers the application of machine learning to fungal effector prediction. EffectorP improves fungal effector prediction from secretomes based on a robust signal of sequence-derived properties, achieving sensitivity and specificity of over 80%. Features that discriminate fungal effectors from secreted noneffectors are predominantly sequence length, molecular weight and protein net charge, as well as cysteine, serine and tryptophan content. We demonstrate that EffectorP is powerful when combined with in planta expression data for predicting high-priority effector candidates. EffectorP is the first prediction program for fungal effectors based on machine learning. Our findings will facilitate functional fungal effector studies and improve our understanding of effectors in plant-pathogen interactions. EffectorP is available at http://effectorp.csiro.au.

Fine mapping of a large-effect QTL conferring Fusarium crown rot resistance on the long arm of chromosome 3B in hexaploid wheat.

BACKGROUND: Fusarium crown rot (FCR) is a major cereal disease in semi-arid areas worldwide. Of the various QTL reported, the one on chromosome arm 3BL (Qcrs.cpi-3B) has the largest effect that can be consistently detected in different genetic backgrounds. Nine sets of near isogenic lines (NILs) for this locus were made available in a previous study. To identify markers that could be reliably used in tagging the Qcrs.cpi-3B locus, a NIL-derived population consisting of 774 F10 lines were generated and exploited to assess markers selected from the existing linkage map and generated from sequences of the 3B pseudomolecule. RESULTS: This is the first report on fine mapping a QTL conferring FCR resistance in wheat. By three rounds of linkage mapping using the NILs and the NIL-derived population, the Qcrs.cpi-3B locus was mapped to an interval of 0.7 cM covering a physical distance of about 1.5 Mb. Seven markers co-segregating with the locus were developed. This interval contains a total of 63 gene-coding sequences based on the 3B pseudomolecule, and six of them were known to encode disease resistance proteins. Several of the genes in this interval were among those responsive to FCR infection detected in an earlier study. CONCLUSIONS: The accurate localization of the Qcrs.cpi-3B locus and the development of the markers co-segregating with it should facilitate the incorporation of this large-effect QTL conferring FCR resistance into breeding programs as well as the cloning of the gene(s) underlying the QTL.

Tissue-specific transcriptome analysis within the maturing sugarcane stalk reveals spatial regulation in the expression of cellulose synthase and sucrose transporter gene families.

Sugarcane (Saccharum spp. hybrids) accumulates high concentrations of sucrose in its mature stalk and a considerable portion of carbohydrate metabolism is also devoted to cell wall synthesis and fibre production. We examined tissue-specific expression patterns to explore the spatial deployment of pathways responsible for sucrose accumulation and fibre synthesis within the stalk. We performed expression profiling of storage parenchyma, vascular bundles and rind dissected from a maturing stalk internode of sugarcane, identifying ten cellulose synthase subunit genes and examining significant differences in the expression of their corresponding transcripts and those of several sugar transporters. These were correlated with differential expression patterns for transcripts of genes encoding COBRA-like proteins and other cell wall metabolism-related proteins. The sugar transporters genes ShPST2a, ShPST2b and ShSUT4 were significantly up-regulated in storage parenchyma while ShSUT1 was up-regulated in vascular bundles. Two co-ordinately expressed groups of cell wall related transcripts were also identified. One group, associated with primary cell wall synthesis (ShCesA1, ShCesA7, ShCesA9 and Shbk2l3), was up-regulated in parenchyma. The other group, associated with secondary cell wall synthesis (ShCesA10, ShCesA11, ShCesA12 and Shbk-2), was up-regulated in rind. In transformed sugarcane plants, the ShCesA7 promoter conferred stable expression of green fluorescent protein preferentially in the storage parenchyma of the maturing stalk internode. Our results indicate that there is spatial separation for elevated expression of these important targets in both sucrose accumulation and cell wall synthesis, allowing for increased clarity in our understanding of sucrose transport and fibre synthesis in sugarcane.

A γ-lactamase from cereal infecting Fusarium spp. catalyses the first step in the degradation of the benzoxazolinone class of phytoalexins.

The benzoxazolinone class of phytoalexins are released by wheat, maize, rye and other agriculturally important species in the Poaceae family upon pathogen attack. Benzoxazolinones show antimicrobial effects on plant pathogens, but certain fungi have evolved mechanisms to actively detoxify these compounds which may contribute to the virulence of the pathogens. In many Fusarium spp. a cluster of genes is thought to be involved in the detoxification of benzoxazolinones. However, only one enzyme encoded in the cluster has been unequivocally assigned a role in this process. The first step in the detoxification of benzoxazolinones in Fusarium spp. involves the hydrolysis of a cyclic ester bond. This reaction is encoded by the FDB1 locus in F. verticillioides but the underlying gene is yet to be cloned. We previously proposed that FDB1 encodes a γ-lactamase, and here direct evidence for this is presented. Expression analyses in the important wheat pathogen F. pseudograminearum demonstrated that amongst the three predicted γ-lactamase genes only the one designated as FDB1, part of the proposed benzoxazolinone cluster in F. pseudograminearum, was strongly responsive to exogenous benzoxazolinone application. Analysis of independent F. pseudograminearum and F. graminearum FDB1 gene deletion mutants, as well as biochemical assays, demonstrated that the γ-lactamase enzyme, encoded by FDB1, catalyses the first step in detoxification of benzoxazolinones. Overall, our results support the notion that Fusarium pathogens that cause crown rot and head blight on wheat have adopted strategies to overcome host-derived chemical defences.

Investigating the Association between Flowering Time and Defense in the Arabidopsis thaliana-Fusarium oxysporum Interaction.

Plants respond to pathogens either by investing more resources into immunity which is costly to development, or by accelerating reproductive processes such as flowering time to ensure reproduction occurs before the plant succumbs to disease. In this study we explored the link between flowering time and pathogen defense using the interaction between Arabidopsis thaliana and the root infecting fungal pathogen Fusarium oxysporum. We report that F. oxysporum infection accelerates flowering time and regulates transcription of a number of floral integrator genes, including FLOWERING LOCUS C (FLC), FLOWERING LOCUS T (FT) and GIGANTEA (GI). Furthermore, we observed a positive correlation between late flowering and resistance to F. oxysporum in A. thaliana natural ecotypes. Late-flowering gi and autonomous pathway mutants also exhibited enhanced resistance to F. oxysporum, supporting the association between flowering time and defense. However, epistasis analysis showed that accelerating flowering time by deletion of FLC in fve-3 or fpa-7 mutants did not alter disease resistance, suggesting that the effect of autonomous pathway on disease resistance occurs independently from flowering time. Indeed, RNA-seq analyses suggest that fve-3 mediated resistance to F. oxysporum is most likely a result of altered defense-associated gene transcription. Together, our results indicate that the association between flowering time and pathogen defense is complex and can involve both pleiotropic and direct effects.

Genome-Wide Analysis in Three Fusarium Pathogens Identifies Rapidly Evolving Chromosomes and Genes Associated with Pathogenicity.

Pathogens and hosts are in an ongoing arms race and genes involved in host-pathogen interactions are likely to undergo diversifying selection. Fusarium plant pathogens have evolved diverse infection strategies, but how they interact with their hosts in the biotrophic infection stage remains puzzling. To address this, we analyzed the genomes of three Fusarium plant pathogens for genes that are under diversifying selection. We found a two-speed genome structure both on the chromosome and gene group level. Diversifying selection acts strongly on the dispensable chromosomes in Fusarium oxysporum f. sp. lycopersici and on distinct core chromosome regions in Fusarium graminearum, all of which have associations with virulence. Members of two gene groups evolve rapidly, namely those that encode proteins with an N-terminal [SG]-P-C-[KR]-P sequence motif and proteins that are conserved predominantly in pathogens. Specifically, 29 F. graminearum genes are rapidly evolving, in planta induced and encode secreted proteins, strongly pointing toward effector function. In summary, diversifying selection in Fusarium is strongly reflected as genomic footprints and can be used to predict a small gene set likely to be involved in host-pathogen interactions for experimental verification.

Fusarium oxysporum triggers tissue-specific transcriptional reprogramming in Arabidopsis thaliana.

Some of the most devastating agricultural diseases are caused by root-infecting pathogens, yet the majority of studies on these interactions to date have focused on the host responses of aerial tissues rather than those belowground. Fusarium oxysporum is a root-infecting pathogen that causes wilt disease on several plant species including Arabidopsis thaliana. To investigate and compare transcriptional changes triggered by F. oxysporum in different Arabidopsis tissues, we infected soil-grown plants with F. oxysporum and subjected root and leaf tissue harvested at early and late timepoints to RNA-seq analyses. At least half of the genes induced or repressed by F. oxysporum showed tissue-specific regulation. Regulators of auxin and ABA signalling, mannose binding lectins and peroxidases showed strong differential expression in root tissue. We demonstrate that ARF2 and PRX33, two genes regulated in the roots, promote susceptibility to F. oxysporum. In the leaves, defensins and genes associated with the response to auxin, cold and senescence were strongly regulated while jasmonate biosynthesis and signalling genes were induced throughout the plant.

Brachypodium as an emerging model for cereal-pathogen interactions.

BACKGROUND: Cereal diseases cause tens of billions of dollars of losses annually and have devastating humanitarian consequences in the developing world. Increased understanding of the molecular basis of cereal host-pathogen interactions should facilitate development of novel resistance strategies. However, achieving this in most cereals can be challenging due to large and complex genomes, long generation times and large plant size, as well as quarantine and intellectual property issues that may constrain the development and use of community resources. Brachypodium distachyon (brachypodium) with its small, diploid and sequenced genome, short generation time, high transformability and rapidly expanding community resources is emerging as a tractable cereal model. SCOPE: Recent research reviewed here has demonstrated that brachypodium is either susceptible or partially susceptible to many of the major cereal pathogens. Thus, the study of brachypodium-pathogen interactions appears to hold great potential to improve understanding of cereal disease resistance, and to guide approaches to enhance this resistance. This paper reviews brachypodium experimental pathosystems for the study of fungal, bacterial and viral cereal pathogens; the current status of the use of brachypodium for functional analysis of cereal disease resistance; and comparative genomic approaches undertaken using brachypodium to assist characterization of cereal resistance genes. Additionally, it explores future prospects for brachypodium as a model to study cereal-pathogen interactions. CONCLUSIONS: The study of brachypodium-pathogen interactions appears to be a productive strategy for understanding mechanisms of disease resistance in cereal species. Knowledge obtained from this model interaction has strong potential to be exploited for crop improvement.

Degradation of the benzoxazolinone class of phytoalexins is important for virulence of Fusarium pseudograminearum towards wheat.

Wheat, maize, rye and certain other agriculturally important species in the Poaceae family produce the benzoxazolinone class of phytoalexins on pest and pathogen attack. Benzoxazolinones can inhibit the growth of pathogens. However, certain fungi can actively detoxify these compounds. Despite this, a clear link between the ability to detoxify benzoxazolinones and pathogen virulence has not been shown. Here, through comparative genome analysis of several Fusarium species, we have identified a conserved genomic region around the FDB2 gene encoding an N-malonyltransferase enzyme known to be involved in benzoxazolinone degradation in the maize pathogen Fusarium verticillioides. Expression analyses demonstrated that a cluster of nine genes was responsive to exogenous benzoxazolinone in the important wheat pathogen Fusarium pseudograminearum. The analysis of independent F. pseudograminearum FDB2 knockouts and complementation of the knockout with FDB2 homologues from F. graminearum and F. verticillioides confirmed that the N-malonyltransferase enzyme encoded by this gene is central to the detoxification of benzoxazolinones, and that Fdb2 contributes quantitatively to virulence towards wheat in head blight inoculation assays. This contrasts with previous observations in F. verticillioides, where no effect of FDB2 mutations on pathogen virulence towards maize was observed. Overall, our results demonstrate that the detoxification of benzoxazolinones is a strategy adopted by wheat-infecting F. pseudograminearum to overcome host-derived chemical defences.

An assessment of heavy ion irradiation mutagenesis for reverse genetics in wheat (Triticum aestivum L.).

Reverse genetic techniques harnessing mutational approaches are powerful tools that can provide substantial insight into gene function in plants. However, as compared to diploid species, reverse genetic analyses in polyploid plants such as bread wheat can present substantial challenges associated with high levels of sequence and functional similarity amongst homoeologous loci. We previously developed a high-throughput method to identify deletions of genes within a physically mutagenized wheat population. Here we describe our efforts to combine multiple homoeologous deletions of three candidate disease susceptibility genes (TaWRKY11, TaPFT1 and TaPLDß1). We were able to produce lines featuring homozygous deletions at two of the three homoeoloci for all genes, but this was dependent on the individual mutants used in crossing. Intriguingly, despite extensive efforts, viable lines possessing homozygous deletions at all three homoeoloci could not be produced for any of the candidate genes. To investigate deletion size as a possible reason for this phenomenon, we developed an amplicon sequencing approach based on synteny to Brachypodium distachyon to assess the size of the deletions removing one candidate gene (TaPFT1) in our mutants. These analyses revealed that genomic deletions removing the locus are relatively large, resulting in the loss of multiple additional genes. The implications of this work for the use of heavy ion mutagenesis for reverse genetic analyses in wheat are discussed.

Transcriptome and allele specificity associated with a 3BL locus for Fusarium crown rot resistance in bread wheat.

Fusarium pathogens cause two major diseases in cereals, Fusarium crown rot (FCR) and head blight (FHB). A large-effect locus conferring resistance to FCR disease was previously located to chromosome arm 3BL (designated as Qcrs-3B) and several independent sets of near isogenic lines (NILs) have been developed for this locus. In this study, five sets of the NILs were used to examine transcriptional changes associated with the Qcrs-3B locus and to identify genes linked to the resistance locus as a step towards the isolation of the causative gene(s). Of the differentially expressed genes (DEGs) detected between the NILs, 12.7% was located on the single chromosome 3B. Of the expressed genes containing SNP (SNP-EGs) detected, 23.5% was mapped to this chromosome. Several of the DEGs and SNP-EGs are known to be involved in host-pathogen interactions, and a large number of the DEGs were among those detected for FHB in previous studies. Of the DEGs detected, 22 were mapped in the Qcrs-3B interval and they included eight which were detected in the resistant isolines only. The enrichment of DEG, and not necessarily those containing SNPs between the resistant and susceptible isolines, around the Qcrs-3B locus is suggestive of local regulation of this region by the resistance allele. Functions for 13 of these DEGs are known. Of the SNP-EGs, 28 were mapped in the Qcrs-3B interval and biological functions for 16 of them are known. These results provide insights into responses regulated by the 3BL locus and identify a tractable number of target genes for fine mapping and functional testing to identify the causative gene(s) at this QTL.