Multiple toxins and a protease contribute to the aphid-killing ability of Pseudomonas fluorescens PpR24.

Aphids are globally important pests causing damage to a broad range of crops. Due to insecticide resistance, there is an urgent need to develop alternative control strategies. In our previous work, we found Pseudomonas fluorescens PpR24 can orally infect and kill the insecticide-resistant green-peach aphid (Myzus persicae). However, the genetic basis of the insecticidal capability of PpR24 remains unclear. Genome sequencing of PpR24 confirmed the presence of various insecticidal toxins such as Tc (toxin complexes), Rhs (rearrangement hotspot) elements, and other insect-killing proteases. Upon aphids infection with PpR24, RNA-Seq analysis revealed 193 aphid genes were differentially expressed with down-regulation of 16 detoxification genes. In addition, 1325 PpR24 genes (542 were upregulated and 783 downregulated) were subject to differential expression, including genes responsible for secondary metabolite biosynthesis, the iron-restriction response, oxidative stress resistance, and virulence factors. Single and double deletion of candidate virulence genes encoding a secreted protease (AprX) and four toxin components (two TcA-like; one TcB-like; one TcC-like insecticidal toxins) showed that all five genes contribute significantly to aphid killing, particularly AprX. This comprehensive host-pathogen transcriptomic analysis provides novel insight into the molecular basis of bacteria-mediated aphid mortality and the potential of PpR24 as an effective biocontrol agent.

The complete genome assemblies of 19 insect pests of worldwide importance to agriculture.

There are many insect pests worldwide that damage agricultural crop and reduce yield either by direct feeding or by the transmission of plant diseases. To date, control of pest insects has been achieved largely by applying synthetic insecticides. However, insecticide use can be seriously impacted by legislation that limits their use or by the evolution of resistance in the target pest. Thus, there is a move towards less use of insecticides and increased adoption of integrated pest management strategies using a wide range of non-chemical and chemical control methods. For good pest control there is a need to understand the mode of action and selectivity of insecticides, the life cycles of the pests and their biology and behaviours, all of which can benefit from good quality genome data. Here we present the complete assembled (chromosome level) genomes (incl. mtDNA) of 19 insect pests, Agriotes lineatus (click beetle/wireworm), Aphis gossypii (melon/cotton aphid), Bemisia tabaci (cotton whitefly), Brassicogethes aeneus (pollen beetle), Ceutorhynchus obstrictus (seedpod weevil), Chilo suppressalis (striped rice stem borer), Chrysodeixis includens (soybean looper), Diabrotica balteata (cucumber beetle), Diatraea saccharalis (sugar cane borer), Nezara viridula (green stink bug), Nilaparvata lugens (brown plant hopper), Phaedon cochleariae (mustard beetle), Phyllotreta striolata (striped flea beetle), Psylliodes chrysocephala (cabbage stem flea beetle), Spodoptera exigua (beet army worm), Spodoptera littoralis (cotton leaf worm), Diabrotica virgifera (western corn root worm), Euschistus heros (brown stink bug) and Phyllotreta cruciferae (crucifer flea beetle). For the first 15 of these we also present the annotation of genes encoding potential xenobiotic detoxification enzymes. This public resource will aid in the elucidation and monitoring of resistance mechanisms, the development of highly selective chemistry and potential techniques to disrupt behaviour in a way that limits the effect of the pests.

Basic LEUCINE ZIPPER TRANSCRIPTION FACTOR67 Transactivates DELAY OF GERMINATION1 to Establish Primary Seed Dormancy in Arabidopsis.

Seed dormancy governs the timing of germination, one of the most important developmental transitions in a plant's life cycle. The DELAY OF GERMINATION1 (DOG1) gene is a key regulator of seed dormancy and a major quantitative trait locus in Arabidopsis (Arabidopsis thaliana). DOG1 expression is under tight developmental and environmental regulation, but the transcription factors involved are not known. Here we show that basic LEUCINE ZIPPER TRANSCRIPTION FACTOR67 (bZIP67) acts downstream of the central regulator of seed development, LEAFY COTYLEDON1, to transactivate DOG1 during maturation and help to establish primary dormancy. We show that bZIP67 overexpression enhances dormancy and that bZIP67 protein (but not transcript) abundance is increased in seeds matured in cool conditions, providing a mechanism to explain how temperature regulates DOG1 expression. We also show that natural allelic variation in the DOG1 promoter affects bZIP67-dependent transactivation, providing a mechanism to explain ecotypic differences in seed dormancy that are controlled by the DOG1 locus.

PHI-base: the pathogen-host interactions database.

The pathogen-host interactions database (PHI-base) is available at www.phi-base.org. PHI-base contains expertly curated molecular and biological information on genes proven to affect the outcome of pathogen-host interactions reported in peer reviewed research articles. PHI-base also curates literature describing specific gene alterations that did not affect the disease interaction phenotype, in order to provide complete datasets for comparative purposes. Viruses are not included, due to their extensive coverage in other databases. In this article, we describe the increased data content of PHI-base, plus new database features and further integration with complementary databases. The release of PHI-base version 4.8 (September 2019) contains 3454 manually curated references, and provides information on 6780 genes from 268 pathogens, tested on 210 hosts in 13,801 interactions. Prokaryotic and eukaryotic pathogens are represented in almost equal numbers. Host species consist of approximately 60% plants (split 50:50 between cereal and non-cereal plants), and 40% other species of medical and/or environmental importance. The information available on pathogen effectors has risen by more than a third, and the entries for pathogens that infect crop species of global importance has dramatically increased in this release. We also briefly describe the future direction of the PHI-base project, and some existing problems with the PHI-base curation process.

Epigenetic patterns within the haplotype phased fig (Ficus carica L.) genome.

Due to DNA heterozygosity and repeat content, assembly of non-model plant genomes is challenging. Herein, we report a high-quality genome reference of one of the oldest known domesticated species, fig (Ficus carica L.), using Pacific Biosciences single-molecule, real-time sequencing. The fig genome is ~333 Mbp in size, of which 80% has been anchored to 13 chromosomes. Genome-wide analysis of N6 -methyladenine and N4 -methylcytosine revealed high methylation levels in both genes and transposable elements, and a prevalence of methylated over non-methylated genes. Furthermore, the characterization of N6 -methyladenine sites led to the identification of ANHGA, a species-specific motif, which is prevalent for both genes and transposable elements. Finally, exploiting the contiguity of the 13 pseudomolecules, we identified 13 putative centromeric regions. The high-quality reference genome and the characterization of methylation profiles, provides an important resource for both fig breeding and for fundamental research into the relationship between epigenetic changes and phenotype, using fig as a model species.

KnetMiner: a comprehensive approach for supporting evidence-based gene discovery and complex trait analysis across species.

The generation of new ideas and scientific hypotheses is often the result of extensive literature and database searches, but, with the growing wealth of public and private knowledge, the process of searching diverse and interconnected data to generate new insights into genes, gene networks, traits and diseases is becoming both more complex and more time-consuming. To guide this technically challenging data integration task and to make gene discovery and hypotheses generation easier for researchers, we have developed a comprehensive software package called KnetMiner which is open-source and containerized for easy use. KnetMiner is an integrated, intelligent, interactive gene and gene network discovery platform that supports scientists explore and understand the biological stories of complex traits and diseases across species. It features fast algorithms for generating rich interactive gene networks and prioritizing candidate genes based on knowledge mining approaches. KnetMiner is used in many plant science institutions and has been adopted by several plant breeding organizations to accelerate gene discovery. The software is generic and customizable and can therefore be readily applied to new species and data types; for example, it has been applied to pest insects and fungal pathogens; and most recently repurposed to support COVID-19 research. Here, we give an overview of the main approaches behind KnetMiner and we report plant-centric case studies for identifying genes, gene networks and trait relationships in Triticum aestivum (bread wheat), as well as, an evidence-based approach to rank candidate genes under a large Arabidopsis thaliana QTL. KnetMiner is available at: https://knetminer.org.

Genome Wide Analysis of Fatty Acid Desaturation and Its Response to Temperature.

Plants modify the polyunsaturated fatty acid content of their membrane and storage lipids in order to adapt to changes in temperature. In developing seeds, this response is largely controlled by the activities of the microsomal ω-6 and ω-3 fatty acid desaturases, FAD2 and FAD3. Although temperature regulation of desaturation has been studied at the molecular and biochemical levels, the genetic control of this trait is poorly understood. Here, we have characterized the response of Arabidopsis (Arabidopsis thaliana) seed lipids to variation in ambient temperature and found that heat inhibits both ω-6 and ω-3 desaturation in phosphatidylcholine, leading to a proportional change in triacylglycerol composition. Analysis of the 19 parental accessions of the multiparent advanced generation intercross (MAGIC) population showed that significant natural variation exists in the temperature responsiveness of ω-6 desaturation. A combination of quantitative trait locus (QTL) analysis and genome-wide association studies (GWAS) using the MAGIC population suggests that ω-6 desaturation is largely controlled by cis-acting sequence variants in the FAD2 5' untranslated region intron that determine the expression level of the gene. However, the temperature responsiveness of ω-6 desaturation is controlled by a separate QTL on chromosome 2. The identity of this locus is unknown, but genome-wide association studies identified potentially causal sequence variants within ∼40 genes in an ∼450-kb region of the QTL.

Insecticide resistance mediated by an exon skipping event.

Many genes increase coding capacity by alternate exon usage. The gene encoding the insect nicotinic acetylcholine receptor (nAChR) α6 subunit, target of the bio-insecticide spinosad, is one example of this and expands protein diversity via alternative splicing of mutually exclusive exons. Here, we show that spinosad resistance in the tomato leaf miner, Tuta absoluta is associated with aberrant regulation of splicing of Taα6 resulting in a novel form of insecticide resistance mediated by exon skipping. Sequencing of the α6 subunit cDNA from spinosad selected and unselected strains of T. absoluta revealed all Taα6 transcripts of the selected strain were devoid of exon 3, with comparison of genomic DNA and mRNA revealing this is a result of exon skipping. Exon skipping cosegregated with spinosad resistance in survival bioassays, and functional characterization of this alteration using modified human nAChR α7, a model of insect α6, demonstrated that exon 3 is essential for receptor function and hence spinosad sensitivity. DNA and RNA sequencing analyses suggested that exon skipping did not result from genetic alterations in intronic or exonic cis-regulatory elements, but rather was associated with a single epigenetic modification downstream of exon 3a, and quantitative changes in the expression of trans-acting proteins that have known roles in the regulation of alternative splicing. Our results demonstrate that the intrinsic capacity of the α6 gene to generate transcript diversity via alternative splicing can be readily exploited during the evolution of resistance and identifies exon skipping as a molecular alteration conferring insecticide resistance.

Transcriptome and metabolite profiling of the infection cycle of Zymoseptoria tritici on wheat reveals a biphasic interaction with plant immunity involving differential pathogen chromosomal contributions and a variation on the hemibiotrophic lifestyle definition.

The hemibiotrophic fungus Zymoseptoria tritici causes Septoria tritici blotch disease of wheat (Triticum aestivum). Pathogen reproduction on wheat occurs without cell penetration, suggesting that dynamic and intimate intercellular communication occurs between fungus and plant throughout the disease cycle. We used deep RNA sequencing and metabolomics to investigate the physiology of plant and pathogen throughout an asexual reproductive cycle of Z. tritici on wheat leaves. Over 3,000 pathogen genes, more than 7,000 wheat genes, and more than 300 metabolites were differentially regulated. Intriguingly, individual fungal chromosomes contributed unequally to the overall gene expression changes. Early transcriptional down-regulation of putative host defense genes was detected in inoculated leaves. There was little evidence for fungal nutrient acquisition from the plant throughout symptomless colonization by Z. tritici, which may instead be utilizing lipid and fatty acid stores for growth. However, the fungus then subsequently manipulated specific plant carbohydrates, including fructan metabolites, during the switch to necrotrophic growth and reproduction. This switch coincided with increased expression of jasmonic acid biosynthesis genes and large-scale activation of other plant defense responses. Fungal genes encoding putative secondary metabolite clusters and secreted effector proteins were identified with distinct infection phase-specific expression patterns, although functional analysis suggested that many have overlapping/redundant functions in virulence. The pathogenic lifestyle of Z. tritici on wheat revealed through this study, involving initial defense suppression by a slow-growing extracellular and nutritionally limited pathogen followed by defense (hyper) activation during reproduction, reveals a subtle modification of the conceptual definition of hemibiotrophic plant infection.

Whole-genome analysis of Fusarium graminearum insertional mutants identifies virulence associated genes and unmasks untagged chromosomal deletions.

BACKGROUND: Identifying pathogen virulence genes required to cause disease is crucial to understand the mechanisms underlying the pathogenic process. Plasmid insertion mutagenesis of fungal protoplasts is frequently used for this purpose in filamentous ascomycetes. Post transformation, the mutant population is screened for loss of virulence to a specific plant or animal host. Identifying the insertion event has previously met with varying degrees of success, from a cleanly disrupted gene with minimal deletion of nucleotides at the insertion point to multiple-copy insertion events and large deletions of chromosomal regions. Currently, extensive mutant collections exist in laboratories globally where it was hitherto impossible to identify all the affected genes. RESULTS: We used a whole-genome sequencing (WGS) approach using Illumina HiSeq 2000 technology to investigate DNA tag insertion points and chromosomal deletion events in mutagenised, reduced virulence F. graminearum isolates identified in disease tests on wheat (Triticum aestivum). We developed the FindInsertSeq workflow to localise the DNA tag insertions to the nucleotide level. The workflow was tested using four mutants showing evidence of single and multi-copy insertions in DNA blot analysis. FindInsertSeq was able to identify both single and multi-copy concatenation insertion sites. By comparing sequencing coverage, unexpected molecular recombination events such as large tagged and untagged chromosomal deletions, and DNA amplification were observed in three of the analysed mutants. A random data sampling approach revealed the minimum genome coverage required to survey the F. graminearum genome for alterations. CONCLUSIONS: This study demonstrates that whole-genome re-sequencing to 22x fold genome coverage is an efficient tool to characterise single and multi-copy insertion mutants in the filamentous ascomycete Fusarium graminearum. In some cases insertion events are accompanied with large untagged chromosomal deletions while in other cases a straight-forward insertion event could be confirmed. The FindInsertSeq analysis workflow presented in this study enables researchers to efficiently characterise insertion and deletion mutants.

The completed genome sequence of the pathogenic ascomycete fungus Fusarium graminearum.

BACKGROUND: Accurate genome assembly and gene model annotation are critical for comparative species and gene functional analyses. Here we present the completed genome sequence and annotation of the reference strain PH-1 of Fusarium graminearum, the causal agent of head scab disease of small grain cereals which threatens global food security. Completion was achieved by combining (a) the BROAD Sanger sequenced draft, with (b) the gene predictions from Munich Information Services for Protein Sequences (MIPS) v3.2, with (c) de novo whole-genome shotgun re-sequencing, (d) re-annotation of the gene models using RNA-seq evidence and Fgenesh, Snap, GeneMark and Augustus prediction algorithms, followed by (e) manual curation. RESULTS: We have comprehensively completed the genomic 36,563,796 bp sequence by replacing unknown bases, placing supercontigs within their correct loci, correcting assembly errors, and inserting new sequences which include for the first time complete AT rich sequences such as centromere sequences, subtelomeric regions and the telomeres. Each of the four F. graminearium chromosomes was found to be submetacentric with respect to centromere positioning. The position of a potential neocentromere was also defined. A preferentially higher frequency of genetic recombination was observed at the end of the longer arm of each chromosome. Within the genome 1529 gene models have been modified and 412 new gene models predicted, with a total gene call of 14,164. The re-annotation impacts upon 69 entries held within the Pathogen-Host Interactions database (PHI-base) which stores information on genes for which mutant phenotypes in pathogen-host interactions have been experimentally tested, of which 59 are putative transcription factors, 8 kinases, 1 ATP citrate lyase (ACL1), and 1 syntaxin-like SNARE gene (GzSYN1). Although the completed F. graminearum contains very few transposon sequences, a previously unrecognised and potentially active gypsy-type long-terminal-repeat (LTR) retrotransposon was identified. In addition, each of the sub-telomeres and centromeres contained either a LTR or MarCry-1_FO element. The full content of the proposed ancient chromosome fusion sites has also been revealed and investigated. Regions with high recombination previously noted to be rich in secretome encoding genes were also found to be rich in tRNA sequences. This study has identified 741 F. graminearum species specific genes and provides the first complete genome assembly for a Sordariomycetes species. CONCLUSIONS: This fully completed F. graminearum PH-1 genome and manually curated annotation, available at Ensembl Fungi, provides the optimum resource to perform interspecies comparative analyses and gene function studies.

Ondex Web: web-based visualization and exploration of heterogeneous biological networks.

SUMMARY: Ondex Web is a new web-based implementation of the network visualization and exploration tools from the Ondex data integration platform. New features such as context-sensitive menus and annotation tools provide users with intuitive ways to explore and manipulate the appearance of heterogeneous biological networks. Ondex Web is open source, written in Java and can be easily embedded into Web sites as an applet. Ondex Web supports loading data from a variety of network formats, such as XGMML, NWB, Pajek and OXL. AVAILABILITY AND IMPLEMENTATION: http://ondex.rothamsted.ac.uk/OndexWeb.

Secondary cell wall composition and candidate gene expression in developing willow (Salix purpurea) stems.

The properties of the secondary cell wall (SCW) in willow largely determine the suitability of willow biomass feedstock for potential bioenergy and biofuel applications. SCW development has been little studied in willow and it is not known how willow compares with model species, particularly the closely related genus Populus. To address this and relate SCW synthesis to candidate genes in willow, a tractable bud culture-derived system was developed in Salix purpurea, and cell wall composition and RNA-Seq transcriptome were followed in stems during early development. A large increase in SCW deposition in the period 0-2 weeks after transfer to soil was characterised by a big increase in xylan content, but no change in the frequency of substitution of xylan with glucuronic acid, and increased abundance of putative transcripts for synthesis of SCW cellulose, xylan and lignin. Histochemical staining and immunolabeling revealed that increased deposition of lignin and xylan was associated with xylem, xylem fibre cells and phloem fibre cells. Transcripts orthologous to those encoding xylan synthase components IRX9 and IRX10 and xylan glucuronyl transferase GUX1 in Arabidopsis were co-expressed, and showed the same spatial pattern of expression revealed by in situ hybridisation at four developmental stages, with abundant expression in proto-xylem, xylem fibre and ray parenchyma cells and some expression in phloem fibre cells. The results show a close similarity with SCW development in Populus species, but also give novel information on the relationship between spatial and temporal variation in xylan-related transcripts and xylan composition.

Proanthocyanidin biosynthesis in the developing wheat seed coat investigated by chemical and RNA-Seq analysis.

The composition of proanthocyanidins in the testa (seed coat) of bread wheat was analyzed by thiolysis of PA oligomers from developing grain and found to consist of (+)-catechin monomers, with a small amount of (+)-gallocatechin. The average chain length of soluble PA stayed relatively constant between 10 and 20 days post-anthesis, whereas that of unextractable PA increased over the same period, suggesting that increases in chain length might account for the insolubility of PAs from mature wheat grain. We carried out RNA-Seq followed by differential expression analysis from dissected tissues of developing grain from red- and white-grained near-isogenic lines differing in the presence of an active R gene that encodes a MYB transcription factor involved in control of PA biosynthesis. In addition to genes already identified encoding chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, and dihydroxyflavonoid 4-reductase, we showed that wheat genes encoding phenylalanine ammonia lyase, flavonoid 3',5'-hydroxylase, leucoanthocyanidin reductase, and a glutathione S-transferase (the orthologue of maize Bronze-2) were more highly expressed in the red NIL. We also identified candidate orthologues of other catalytic and regulatory components of flavonoid biosynthesis in wheat.

Novel molecules and target genes for vegetative heat tolerance in wheat.

To prevent yield losses caused by climate change, it is important to identify naturally tolerant genotypes with traits and related pathways that can be targeted for crop improvement. Here we report on the characterization of contrasting vegetative heat tolerance in two UK bread wheat varieties. Under chronic heat stress, the heat-tolerant cultivar Cadenza produced an excessive number of tillers which translated into more spikes and higher grain yield compared to heat-sensitive Paragon. RNAseq and metabolomics analyses revealed that over 5000 genotype-specific genes were differentially expressed, including photosynthesis-related genes, which might explain the observed ability of Cadenza to maintain photosynthetic rate under heat stress. Around 400 genes showed a similar heat-response in both genotypes. Only 71 genes showed a genotype × temperature interaction. As well as known heat-responsive genes such as heat shock proteins (HSPs), several genes that have not been previously linked to the heat response, particularly in wheat, have been identified, including dehydrins, ankyrin-repeat protein-encoding genes, and lipases. Contrary to primary metabolites, secondary metabolites showed a highly differentiated heat response and genotypic differences. These included benzoxazinoid (DIBOA, DIMBOA), and phenylpropanoids and flavonoids with known radical scavenging capacity, which was assessed via the DPPH assay. The most highly heat-induced metabolite was (glycosylated) propanediol, which is widely used in industry as an anti-freeze. To our knowledge, this is the first report on its response to stress in plants. The identified metabolites and candidate genes provide novel targets for the development of heat-tolerant wheat.

Assessing the functional coherence of modules found in multiple-evidence networks from Arabidopsis.

BACKGROUND: Combining multiple evidence-types from different information sources has the potential to reveal new relationships in biological systems. The integrated information can be represented as a relationship network, and clustering the network can suggest possible functional modules. The value of such modules for gaining insight into the underlying biological processes depends on their functional coherence. The challenges that we wish to address are to define and quantify the functional coherence of modules in relationship networks, so that they can be used to infer function of as yet unannotated proteins, to discover previously unknown roles of proteins in diseases as well as for better understanding of the regulation and interrelationship between different elements of complex biological systems. RESULTS: We have defined the functional coherence of modules with respect to the Gene Ontology (GO) by considering two complementary aspects: (i) the fragmentation of the GO functional categories into the different modules and (ii) the most representative functions of the modules. We have proposed a set of metrics to evaluate these two aspects and demonstrated their utility in Arabidopsis thaliana. We selected 2355 proteins for which experimentally established protein-protein interaction (PPI) data were available. From these we have constructed five relationship networks, four based on single types of data: PPI, co-expression, co-occurrence of protein names in scientific literature abstracts and sequence similarity and a fifth one combining these four evidence types. The ability of these networks to suggest biologically meaningful grouping of proteins was explored by applying Markov clustering and then by measuring the functional coherence of the clusters. CONCLUSIONS: Relationship networks integrating multiple evidence-types are biologically informative and allow more proteins to be assigned to a putative functional module. Using additional evidence types concentrates the functional annotations in a smaller number of modules without unduly compromising their consistency. These results indicate that integration of more data sources improves the ability to uncover functional association between proteins, both by allowing more proteins to be linked and producing a network where modular structure more closely reflects the hierarchy in the gene ontology.

Gene-based mapping of trehalose biosynthetic pathway genes reveals association with source- and sink-related yield traits in a spring wheat panel.

Trehalose 6-phosphate (T6P) signalling regulates carbon use and allocation and is a target to improve crop yields. However, the specific contributions of trehalose phosphate synthase (TPS) and trehalose phosphate phosphatase (TPP) genes to source- and sink-related traits remain largely unknown. We used enrichment capture sequencing on TPS and TPP genes to estimate and partition the genetic variation of yield-related traits in a spring wheat (Triticum aestivum) breeding panel specifically built to capture the diversity across the 75,000 CIMMYT wheat cultivar collection. Twelve phenotypes were correlated to variation in TPS and TPP genes including plant height and biomass (source), spikelets per spike, spike growth and grain filling traits (sink) which showed indications of both positive and negative gene selection. Individual genes explained proportions of heritability for biomass and grain-related traits. Three TPS1 homologues were particularly significant for trait variation. Epistatic interactions were found within and between the TPS and TPP gene families for both plant height and grain-related traits. Gene-based prediction improved predictive ability for grain weight when gene effects were combined with the whole-genome markers. Our study has generated a wealth of information on natural variation of TPS and TPP genes related to yield potential which confirms the role for T6P in resource allocation and in affecting traits such as grain number and size confirming other studies which now opens up the possibility of harnessing natural genetic variation more widely to better understand the contribution of native genes to yield traits for incorporation into breeding programmes.

The Wheat GENIE3 Network Provides Biologically-Relevant Information in Polyploid Wheat.

Gene regulatory networks are powerful tools which facilitate hypothesis generation and candidate gene discovery. However, the extent to which the network predictions are biologically relevant is often unclear. Recently a GENIE3 network which predicted targets of wheat transcription factors was produced. Here we used an independent RNA-Seq dataset to test the predictions of the wheat GENIE3 network for the senescence-regulating transcription factor NAM-A1 (TraesCS6A02G108300). We re-analyzed the RNA-Seq data against the RefSeqv1.0 genome and identified a set of differentially expressed genes (DEGs) between the wild-type and nam-a1 mutant which recapitulated the known role of NAM-A1 in senescence and nutrient remobilisation. We found that the GENIE3-predicted target genes of NAM-A1 overlap significantly with the DEGs, more than would be expected by chance. Based on high levels of overlap between GENIE3-predicted target genes and the DEGs, we identified candidate senescence regulators. We then explored genome-wide trends in the network related to polyploidy and found that only homeologous transcription factors are likely to share predicted targets in common. However, homeologs which vary in expression levels across tissues are less likely to share predicted targets than those that do not, suggesting that they may be more likely to act in distinct pathways. This work demonstrates that the wheat GENIE3 network can provide biologically-relevant predictions of transcription factor targets, which can be used for candidate gene prediction and for global analyses of transcription factor function. The GENIE3 network has now been integrated into the KnetMiner web application, facilitating its use in future studies.

A roadmap for gene functional characterisation in crops with large genomes: Lessons from polyploid wheat.

Understanding the function of genes within staple crops will accelerate crop improvement by allowing targeted breeding approaches. Despite their importance, a lack of genomic information and resources has hindered the functional characterisation of genes in major crops. The recent release of high-quality reference sequences for these crops underpins a suite of genetic and genomic resources that support basic research and breeding. For wheat, these include gene model annotations, expression atlases and gene networks that provide information about putative function. Sequenced mutant populations, improved transformation protocols and structured natural populations provide rapid methods to study gene function directly. We highlight a case study exemplifying how to integrate these resources. This review provides a helpful guide for plant scientists, especially those expanding into crop research, to capitalise on the discoveries made in Arabidopsis and other plants. This will accelerate the improvement of crops of vital importance for food and nutrition security.