Identification of quantitative trait loci associated with leaf rust resistance in rye by precision mapping.

BACKGROUND: Leaf rust (LR) is among the most destructive fungal diseases of rye (Secale cereale L.). Despite intensive research using various analytical and methodological approaches, such as quantitative trait locus (QTL) mapping, candidate gene expression analysis, and transcriptome sequencing, the genetic basis of the rye immune response to LR remains unclear. RESULTS: A genome-wide association study was employed to detect QTLs controlling the immune response to LR of rye. A mapping population, G38A, was constructed by crossing two inbred lines: 723 (susceptible to LR) and JKI-NIL-Pr3 (a donor of the LR resistance gene Pr3). For genotyping, SNP-DArT and silico-DArT markers were used. Resistance phenotyping was conducted by visual assessment of the infection severity in detached leaf segments inoculated with two isolates of Puccinia recondita f. sp. secalis, namely, 60/17/2.1 (isolate S) in the main experiment and 86/n/2.1_5x (isolate N) in the validation experiment, at 10 and 17 days post-infection (dpi), respectively. In total, 42,773 SNP-DArT and 105,866 silico-DArT markers were included in the main analysis including isolate S, of which 129 and 140 SNP-DArTs and 767 and 776 silico-DArTs were significantly associated (p ≤ 0.001; - log10(p) ≥ 3.0) with the immune response to LR at 10 and 17 dpi, respectively. Most significant markers were mapped to chromosome 1R. The number of common markers from both systems and at both time points occupying common chromosomal positions was 37, of which 21 were positioned in genes, comprising 18 markers located in exons and three in introns. This gene pool included genes encoding proteins with a known function in response to LR (e.g., a NBS-LRR disease resistance protein-like protein and carboxyl-terminal peptidase). CONCLUSION: This study has expanded and supplemented existing knowledge of the genetic basis of rye resistance to LR by (1) detecting two QTLs associated with the LR immune response of rye, of which one located on the long arm of chromosome 1R is newly detected, (2) assigning hundreds of markers significantly associated with the immune response to LR to genes in the 'Lo7' genome, and (3) predicting the potential translational effects of polymorphisms of SNP-DArT markers located within protein-coding genes.

Genotype-Specific Expression of Selected Candidate Genes Conferring Resistance to Leaf Rust of Rye (Secale cereale L.).

Leaf rust (LR) caused by Puccinia recondita f. sp. secalis (Prs) is a highly destructive disease in rye. However, the genetic mechanisms underlying the rye immune response to this disease remain relatively uncharacterised. In this study, we analysed the expression of four genes in 12 rye inbred lines inoculated with Prs at 20 and 36 h post-treatment (hpt): DXS (1-deoxy-D-xylulose 5-phosphate synthase), Glu (ß-1,3-glucanase), GT (UDP-glycosyltransferase) and PR-1 (pathogenesis-related protein 1). The RT-qPCR analysis revealed the upregulated expression of the four genes in response to Prs in all inbred lines and at both time-points. The gene expression data were supported by microscopic and macroscopic examinations, which revealed that eight lines were susceptible to LR and four lines were highly resistant to LR. A relationship between the infection profiles and the expression of the analysed genes was observed: in the resistant lines, the expression level fold changes were usually higher at 20 hpt than at 36 hpt, while the opposite trend was observed in the susceptible lines. The study results indicate that DXS, Glu, GT and PR-1 may encode proteins crucial for the rye defence response to the LR pathogen.

Identification and characteristics of wheat Lr orthologs in three rye inbred lines.

The genetic background of the immune response of rye to leaf rust (LR), although extensively studied, is still not well understood. The recent publication of the genome of rye line Lo7 and the development of efficient transcriptomic methods has aided the search for genes that confer resistance to this disease. In this study, we investigated the potential role of rye orthologs of wheat Lr genes (Lr1, Lr10, Lr21, Lr22a, and RGA2/T10rga2-1A) in the LR seedling-stage resistance of inbred rye lines D33, D39, and L318. Bioinformatics analysis uncovered numerous Lr orthologs in the Lo7 genome, namely, 14 ScLr1, 15 ScRga2, and 2 ScLr21 paralogs, and 1 each of ScLr10 and ScLr22a genes. The paralogs of ScLr1, ScRga2, and ScLr21 were structurally different from one another and their wheat counterparts. According to an RNA sequencing analysis, only four wheat Lr gene orthologs identified in the Lo7 genome (ScLr1_3, ScLr1_4, ScLr1_8, and ScRga2_6) were differentially expressed; all four were downregulated after infection with compatible or incompatible isolates of Puccinia recondita f. sp. secalis (Prs). Using a more precise tool, RT-qPCR, we found that two genes were upregulated at 20 h post-infection, namely, ScLr1_4 and ScLr1_8 in lines D33 and D39, respectively, both of which have been found to be resistant to LR under field conditions and after treatment with a semi-compatible Prs strain. We were unable to discern any universal pattern of gene expression after Prs infection; on the contrary, all detected relationships were plant genotype-, Prs isolate-, or time-specific. Nevertheless, at least some Lr orthologs in rye (namely, ScLr1_3 ScLr1_4, ScLr1_8, and ScRga2_6), even though mainly downregulated, may play an important role in the response of rye to LR.

Identification of candidate genes responsible for chasmogamy in wheat.

BACKGROUND: The flowering biology of wheat plants favours self-pollination which causes obstacles in wheat hybrid breeding. Wheat flowers can be divided into two groups, the first one is characterized by flowering and pollination within closed flowers (cleistogamy), while the second one possesses the ability to open flowers during processes mentioned above (chasmogamy). The swelling of lodicules is involved in the flowering of cereals and among others their morphology, calcium and potassium content differentiate between cleistogamic and non-cleistogamous flowers. A better understanding of the chasmogamy mechanism can lead to the development of tools for selection of plants with the desired outcrossing rate. To learn more, the sequencing of transcriptomes (RNA-Seq) and Representational Difference Analysis products (RDA-Seq) were performed to investigate the global transcriptomes of wheat lodicules in two highly chasmogamous (HCH, Piko and Poezja) and two low chasmogamous (LCH, Euforia and KWS Dacanto) varieties at two developmental stages-pre-flowering and early flowering. RESULTS: The differentially expressed genes were enriched in five, main pathways: "metabolism", "organismal systems", "genetic information processing", "cellular processes" and "environmental information processing", respectively. Important genes with opposite patterns of regulation between the HCH and LCH lines have been associated with the lodicule development i.e. expression levels of MADS16 and MADS58 genes may be responsible for quantitative differences in chasmogamy level in wheat. CONCLUSIONS: We conclude that the results provide a new insight into lodicules involvement in the wheat flowering process. This study generated important genomic information to support the exploitation of the chasmogamy in wheat hybrid breeding programs.

Genome-Wide Association Analysis for Hybrid Breeding in Wheat.

Disclosure of markers that are significantly associated with plant traits can help develop new varieties with desirable properties. This study determined the genome-wide associations based on DArTseq markers for six agronomic traits assessed in eight environments for wheat. Moreover, the association study for heterosis and analysis of the effects of markers grouped by linkage disequilibrium were performed based on mean values over all experiments. All results were validated using data from post-registration trials. GWAS revealed 1273 single nucleotide polymorphisms with biologically significant effects. Most polymorphisms were predicted to be modifiers of protein translation, with only two having a more pronounced effect. Markers significantly associated with the considered set of features were clustered within chromosomes based on linkage disequilibrium in 327 LD blocks. A GWAS for heterosis revealed 1261 markers with significant effects.

Multiplexed SSR and agronomic data used in an investigation of obsolete diversity of rye (Secale cereale L.).

Rye (Secale cereale L.) is one of the most important cereal crops cultivated in the world due to its ability to produce high yields even when grown under environmental stress conditions. About 27,000 Secale accessions have been collected and preserved in 70 gene banks worldwide. Although the germplasm represents a great source of genetic diversity, the molecular characteristics refers only to the part of them. Here, we present data obtained by the Simple Sequence Repeat markers (SSR) analysis of 100 rye accessions preserved in the gene bank of the Polish Academy of Sciences Botanical Garden - Center for Biological Diversity Conservation in Powsin. Additionally, the data presented in this article refers to evaluation of agronomoic traits and weather conditions measured for 14 years. The data was used in the research article "Investigation of obsolete diversity of rye (Secale cereale L.) using multiplexed SSR fingerprinting and evaluation of agronomic traits" [1].

Identification of Rf Genes in Hexaploid Wheat (Triticumaestivum L.) by RNA-Seq and Paralog Analyses.

Among the natural mechanisms used for wheat hybrid breeding, the most desirable is the system combining the cytoplasmic male sterility (cms) of the female parent with the fertility-restoring genes (Rf) of the male parent. The objective of this study was to identify Rf candidate genes in the wheat genome on the basis of transcriptome sequencing (RNA-seq) and paralog analysis data. Total RNA was isolated from the anthers of two fertility-restorer (Primépi and Patras) and two non-restorer (Astoria and Grana) varieties at the tetrad and late uninucleate microspore stages. Of 36,912 differentially expressed genes (DEGs), 21 encoding domains in known fertility-restoring proteins were selected. To enrich the pool of Rf candidates, 52 paralogs (PAGs) of the 21 selected DEGs were included in the analyses. The expression profiles of most of the DEGs and PAGs determined bioinformatically were as expected (i.e., they were overexpressed in at least one fertility-restorer variety). However, these results were only partially consistent with the quantitative real-time PCR data. The DEG and PAG promoters included cis-regulatory elements common among PPR-encoding genes. On the basis of the obtained results, we designated seven genes as Rf candidate genes, six of which were identified for the first time in this study.

Barley Seeds miRNome Stability during Long-Term Storage and Aging.

Seed aging is a complex biological process that has been attracting scientists' attention for many years. High-throughput small RNA sequencing was applied to examine microRNAs contribution in barley seeds senescence. Unique samples of seeds that, despite having the same genetic makeup, differed in viability after over 45 years of storage in a dry state were investigated. In total, 61 known and 81 novel miRNA were identified in dry seeds. The highest level of expression was found in four conserved miRNA families, i.e., miR159, miR156, miR166, and miR168. However, the most astonishing result was the lack of significant differences in the level of almost all miRNAs in seed samples with significantly different viability. This result reveals that miRNAs in dry seeds are extremely stable. This is also the first identified RNA fraction that is not deteriorating along with the loss of seed viability. Moreover, the novel miRNA hvu-new41, with higher expression in seeds with the lowest viability as detected by RT-qPCR, has the potential to become an indicator of the decreasing viability of seeds during storage in a dry state.

The Roots of Rye (Secale cereale L.) Are Capable of Synthesizing Benzoxazinoids.

According to current opinion, the first step of benzoxazinoids (BXs) synthesis, that is, the conversion of indole-3-glycerol phosphate to indole, occurs exclusively in the photosynthesising parts of plants. However, the results of our previous work and some other studies suggest that this process may also occur in the roots. In this study, we provide evidence that the first step of BXs synthesis does indeed occur in the roots of rye seedlings. We detected ScBx1 transcripts, BX1 enzyme, and six BXs (2-hydroxy-1,4-benzoxazin-3-one, 2,4-dihydroxy-1,4-benzoxazin-3-one, (2R)-2-O-ß-d-glucopyranosyl-4-hydroxy-(2H)-1,4-benzoxazin-3(4H)-one glucoside, 2,4-dihydroxy- 7-methoxy-1,4-benzoxazin-3-one, 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one glucoside, and 6-methoxy-2-benzoxazolinone) in the roots developed from seeds deprived of the coleoptile at 2 days after sowing (i.e., roots without contact with aerial parts). In roots regenerated in vitro, both ScBx1 transcripts and BX1 enzyme were detected at a low but still measurable levels. Thus, BXs are able to be synthesised in both the roots and above-ground parts of rye plants.

Spikelet movements, anther extrusion and pollen production in wheat cultivars with contrasting tendencies to cleistogamy.

BACKGROUND: Cleistogamic flowers are a main barrier in pollen dispersal for cross-pollination necessary in wheat hybrid breeding. The aim of our study was to gain new knowledge on the biology of wheat flowering, in particular on the differences between the cleisto- and chasmogamic forms which has certainly cognitive significance, but it can also be used in practice when seeking a female and male ideotypes for cross hybridization. RESULTS: We characterized the most significant features defining the flowering specificity in two wheat cultivars with contrasting tendency to cleistogamy: Piko (chasmogamous) and Dacanto (cleistogamous). In the field observations we assessed diurnal pattern of anther extrusion and anther extrusion capacity. For the first time we adapted the time lapse method for measuring kinetics of the spikelet movement and 3-D image correlation technique for the non-invasive measurements of potential deformations of the spikelet lemmas. We found that the two cultivars differ in the potential of pollen dispersion for-cross-pollination and in the spikelet kinetics. We also described some anatomical traits that can have potential functional role in floret opening. None of the cultivars showed any symptoms of lemma surface deformation. CONCLUSIONS: The cleistogamic and chasmogamic wheat cultivars differ significantly in the potential for pollen dispersion for cross-pollination, which is mainly related to anther extrusion capacity. Although none of these features differentiated the cultivars clearly, we assume, based on spikelet kinetics and the lack of lemmas surface deformation, that the water transport and turgor of cells is essential for the floret opening and anther extrusion in wheat. The search for parental ideotype should be supported by marker assisted selection, e.g. based of polymorphisms in genes related to aquaporin biosynthesis.

Evaluation of genetic structure in European wheat cultivars and advanced breeding lines using high-density genotyping-by-sequencing approach.

BACKGROUND: The genetic diversity and gene pool characteristics must be clarified for efficient genome-wide association studies, genomic selection, and hybrid breeding. The aim of this study was to evaluate the genetic structure of 509 wheat accessions representing registered varieties and advanced breeding lines via the high-density genotyping-by-sequencing approach. RESULTS: More than 30% of 13,499 SNP markers representing 2162 clusters were mapped to genes, whereas 22.50% of 26,369 silicoDArT markers overlapped with coding sequences and were linked in 3527 blocks. Regarding hexaploidy, perfect sequence matches following BLAST searches were not sufficient for the unequivocal mapping to unique loci. Moreover, allelic variations in homeologous loci interfered with heterozygosity calculations for some markers. Analyses of the major genetic changes over the last 27 years revealed the selection pressure on orthologs of the gibberellin biosynthesis-related GA2 gene and the senescence-associated SAG12 gene. A core collection representing the wheat population was generated for preserving germplasm and optimizing breeding programs. CONCLUSIONS: Our results confirmed considerable differences among wheat subgenomes A, B and D, with D characterized by the lowest diversity but the highest LD. They revealed genomic regions that have been targeted by breeding.

Investigation of obsolete diversity of rye (Secale cereale L.) using multiplexed SSR fingerprinting and evaluation of agronomic traits.

Common rye (Secale cereale L.) is one of the most important cereals in Europe. Nevertheless, its germplasm collections are among the least numerous compared with cereals. There are only about 27,000 Secale accessions in 70 gene banks around the world. Despite extensive research on the molecular characterization of genetic resources, only a fraction of this collection has been described. The main objective of the presented study was to perform genotypic and phenotypic characterization of an obsolete gene pool represented by 100 accessions originated from 28 countries around the world and preserved in the gene bank of the Polish Academy of Sciences Botanical Garden - Center for Biological Diversity Conservation in Powsin. Genetic analysis using simple sequence repeat markers showed that the obsolete gene pool is relatively large. This indicates that different sources of variability were used in breeding programs. However, the genetic variation is in no way related to the place of origin. Despite the great differences in the genetic make-up, the collection showed a broadly common phenotype. This could result in a low level of interest among breeders in the stored germplasm, undervalued as a source of important but not easily observable traits, e.g., high disease resistance, which was found in some accessions.

Changes in benzoxazinoid contents and the expression of the associated genes in rye (Secale cereale L.) due to brown rust and the inoculation procedure.

Benzoxazinoids (BXs) are secondary metabolites with diverse functions, but are primarily involved in protecting plants, mainly from the family Poaceae, against insects and fungal pathogens. Rye is a cereal crop that is highly resistant to biotic stresses. However, its susceptibility to brown rust caused by Puccinia recondita f. sp. secalis (Prs) is still a major problem affecting its commercial production. Additionally, the genetic and metabolic factors related to this disease remain poorly characterized. In this study, we investigated whether and to what extent the brown rust infection and the inoculation procedure affect the contents of specific BXs (HBOA, GDIBOA, DIBOA, GDIMBOA, DIMBOA, and MBOA) and the expression of genes related to BX (ScBx1-5, ScIgl, and Scglu). We revealed that treatments with water and a urediniospore suspension usually downregulate gene expression levels. Moreover, HBOA and DIBOA contents decreased, whereas the contents of the remaining metabolites increased. Specifically, the MBOA content increased more after the mock treatment than after the Prs treatment, whereas the increase in GDIBOA and GDIMBOA levels was usually due to the Prs infection, especially at two of the most critical time-points, 17 and 24 h post-treatment. Therefore, GDIBOA and GDIMBOA are glucosides that are important components of rye defence responses to brown rust. Furthermore, along with MBOA, they protect rye against the stress associated with the inoculation procedure used in this study.

Genes ScBx1 and ScIgl-Competitors or Cooperators?

Two genes, Bx1 and Igl, both encoding indole-3-glycerol phosphate lyase (IGL), are believed to control the conversion of indole-3-glycerol phosphate (IGP) to indole. The first of these has generally been supposed to be regulated developmentally, being expressed at early stages of plant development with the indole being used in the benzoxazinoid (BX) biosynthesis pathway. In contrast, it has been proposed that the second one is regulated by stresses and that the associated free indole is secreted as a volatile. However, our previous results contradicted this. In the present study, we show that the ScIgl gene takes over the role of ScBx1 at later developmental stages, between the 42nd and 70th days after germination. In the majority of plants with silenced ScBx1 expression, ScIgl was either expressed at a significantly higher level than ScBx1 or it was the only gene with detectable expression. Therefore, we postulate that the synthesis of indole used in BX biosynthesis in rye is controlled by both ScBx1 and ScIgl, which are both regulated developmentally and by stresses. In silico and in vivo analyses of the promoter sequences further confirmed our hypothesis that the roles and modes of regulation of the ScBx1 and ScIgl genes are similar.

QTL mapping for benzoxazinoid content, preharvest sprouting, α-amylase activity, and leaf rust resistance in rye (Secale cereale L.).

Mapping population of recombinant inbred lines (RILs) representing 541 × Ot1-3 cross exhibited wide variations of benzoxazinoid (BX) content in leaves and roots, brown rust resistance, α-amylase activity in the grain, and resistance to preharvest sprouting. QTL mapping of major BX species using a DArT-based map revealed a complex genetic architecture underlying the production of these main secondary metabolites engaged in stress and allelopathy responses. The synthesis of BX in leaves and roots was found to be regulated by different QTL. The QTL for the BX content, rust resistance, α-amylase activity, and preharvest sprouting partially overlapped; this points to their common genetic regulation by a definite subset of genes. Only one QTL for BX located on chromosome 7R coincided with the loci of the ScBx genes, which were mapped as two clusters on chromosomes 5RS (Bx3-Bx5) and 7R (Bx1-Bx2). The QTL common for several BX species, rust resistance, preharvest sprouting, and α-amylase activity are interesting objects for further exploration aimed at developing common markers for these important agronomic traits.

Identification of Single Nucleotide Polymorphisms Associated with Brown Rust Resistance, α-Amylase Activity and Pre-harvest Sprouting in Rye (Secale cereale L.).

Rye is a crop with relatively high resistance to biotic and abiotic stresses. However, the resistance to brown rust (Puccinia recondita f. sp. secalis) and pre-harvest sprouting are still not satisfactory. High α-amylase activity is also among the main disadvantages of this species. Therefore, effective tools, e.g. molecular markers, allowing precise and environmentally independent selection of favourable alleles are desirable. In the present study, two kinds of association mapping-genome-wide association mapping (GWAM) based on sequences of DArTSeq markers and candidate gene association mapping (CGAM) based on sequences of ScBx genes-were chosen for development of molecular markers fulfilling these criteria. The analysed population consisted of 149 diverse inbred lines (DILs). Altogether, 67 and 11 single nucleotide polymorphisms (SNPs) identified in, respectively, GWAM and CGAM, were significantly associated with the investigated traits: 2 SNPs with resistance to brown rust, 71 SNPs with resistance to pre-harvest sprouting and 5 SNPs with α-amylase activity in the grain. Fifteen SNPs were stable across all environments. The highest number (13) of environmentally stable SNPs was associated with pre-harvest sprouting resistance. The test employing the Kompetitive Allele Specific PCR method proved the versatility of four markers identified in both GWAM and CGAM.

Identification and VIGS-based characterization of Bx1 ortholog in rye (Secale cereale L.).

The first step of the benzoxazinoid (BX) synthesis pathway is catalyzed by an enzyme with indole-3-glycerol phosphate lyase activity encoded by 3 genes, Bx1, TSA and Igl. A gene highly homologous to maize and wheat Bx1 has been identified in rye. The goal of the study was to analyze the gene and to experimentally verify its role in the rye BX biosynthesis pathway as a rye ortholog of the Bx1 gene. Expression of the gene showed peak values 3 days after imbibition (dai) and at 21 dai it was undetectable. Changes of the BX content in leaves were highly correlated with the expression pattern until 21 dai. In plants older than 21 dai despite the undetectable expression of the analyzed gene there was still low accumulation of BXs. Function of the gene was verified by correlating its native expression and virus-induced silencing with BX accumulation. Barley stripe mosaic virus (BSMV)-based vectors were used to induce transcriptional (TGS) and posttranscriptional (PTGS) silencing of the analyzed gene. Both strategies (PTGS and TGS) significantly reduced the transcript level of the analyzed gene, and this was highly correlated with lowered BX content. Inoculation with virus-based vectors specifically induced expression of the analyzed gene, indicating up-regulation by biotic stressors. This is the first report of using the BSMV-based system for functional analysis of rye gene. The findings prove that the analyzed gene is a rye ortholog of the Bx1 gene. Its expression is developmentally regulated and is strongly induced by biotic stress. Stable accumulation of BXs in plants older than 21 dai associated with undetectable expression of ScBx1 indicates that the function of the ScBx1 in the BX biosynthesis is redundant with another gene. We anticipate that the unknown gene is a putative ortholog of the Igl, which still remains to be identified in rye.

ScBx gene based association analysis of hydroxamate content in rye (Secale cereale L.).

Hydroxamates (HX) are major secondary metabolites synthesized by rye and are responsible for some of the unique properties of this cereal, including good tolerance of biotic and abiotic stresses and allelopathy. Recently, five genes encoding enzymes taking part in HX biosynthesis have been sequenced and characterized, which was the starting point to undertake the present study. Association analysis of the content of six HX-HBOA (2-hydroxy-1,4-benzoxazin-3-one), GDIBOA (2,4-dihydroxy-1,4-benzoxazin-3(4H)-one glucoside), DIBOA (2,4-dihydroxy-1,4-benzoxazin-3(4H)-one), GDIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3(4H)-one glucoside), DIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3(4H)-one) and MBOA (6-methoxy-benzoxazolin-2(3H)-one) in the above-ground parts of plants and roots was performed on a population consisting of 102 and 121 diverse inbred lines, in 2013 and 2014, respectively. Altogether, 48 single nucleotide polymorphisms (SNPs) were found to be associated with the content of at least one HX: 20 SNPs were associated with HX synthesized in the above-ground parts of rye plants (AG-SNP), and 28 were associated with HX synthesized in the roots (R-SNP). The highest number of SNPs was present in genes ScBx1 (9) and ScBx5 (11). The majority of SNPs were affected by environmental factors, except for two: ScBx4_1702 associated with GDIBOA and MBOA contents, and ScBx5_1105 associated with HBOA content in roots.