Detection of the quantitative trait loci for α-amylase activity on a high-density genetic map of rye and comparison of their localization to loci controlling preharvest sprouting and earliness

Construction of a high-density genetic map and mapping of a spike length locus for rye

Genetic maps provide the foundation for QTL mapping of important traits of crops. As a valuable food and forage crop, rye (Secale cereale L., RR) is also one of the tertiary gene sources of wheat, especially wild rye, Secale cereale subsp. segetale, possessing remarkable stress tolerance, tillering capacity and numerous valuable traits. In this study, based on the technique of specific-locus amplified fragment sequencing (SLAF-seq), a high-density single nucleotide polymorphism (SNP) linkage map of the cross-pollinated (CP) hybrid population crossed by S. cereale L (female parent) and S. cereale subsp. segetale (male parent) was successfully constructed. Following preprocessing, the number of 1035.11 M reads were collected and 2425800 SNP were obtained, of which 409134 SNP were polymorphic. According to the screening process, 9811 SNP markers suitable for constructing linkage groups (LGs) were selected. Subsequently, all of the markers with MLOD values lower than 3 were filtered out. Finally, an integrated map was constructed with 4443 markers, including 1931 female mapping markers and 3006 male mapping markers. A major quantitative trait locus (QTL) linked with spike length (SL) was discovered at 73.882 cM on LG4, which explained 25.29% of phenotypic variation. Meanwhile two candidate genes for SL, ScWN4R01G329300 and ScWN4R01G329600, were detected. This research presents the first high-quality genetic map of rye, providing a substantial number of SNP marker loci that can be applied to marker-assisted breeding. Additionally, the finding could help to use SLAF marker mapping to identify certain QTL contributing to important agronomic traits. The QTL and the candidate genes identified through the high-density genetic map above may provide diverse potential gene resources for the genetic improvement of rye.

How Machine Learning Methods Helped Find Putative Rye Wax Genes Among GBS Data

The standard approach to genetic mapping was supplemented by machine learning (ML) to establish the location of the rye gene associated with epicuticular wax formation (glaucous phenotype). Over 180 plants of the biparental F2 population were genotyped with the DArTseq (sequencing-based diversity array technology). A maximum likelihood (MLH) algorithm (JoinMap 5.0) and three ML algorithms: logistic regression (LR), random forest and extreme gradient boosted trees (XGBoost), were used to select markers closely linked to the gene encoding wax layer. The allele conditioning the nonglaucous appearance of plants, derived from the cultivar Karlikovaja Zelenostebelnaja, was mapped at the chromosome 2R, which is the first report on this localization. The DNA sequence of DArT-Silico 3585843, closely linked to wax segregation detected by using ML methods, was indicated as one of the candidates controlling the studied trait. The putative gene encodes the ABCG11 transporter.

The GAMYB gene in rye: sequence, polymorphisms, map location, allele-specific markers, and relationship with α-amylase activity

BACKGROUND

Transcription factor (TF) GAMYB, belonging to MYB family (named after the gene of the avian myeloblastosis virus) is a master gibberellin (GA)-induced regulatory protein that is crucial for development and germination of cereal grain and involved in anther formation. It activates many genes including high-molecular-weight glutenin and α-amylase gene families. This study presents the first attempt to characterize the rye gene encoding GAMYB in relation to its sequence, polymorphisms, and phenotypic effects.

RESULTS

ScGAMYB was mapped on rye chromosome 3R using high-density Diversity Arrays Technology (DArT)/DArTseq-based maps developed in three mapping populations. The ScGAMYB sequences were identified in RNA-seq libraries of four rye inbred lines. The transcriptome used for the search contained almost 151,000 transcripts with a median contig length of 500 nt. The average amount of total base raw data was approximately 9 GB. Comparative analysis of the ScGAMYB sequence revealed its high level of homology to wheat and barley orthologues. Single nucleotide polymorphisms (SNPs) detected among rye inbred lines allowed the development of allele specific-PCR (AS-PCR) markers for ScGAMYB that might be used to detect this gene in wide genetic stocks of rye and triticale. Segregation of the ScGAMYB alleles showed significant relationship with α-amylase activity (AMY).

CONCLUSIONS

The research showed the strong similarity of rye GAMYB sequence to its orthologues in other Graminae and confirmed the position in the genome consistent with the collinearity rule of cereal genomes. Concurrently, the ScGAMYB coding sequence (cds) showed stronger variability (24 SNPs) compared to the analogous region of wheat (5 SNPs) and barley (7 SNPs). The moderate regulatory effect of ScGAMYB on AMY was confirmed, therefore, ScGAMYB was identified as a candidate gene for partial control of α-amylase production in rye grain. The predicted structural protein change in the aa region 362-372, caused by a single SNP (C/G) at the 1100 position in ScGAMYB cds and single aa sequence change (S/C) at the 367 position, is the likely cause of the differences in the effectiveness of ScGAMYB regulatory function associated with AMY. The development of sequence-based, allele-specific (AS) PCR markers could be useful in research and application.

Identification and mapping of a new recessive dwarfing gene dw9 on the 6RL rye chromosome and its phenotypic effects

The introduction of high-yielding semi-dwarf varieties of wheat into cultivation has led to a "green revolution." This has required intensive research into various sources of dwarfism in wheat. However, there has been very little advancement in research on dwarfing genes in rye in comparison to wheat or barley. So far, three dominant dwarfing genes (Ddw1, Ddw3, and Ddw4) and three recessive genes (ct1, ct2, and np) have been characterized and precisely mapped in rye. There is no complete catalog of dwarfing genes available in rye. This paper presents an identification of the source of dwarfism and preliminary characterization of the new recessive gene dw9 from the BK-1 line. The gene was mapped on the long arm of the 6R chromosome and belongs to the GA-insensitive group. The initial characterization of the influence of this gene on morphological traits shows that it significantly affects the decrease of yielding trait parameters. A full evaluation can be performed after detailed breeding studies.

Putative candidate genes responsible for leaf rolling in rye (Secale cereale L.)

BACKGROUND

Rolling of leaves (RL) is a phenomenon commonly found in grasses. Morphology of the leaf is an important agronomic trait in field crops especially in rice; therefore, majority of the rice breeders are interested in RL. There are only few studies with respect to RL of wheat and barley; however, the information regarding the genetic base of RL with respect to the shape of leaf in rye is lacking. To the best of our knowledge, this is the first study on the localization of loci controlling RL on high density consensus genetic map of rye.

RESULTS

Genotypic analysis led to the identification of 43 quantitative trait loci (QTLs) for RL, grouped into 28 intervals, which confirms the multigenic base of the trait stated for wheat and rice. Four stable QTLs were located on chromosomes 3R, 5R, and 7R. Co-localization of QTL for RL and for different morphological, biochemical and physiological traits may suggests pleiotropic effects of some QTLs. QTLs for RL were associated with QTLs for such morphological traits as: grain number and weight, spike number per plant, compactness of spike, and plant height. Two QTLs for RL were found to coincide with QTLs for drought tolerance (4R, 7R), two with QTLs for heading earliness (2R, 7R), one with α-amylase activity QTL (7R) and three for pre-harvest sprouting QTL (1R, 4R, 7R). The set of molecular markers strongly linked to RL was selected, and the putative candidate genes controlling the process of RL were identified. Twelve QTLs are considered as linked to candidate genes on the base of DArT sequences alignment, which is a new information for rye.

CONCLUSIONS

Our results expand the knowledge about the network of QTLs for different morphological, biochemical and physiological traits and can be a starting point to studies on particular genes controlling RL and other important agronomic traits (yield, earliness, pre-harvest sprouting, reaction to water deficit) and to appoint markers useful in marker assisted selection (MAS). A better knowledge of the rye genome and genes could both facilitate rye improvement itself and increase the efficiency of utilizing rye genes in wheat breeding.

Identification of a novel, dominant dwarfing gene (Ddw4) and its effect on morphological traits of rye

Shortening rye stems to improve lodging resistance is among the major tasks awaiting breeders of this cereal. The most straightforward way to achieve this goal is the implementation of a dominant dwarfing gene into high yielding cultivars. The choice of dominant dwarfing genes in rye is limited to Ddw1 and Ddw3 loci, which are well characterized with respect to map position and tightly linked molecular markers on the long arms of chromosomes 5RL and 1RL, respectively. This paper reports on the identification and preliminary characterization of a novel dominant dwarfing gene, Ddw4, from line S44. This was mapped within the centromeric region of chromosome 3R. The Ddw4 gene is sensitive to exogenous gibberellin. Its introduction into the rye populational cultivar Dańkowskie Amber decreased plant height by c. 54% without any negative effects on spike length and number of kernels per spike. Further genetic studies are needed to determine the perspectives for application of the newly detected dwarfing gene into breeding programs for short-stem rye.

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.

QTL mapping and comparative genome analysis of agronomic traits including grain yield in winter rye

Genetic diversity in elite rye germplasm as well as F _2:3 testcross design enables fast QTL mapping to approach genes controlling grain yield, grain weight, tiller number and heading date in rye hybrids. Winter rye (Secale cereale L.) is a multipurpose cereal crop closely related to wheat, which offers the opportunity for a sustainable production of food and feed and which continues to emerge as a renewable energy source for the production of bioethanol and biomethane. Rye contributes to increase agricultural crop species diversity particularly in Central and Eastern Europe. In contrast to other small grain cereals, knowledge on the genetic architecture of complex inherited, agronomic important traits is yet limited for the outbreeding rye. We have performed a QTL analysis based on a F_2:3 design and testcross performance of 258 experimental hybrids in multi-environmental field trials. A genetic linkage map covering 964.9 cM based on SSR, conserved-orthologous set (COS), and mixed-phase dominant DArT markers allowed to describe 22 QTL with significant effects for grain yield, heading date, tiller number, and thousand grain weight across seven environments. Using rye COS markers, orthologous segments for these traits have been identified in the rice genome, which carry cloned and functionally characterized rice genes. The initial genome scan described here together with the existing knowledge on candidate genes provides the basis for subsequent analyses of the genetic and molecular mechanisms underlying agronomic important traits in rye.

Comparative analysis of genetic architectures for nine developmental traits of rye

Genetic architectures of plant height, stem thickness, spike length, awn length, heading date, thousand-kernel weight, kernel length, leaf area and chlorophyll content were aligned on the DArT-based high-density map of the 541 × Ot1–3 RILs population of rye using the genes interaction assorting by divergent selection (GIABDS) method. Complex sets of QTL for particular traits contained 1–5 loci of the epistatic D class and 10–28 loci of the hypostatic, mostly R and E classes controlling traits variation through D–E or D–R types of two-loci interactions. QTL were distributed on each of the seven rye chromosomes in unique positions or as a coinciding loci for 2–8 traits. Detection of considerable numbers of the reversed (D′, E′ and R′) classes of QTL might be attributed to the transgression effects observed for most of the studied traits. First examples of E* and F QTL classes, defined in the model, are reported for awn length, leaf area, thousand-kernel weight and kernel length. The results of this study extend experimental data to 11 quantitative traits (together with pre-harvest sprouting and alpha-amylase activity) for which genetic architectures fit the model of mechanism underlying alleles distribution within tails of bi-parental populations. They are also a valuable starting point for map-based search of genes underlying detected QTL and for planning advanced marker-assisted multi-trait breeding strategies.

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.

DArT Markers Effectively Target Gene Space in the Rye Genome

Large genome size and complexity hamper considerably the genomics research in relevant species. Rye (Secale cereale L.) has one of the largest genomes among cereal crops and repetitive sequences account for over 90% of its length. Diversity Arrays Technology is a high-throughput genotyping method, in which a preferential sampling of gene-rich regions is achieved through the use of methylation sensitive restriction enzymes. We obtained sequences of 6,177 rye DArT markers and following a redundancy analysis assembled them into 3,737 non-redundant sequences, which were then used in homology searches against five Pooideae sequence sets. In total 515 DArT sequences could be incorporated into publicly available rye genome zippers providing a starting point for the integration of DArT- and transcript-based genomics resources in rye. Using Blast2Go pipeline we attributed putative gene functions to 1101 (29.4%) of the non-redundant DArT marker sequences, including 132 sequences with putative disease resistance-related functions, which were found to be preferentially located in the 4RL and 6RL chromosomes. Comparative analysis based on the DArT sequences revealed obvious inconsistencies between two recently published high density consensus maps of rye. Furthermore we demonstrated that DArT marker sequences can be a source of SSR polymorphisms. Obtained data demonstrate that DArT markers effectively target gene space in the large, complex, and repetitive rye genome. Through the annotation of putative gene functions and the alignment of DArT sequences relative to reference genomes we obtained information, that will complement the results of the studies, where DArT genotyping was deployed, by simplifying the gene ontology and microcolinearity based identification of candidate genes.

Bidirectional selective genotyping approach for the identification of quantitative trait loci controlling earliness per se in winter rye (Secale cereale L.)

The genes controlling earliness of plants include genes responsible for vernalisation (Vrn) and photoperiod (Ppd), and those that are not entirely associated with a response to temperature or light. The last group of loci is known as earliness per se (Eps). Eps genes have been most commonly reported in the scientific literature as quantitative trait loci (QTL). The objective of this study was to use a bidirectional selective genotyping (BSG) method for the identification of loci controlling Eps, assessed at the heading stage in a newly developed rye population of recombinant inbred lines. We identified four linkage groups representing fragments of chromosomes, 1R, 5R, 6R and 7R, and containing, in total, 155 DArT markers significantly associated with earliness in rye. The results are discussed in the context of previously published QTL analyses.

Genetic mapping of the ScHd1 gene in rye and an assessment of its relationship with earliness per se and plant morphology

A fragment of the ScHd1 gene derived from eight inbred lines was sequenced and showed homology to other Hd1 genes from different cereals. Sequences were analysed with respect to the presence of a single-nucleotide polymorphism (SNP) difference. A C-T transition at position 312 of the consensus sequence was found, which distinguished two lines from the remaining six. The deduced amino acid sequence revealed a high identity (93%) to a Hd1-like protein from wheat. The identified mutation allowed the localisation of ScHd1 on a genetic map of rye (6RS). A small, statistically significant linkage between ScHd1 and earliness per se (eps) and some morphological traits was also established. The chromosomal region, including the S76 allele for the ScHd1 gene was linked to earlier heading, elongated spikes, a greater number of spikelets per spike and an increased weight of 1000 kernels.

Comparative QTL analysis of early short-time drought tolerance in Polish fodder and malting spring barleys

An effective approach for the further evolution of QTL markers, may be to create mapping populations for locally adapted gene pools, and to phenotype the studied trait under local conditions. Mapping populations of Polish fodder and malting spring barleys (Hordeum vulgare L.) were used to analyze traits describing short-time drought response at the seedlings stage. High-throughput genotyping (Diversity Array Technology (DArT) markers) and phenotyping techniques were used. The results showed high genetic diversity of the studied populations which allowed the creation of high-density linkage maps. There was also high diversity in the physiological responses of the barleys. Quantitative trait locus (QTL) analysis revealed 18 QTLs for nine physiological traits on all chromosomes except 1H in malting barley and 15 QTLs for five physiological traits on chromosomes 2H, 4H, 5H and 6H in fodder barley. Chromosomes 4H and 5H contained QTLs which explained most of the observed phenotypic variations in both populations. There was a major QTL for net photosynthetic rate in the malting barley located on chromosome 5H and two major QTLs for overall photochemical performance (PI) located on 5H and 7H. One major QTL related to photochemical quenching of chlorophyll fluorescence was located on chromosome 4H in fodder barley. Three QTL regions were common to both mapping populations but the corresponding regions explained different drought-induced traits. One region was for QTLs related to PSII photosynthetic activity stress index in malting barley, and the corresponding region in fodder barley was related to the water content stress index. These results are in accordance with previous studies which showed that different traits were responsible for drought tolerance variations in fodder and malting barleys.

Genetic architecture of complex agronomic traits examined in two testcross populations of rye (Secale cereale L.)

Background

Rye is an important European crop used for food, feed, and bioenergy. Several quality and yield-related traits are of agronomic relevance for rye breeding programs. Profound knowledge of the genetic architecture of these traits is needed to successfully implement marker-assisted selection programs. Nevertheless, little is known on quantitative loci underlying important agronomic traits in rye.

Results

We used 440 F_3:4 inbred lines from two biparental populations (Pop-A, Pop-B) fingerprinted with about 800 to 900 SNP, SSR and/or DArT markers and outcrossed them to a tester for phenotyping. The resulting hybrids and their parents were evaluated for grain yield, single-ear weight, test weight, plant height, thousand-kernel weight, falling number, protein, starch, soluble and total pentosan contents in up to ten environments in Central Europe. The quality of the phenotypic data was high reflected by moderate to high heritability estimates. QTL analyses revealed a total of 31 QTL for Pop-A and 52 for Pop-B. QTL x environment interactions were significant (P < 0.01) in most cases but variance of QTL main effect was more prominent.

Conclusions

QTL mapping was successfully applied based on two segregating rye populations. QTL underlying grain yield and several quality traits had small effects. In contrast, thousand-kernel weight, test weight, falling number and starch content were affected by several major QTL with a high frequency of occurrence in cross validation. These QTL explaining a large proportion of the genotypic variance can be exploited in marker-assisted selection programs and are candidates for further genetic dissection.

Genetic architecture of complex agronomic traits examined in two testcross populations of rye (Secale cereale L.)

BACKGROUND

Rye is an important European crop used for food, feed, and bioenergy. Several quality and yield-related traits are of agronomic relevance for rye breeding programs. Profound knowledge of the genetic architecture of these traits is needed to successfully implement marker-assisted selection programs. Nevertheless, little is known on quantitative loci underlying important agronomic traits in rye.

RESULTS

We used 440 F(3:4) inbred lines from two biparental populations (Pop-A, Pop-B) fingerprinted with about 800 to 900 SNP, SSR and/or DArT markers and outcrossed them to a tester for phenotyping. The resulting hybrids and their parents were evaluated for grain yield, single-ear weight, test weight, plant height, thousand-kernel weight, falling number, protein, starch, soluble and total pentosan contents in up to ten environments in Central Europe. The quality of the phenotypic data was high reflected by moderate to high heritability estimates. QTL analyses revealed a total of 31 QTL for Pop-A and 52 for Pop-B. QTL x environment interactions were significant (P < 0.01) in most cases but variance of QTL main effect was more prominent.

CONCLUSIONS

QTL mapping was successfully applied based on two segregating rye populations. QTL underlying grain yield and several quality traits had small effects. In contrast, thousand-kernel weight, test weight, falling number and starch content were affected by several major QTL with a high frequency of occurrence in cross validation. These QTL explaining a large proportion of the genotypic variance can be exploited in marker-assisted selection programs and are candidates for further genetic dissection.

A high density consensus map of rye (Secale cereale L.) based on DArT markers

Background

Rye (Secale cereale L.) is an economically important crop, exhibiting unique features such as outstanding resistance to biotic and abiotic stresses and high nutrient use efficiency. This species presents a challenge to geneticists and breeders due to its large genome containing a high proportion of repetitive sequences, self incompatibility, severe inbreeding depression and tissue culture recalcitrance. The genomic resources currently available for rye are underdeveloped in comparison with other crops of similar economic importance. The aim of this study was to create a highly saturated, multilocus linkage map of rye via consensus mapping, based on Diversity Arrays Technology (DArT) markers.

Methodology/Principal Findings

Recombinant inbred lines (RILs) from 5 populations (564 in total) were genotyped using DArT markers and subjected to linkage analysis using Join Map 4.0 and Multipoint Consensus 2.2 software. A consensus map was constructed using a total of 9703 segregating markers. The average chromosome map length ranged from 199.9 cM (2R) to 251.4 cM (4R) and the average map density was 1.1 cM. The integrated map comprised 4048 loci with the number of markers per chromosome ranging from 454 for 7R to 805 for 4R. In comparison with previously published studies on rye, this represents an eight-fold increase in the number of loci placed on a consensus map and a more than two-fold increase in the number of genetically mapped DArT markers.

Conclusions/Significance

Through the careful choice of marker type, mapping populations and the use of software packages implementing powerful algorithms for map order optimization, we produced a valuable resource for rye and triticale genomics and breeding, which provides an excellent starting point for more in-depth studies on rye genome organization.