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.

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.

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.

MiR-93/miR-375: Diagnostic Potential, Aggressiveness Correlation and Common Target Genes in Prostate Cancer.

Dysregulation of miRNAs has a fundamental role in the initiation, development and progression of prostate cancer (PCa). The potential of miRNA in gene therapy and diagnostic applications is well documented. To further improve miRNAs' ability to distinguish between PCa and benign prostatic hyperplasia (BPH) patients, nine miRNA (-21, -27b, -93, -141, -205, -221, -182, -375 and let-7a) with the highest reported differentiation power were chosen and for the first time used in comparative studies of serum and prostate tissue samples. Spearman correlations and response operating characteristic (ROC) analyses were applied to assess the capability of the miRNAs present in serum to discriminate between PCa and BPH patients. The present study clearly demonstrates that miR-93 and miR-375 could be taken into consideration as single blood-based non-invasive molecules to distinguish PCa from BPH patients. We indicate that these two miRNAs have six common, PCa-related, target genes (CCND2, MAP3K2, MXI1, PAFAH1B1, YOD1, ZFYVE26) that share the molecular function of protein binding (GO:0005515 term). A high diagnostic value of the new serum derived miR-182 (AUC = 0.881, 95% confidence interval, CI = 0.816-0.946, p < 0.0001, sensitivity and specificity were 85% and 79%, respectively) is also described.

Local and systemic regulation of PSII efficiency in triticale infected by the hemibiotrophic pathogen Microdochium nivale.

Microdochium nivale is a fungal pathogen that causes yield losses of cereals during winter. Cold hardening under light conditions induces genotype-dependent resistance of a plant to infection. We aim to show how photosystem II (PSII) regulation contributes to plant resistance. Using mapping population of triticale doubled haploid lines, three M. nivale strains and different infection assays, we demonstrate that plants that maintain a higher maximum quantum efficiency of PSII show less leaf damage upon infection. The fungus can establish necrotrophic or biotrophic interactions with susceptible or resistant genotypes, respectively. It is suggested that local inhibition of photosynthesis during the infection of sensitive genotypes is not balanced by a supply of energy from the tissue surrounding the infected cells as efficiently as in resistant genotypes. Thus, defence is limited, which in turn results in extensive necrotic damage. Quantitative trait loci regions, involved in the control of both PSII functioning and resistance, were located on chromosomes 4 and 6, similar to a wide range of PSII- and resistance-related genes. A meta-analysis of microarray experiments showed that the expression of genes involved in the repair and de novo assembly of PSII was maintained at a stable level. However, to establish a favourable energy balance for defence, genes encoding PSII proteins resistant to oxidative degradation were downregulated to compensate for the upregulation of defence-related pathways. Finally, we demonstrate that the structural and functional integrity of the plant is a factor required to meet the energy demand of infected cells, photosynthesis-dependent systemic signalling and defence responses.

Changes in fenugreek transcriptome induced by methyl jasmonate and steroid precursors revealed by RNA-Seq.

Changes in fenugreek transcriptome related to enhanced production of steroids were induced by methyl jasmonate, cholesterol and squalene, and recorded using RNA-seq. A total of 112,850 unigenes were obtained after de novo assembling of next generation sequencing data, and used for functional annotations. In steroidal saponins pathway, transcripts involved in mevalonate, terpenoid backbone and plant sterol synthesis were annotated. Overexpression of several transcripts from phytosterol biosynthesis pathway was confirmed by quantitative RT-PCR. In diosgenin biosynthesis pathway, fatty acid ω-hydroxylase (CYP86A2) and steroid 22-alpha-hydroxylase (CYP90B1) genes were annotated in all induced transcriptomes. Moreover, direct sequencing confirmed increased levels of CYP90B1, unspecific monooxygenase and 26-hydroxylase genes in plants with elevated level of diosgenin. New unigenes corresponding to enzymes involved in biosynthesis of diosgenin from cycloartenol via cholesterol were obtained and the role of CYP72A family in steroidal saponin biosynthesis was proposed. Additional support for biosynthetic pathway from cycloartenol to diosgenin was provided.

QTLs for cell wall-bound phenolics in relation to the photosynthetic apparatus activity and leaf water status under drought stress at different growth stages of triticale.

The present study aimed at identifying the regions of triticale genome responsible for cell wall saturation with phenolic compounds under drought stress during vegetative and generative growth. Moreover, the loci determining the activity of the photosynthetic apparatus, leaf water content (LWC) and osmotic potential (Ψ o) were identified, as leaf hydration and functioning of the photosynthetic apparatus under drought are associated with the content of cell wall-bound phenolics (CWPh). Compared with LWC and Ψ o, CWPh fluctuations were more strongly associated with changes in chlorophyll fluorescence. At the vegetative stage, CWPh fluctuations were due to the activity of three loci, of which only QCWPh.4B was also related to changes in F v/F m and ABS/CSm. In the other QTLs (QCWPh.6R.2 and QCWPh.6R.3), the genes of these loci determined also the changes in majority of chlorophyll fluorescence parameters. At the generative stage, the changes in CWPh in loci QCWPh.4B, QCWPh.3R and QCWPh.6R.1 corresponded to those in DIo/CSm. The locus QCWPh.6R.3, active at V stage, controlled majority of chlorophyll fluorescence parameters. This is the first study on mapping quantitative traits in triticale plants exposed to drought at different stages of development, and the first to present the loci for cell wall-bound phenolics.

Next-generation sequencing of representational difference analysis products for identification of genes involved in diosgenin biosynthesis in fenugreek (Trigonella foenum-graecum).

MAIN CONCLUSION: Representational difference analysis of cDNA was performed and differential products were sequenced and annotated. Candidate genes involved in biosynthesis of diosgenin in fenugreek were identified. Detailed mechanism of diosgenin synthesis was proposed. Fenugreek (Trigonella foenum-graecum L.) is a valuable medicinal and crop plant. It belongs to Fabaceae family and has a unique potential to synthesize valuable steroidal saponins, e.g., diosgenin. Elicitation (methyl jasmonate) and precursor feeding (cholesterol and squalene) were used to enhance the content of sterols and steroidal sapogenins in in vitro grown plants for representational difference analysis of cDNA (cDNA-RDA). To identify candidate genes involved in diosgenin biosynthesis, differential, factor-specific libraries were subject to the next-generation sequencing. Approximately 9.9 million reads were obtained, trimmed, and assembled into 31,491 unigenes with an average length of 291 bp. Then, functional annotation and gene ontogeny enrichment analysis was performed by aligning all-unigenes with public databases. Within the transcripts related to sterol and steroidal saponin biosynthesis, we discovered novel candidate genes of diosgenin biosynthesis and validated their expression using quantitative RT-PCR analysis. Based on these findings, we supported the idea that diosgenin is biosynthesized from cycloartenol via cholesterol. This is the first report on the next-generation sequencing of cDNA-RDA products. Analysis of the transcriptomes enriched in low copy sequences contributed substantially to our understanding of the biochemical pathways of steroid synthesis in fenugreek.

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

KEY MESSAGE: 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.

Genome-wide association mapping in elite winter wheat breeding for yield improvement.

Increased grain yield (GY) is the primary breeding target of wheat breeders. We performed the genome-wide association study (GWAS) on 168 elite winter wheat lines from an ongoing breeding program to identify the main determinants of grain yield. Sequencing of Diversity Array Technology fragments (DArTseq) resulted in 19,350 single-nucleotide polymorphism (SNP) and presence-absence variation (PAV) markers. We identified 15 main genomic regions located in ten wheat chromosomes (1B, 2B, 2D, 3A, 3D, 5A, 5B, 6A, 6B, and 7B) that explained from 7.9 to 20.3% of the variation in grain yield and 13.3% of the yield stability. Loci identified in the reduced genepool are important for wheat improvement using marker-assisted selection. We found marker-trait associations between three genes involved in starch biosynthesis and grain yield. Two starch synthase genes (TraesCS2B03G1238800 and TraesCS2D03G1048800) and a sucrose synthase gene (TraesCS3D03G0024300) were found in regions of QGy.rut-2B.2, QGy.rut-2D.1, and QGy.rut-3D, respectively. These loci and other significantly associated SNP markers found in this study can be used for pyramiding favorable alleles in high-yielding varieties or to improve the accuracy of prediction in genomic selection.

Genetic mapping of adult-plant resistance genes to powdery mildew in triticale.

Triticale is a cereal of high economic importance; however, along with the increase in the area of this cereal, it is more often infected by the fungal pathogen Blumeria graminis, which causes powdery mildew. The rapid development of molecular biology techniques, in particular methods based on molecular markers may be an important tool used in modern plant breeding. Development of genetic maps, location of the QTLs defining the region of the genome associated with resistance and selection of markers linked to particular trait can be used to select resistant genotypes as well as to pyramidize several resistance genes in one variety. In this paper, we present a new, high-density genetic map of triticale doubled haploids (DH) population "Grenado" × "Zorro" composed of DArT, silicoDArT, and SNP markers. Composite interval mapping method was used to detect eight QTL regions associated with the area under disease progress curve (AUDPC) and 15 regions with the average value of powdery mildew infection (avPM) based on observation conducted in 3-year period in three different locations across the Poland. Two regions on rye chromosome 4R, and single loci on 5R and 6R were reported for the first time as regions associated with powdery mildew resistance. Among all QTLs, 14 candidate genes were identified coded cyclin-dependent kinase, serine/threonine-protein kinase-like protein as well as AMEIOTIC 1 homolog DYAD-like protein, DETOXIFICATION 16-like protein, and putative disease resistance protein RGA3. Three of identified candidate genes were found among newly described QTL regions associated with powdery mildew resistance in triticale.

Representational Difference Analysis of Transcripts Involved in Jervine Biosynthesis.

Veratrum-type steroidal alkaloids (VSA) are the major bioactive ingredients that strongly determine the pharmacological activities of Veratrum nigrum. Biosynthesis of VSA at the molecular and genetic levels is not well understood. Next-generation sequencing of representational difference analysis (RDA) products after elicitation and precursor feeding was applied to identify candidate genes involved in VSA biosynthesis. A total of 12,048 contigs with a median length of 280 bases were received in three RDA libraries obtained after application of methyl jasmonate, squalene and cholesterol. The comparative analysis of annotated sequences was effective in identifying candidate genes. GABAT2 transaminase and hydroxylases active at C-22, C-26, C-11, and C-16 positions in late stages of jervine biosynthesis were selected. Moreover, genes coding pyrroline-5-carboxylate reductase and enzymes from the short-chain dehydrogenases/reductases family (SDR) associated with the reduction reactions of the VSA biosynthesis process were proposed. The data collected contribute to better understanding of jervine biosynthesis and may accelerate implementation of biotechnological methods of VSA biosynthesis.

Identification of candidate genes involved in steroidal alkaloids biosynthesis in organ-specific transcriptomes of Veratrum nigrum L.

Veratrum nigrum is protected plant of Melanthiaceae family, able to synthetize unique steroidal alkaloids important for pharmacy. Transcriptomes from leaves, stems and rhizomes of in vitro maintained V. nigrum plants were sequenced and annotated for genes and markers discovery. Sequencing of samples derived from the different organs resulted in a total of 108,511 contigs with a mean length of 596 bp. Transcripts derived from leaf and stalk were annotated at 28%, and 38% in Nr nucleotide database, respectively. The sequencing revealed 949 unigenes related with lipid metabolism, including 73 transcripts involved in steroids and genus-specific steroid alkaloids biosynthesis. Additionally, 3203 candidate SSRs markers we identified in unigenes with average density of one SSR locus every 6.2 kb sequence. Unraveling of biochemical machinery of the pathway responsible for steroidal alkaloids will open possibility to design and optimize biotechnological process. The transcriptomic data provide valuable resources for biochemical, molecular genetics, comparative transcriptomics, functional genomics, ecological and evolutionary studies of V. nigrum.

The effects of cold acclimation on photosynthetic apparatus and the expression of COR14b in four genotypes of barley (Hordeum vulgare) contrasting in their tolerance to freezing and high-light treatment in cold conditions.

BACKGROUND AND AIMS: Cold acclimation modifies the balance of the energy absorbed and metabolized in the dark processes of photosynthesis, which may affect the expression of cold-regulated (COR) genes. At the same time, a gradual acclimation to the relatively high light conditions is observed, thereby minimizing the potential for photo-oxidative damage. As a result, the resistance to photoinhibition in the cold has often been identified as a trait closely related to freezing tolerance. Using four barley genotypes that differentially express both traits, the effect of cold acclimation on freezing tolerance and high-light tolerance was studied together with the expression of COR14b, one of the best-characterized barley COR genes. METHODS: Plants were cold acclimated for 2 weeks at 2 degrees C. Freezing tolerance was studied by means of electrolyte leakage. Changes in photosynthetic apparatus and high-light tolerance were monitored by means of chlorophyll fluorescence. Accumulation of COR14b and some proteins important in photosynthetic acclimation to cold were studied with western analysis. COR14b transcript accumulation during cold acclimation was assessed with real-time PCR. KEY RESULTS: Cold acclimation increased both freezing tolerance and high-light tolerance, especially when plants were treated with high light after non-lethal freezing. In all plants, cold acclimation triggered the increase in photosynthetic capacity during high-light treatment. In two plants that were characterized by higher high-light tolerance but lower freezing tolerance, higher accumulation of COR14b transcript and protein was observed after 7 d and 14 d of cold acclimation, while a higher transient induction of COR14b expression was observed in freezing-tolerant plants during the first day of cold acclimation. High-light tolerant plants were also characterized with a higher level of PsbS accumulation and more efficient dissipation of excess light energy. CONCLUSIONS: Accumulation of COR14b in barley seems to be important for resistance to combined freezing and high-light tolerance, but not for freezing tolerance per se.

A new diagnostic SSR marker for selection of the Rym4/Rym5 locus in barley breeding.

Genomic sequence AY661558, representing a part of the BAC contig of the Rym4/Rym5 locus conferring resistance to the barley yellow mosaic virus complex (BaMMV/BaYMV), was exploited in order to develop SSR markers for practical barley breeding. Out of 57 SSR motifs found within this sequence, primers were designed and tested for the 5 SSRs with the highest repeat length. The polymorphic SSR marker QLB1 co-segregated with rym4 and rym5 phenotypes in respective high-resolution mapping populations developed for the construction of the original BAC contig. The primers targeted 2 sites located 756 bp and 5173 bp downstream of the translation initiation factor 4E (Hv-eIF4E). Physical linkage of the QLB1 marker to the Rym4/Rym5 locus was confirmed experimentally on Morex BAC 519J14, a seed BAC of Hv-eIF4E, and BAC 801A11, which is located proximally to Hv-eIF4E. QLB1 revealed 7 alleles in a set of 100 winter barley lines and cultivars. Five alleles were found within 673 advanced breeding lines derived from applied Polish winter barley breeding programmes, which corresponds to a PIC value of 0.684. No recombinants between Rym4/5 and QLB1 were detected, suggesting that QLB1 can be used efficiently in marker-assisted selection of the Hv-eIF4E-mediated bymovirus resistance.

Can telomere shortening be the main indicator of non-viable fetus elimination?

BACKGROUND: Telomeres are transcriptionally inactive genomic areas, which, if shortened, are associated with pathological processes, unsuccessful fertilization, aging, and death. Telomere dysfunction has also been linked to chromosomal rearrangements and genomic instability. The role of telomeres in postnatal life has been extensively studied and discussed both in physiological as well as in pathological processes. However, the role of telomere length in prenatal development is still poorly understood, and mainly concerns the preimplantation stage. The aim of this study was to estimate relative telomere length in spontaneously eliminated human embryos between 5th and 12th week of gestation. RESULTS: Relative telomere length was measured from total genomic DNA using a real-time polymerase chain reaction approach. In this study, we examined relative telomere length in 80 spontaneously eliminated embryos and in 25 embryos eliminated due to induced abortions. Relative telomere length in spontaneous abortions was significantly lower (P = 0.000001) compared to the induced abortions. Spontaneous abortions with aneuploid anomalies (monosomy X, trisomy 21, trisomy 16 and triploidy) were characterized by shorter telomeres, compared to spontaneous abortions, subgroup with euploid (46,XN) karyotype. CONCLUSION: Spontaneously lost pregnancies are characterized by shortened telomeres, especially in embryos with aneuploidies. We hypothesize that the shortening of telomeres is involved in the processes leading to spontaneous abortions.

Combination of three miRNA (miR-141, miR-21, and miR-375) as potential diagnostic tool for prostate cancer recognition.

PURPOSE: Prostate cancer (PCa) is a common tumor disease in western countries and a leading cause of cancer-driven mortality in men. Current methods for prostate cancer detection, like prostate-specific antigen screening, lead to significant overtreatment. The purpose of the study was to analyze circulating microRNAs in serum as non-invasive biomarkers in patients with diagnosis of prostate cancer and healthy individuals. METHODS: This preliminary study included a population of 20 patients with mean age of 68.6 years and mean PSA of 21.3 ng/ml. Eight healthy patients were used as control. MiRNAs were quantified in the total RNA fraction extracted from serum and levels of five microRNAs (miR-106b, miR-141, miR-21, mir-34a, and miR-375) were quantified by RT-qPCR. Statistical analyses evaluated correlation between clinicopathological data and miRNAs expression levels. RESULTS: Relative expression ratios of miR-106b, miR-141-3p, miR-21, and miR-375 were significantly increased (1.8-, 1.9-, 2.4-, and 2.6-fold, respectively) in the PCa group compared to healthy control. Using receiver operating characteristics, the highest area under the curve equal to 0.906 was obtained for miR-357 and indicates a very good diagnostic properties of this biomarker. We found expression level of mir-34a not related with PCa. CONCLUSIONS: Our results support previous findings on the possibility of discriminating prostate cancer patients from healthy controls by detecting miRNA (miR-141-3p, miR-21, and miR-375). Further insights into miRNA abundance and characteristics are necessary to validate the panel of miRNA as surrogate markers in diagnosis of prostate cancer.

The application of GBS markers for extending the dense genetic map of rye (Secale cereale L.) and the localization of the Rfc1 gene restoring male fertility in plants with the C source of sterility-inducing cytoplasm.

Genotyping by sequencing (GBS) is an efficient method of genotyping in numerous plant species. One of the crucial steps toward the application of GBS markers in crop improvement is anchoring them on particular chromosomes. In rye (Secale cereale L.), chromosomal localization of GBS markers has not yet been reported. In this paper, the application of GBS markers generated by the DArTseq platform for extending the high-density map of rye is presented. Additionally, their application is used for the localization of the Rfc1 gene that restores male fertility in plants with the C source of sterility-inducing cytoplasm. The total number of markers anchored on the current version of the map is 19,081, of which 18,132 were obtained from the DArTseq platform. Numerous markers co-segregated within the studied mapping population, so, finally, only 3397 unique positions were located on the map of all seven rye chromosomes. The total length of the map is 1593 cM and the average distance between markers is 0.47 cM. In spite of the resolution of the map being not very high, it should be a useful tool for further studies of the Secale cereale genome because of the presence on this map of numerous GBS markers anchored for the first time on rye chromosomes. The Rfc1 gene was located on high-density maps of the long arm of the 4R chromosome obtained for two mapping populations. Genetic maps were composed of DArT, DArTseq, and PCR-based markers. Consistent mapping results were obtained and DArTs tightly linked to the Rfc1 gene were successfully applied for the development of six new PCR-based markers useful in marker-assisted selection.