Multi-locus genome-wide association study of fusarium head blight in relation to days to anthesis and plant height in a spring wheat association panel

Fusarium head blight (FHB) is a highly destructive fungal disease of wheat to which host resistance is quantitatively inherited and largely influenced by the environment. Resistance to FHB has been associated with taller height and later maturity; however, a further understanding of these relationships is needed. An association mapping panel (AMP) of 192 predominantly Canadian spring wheat was genotyped with the wheat 90K single-nucleotide polymorphism (SNP) array. The AMP was assessed for FHB incidence (INC), severity (SEV) and index (IND), days to anthesis (DTA), and plant height (PLHT) between 2015 and 2017 at three Canadian FHB-inoculated nurseries. Seven multi-environment trial (MET) datasets were deployed in a genome-wide association study (GWAS) using a single-locus mixed linear model (MLM) and a multi-locus random SNP-effect mixed linear model (mrMLM). MLM detected four quantitative trait nucleotides (QTNs) for INC on chromosomes 2D and 3D and for SEV and IND on chromosome 3B. Further, mrMLM identified 291 QTNs: 50 (INC), 72 (SEV), 90 (IND), 41 (DTA), and 38 (PLHT). At two or more environments, 17 QTNs for FHB, DTA, and PLHT were detected. Of these 17, 12 QTNs were pleiotropic for FHB traits, DTA, and PLHT on chromosomes 1A, 1D, 2D, 3B, 5A, 6B, 7A, and 7B; two QTNs for DTA were detected on chromosomes 1B and 7A; and three PLHT QTNs were located on chromosomes 4B and 6B. The 1B DTA QTN and the three pleiotropic QTNs on chromosomes 1A, 3B, and 6B are potentially identical to corresponding quantitative trait loci (QTLs) in durum wheat. Further, the 3B pleiotropic QTN for FHB INC, SEV, and IND co-locates with TraesCS3B02G024900 within the Fhb1 region on chromosome 3B and is ~3 Mb from a cloned Fhb1 candidate gene TaHRC. While the PLHT QTN on chromosome 6B is putatively novel, the 1B DTA QTN co-locates with a disease resistance protein located ~10 Mb from a Flowering Locus T1-like gene TaFT3-B1, and the 7A DTA QTN is ~5 Mb away from a maturity QTL QMat.dms-7A.3 of another study. GWAS and QTN candidate genes enabled the characterization of FHB resistance in relation to DTA and PLHT. This approach should eventually generate additional and reliable trait-specific markers for breeding selection, in addition to providing useful information for FHB trait discovery.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Conditional Mapping Identified Quantitative Trait Loci for Grain Protein Concentration Expressing Independently of Grain Yield in Canadian Durum Wheat

Grain protein concentration (GPC) is an important trait in durum cultivar development as a major determinant of the nutritional value of grain and end-use product quality. However, it is challenging to simultaneously select both GPC and grain yield (GY) due to the negative correlation between them. To characterize quantitative trait loci (QTL) for GPC and understand the genetic relationship between GPC and GY in Canadian durum wheat, we performed both traditional and conditional QTL mapping using a doubled haploid (DH) population of 162 lines derived from Pelissier × Strongfield. The population was grown in the field over 5 years and GPC was measured. QTL contributing to GPC were detected on chromosome 1B, 2B, 3A, 5B, 7A, and 7B using traditional mapping. One major QTL on 3A (QGpc.spa-3A.3) was consistently detected over 3 years accounting for 9.4-18.1% of the phenotypic variance, with the favorable allele derived from Pelissier. Another major QTL on 7A (QGpc.spa-7A) detected in 3 years explained 6.9-14.8% of the phenotypic variance, with the beneficial allele derived from Strongfield. Comparison of the QTL described here with the results previously reported led to the identification of one novel major QTL on 3A (QGpc.spa-3A.3) and five novel minor QTL on 1B, 2B and 3A. Four QTL were common between traditional and conditional mapping, with QGpc.spa-3A.3 and QGpc.spa-7A detected in multiple environments. The QTL identified by conditional mapping were independent or partially independent of GY, making them of great importance for development of high GPC and high yielding durum.

High Density Mapping of Quantitative Trait Loci Conferring Gluten Strength in Canadian Durum Wheat

Gluten strength is one of the factors that determine the end-use quality of durum wheat and is an important breeding target for this crop. To characterize the quantitative trait loci (QTL) controlling gluten strength in Canadian durum wheat cultivars, a population of 162 doubled haploid (DH) lines segregating for gluten strength and derived from cv. Pelissier × cv. Strongfield was used in this study. The DH lines, parents, and controls were grown in 3 years and two seeding dates in each year and gluten strength of grain samples was measured by sodium dodecyl sulfate (SDS)-sedimentation volume (SV). With a genetic map created by genotyping the DH lines using the Illumina Infinium iSelect Wheat 90K SNP (single nucleotide polymorphism) chip, QTL contributing to gluten strength were detected on chromosome 1A, 1B, 2B, and 3A. Two major and stable QTL detected on chromosome 1A (QGlu.spa-1A) and 1B (QGlu.spa-1B.1) explaining 13.7-18.7% and 25.4-40.1% of the gluten strength variability respectively were consistently detected over 3 years, with the trait increasing alleles derived from Strongfield. Putative candidate genes underlying the major QTL were identified. Two novel minor QTL (QGlu.spa-3A.1 and QGlu.spa-3A.2) with the trait increasing allele derived from Pelissier were mapped on chromosome 3A explaining up to 8.9% of the phenotypic variance; another three minor QTL (QGlu.spa-2B.1, QGlu.spa-2B.2, and QGlu.spa-2B.3) located on chromosome 2B explained up to 8.7% of the phenotypic variance with the trait increasing allele derived from Pelissier. QGlu.spa-2B.1 is a new QTL and has not been reported in the literature. Multi-environment analysis revealed genetic (QTL) × environment interaction due to the difference of effect in magnitude rather than the direction of the QTL. Eleven pairs of digenic epistatic QTL were identified, with an epistatic effect between the two major QTL of QGlu.spa-1A and QGlu.spa-1B.1 detected in four out of six environments. The peak SNPs and SNPs flanking the QTL interval of QGlu.spa-1A and QGlu.spa-1B.1 were converted to Kompetitive Allele Specific PCR (KASP) markers, which can be deployed in marker-assisted breeding to increase the efficiency and accuracy of phenotypic selection for gluten strength in durum wheat. The QTL that were expressed consistently across environments are of great importance to maintain the gluten strength of Canadian durum wheat to current market standards during the genetic improvement.

Genetic analysis of resistance to stripe rust in durum wheat (Triticum turgidum L. var. durum)

Stripe rust, caused by the fungal pathogen Puccinia striiformis Westend. f. sp. tritici Eriks, is an important disease of bread wheat (Triticum aestivum L.) worldwide and there is an indication that it may also become a serious disease of durum wheat (T. turgidum L. var. durum). Therefore, we investigated the genetic architecture underlying resistance to stripe rust in adapted durum wheat germplasm. Wheat infection assays were conducted under controlled conditions in Canada and under field conditions in Mexico. Disease assessments were performed on a population of 155 doubled haploid (DH) lines derived from the cross of Kofa (susceptible) and W9262-260D3 (moderately resistant) and on a breeding panel that consisted of 92 diverse cultivars and breeding lines. Both populations were genotyped using the 90K single-nucleotide polymorphism (SNP) iSelect assay. In the DH population, QTL for stripe rust resistance were identified on chromosome 7B (LOD 6.87-11.47) and chromosome 5B (LOD 3.88-9.17). The QTL for stripe rust resistance on chromosome 7B was supported in the breeding panel. Both QTL were anchored to the genome sequence of wild emmer wheat, which identified gene candidates involved in disease resistance. Exome capture sequencing identified variation in the candidate genes between Kofa and W9262-260D3. These genetic insights will be useful in durum breeding to enhance resistance to stripe rust.

Genetics and mapping of seedling resistance to Ug99 stem rust in winter wheat cultivar Triumph 64 and differentiation of SrTmp, SrCad, and Sr42

Resistance to Ug99 stem rust in Triumph 64 was conferred by SrTmp on chromosome arm 6DS and was mapped to the same position as SrCad and Sr42 , however, the three genes show functional differences. Stem rust, caused by Puccinia graminis f. sp. tritici (Pgt), is an important disease of wheat that can be controlled by effective stem rust resistance (Sr) genes. The emergence of virulent Pgt races in Africa, namely Ug99 and its variants, has stimulated the search for new Sr genes and genetic characterization of known sources of resistance. Triumph 64 is a winter wheat cultivar that carries gene SrTmp, which confers resistance to Ug99. The goals of this study were to genetically map SrTmp and examine its relationship with other Sr genes occupying a similar chromosome location. A doubled haploid (DH) population from the cross LMPG-6S/Triumph 64 was inoculated with Ug99 at the seedling stage. A single gene conditioning resistance to Ug99 segregated in the population. Genetic mapping with SSR markers placed SrTmp on chromosome arm 6DS in a region similar to SrCad and Sr42. SNP markers developed for SrCad were used to further map SrTmp and were also added to a genetic map of Sr42 using a DH population (LMPG-6S/Norin 40). Three SNP markers that co-segregated with SrTmp also co-segregated with SrCad and Sr42. The SNP markers showed no difference in the map locations of SrTmp, SrCad, and Sr42. Multi-race testing with DH lines from the Triumph 64 and Norin 40 populations and a recombinant inbred-line population from the cross LMPG-6S/AC Cadillac showed that SrTmp, SrCad, and Sr42 confer different spectra of resistance. Markers closely linked to SrTmp are suitable for marker-assisted breeding and germplasm development.

Genetic mapping of SrCad and SNP marker development for marker-assisted selection of Ug99 stem rust resistance in wheat

New SNP markers that can be used for marker-assisted selection and map-based cloning saturate the chromosome region carrying SrCad , a wheat gene that confers resistance to Ug99 stem rust. Wheat stem rust, caused by Puccinia graminis f. sp. tritici, is a devastating disease of wheat worldwide. Development of cultivars with effective resistance has been the primary means to control this disease, but the appearance of new virulent strains such as Ug99 has rendered most wheat varieties vulnerable. The stem rust resistance gene SrCad located on chromosome arm 6DS has provided excellent resistance to various strains of Ug99 in field nurseries conducted in Njoro, Kenya since 2005. Three genetic populations were used to identify SNP markers closely linked to the SrCad locus. Of 220 SNP markers evaluated, 27 were found to be located within a 2 cM region surrounding SrCad. The diagnostic potential of these SNPs was evaluated in a diverse set of 50 wheat lines that were primarily of Canadian origin with known presence or absence of SrCad. Three SNP markers tightly linked proximally to SrCad and one SNP that co-segregated with SrCad were completely predictive of the presence or absence of SrCad. These markers also differentiated SrCad from Sr42 and SrTmp which are also located in the same region of chromosome arm 6DS. These markers should be useful in marker-assisted breeding to develop new wheat varieties containing SrCad-based resistance to Ug99 stem rust.

A Systems Genetics Approach Identifies Gene Regulatory Networks Associated with Fatty Acid Composition in Brassica rapa Seed

Fatty acids in seeds affect seed germination and seedling vigor, and fatty acid composition determines the quality of seed oil. In this study, quantitative trait locus (QTL) mapping of fatty acid and transcript abundance was integrated with gene network analysis to unravel the genetic regulation of seed fatty acid composition in a Brassica rapa doubled haploid population from a cross between a yellow sarson oil type and a black-seeded pak choi. The distribution of major QTLs for fatty acids showed a relationship with the fatty acid types: linkage group A03 for monounsaturated fatty acids, A04 for saturated fatty acids, and A05 for polyunsaturated fatty acids. Using a genetical genomics approach, expression quantitative trait locus (eQTL) hotspots were found at major fatty acid QTLs on linkage groups A03, A04, A05, and A09. An eQTL-guided gene coexpression network of lipid metabolism-related genes showed major hubs at the genes BrPLA2-ALPHA, BrWD-40, a number of seed storage protein genes, and the transcription factor BrMD-2, suggesting essential roles for these genes in lipid metabolism. Three subnetworks were extracted for the economically important and most abundant fatty acids erucic, oleic, linoleic, and linolenic acids. Network analysis, combined with comparison of the genome positions of cis- or trans-eQTLs with fatty acid QTLs, allowed the identification of candidate genes for genetic regulation of these fatty acids. The generated insights in the genetic architecture of fatty acid composition and the underlying complex gene regulatory networks in B. rapa seeds are discussed.

Localization of DIR1 at the tissue, cellular and subcellular levels during Systemic Acquired Resistance in Arabidopsis using DIR1:GUS and DIR1:EGFP reporters

BACKGROUND

Systemic Acquired Resistance (SAR) is an induced resistance response to pathogens, characterized by the translocation of a long-distance signal from induced leaves to distant tissues to prime them for increased resistance to future infection. DEFECTIVE in INDUCED RESISTANCE 1 (DIR1) has been hypothesized to chaperone a small signaling molecule to distant tissues during SAR in Arabidopsis.

RESULTS

DIR1 promoter:DIR1-GUS/dir1-1 lines were constructed to examine DIR1 expression. DIR1 is expressed in seedlings, flowers and ubiquitously in untreated or mock-inoculated mature leaf cells, including phloem sieve elements and companion cells. Inoculation of leaves with SAR-inducing avirulent or virulent Pseudomonas syringae pv tomato (Pst) resulted in Type III Secretion System-dependent suppression of DIR1 expression in leaf cells. Transient expression of fluorescent fusion proteins in tobacco and intercellular washing fluid experiments indicated that DIR1's ER signal sequence targets it for secretion to the cell wall. However, DIR1 expressed without a signal sequence rescued the dir1-1 SAR defect, suggesting that a cytosolic pool of DIR1 is important for the SAR response.

CONCLUSIONS

Although expression of DIR1 decreases during SAR induction, the protein localizes to all living cell types of the vasculature, including companion cells and sieve elements, and therefore DIR1 is well situated to participate in long-distance signaling during SAR.

Towards unambiguous transcript mapping in the allotetraploid Brassica napusThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

The architecture of the Brassica napus genome is marked by its evolutionary origins. The genome of B. napus was formed from the hybridization of two closely related diploid Brassica species, both of which evolved from an hexaploid ancestor. The extensive whole genome duplication events in its near and distant past result in the allotetraploid genome of B. napus maintaining multiple copies of most genes, which predicts a highly complex and redundant transcriptome that can confound any expression analyses. A stringent assembly of 142 399 B. napus expressed sequence tags allowed the development of a well-differentiated set of reference transcripts, which were used as a foundation to assess the efficacy of available tools for identifying and distinguishing transcripts in B. napus ; including microarray hybridization and 3′ anchored sequence tag capture. Microarray platforms cannot distinguish transcripts derived from the two progenitors or close homologues, although observed differential expression appeared to be biased towards unique transcripts. The use of 3′ capture enhanced the ability to unambiguously identify homologues within the B. napus transcriptome but was limited by tag length. The ability to comprehensively catalogue gene expression in polyploid species could be transformed by the application of cost-efficient next generation sequencing technologies that will capture millions of long sequence tags.

A novel role for protein farnesylation in plant innate immunity

Plants utilize tightly regulated mechanisms to defend themselves against pathogens. Initial recognition results in activation of specific Resistance (R) proteins that trigger downstream immune responses, in which the signaling networks remain largely unknown. A point mutation in SUPPRESSOR OF NPR1 CONSTITUTIVE1 (SNC1), a RESISTANCE TO PERONOSPORA PARASITICA4 R gene homolog, renders plants constitutively resistant to virulent pathogens. Genetic suppressors of snc1 may carry mutations in genes encoding novel signaling components downstream of activated R proteins. One such suppressor was identified as a novel loss-of-function allele of ENHANCED RESPONSE TO ABSCISIC ACID1 (ERA1), which encodes the beta-subunit of protein farnesyltransferase. Protein farnesylation involves attachment of C15-prenyl residues to the carboxyl termini of specific target proteins. Mutant era1 plants display enhanced susceptibility to virulent bacterial and oomycete pathogens, implying a role for farnesylation in basal defense. In addition to its role in snc1-mediated resistance, era1 affects several other R-protein-mediated resistance responses against bacteria and oomycetes. ERA1 acts partly independent of abscisic acid and additively with the resistance regulator NON-EXPRESSOR OF PR GENES1 in the signaling network. Defects in geranylgeranyl transferase I, a protein modification similar to farnesylation, do not affect resistance responses, indicating that farnesylation is most likely specifically required in plant defense signaling. Taken together, we present a novel role for farnesyltransferase in plant-pathogen interactions, suggesting the importance of protein farnesylation, which contributes to the specificity and efficacy of signal transduction events.

Development of a Brassica seed cDNA microarray

Brassica species represent several important crops including canola (Brassica napus). Understanding of genetic elements that contribute to seed-associated functions will impact future improvements in the canola crop. Brassica species share a very close taxonomic and molecular relationship with Arabidopsis thaliana. However, there are several subtle but distinct seed-associated agronomic characteristics that differ among the oil seed crop species. To address these, we have generated 67,535 ESTs predominately from Brassica seeds, analyzed these sequences, and identified 10,642 unigenes for the preparation of a targeted seed cDNA array. A set of 10,642 PCR primer pairs was designed and corresponding amplicons were produced for spotting, along with relevant controls. Critical quality control tests produced satisfactory results for use of this microarray in biological experiments. The microarray was also tested with specific RNA targets from embryos, germinating seeds, and leaf tissues. The hybridizations, signal intensities, and overall quality of these slides were consistent and reproducible. Additionally, there are 429 ESTs represented on the array that show no homology with any A. thaliana annotated gene or any gene in the Brassica genome databases or other plant databases; however, all of these probes hybridized to B. napus transcripts, indicating that the array also will be useful in defining expression patterns for genes so far unique to Brassica species.

Characterization of an AGAMOUS homologue from the conifer black spruce (Picea mariana) that produces floral homeotic conversions when expressed in Arabidopsis

Advances in elucidating the molecular processes controlling flower initiation and development have provided unique opportunities to investigate the developmental genetics of non-flowering plants. In addition to providing insights into the evolutionary aspects of seed plants, identification of genes regulating reproductive organ development in gymnosperms could help determine the level of homology with current models of flower induction and floral organ identity. Based upon this, we have searched for putative developmental regulators in conifers with amino acid sequence homology to MADS-box genes. PCR cloning using degenerate primers targeted to the MADS-box domain revealed the presence of over 27 MADS-box genes within black spruce (Picea mariana), including several with extensive homology to either AP1 or AGAMOUS, both known to regulate flower development in Arabidopsis. This indicates that like angiosperms, conifers contain a large and diverse MADS-box gene family that probably includes regulators of reproductive organ development. Confirmation of this was provided by the characterization of an AGAMOUS-like cDNA clone called SAG1, whose conservation of intron position and tissue-specific expression within reproductive organs indicate that it is a homologue of AGAMOUS. Functional homology with AGAMOUS was demonstrated by the ability of SAG1 to produce homeotic conversions of sepals to carpels and petals to stamens when ectopically expressed in transgenic Arabidopsis. This suggests that some of the genetic pathways controlling flower and cone development are homologous, and antedate the 300-million-year-old divergence of angiosperms and gymnosperms.

Conserved and novel regulators of the plant cell cycle

Cell division is highly regulated, both spatially and temporally, during plant development. Recent evidence implicates cyclin-dependent kinases (cdks) and their associated proteins as the principal temporal regulators of cell division. It is now known that plants contain an extended family of cdks, some of which appear to be unique to this group. Positive rate-limiting regulators of cell proliferation and growth include mitotic or B-type cyclins whose transcription is restricted to the G2 and M phases. Current research suggests that MYB-related transcription factors may be responsible for this restriction. Cdk-interacting proteins, such as cdk inhibitors and suc1 homologues, have been isolated using yeast two-hybrid approaches.

A nimA-like protein kinase transcript is highly expressed in meristems of Antirrhinum majus

In plants, cell proliferation occurs mostly within meristems but a significant amount also occurs at other well-defined sites during specific stages of development. We have developed molecular markers to follow the location and progress of cell division within multicellular plant organs and thereby gain some insight into how cell division might be regulated during morphogenesis. As in other eukaryotes, cell division in plants is regulated by a highly conserved set of protein kinases and phosphatases. The molecular information available on these molecules from other eukaryotes has allowed the design of strategies by which plant homologues can be isolated. In this report we describe the identification of a nimA-like gene from Antirrhinum majus and describe the pattern in which its transcript is expressed. Comparison of the pattern of AmnimA gene expression with that of genes which are expressed in a cell cycle-dependent manner suggests that this gene is expressed in actively dividing tissues but expression is not specific to any particular phase of the cell cycle nor specific to any particular tissue type.

Interactions of growth regulators, explant age, and culture environment controlling organogenesis from lettuce cotyledons in vitro

Organogenesis in cultured Lactuca sativa L. cv. Grand Rapids cotyledons was controlled by growth regulators and influenced by explant age as well as culture environment. Control explants formed roots with a low frequency. Exposure to indoleacetic acid (IAA) greatly stimulated root initiation. Kinetin (K) applied singly or in combination with IAA induced bud formation. Benzylaminopurine (BA) also caused bud initiation. Root and shoot initiation was optimal with 4-day-old explants. Organogenesis declined with older explants. More buds formed in light than in darkness. Other environmental variables tested had no marked effect on K-induced bud initiation. Control explants on medium containing no growth regulator never formed buds.