Background {#Sec1}
==========

Leaf rolling (RL) is a typical response of a plant during water deficit that is observed in various field crops such as rice, maize, wheat, and sorghum. It decreases transpiration by decreasing the effective leaf area, and thus is a potentially useful drought tolerance mechanism in dry areas \[[@CR1]\]. Although RL is a phenomenon commonly found in grasses, it has attracted much attention from rice researchers and breeders \[[@CR2]\]. Only a few studies reported on RL of wheat \[[@CR3]--[@CR5]\] and barley \[[@CR6]\], but studies regarding RL of rye are lacking.

Flag RL function is an important drought tolerance mechanism, enabling the plant to conserve water by decreasing transpiration during water stress and reduce leaf temperature \[[@CR4]\]. Some Mediterranean grasses decrease transpiration as much as 46 to 63% by RL. In many species, RL does not occur until the water content in the leaf decreases to lethal levels (Parker 1968, after \[[@CR7]\]).

In rice, RL is classified as abaxial leaf roll (both sides of the leaf roll inward along the vein) and adaxial leaf roll (both sides of the leaf roll outward along the vein) according to the direction of RL \[[@CR8]\]. Analysis of some of the leaf development mutants of *Arabidopsis thaliana* and maize has shown that some mutations with respect to RL are related to the development of the leaf along the adaxial--abaxial axis. Establishment of leaf polarity and cell differentiation affecting leaf shape are controlled by both transcription factors and small RNAs \[[@CR2], [@CR9]\]. In higher plants, two types of cells are involved in RL: bulliform and hypodermis cells. Bulliform cells, which are located in the upper epidermis of the leaf near the midrib or vascular bundles of leaves, cause rolling in some *Gramineae* species such as rice, maize, wheat, and *Sorghum* spp. \[[@CR1]\].

Due to its importance, many studies have been performed to characterize the genes controlling RL in rice. To date, no less than 17 rice mutants with rolled leaves have been characterized \[[@CR10]\], no fewer than 70 genes/QTLs for the rolled leaf trait have been mapped or cloned throughout the rice genome \[[@CR11]\], and at least 28 differentially expressed proteins related to RL traits have been isolated and identified \[[@CR12]\]. A study on durum × wild emmer wheat recombinant inbred line population has reported 11 significant QTLs associated with flag RL mechanism \[[@CR4]\]. In this study, we aimed to localize the loci controlling RL on high density consensus genetic map of rye. We also aimed to detect QTL co-localization of RL and other agronomic traits.

Methods {#Sec2}
=======

Plant material and genetic map {#Sec3}
------------------------------

In this study, we used a population of recombinant inbred lines (RILs) namely, RIL-M, which is a cross between S120 and S76 lines. S120 and S76 were developed within the commercial breeding programs conducted at Danko Plant Breeding Ltd. (Choryń, Poland) and are partially related but are genetically different \[[@CR13]\] with respect to the following: α-amylase activity (AA), preharvest sprouting (PHS), heading and flowering time (HE), and different morphological traits.

Mapping population consisted of 143 genotypes of RIL-F~8~ generation. Consensus genetic map for the RIL-M (with consideration of the data from four rye populations) was created using the Multipoint Consensus 2.2 software package \[[@CR14]\]. Detailed information on RIL-M mapping population and algorithms used to release genetic map are provided by Milczarski et al. \[[@CR15]\].

All seven linkage groups accounted for 1318 markers (1256 DArTs, i.e. markers detected by Diversity Array Technology and 62 PCR-based loci). Individual chromosomes included from 117 (5R) to 257 (6R) loci and spanned the distance of 128 cM (5R)--251 cM (7R). The whole length of map was 1355 cM, and average distance varied from 0.7 cM to 1.9 cM.

Phenotype analyses {#Sec4}
------------------

All the experiments in this study were conducted at West Pomeranian University of Technology, Szczecin (53.45°N, 14.53°E). The RIL-M population (S7--S12 generations) consisted of 143 genotypes and were analyzed during six vegetation seasons (years: 2008, 2009, 2010, 2012, 2013, and 2017). In the years 2008--2010, genotypes of mapping population were planted and analyzed in duplicates. Each RIL was represented by 7--8 plants grown in the field or in pots under natural or near natural conditions.

Due to inbreeding depression of numerous lines, individuals representing each genotype were first germinated in a glasshouse (15th--25th September) and then (1st--20th October) vital seedlings were planted manually in the field---each line into one row. Finally, 7--8 individuals were grown in each row (length of rows was 100 cm, dimension between rows was 17.5 cm). The order of lines grown in the field was random and different in each year of study for the following replicates: RL08--1, RL09--1, RL09--2, RL10--1, RL10--2, RL12, and RL17. Variants RL09--1 and RL09--2 as well as RL10--1 and RL10--2 were grown in different plots.

Genotypes of two variants (replicates RL08--2 and RL13) were grown in the buckets of 10 dm^3^ volume filled with an equal mass (1440 g) of soil and sand mixture (1:1 *v*/v), seven plants of each genotype together in one bucket. Seeds were sown in January; initially plants were grown at a temperature of 15 °C for first 3 weeks, followed by gradually decreasing the temperature until it approached atmospheric conditions. Further vegetation proceeded under natural conditions of the winter--summer period (February--August), with natural daylight duration.

Our plant material was characterized by adaxial rolling leaves. A visual score of the degree of RL was made on the whole plants at the tillering stage, using a 5-point scale: 1---no signs of rolling, 1.5---the first evidence of rolling, 2---slightly rolled leaves, 2.5---strongly rolled leaves, and 3---completely rolled leaves (a closed cylinder). Figure [1](#Fig1){ref-type="fig"} shows the parental lines with extreme rate of traits. The assessment was performed after rainfall or watering the plants in buckets, to avoid the effect of leaves rolling under the influence of water deficit.Fig. 1Parental inbred lines of rye mapping population RIL-M, differing in terms of rolling leaves. **a** Paternal line S76 with rolled leaves, maternal line S120 with straight leaves, (**b**) line S76, and (**c**) line S120

Statistical analysis {#Sec5}
--------------------

QTL mapping was performed using the method of composite interval mapping (CIM) with Windows QTL Cartographer 2.5 software \[[@CR8]\]. The step size chosen for all QTLs was 2 cM. Significant thresholds to declare the presence of a QTL were estimated from 1000 permutations of the data at *P* \< 0.05. In addition, nonparametric Fisher's test was conducted to point markers significantly connected with rolling leaves segregations.

Variation parameters and correlation coefficient between replicates was established with STATISTICA 12.0 software package (Stat-Soft, Inc., USA, <http://www.statsoft.com>).

Results of CIM were compared with previous QTL analyses for AA, PHS, HE \[[@CR16]\], with recalculated data for consensus map), morphological traits \[[@CR17]\], spike compactness, and drought index of morphological traits \[[@CR18]\] for the same genetic map of RIL-M population.

Results {#Sec6}
=======

Phenotypic variation and correlation analyses {#Sec7}
---------------------------------------------

The mapping population RIL-M was characterized by an average coefficient of variation of 32.15--44.51 in terms of RL, depending on the replication (Table [1](#Tab1){ref-type="table"}). Mean values of the trait ranged from 1.74 to 2.11.Table 1Parameters of phenotypic variation and correlation coefficient (r) of leaf rolling (RL) in rye mapping population RIL-M, assessed in nine experimental replicationsExperiment replicationNumber of genotypes evaluated NRL parameter AvVarianceSDCoefficient of variationSESkewnessKurtosisrMinMaxRL08--11391.920.710.8443.900.070.15−1.580.36\*0.60\*RL08-2901.760.590.7743.800.080.45−1.170.250.60\*RL09-11421.830.610.7842.640.070.31−1.300.29\*0.47\*RL09-21401.740.510.7140.930.060.42−0.950.190.40\*RL10-11381.820.600.7842.660.070.33−1.270.260.60\*RL10-21311.750.610.7844.510.070.47−1.200.250.48\*RL121281.940.510.7136.780.060.19−1.140.38\*0.56\*RL131272.080.520.7234.590.06−0.12−1.170.290.46\*RL171352.110.460.6832.150.06−0.17−1.120.29\*0.60\*\*significant at *P* = 0.01

Significant positive correlations were found between results of the most pairs of replicates, and there was no correlation found in four cases (Additional file [1](#MOESM1){ref-type="media"}). Assessment of two variants in the same vegetative season showed a correlation coefficient of 0.60, 0.35, and 0.48 in the year 2008, 2009, and 2010, respectively. The weakest correlation (0.19) was noted between RL08--1 and RL09--2 and the strongest (0.60) correlation was found between RL08--2 and three other replications: RL08--1, RL10, RL17.

QTL analyses {#Sec8}
------------

Putative QTLs for RL were detected in each season except the year 2013. A total of 43 QTLs grouped in 28 intervals on all 7 chromosomes (2--8 *per* chromosome) were identified (Fig. [2](#Fig2){ref-type="fig"}). There were seven groups of coinciding QTLs found on chromosomes 3R, 4R, 5R, and 7R; four of them were detected 3--5 times (Fig. [2](#Fig2){ref-type="fig"}, Table [2](#Tab2){ref-type="table"}).Fig. 2QTLs for leaf rolling (RL) detected on genetic map of rye population RIL-M. On the left side of chromosomes---markers nearest to the QTL peak, underlined---markers pointed as the nearest the QTL peak at least twice. On the right side of chromosomes---white rectangle---QTL detected once out of nine replicates, gray rectangle---QTL detected at least twice, black rectangle---QTL for other traits mapped previously in this population. Abbreviations: QTLs, quantitative trait loci; TGW, thousand grain weight; PHS, preharvest sprouting; HE -- heading earliness; GNPS, grain number *per* spike; GNPS_DI, grain number *per* spike drought index; SCT, spike compactness; SNPP_DI, spike number per plant drought index; PH, plant height; AA, α-amylase activity; GW, grain weight per spikeTable 2QTLs for leaf rolling (RL) in rye mapping population RIL-M, detected at least twice, out of nine replicationsExperiment replicationQTL positionLODar^2^Marker nearest to the QTL peakChromosomePeak \[cM\]Interval \[cM\]RL10--13R9491--962.09−0.36.73XrPt507717RL173R9592--983.55−0.2411.4XrPt507717RL123R9894--1012.25−0.228.17XrPt507473RL08--13R9996--1004.48−0.3512.12XrPt507473RL08--23R9995--1013.12−0.3113XrPt507473RL09--14R9087--963.21−0.5111.73XrPt508114RL174R9191--923.61−0.3311.93XrPt507139RL08--24R164162--1643.58−0.5115.39XrPt507812RL09--24R164162--1652.59−0.387.75XrPt507812RL175R65--94.51−0.2411.91XrPt505767RL08--15R98--127.34−0.421XrPt505693RL10--15R127--122.73−0.238.09XrPt398579RL09--15R8381--845.38−0.3518.59Xpr665L1050RL10--15R8380--842.94−0.2910.06Xpr665L1050RL08--15R8482--844.01−0.319.85Xpr665L1050RL127R6964--762.58−0.217.95XrPt402607RL177R7674--843.46−0.329.17Xopo7RL09--27R213212--2153.67−0.269.41XrPt400276RL177R215210--2163.39−0.2510.18XrPt400276RL127R221210--2214.24−0.2713.52XrPt507004RL08--27R225215--2263.64−0.3517.84XrPt505931a -- additive effect of the maternal allele, r^2^ -- variance explained by a QTL

The QTLs detected were found to be responsible for 6.61--21.00% of the variation in the trait. Eighteen loci were characterized with coefficient of variation more than 10%. LOD values ranged from 2.09 to 7.34. Twenty-two QTLs peaks achieved LOD value above 3.0. The highest LOD was achieved by locus from 5RS. Absolute value of parameter *a* ranged from 0.21 to 0.51. Most of the alleles causing RL (31; 16 if consider common localization) were contributed by S76; however, 12 alleles responsible for RL originated from S120.

Thirty-five markers were indicated as the nearest to the QTL peaks: 31 DArTs, 2 random amplified polymorphic DNAs (RAPDs), and 2 sequence characterized amplified regions (SCARs); among them, three loci were revealed twice (XrPt507717, XrPt507812, and XrPt400276), and two loci were revealed three times (XrPt507473 and Xpr665L1050).

Figure [2](#Fig2){ref-type="fig"} shows the location of all QTLs. If any locus was found in different replicates, then we presented those only as a single gray rectangle.

In addition to CIM, the nonparametric, Fisher's (*F*) test was used to detect markers significantly linked to the phenotype of rolled leaves. A total of 155 loci from 6 chromosomes were pointed as linked to the trait at least once, most of them (46) from chromosome 7R. Chromosome 1R had no representation; 2R was represented by one marker. By limiting the pool of markers to those detected at least thrice, a set of 128 loci was obtained (Table [3](#Tab3){ref-type="table"}). Only 15 out of 35 loci nearest to the QTL peaks were in this group. Two strongly linked markers from 7R were revealed by *F* test in all 9 replications of the experiment. The other 26 markers were detected 8 times; among them, 1 marker from 5R, 3 from 3R, and the rest from 7R (Table [3](#Tab3){ref-type="table"}). Next, 26 loci were found to be significant in 7 replications, 19 markers were pointed 6 times, 40 markers were pointed 5 times, 7 markers were pointed 4 times, 8 markers were pointed 3 times, and 10 markers were pointed twice. All intervals containing markers detected in *F* test included QTL designated by CIM method.Table 3Markers linked to leaf rolling (RL), significant in Fisher's test at *P* = 0.01, at least in three replicatesMarker12345Marker12345Marker12345XrPt4107637R22193229XrPt4023775R596257XrPt3905855R12143115XrPt4108307R22296229XrPt4024057R24066187XrPt3988005R12170205XrPt3461847R22497208XrPt4108183R9275197[XrPt400297]{.ul}5R9850135XrPt3892587R22296248**XrPt505767**5R6103277XrPt4008483R8460155**[XrPt390593]{.ul}**7R240116248XrPt5061473R8972177XrPt4012283R8452145XrPt3907507R251110218XrPt5062607R24069177[XrPt401454]{.ul}5R1769175**XrPt398579**5R12103298XrPt5071805R12112347XrPt4014813R10843135XrPt3997757R23368148XrPt5075005R5103277XrPt4022455R1676205XrPt4002527R22490178**[XrPt507717]{.ul}**3R9482187[XrPt402531]{.ul}5R12547165XrPt4003193R8963138XrPt5079485R083237**XrPt410924**5R1592225XrPt4007937R23368148XrPt5080617R7268147XrPt4111095R11357165XrPt4012007R227114198**Xpr665L1050**5R8485226XrPt4112525R1776205XrPt4013727R23375158XrPt3469083R7856146XrPt4113205R11956165**[XrPt401480]{.ul}**7R22794168XrPt3475747R6571166XrPt5052197R6547175XrPt4018287R247121248XrPt3479985R1492226**XrPt505693**5R889295XrPt4022627R247109218[XrPt398502]{.ul}3R10069176XrPt5057215R11357165XrPt4108847R23368148XrPt3985197R6571166XrPt5060015R1776205XrPt5052153R9167138[XrPt401081]{.ul}3R10069176XrPt5068215R10248135XrPt5055237R23368148XrPt4017957R6571166[XrPt506905]{.ul}3R7951145XrPt5058647R22497208XrPt4025893R7657156XrPt5073693R8452145**XrPt505931**7R22497208**[XrPt402607]{.ul}**7R6874176XrPt5079265R9850135XrPt5064947R22296208[XrPt410783]{.ul}5R1292226XrPt5079535R12272205XrPt5066073R9167138XrPt4115225R12157146XrPt5089055R9850135XrPt5067647R22495198XrPt5055933R8977196XrPt5089255R782245XrPt5077547R22296218XrPt5068743R10348146**XrPt509159**3R7962175XrPt5079367R22291218**XrPt507004**7R22081186Xscm1415R7953125XrPt5085597R7286168XrPt5074623R7657156Xscsz877L9505R12159165XrPt5088377R22296208**[XrPt507473]{.ul}**3R9995256XrPt3465835R10836114**Xopo7**7R7682187[XrPt508197]{.ul}5R598276XrPt3907417R20848144Xpr57L4705R10168157XrPt5096473R7856146XrPt3996547R20944124XrPt1172523R8972177Amy3.25R12850135XrPt4110205R12155154XrPt1197185R8451117Xpr665L4303R8141105XrPt4111845R1859164XrPt1209903R8972177XrPt3467555R1776205XrPt5054377R20850134XrPt3454393R8972177XrPt3468925R1776205XrPt5090515R10643124XrPt3480933R8972177XrPt3469465R1776205XrPt3467796R5935163XrPt3892615R12108327XrPt3469805R1776205XrPt3471146R6333163XrPt3895857R21681177XrPt3470725R1776205XrPt3897117R20738133XrPt3898955R5103277XrPt3472125R1679215XrPt3904426R5935163XrPt3899597R21096237XrPt3474545R1776205XrPt4005097R20734113XrPt3903625R696257XrPt3478095R1776205XrPt4007327R20238133XrPt3986273R8972177XrPt3894543R8153145XrPt4015237R20748123**XrPt400276**7R21496237XrPt3897595R10738125XrPt5091767R20748123XrPt4017547R7355127XrPt3898153R8153145Bolded---markers nearest to the QTL peak, underlined---markers with predicted function (Table [4](#Tab4){ref-type="table"}), 1---chromosome, 2---position \[cM\], 3---sum of *F* statistic from all replicates, 4---maximal value of *F* statistic, 5---number of replicates, in which marker was pointed as significant

Markers of RL putative homologs {#Sec9}
-------------------------------

The set of DArT sequences \[[@CR19]\] was screened to find sequences of markers linked to RL. A total of 67 found sequences were directed to BLAST sequence analysis in NCBI nucleotide collection database (<https://www.ncbi.nlm.nih.gov/>), using megaBLAST algorithm. Hits exceeding a score of 200 or identity of 95% and those with known (predicted) identity were chosen. Finally, 12 records met the set criteria (Table [4](#Tab4){ref-type="table"}); one of them was matched to *Secale cereale* cds, one to *Triticum aestivum* genomic sequence, and the rest to *Aegilops tauschii* mRNA.Table 4DArTs linked to leaf rolling (RL) with known sequences, annotated in NCBIMarkerNCBI annotationMarkerDArT sequence length (bp)Chrom.DescriptionSpeciesMax scoreTotal scoreQuery coverE valueIdentityAccessionXrPt5069056653Rsubtilisin-like protease SBT2.2 (LOC109761160), mRNA*Aegilops tauschii subsp. tauschii*97978%6.00E-1695%XM_020319955.1**XrPt507717**6783Rprotein LIKE COV 2 (LOC109734806), mRNA*Aegilops tauschii* subsp. *tauschii*27927924%5.00E-7196%XM_020293989.1**XrPt507473**5993RFhb1 region genomic sequence*Triticum aestivum cv Sumai 3*60605%7.00E-0597%KX907434.1XrPt401081 XrPt3985025453Rjasmonate O-methyltransferase-like (LOC109744228), mRNA*Aegilops tauschii* subsp. *tauschii*35969286%5.00E-9595%XM_020303320.1XrPt5081977205RRGA1-G gene, complete cds*Secale cereale*22035819%4.00E-5396%KT725818.1XrPt4107834635Rcancer-related nucleoside-triphosphatase (LOC109750338), mRNA*Aegilops tauschii* subsp. *tauschii*52762475%1.00E-14594%XM_020309304.1XrPt4014542185RLRR receptor-like serine/threonine-protein kinase FEI 1 (LOC109774869), mRNA*Aegilops tauschii* subsp. *tauschii*575713%3.00E-04100%XM_020333663.1XrPt4002977735Rserine/threonine-protein kinase AFC1-like (LOC109779169), transcript variant X4, mRNA*Aegilops tauschii* subsp. *tauschii*22828219%2.00E-55100%XM_020337784.1XrPt4025314225Rtranscription factor bHLH79-like (LOC109737374), mRNA*Aegilops tauschii* subsp. *tauschii*999913%1.00E-1698%XM_020296517.1**XrPt402607**5517Rpolyadenylate-binding protein-interacting protein 7-like (LOC109738116), transcript variant X4, mRNA*Aegilops tauschii* subsp. *tauschii*47047060%2.00E-12893%XM_020297211.1**XrPt401480**5617Rtubby-like F-box protein 12 (LOC109763585), transcript variant X3, mRNA*Aegilops tauschii* subsp. *tauschii*60160165%8.00E-16896%XM_020322436.1XrPt3905934517Rvegetative cell wall protein gp1-like (LOC109748681), transcript variant X4, mRNA*Aegilops tauschii* subsp. *tauschii*57576%6.00E-04100%XM_020307705.1Bolded---markers nearest to the QTL peak

Co-localization of QTLs {#Sec10}
-----------------------

Some QTL intervals for RL overlapped partially or completely with QTLs for other agronomic traits analyzed previously in this population. There were nine QTLs overlapping most precisely, which means that they had the same markers nearest to the QTL peak and nearly the same intervals (Fig. [2](#Fig2){ref-type="fig"}, Table [5](#Tab5){ref-type="table"}).Table 5QTLs for leaf rolling (RL) in rye mapping population RIL-M, coinciding with QTLs for other traits, mapped previously in this populationChrom.Coinciding QTLsMarker nearest to the QTL peakQTL peak \[cM\]QTL interval \[cM\]LODar^2^1RRL08--1Xpr57L19009896--1012.70.338.61TGW099898--1023.03−1.949.421RRL12XrPt507993181180--1822.80.248.57PHS08181178--1833.02−6.058.372RHE08Xrpt507900202200--2063.330.558.74RL08--1204202--2083.1−0.319.634RGNPS08XrPt5081148987--913.056.7710.8PHS108985--1062.69−8.649.2RL09--19087--963.21− 0.5111.734RGNPS08XrPt401039207206--2084.47−3.9412.54RL10--2207205--2083.36−0.3610.77GNPS_DI207205--2102.7312.7615.455RSCT09Xpr665L10508173--842.8−0.112.02RL09--18181--845.38−0.3518.59RL10--18380--842.94−0.2910.06RL-08-18482--844.01−0.319.857RRL17Xopo77674--843.46−0.329.17PHS077769--872.7413.7618.49HE098474--952.30.8316.977RSNPP_DIXrPt507004219214--2281.85−8.389.9PH08220216--2222.764.158.31RL12221210--2214.24−0.2713.527RAA09XrPt401480227226--2388.251.2122.85RL10--1227227--2283.68−0.2811.51PHS08N228227--2385.179.8916.64GW09229224--2292.11−0.096.62a -- additive effect of the maternal allele, r^2^ -- variance explained by a QTL*TGW* thousand grain weight, *PHS* preharvest sprouting, *HE* heading earliness, *GNPS* grain number *per* spike, *GNPS_DI* grain number *per* spike drought index, *SCT* spike compactness, *SNPP_DI* spike number *per* plant drought index, *PH* plant hight, *AA* α-amylase activity, *GW* grain weight *per* spike

QTLs for RL were associated with QTLs for morphological traits such as grain number and weight, spike number per plant, compactness of spike, and plant height. Two QTLs of this set were loci controlling drought index (GNPS_DI and SNPP_DI), that is, those that caused differences in trait expression in normal and drought conditions. QTL for RL were also co-localized with two QTLs for heading earliness (2R, 7R), one with α-amylase activity QTL (7R) and three for pre-harvest sprouting QTL (1R, 4R, and 7R).

Distal part of chromosome 7RL is a special region with numerous overlapping QTLs: five controlling RL and five engaged in the expression of other traits (Fig. [2](#Fig2){ref-type="fig"}, Table [5](#Tab5){ref-type="table"}).

Discussion {#Sec11}
==========

RL is a phenomenon commonly found in grasses and quite commonly described for field crops, but it has not been studied in rye. As leaf morphology is an important agronomic trait in the breeding of rice \[[@CR20]\], the mechanism of RL was of interest, especially to rice researchers. There is no information about the genetic base of molecular mechanisms involved in leaf shape in rye, except one report regarding QTLs controlling leaf area \[[@CR21]\].

A study similar to ours has been conducted on tetraploid wheat \[[@CR4]\]; however, we did not study RL as the reaction to water deficit. The research on drought resistance in durum wheat × wild emmer wheat recombinant inbred line population allowed to detect 11 significant QTLs associated with flag RL, mapped on chromosomes: 1A, 2A, 2B, 4B, 5A, 5B, 6A, 6B, 7A, and 7B \[[@CR4]\]. Three of these QTLs were found to be environment responsive.

We detected 43 QTLs, grouped into 28 intervals, which confirms the multigenic base of the trait stated for wheat and rice. Studies on rice showed that no fewer than 70 genes/QTLs for RL have been mapped or cloned till now \[[@CR11]\]. Four QTLs for RL were stable in different environments. Due to many QTLs for RL were detected once, their genotype×environment interaction (GEI) could be inferred. GEI is a common characteristic for quantitative traits. For possible breeding purposes (like marker Assisted Selection - MAS), QTLs that are more environment-specific should be treated with caution and more attention should be paid to the repetitive QTLs. However, QTLs of varying manifestation, dependent on the environmental influence are also interesting for cognitive purposes.

In this study, we focused on the analysis of RL per se, without linking this trait with the response to drought. However, co-localization of two QTLs for RL with QTLs for drought index of grain number per spike (4R) and spike number per plant (7R) revealed in other experiments \[[@CR18]\] indicates the relationship of detected loci with adaptive mechanisms to drought-related stress conditions.

Furthermore, 5 out of 11 QTLs for RL mapped on tetraploid wheat were co-localized with QTLs associated with plant productivity \[[@CR4]\]. RL was also found to be associated with plant height in four regions (2B, 4B, 6A, and 7A) and with heading earliness (days from planting to heading) in two intervals on chromosomes 4B and 5A \[[@CR4]\]. We also found nine QTLs co-localized with QTLs for other agronomic traits mapped in the same population \[[@CR17]\], such as grain number and weight, spike number per plant, compactness of spike and plant height. QTLs for RL were also found to be co-localized with two QTLs for heading earliness (2R and 7R), one with α-amylase activity QTL (7R), and three with preharvest sprouting QTLs (1R, 4R, and 7R) \[[@CR16]\].

Additional confirmation of the association between loci engaged in controlling RL and QTLs responsible for different agronomic traits is the result of comparing markers linked with RL and markers for nine features studied in the other, unrelated rye RILs' mapping population 541 × Ot1--3 \[[@CR21]\]. There were loci for plant height, stem thickness, spike length, awn length, heading date, thousand grain weight, grain length, leaf area, and chlorophyll content localized on the DArT-based high-density map of this population.

Each of the nine traits was characterized by some markers common with these, linked to RL in our population. There were 45 such markers distributed throughout the five chromosomes (3R, 4R, 5R, 6R, and 7R). Among them three DArTs from 7R were common for RL and leaf size (XrPt505931, XrPt389959, and XrPt400276). Seven DArTs were linked to RL and chlorophyll content: three from 5R (XrPt389759, XrPt346583, and XrPt505721) and four from 7R chromosome (XrPt402607, XrPt398519, XrPt347574, and XrPt401795). All these relationships suggest very strong linkages and/or pleiotropic effects of many genes, which remains in agreement with previous results for rye \[[@CR16], [@CR17], [@CR21]\], wheat \[[@CR4]\], and rice \[[@CR22]\].

All DArTs significantly linked to RL were subjected to screening the DArT sequences database, followed by NCBI database blasting, in order to find the homologs. A total of 12 records with a known identity were found; majority of rye DArTs sequences were most similar to *Aegilops tauschii* mRNAs. Only one of them was matched to *Secale cereale* cds of resistant gene analog (putative disease resistance gene), and one to *Triticum aestivum* genomic sequence---also connected with resistance, in this case to fusarium head blight (FHB).

Sequence of DArT XrPt506905 from 3R showed similarity to predicted gene, namely, subtilisin-like protease. Subtilisin-like proteases (subtilases) are serine proteases and constitute the largest group of peptidases. Although several subtilases have been identified in plants (e.g., about 60 subtilase genes are known in *Oryza sativa* and *Arabidopsis thaliana*), most of their functions in plants remain unknown \[[@CR23]\] (and bibliography therein). It is likely that subtilases contribute significantly to the developmental processes and signaling cascades in plants (Rautengarten et al. 2005, after \[[@CR23]\]). For instance, the loss-of-function mutation in *ALE1* leads to abnormal leaf shape \[[@CR24]\]. Marker XrPt506905, which is a predicted gene for subtilisin-like protein, in addition to the linkage with RL also showed a relationship with awn and grain length \[[@CR21]\].

DArT XrPt507717, nearest to the peak of a QTL for RL from 3R, revealed similarity to protein-like COV2 mRNA. The role of COV2, inferred from the sequence or structural similarity to COV1 is stem vascular tissue pattern formation. COV1 is predicted to be an integral membrane protein that may be involved in the perception or transport of a signaling molecule that negatively regulates the differentiation of vascular tissue in the developing stem of *Arabidopsis* \[[@CR25]\]. Marker XrPt507717, in addition to the linkage with RL also showed a relationship with awn length \[[@CR21]\].

Two DArTs, XrPt401081 and XrPt398502 from 3R, has sequences homologous to jasmonate *O*-methyltransferase. This enzyme catalyzes the methylation of jasmonate into methyljasmonate, a plant volatile that acts as an important cellular regulator mediating diverse developmental processes and defense responses (<http://www.uniprot.org/uniprot/Q9AR07>). It is involved in the pathway of oxylipin biosynthesis, which is a part of lipid metabolism. To this end, 28 differentially expressed proteins related to rolled leaf traits were isolated and identified. Some of the proteins and genes detected are involved in lipid metabolism, which is related to the development of bulliform cells, such as phosphoinositide phospholipase C, Mgll, and At4g26790 \[[@CR12]\].

Sequence of next DArT from 3R, XrPt507473, proved to be similar to *Fhb1*, a major FHB-resistant gene. *Fhb1* was fine mapped on the distal segment of chromosome 3BS of spring wheat (*Triticum aestivum* L.). One of the recent studies has reported that wheat *Fhb1* encodes a chimeric lectin with agglutinin domains and a pore-forming toxin-like domain conferring resistance to FHB \[[@CR26]\]. Marker XrPt507473, in addition to the linkage with RL, also showed a relationship with awn length and stem thickness \[[@CR21]\].

XrPt401454 from 5R showed homology to LRR receptor-like serine/threonine-protein kinase FEI11 mRNA. FEI1 is involved in the signaling pathway that regulates cell wall function, including cellulose biosynthesis, likely via an 1-aminocyclopropane-1-carboxylic acid (ACC)-mediated signal; a precursor of ethylene (<http://www.uniprot.org/uniprot/C0LGF4>). To date, 13 genes associated with rice RL have been isolated or cloned. The cytological mechanism of RL has been found to be largely related to the abnormal development of bulliform cells. NRL1 encodes cellulose synthase and plays a positive role in the regulation of bulliform cell development. In mutant rice plants that lack this gene, shrinkage is found in the area of the bulliform cells, thereby causing inward rolling of rice leaves \[[@CR12]\].

Sequence of XrPt402531 from 5R demonstrated high similarity to predicted *A. tauschii* transcription factor bHLH79 (basic helix-loop-helix protein 79). The function of this factor is unknown; however, several transcription factors are known to be engaged in the establishment of abaxial/adaxial leaf polarity. For example, mutation in SLL1/RL9, a member of the KANADI family, encoding a transcription factor \[[@CR27]\], leads to the failure of programmed cell death of abaxial mesophyll cells and the suppression of the differentiation of the abaxial cells, and finally to generate adaxially rolled leaves. ROC5 encodes a protein containing a leucine zipper domain, homologous to GLABRA2 in *Arabidopsis* which results in the negative regulation in the development of the bulliform cells. The number and size of the bulliform cells increased when ROC5 was knocked out, thereby leading to the generation of adaxially rolled leaves, whereas co-suppression of ROC5 resulted in abaxial RL \[[@CR28]\].

Overexpression of a rice gene OsLBD3--7 that encodes a LBD family transcription factor promoted narrow and adaxially rolled leaves by decreasing the size and number of bulliform cells. OsLBD3--7 also upregulated the expression of negative regulators of bulliform cells, which implies that OsLBD3--7 acts as a suppressor of bulliform cell development \[[@CR23]\].

The other rice gene ACL1 encodes an unknown protein with a conserved functional domain; OsZHD1 encodes a domain transcription factor with homologous zinc finger structure. These genes also play a positive role in the regulation of bulliform cell development, and overexpression of these two genes results in an increased number of bulliform cells, thereby causing outward rolling of rice leaves \[[@CR12], [@CR29]\].

XrPt402607 from 7R, marker nearest to the peak of QTL for RL and also marker linked with grain length and chlorophyll content \[[@CR21]\], showed homology to predicted polyadenylate-binding protein-interacting protein 7-like gene. The poly(A) binding proteins (PABP) play an important role in the regulation of translation; however, the role of this particular factor is unknown.

Next marker from 7R, the peak of QTL for RL, DArT XrPt401480, seems to be homologous to the predicted gene encoding tubby-like F-box protein 12. Plants include a large number of tubby-like proteins (TLPs/TULPs). For example, there are 11 members of the tubby gene family in *Arabidopsis* \[[@CR30]\], 14 in rice \[[@CR31]\], 11 in poplar \[[@CR32]\], 4 in wheat \[[@CR33]\], and 8 in sorghum (<http://www.ncbi.nlm.nih.gov>). The existence of multiple TLPs implies their vital function in plants. F-box proteins regulate diverse cellular processes, including cell cycle transition, transcriptional regulation, and signal transduction. Lai et al. \[[@CR30]\] have demonstrated that AtTLP9 interacts with ASK1 (*Arabidopsis* Skp1-like 1). According to them, F-box domain containing plant TLPs acting as transcription regulators should have cellular function activities of F-box proteins in signal transduction. AtTLP9 might participate in the abscisic acid signaling pathway \[[@CR30]\]. The other example of F-box protein regulating plant growth and development include TIR1 acting in response to auxin \[[@CR34]\].

The function of tubby-like F-box protein 12 encoded by genes of *Aegilops tauschii* and *Brachypodium distachyon*, homologous to XrPt401480, is unknown. However, some rice TLPs, especially OsTLP12, were probably involved in the abscisic acid and gibberellin signaling processes. This role might also be attributed to rye TLP12, because the same DArT was pointed as a marker linked to plant height \[[@CR21]\].

XrPt390593 sequence from 7R was similar to vegetative cell wall protein gp1-like mRNA. The nature of cell wall proteins is as varied as the many functions of plant cell walls. Majority of the cell wall proteins are cross-linked into the cell wall and probably have structural functions. If this protein was associated with bulliform and/or hypodermis cells it might have an effect on RL, because these two types of cells are involved in RL in higher plants \[[@CR1]\].

Although the roles described for the aforementioned markers as the potential genes that control RL are likely, their functions and association with RL and other traits should be verified in expression tests and will be studied during further research.

Conclusion {#Sec12}
==========

The first localization of loci controlling leaf rolling (RL) on high density consensus genetic map of rye confirms the multigenic base of the trait stated for wheat and rice. Four stable quantitative QTLs on chromosomes 3R, 5R, and 7R were identified. Coinciding QTLs for RL and agronomic traits (e.g. drought tolerance) suggested pleiotropic effects of genes engaged in leaf rolling control. Four reproducible QTLs may be interesting for breeding purposes. The composite QTL analyze, the knowledge of the sequence of markers strongly linked to RL and the comparison of their homology with sequences data of related species allowed to indicate rye candidate genes controlling RL process.

Additional files
================

 {#Sec13}

Additional file 1:Matrix of coefficient correlation of leaf rolling (RL) in nine experimental variants, assessed in rye mapping population RIL-M. (XLSX 11 kb) Additional file 2:Phenotypic data of rye mapping population RIL-M. (XLSX 66 kb)

AA

:   α-amylase activity

bHLH

:   Basic helix-loop-helix protein

DArT

:   Diversity array technology

GNPS

:   Grain number *per* spike

GNPS_DI

:   Grain number *per* spike drought index

GW

:   Grain weight per spike

HE

:   Heading earliness

MAS

:   Marker assisted selection

PH

:   Plant height

PHS

:   Pre-harvest sprouting

QTLs

:   Quantitative trait loci

RIL

:   Recombinant inbred line

RL

:   Leaf rolling

SCT

:   Spike compactness

SNPP_DI

:   Spike number per plant drought index

TGW

:   Thousand grain weight

TLP

:   Tubby-like proteins

Thanks to Elżbieta Pauker, Henryk Pauker and Jarosław Myśków for technical assistance.

Funding {#FPar1}
=======

The study was financially supported by Ministry of Science and Higher Education, Republic of Poland under grant No 518070173171--01/18. The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials {#FPar2}
==================================

The raw mapping data are available in the article by Milczarski et al. \[[@CR15]\] and the phenotypic data are submitted as the supplementary material in Additional file [2](#MOESM2){ref-type="media"}.

Conceived and designed the experiments: BM. Performed the experiments: BM, NL. Analyzed the data: BM, MG, NL, IC-M. Wrote the paper: BM. Critically revised the paper: SS. Proportions of authors' contribution: 6:1:1:1:1 (BM, MG, NL, IC-M, SS). All authors read and approved the final version of the manuscript.

This study has not directly involved humans and animals. The plant material belonged to the West-Pomeranian University of Technology, Szczecin; the mapping population was bred by the authors. No special permissions and/or licences were necessary.

Not applicable

The authors declare that they have no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
