Comparative mapping of DNA sequences in rye (Secale cereale L.) in relation to the rice genome

A Genome-Wide Association Study Pinpoints Quantitative Trait Genes for Plant Height, Heading Date, Grain Quality, and Yield in Rye (Secale cereale L.)

Rye is the only cross-pollinating Triticeae crop species. Knowledge of rye genes controlling complex-inherited traits is scarce, which, currently, largely disables the genomics assisted introgression of untapped genetic variation from self-incompatible germplasm collections in elite inbred lines for hybrid breeding. We report on the first genome-wide association study (GWAS) in rye based on the phenotypic evaluation of 526 experimental hybrids for plant height, heading date, grain quality, and yield in 2 years and up to 19 environments. We established a cross-validated NIRS calibration model as a fast, effective, and robust analytical method to determine grain quality parameters. We observed phenotypic plasticity in plant height and tiller number as a resource use strategy of rye under drought and identified increased grain arabinoxylan content as a striking phenotype in osmotically stressed rye. We used DArTseq™ as a genotyping-by-sequencing technology to reduce the complexity of the rye genome. We established a novel high-density genetic linkage map that describes the position of almost 19k markers and that allowed us to estimate a low genome-wide LD based on the assessed genetic diversity in elite germplasm. We analyzed the relationship between plant height, heading date, agronomic, as well as grain quality traits, and genotype based on 20k novel single-nucleotide polymorphism markers. In addition, we integrated the DArTseq™ markers in the recently established 'Lo7' reference genome assembly. We identified cross-validated SNPs in 'Lo7' protein-coding genes associated with all traits studied. These include associations of the WUSCHEL-related homeobox transcription factor DWT1 and grain yield, the DELLA protein gene SLR1 and heading date, the Ethylene overproducer 1-like protein gene ETOL1 and thousand-grain weight, protein and starch content, as well as the Lectin receptor kinase SIT2 and plant height. A Leucine-rich repeat receptor protein kinase and a Xyloglucan alpha-1,6-xylosyltransferase count among the cross-validated genes associated with water-extractable arabinoxylan content. This study demonstrates the power of GWAS, hybrid breeding, and the reference genome sequence in rye genetics research to dissect and identify the function of genes shaping genetic diversity in agronomic and grain quality traits of rye. The described links between genetic causes and phenotypic variation will accelerate genomics-enabled rye improvement.

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.

Gene Expression Profiling and Fine Mapping Identifies a Gibberellin 2-Oxidase Gene Co-segregating With the Dominant Dwarfing Gene Ddw1 in Rye (Secale cereale L.)

The gibberellin (GA)-sensitive dwarfing gene Ddw1 provides an opportunity to genetically reduce plant height in rye. Genetic analysis in a population of recombinant inbred lines confirmed a monogenetic dominant inheritance of Ddw1. Significant phenotypic differences in PH between homo- and heterozygotic genotypes indicate an incomplete dominance of Ddw1. De novo transcriptome sequencing of Ddw1 mutant as well as tall genotypes resulted in 113,547 contigs with an average length of 318 bp covering 36.18 Mbp rye DNA. A hierarchical cluster analysis based on individual groups of rye homologs of functionally characterized rice genes controlling morphological or physiological traits including plant height, flowering time, and source activity identified the gene expression profile of stems at the begin of heading to most comprehensively mirror effects of Ddw1. Genome-wide expression profiling identified 186 transcripts differentially expressed between semi-dwarf and tall genotypes in stems. In total, 29 novel markers have been established and mapped to a 27.2 cM segment in the distal part of the long arm of chromosome 5R. Ddw1 could be mapped within a 0.4 cM interval co-segregating with a marker representing the C20-GA2-oxidase gene ScGA2ox12, that is up-regulated in stems of Ddw1 genotypes. The increased expression of ScGA2ox12 observed in semi-dwarf rye as well as structural alterations in transcript sequences associated with the ScGA2ox12 gene implicate, that Ddw1 is a dominant gain-of-function mutant. Integration of the target interval in the wheat reference genome sequence indicated perfect micro-colinearity between the Ddw1 locus and a 831 kb segment on chromosome 5A, which resides inside of a 11.21 Mb interval carrying the GA-sensitive dwarfing gene Rht12 in wheat. The potential of Ddw1 as a breeder's option to improve lodging tolerance in rye is discussed.

Exploring new alleles for frost tolerance in winter rye

Rye genetic resources provide a valuable source of new alleles for the improvement of frost tolerance in rye breeding programs. Frost tolerance is a must-have trait for winter cereal production in northern and continental cropping areas. Genetic resources should harbor promising alleles for the improvement of frost tolerance of winter rye elite lines. For frost tolerance breeding, the identification of quantitative trait loci (QTL) and the choice of optimum genome-based selection methods are essential. We identified genomic regions involved in frost tolerance of winter rye by QTL mapping in a biparental population derived from a highly frost tolerant selection from the Canadian cultivar Puma and the European elite line Lo157. Lines per se and their testcrosses were phenotyped in a controlled freeze test and in multi-location field trials in Russia and Canada. Three QTL on chromosomes 4R, 5R, and 7R were consistently detected across environments. The QTL on 5R is congruent with the genomic region harboring the Frost resistance locus 2 (Fr-2) in Triticeae. The Puma allele at the Fr-R2 locus was found to significantly increase frost tolerance. A comparison of predictive ability obtained from the QTL-based model with different whole-genome prediction models revealed that besides a few large, also small QTL effects contribute to the genomic variance of frost tolerance in rye. Genomic prediction models assigning a high weight to the Fr-R2 locus allow increasing the selection intensity for frost tolerance by genome-based pre-selection of promising candidates.

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.

Fine mapping of the restorer gene Rfp3 from an Iranian primitive rye (Secale cereale L.)

A comparative genetics approach allowed to precisely determine the map position of the restorer gene Rfp3 in rye and revealed that Rfp3 and the restorer gene Rfm1 in barley reside at different positions in a syntenic 4RL/6HS segment. Cytoplasmic male sterility (CMS) is a reliable and striking genetic mechanism for hybrid seed production. Breeding of CMS-based hybrids in cereals requires the use of effective restorer genes as an indispensable pre-requisite. We report on the fine mapping of a restorer gene for the Pampa cytoplasm in winter rye that has been tapped from the Iranian primitive rye population Altevogt 14160. For this purpose, we have mapped 41 gene-derived markers to a 38.8 cM segment in the distal part of the long arm of chromosome 4R, which carries the restorer gene. Male fertility restoration was comprehensively analyzed in progenies of crosses between a male-sterile tester genotype and 21 recombinant as well as six non-recombinant BC₄S₂ lines. This approach allowed us to validate the position of this restorer gene, which we have designated Rfp3, on chromosome 4RL. Rfp3 was mapped within a 2.5 cM interval and cosegregated with the EST-derived marker c28385. The gene-derived conserved ortholog set (COS) markers enabled us to investigate the orthology of restorer genes originating from different genetic resources of rye as well as barley. The observed localization of Rfp3 and Rfm1 in a syntenic 4RL/6HS segment asks for further efforts towards cloning of both restorer genes as an option to study the mechanisms of male sterility and fertility restoration in cereals.

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.

Genetic Map of Triticale Integrating Microsatellite, DArT and SNP Markers

Triticale (×Triticosecale Wittm) is an economically important crop for fodder and biomass production. To facilitate the identification of markers for agronomically important traits and for genetic and genomic characteristics of this species, a new high-density genetic linkage map of triticale was constructed using doubled haploid (DH) population derived from a cross between cultivars 'Hewo' and 'Magnat'. The map consists of 1615 bin markers, that represent 50 simple sequence repeat (SSR), 842 diversity array technology (DArT), and 16888 DArTseq markers mapped onto 20 linkage groups assigned to the A, B, and R genomes of triticale. No markers specific to chromosome 7R were found, instead mosaic linkage group composed of 1880 highly distorted markers (116 bins) from 10 wheat chromosomes was identified. The genetic map covers 4907 cM with a mean distance between two bins of 3.0 cM. Comparative analysis in respect to published maps of wheat, rye and triticale revealed possible deletions in chromosomes 4B, 5A, and 6A, as well as inversion in chromosome 7B. The number of bin markers in each chromosome varied from 24 in chromosome 3R to 147 in chromosome 6R. The length of individual chromosomes ranged between 50.7 cM for chromosome 2R and 386.2 cM for chromosome 7B. A total of 512 (31.7%) bin markers showed significant (P < 0.05) segregation distortion across all chromosomes. The number of 8 the segregation distorted regions (SDRs) were identified on 1A, 7A, 1B, 2B, 7B (2 SDRs), 5R and 6R chromosomes. The high-density genetic map of triticale will facilitate fine mapping of quantitative trait loci, the identification of candidate genes and map-based cloning.

CNMS: The preferred genic markers for comparative genomic, molecular phylogenetic, functional genetic diversity and differential gene regulatory expression analyses in chickpea

The intra/inter-genomic comparative mapping-based phylogenetic footprinting identified 5 paralogous and 656 orthologous genome-wide CNMS markers in the upstream sequences of chickpea genes. These CNMS markers revealed a high-degree of gene-based syntenic relationship between chickpea and Medicago genomes while minimum between chickpea and Vitis genomes. The time of divergence and duplication estimated using CNMS markers highlight the expected phylogenetic relationships between chickpea and six dicot (legume) species as well as occurrence of ancient genome (approximately 53 Mya) with small-scale recent segmental (approximately 10 Mya) duplication events in chickpea. A wider level of functional molecular diversity (14 to 88 percent) and admixed population genetic structure was detected among desi, kabuli and wild genotypes by genic CNMS markers at a genome-wide scale suggesting their utility in large-scale genetic analysis in chickpea. The subfunctionalization at the cis-regulatory element region and TFBS (transcription factor binding site) motif levels in the upstream sequences of CNMS marker-associated orthologous genes than the paralogues was predominant. Functional constraint might have considerable effect on these CNMScontaining regulatory elements controlling consistent orthologous gene expression in dicots. A rapid subfunctionalization based on diverge differential expression of paralogous CNMS marker-associated genes particularly those that underwent recent small-scale segmental duplication events in chickpea was apparent. The differential regulation of expression and subfunctionalization potential of Ultra CNMS marker-associated genes suggest their utility in deciphering the complex gene regulatory function as well as identification and targeted mapping of potential genes/QTLs governing vital agronomic traits in chickpea. The gene-based CNMS markers with desirable inherent genetic attributes like higher degree of comparative genome mapping, functional genetic diversity and differential gene regulatory expression potential can significantly propel the genomics-assisted chickpea crop improvement.

Assessment of Genetic Diversity in Secale cereale Based on SSR Markers

The primary aim of this study was to estimate genetic diversity among Secale cereale L. accessions using 22 previously published simple sequence repeat (SSR) markers. The plant material included 367 rye accessions comprising historical and contemporary cultivars, cultivated materials, landraces, and breeding strains from the Polish breeding company Danko. The studied accessions represented a wide geographical diversity. Several methods were employed to analyze genetic diversity among the Secale cereale L. accessions and to determine population structure: principal coordinate analysis (PCoA), neighbor-joining (NJ), and Bayesian clustering. We also defined a core collection of 25 rye accessions representing over 93 % of SSR alleles. The results of these analyses showed that accessions from the rye gene bank are clearly divergent in comparison with materials received directly from European breeding companies. Our findings suggest also that the genetic pool of current rye cultivars is becoming narrower during breeding processes. The selected panel of SSR markers performed well in detection of genetic diversity patterns and can be recommended for future germplasm characterization studies in rye.

Development of T. aestivum L.-H. californicum alien chromosome lines and assignment of homoeologous groups of Hordeum californicum chromosomes

Hordeum californicum (2n = 2x = 14, HH) is resistant to several wheat diseases and tolerant to lower nitrogen. In this study, a molecular karyotype of H. californicum chromosomes in the Triticum aestivum L. cv. Chinese Spring (CS)-H. californicum amphidiploid (2n = 6x = 56, AABBDDHH) was established. By genomic in situ hybridization (GISH) and multicolor fluorescent in situ hybridization (FISH) using repetitive DNA clones (pTa71, pTa794 and pSc119.2) as probes, the H. californicum chromosomes could be differentiated from each other and from the wheat chromosomes unequivocally. Based on molecular karyotype and marker analyses, 12 wheat-alien chromosome lines, including four disomic addition lines (DAH1, DAH3, DAH5 and DAH6), five telosomic addition lines (MtH7L, MtH1S, MtH1L, DtH6S and DtH6L), one multiple addition line involving H. californicum chromosome H2, one disomic substitution line (DSH4) and one translocation line (TH7S/1BL), were identified from the progenies derived from the crosses of CS-H. californicum amphidiploid with common wheat varieties. A total of 482 EST (expressed sequence tag) or SSR (simple sequence repeat) markers specific for individual H. californicum chromosomes were identified, and 47, 50, 45, 49, 21, 51 and 40 markers were assigned to chromosomes H1, H2, H3, H4, H5, H6 and H7, respectively. According to the chromosome allocation of these markers, chromosomes H2, H3, H4, H5, and H7 of H. californicum have relationship with wheat homoeologous groups 5, 2, 6, 3, and 1, and hence could be designated as 5H(c), 2H(c), 6H(c), 3H(c) and 1H(c), respectively. The chromosomes H1 and H6 were designated as 7H(c) and 4H(c), respectively, by referring to SSR markers located on rye chromosomes.

Homoeologous relationship of rye chromosome arms as detected with wheat PLUG markers

Based on the similarity in gene structure between rice and wheat, the polymerase chain reaction (PCR)-based landmark unique gene (PLUG) system enabled us to design primer sets that amplify wheat genic sequences including introns. From the previously reported wheat PLUG markers, we chose 144 markers that are distributed on different chromosomes and in known chromosomal regions (bins) to obtain rye-specific PCR-based markers. We conducted PCR with the 144 primer sets and the template of the Imperial rye genomic DNA and found that 131 (91.0%) primer sets successfully amplified PCR products. Of the 131 PLUG markers, 110 (76.4%) markers showed rye-specific PCR amplification with or without restriction enzyme digestion. We assigned 79 of the 110 markers to seven rye chromosomes (1R to 7R) using seven wheat-rye (cv. Imperial) chromosome addition and substitution lines: 12 to 1R, 8 to 2R, 11 to 3R, 8 to 4R, 16 to 5R, 12 to 6R, and 12 to 7R. Furthermore, we located their positions on the short or long (L) chromosome arm, using 13 Imperial rye telosomic lines of common wheat (except for 3RL). Referring to the chromosome bin locations of the 79 PLUG markers in wheat, we deduced the syntenic relationships between rye and wheat chromosomes. We also discussed chromosomal rearrangements in the rye genome with reference to the cytologically visible chromosomal gaps.

The perennial ryegrass GenomeZipper: targeted use of genome resources for comparative grass genomics

Whole-genome sequences established for model and major crop species constitute a key resource for advanced genomic research. For outbreeding forage and turf grass species like ryegrasses (Lolium spp.), such resources have yet to be developed. Here, we present a model of the perennial ryegrass (Lolium perenne) genome on the basis of conserved synteny to barley (Hordeum vulgare) and the model grass genome Brachypodium (Brachypodium distachyon) as well as rice (Oryza sativa) and sorghum (Sorghum bicolor). A transcriptome-based genetic linkage map of perennial ryegrass served as a scaffold to establish the chromosomal arrangement of syntenic genes from model grass species. This scaffold revealed a high degree of synteny and macrocollinearity and was then utilized to anchor a collection of perennial ryegrass genes in silico to their predicted genome positions. This resulted in the unambiguous assignment of 3,315 out of 8,876 previously unmapped genes to the respective chromosomes. In total, the GenomeZipper incorporates 4,035 conserved grass gene loci, which were used for the first genome-wide sequence divergence analysis between perennial ryegrass, barley, Brachypodium, rice, and sorghum. The perennial ryegrass GenomeZipper is an ordered, information-rich genome scaffold, facilitating map-based cloning and genome assembly in perennial ryegrass and closely related Poaceae species. It also represents a milestone in describing synteny between perennial ryegrass and fully sequenced model grass genomes, thereby increasing our understanding of genome organization and evolution in the most important temperate forage and turf grass species.

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.

Development and application of EST-based markers specific for chromosome arms of rye (Secale cereale L.)

To develop a set of molecular markers specific for the chromosome arms of rye, a total of 1,098 and 93 primer pairs derived from the expressed sequence tag (EST) sequences distributed on all 21 wheat chromosomes and 7 rye chromosomes, respectively, were initially screened on common wheat 'Chinese Spring' and rye cultivar 'Imperial'. Four hundred and fourteen EST-based markers were specific for the rye genome. Seven disomic chromosome addition lines, 10 telosomic addition lines and 1 translocation line of 'Chinese Spring-Imperial' were confirmed by genomic in situ hybridization and fluorescencein situ hybridization, and used to screen the rye-specific markers. Thirty-one of the 414 markers produced stable specific amplicons in 'Imperial', as well as individual addition lines and were assigned to 13 chromosome arms of rye except for 6RS. Six rye cultivars, wheat cultivar 'Xiaoyan 6' and accessions of 4 wheat relatives were then used to test the specificity of the 31 EST-based markers. To confirm the specificity, 4 wheat-rye derivatives of 'Xiaoyan 6 × German White', with chromosomes 1RS, 2R and 4R, were amplified by some of the EST-based markers. The results indicated that they can effectively be used to detect corresponding rye chromosomes or chromosome arms introgressed into a wheat background, and hence to accelerate the utilization of rye genes in wheat breeding.

Sequence composition and gene content of the short arm of rye (Secale cereale) chromosome 1

Background

The purpose of the study is to elucidate the sequence composition of the short arm of rye chromosome 1 (Secale cereale) with special focus on its gene content, because this portion of the rye genome is an integrated part of several hundreds of bread wheat varieties worldwide.

Methodology/Principal Findings

Multiple Displacement Amplification of 1RS DNA, obtained from flow sorted 1RS chromosomes, using 1RS ditelosomic wheat-rye addition line, and subsequent Roche 454FLX sequencing of this DNA yielded 195,313,589 bp sequence information. This quantity of sequence information resulted in 0.43× sequence coverage of the 1RS chromosome arm, permitting the identification of genes with estimated probability of 95%. A detailed analysis revealed that more than 5% of the 1RS sequence consisted of gene space, identifying at least 3,121 gene loci representing 1,882 different gene functions. Repetitive elements comprised about 72% of the 1RS sequence, Gypsy/Sabrina (13.3%) being the most abundant. More than four thousand simple sequence repeat (SSR) sites mostly located in gene related sequence reads were identified for possible marker development. The existence of chloroplast insertions in 1RS has been verified by identifying chimeric chloroplast-genomic sequence reads. Synteny analysis of 1RS to the full genomes of Oryza sativa and Brachypodium distachyon revealed that about half of the genes of 1RS correspond to the distal end of the short arm of rice chromosome 5 and the proximal region of the long arm of Brachypodium distachyon chromosome 2. Comparison of the gene content of 1RS to 1HS barley chromosome arm revealed high conservation of genes related to chromosome 5 of rice.

Conclusions

The present study revealed the gene content and potential gene functions on this chromosome arm and demonstrated numerous sequence elements like SSRs and gene-related sequences, which can be utilised for future research as well as in breeding of wheat and rye.

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.

From RNA-seq to large-scale genotyping - genomics resources for rye (Secale cereale L.)

BACKGROUND

The improvement of agricultural crops with regard to yield, resistance and environmental adaptation is a perpetual challenge for both breeding and research. Exploration of the genetic potential and implementation of genome-based breeding strategies for efficient rye (Secale cereale L.) cultivar improvement have been hampered by the lack of genome sequence information. To overcome this limitation we sequenced the transcriptomes of five winter rye inbred lines using Roche/454 GS FLX technology.

RESULTS

More than 2.5 million reads were assembled into 115,400 contigs representing a comprehensive rye expressed sequence tag (EST) resource. From sequence comparisons 5,234 single nucleotide polymorphisms (SNPs) were identified to develop the Rye5K high-throughput SNP genotyping array. Performance of the Rye5K SNP array was investigated by genotyping 59 rye inbred lines including the five lines used for sequencing, and five barley, three wheat, and two triticale accessions. A balanced distribution of allele frequencies ranging from 0.1 to 0.9 was observed. Residual heterozygosity of the rye inbred lines varied from 4.0 to 20.4% with higher average heterozygosity in the pollen compared to the seed parent pool.

CONCLUSIONS

The established sequence and molecular marker resources will improve and promote genetic and genomic research as well as genome-based breeding in rye.

High-Resolution Genotyping of Wild Barley Introgression Lines and Fine-Mapping of the Threshability Locus thresh-1 Using the Illumina GoldenGate Assay

Genetically well-characterized mapping populations are a key tool for rapid and precise localization of quantitative trait loci (QTL) and subsequent identification of the underlying genes. In this study, a set of 73 introgression lines (S42ILs) originating from a cross between the spring barley cultivar Scarlett (Hordeum vulgare ssp. vulgare) and the wild barley accession ISR42-8 (H. v. ssp. spontaneum) was subjected to high-resolution genotyping with an Illumina 1536-SNP array. The array enabled a precise localization of the wild barley introgressions in the elite barley background. Based on 636 informative SNPs, the S42IL set represents 87.3% of the wild barley genome, where each line contains on average 3.3% of the donor genome. Furthermore, segregating high-resolution mapping populations (S42IL-HRs) were developed for 70 S42ILs in order to facilitate QTL fine-mapping and cloning. As a case study, we used the developed genetic resources to rapidly identify and fine-map the novel locus thresh-1 on chromosome 1H that controls grain threshability. Here, the recessive wild barley allele confers a difficult to thresh phenotype, suggesting that thresh-1 played an important role during barley domestication. Using a S42IL-HR population, thresh-1 was fine-mapped within a 4.3cM interval that was predicted to contain candidate genes involved in regulation of plant cell wall composition. The set of wild barley introgression lines and derived high-resolution populations are ideal tools to speed up the process of mapping and further dissecting QTL, which ultimately clears the way for isolating the genes behind QTL effects.

Genetic map of triticale compiling DArT, SSR, and AFLP markers

A set of 90 doubled haploid (DH) lines derived from F(1) plants that originated from a cross between × Triticosecale Wittm. 'Saka3006' and ×Triticosecale Wittm. 'Modus', via wide crossing with maize, were used to create a genetic linkage map of triticale. The map has 21 linkage groups assigned to the A, B, and R genomes including 155 simple sequence repeat (SSR), 1385 diversity array technology (DArT), and 28 amplified fragment length polymorphism (AFLP) markers covering 2397 cM with a mean distance between two markers of 4.1 cM. Comparative analysis with wheat consensus maps revealed that triticale chromosomes of the A and B genomes were represented by 15 chromosomes, including combinations of 2AS.2AL#, 2AL#2BL, 6AS.6AL#, and 2BS.6AL# instead of 2A, 2B, and 6A. In respect to published maps of rye, substantial rearrangements were found also for chromosomes 1R, 2R, and 3R of the rye genome. Chromosomes 1R and 2R were truncated and the latter was linked with 3R. A nonhomogeneous distribution of markers across the triticale genome was observed with evident bias (48%) towards the rye genome. This genetic map may serve as a reference linkage map of triticale for efficient studies of structural rearrangements, gene mapping, and marker-assisted selection.

Mapping of sequence-specific markers and loci controlling preharvest sprouting and alpha-amylase activity in rye (Secale cereale L.) on the genetic map of an F2 (S120&times;S76) population

Location of the loci that control preharvest sprouting and alpha-amylase activity in rye was studied based on intercross S120&times;S76, consisting of 110 genotypes of F2 and F3 progenies. The genetic map currently consists of 141 loci distributed in 11 linkage groups, covering a distance of 506.4 cM, and was enriched during this study with 24 sequence-specific markers (7 SCARs, 7 SSRs, and 10 STSs). The extended map was applied for composite interval mapping of the loci controlling preharvest sprouting and &alpha;-amylase activity, revealing 3 significant QTLs for preharvest sprouting, located on chromosomes 3R, 5R and 6R (in 1999), and one QTL for &alpha;-amylase activity found on chromosome 2R (in 2000).

Dissection of QTL effects for root traits using a chromosome arm-specific mapping population in bread wheat

A high-resolution chromosome arm-specific mapping population was used in an attempt to locate/detect gene(s)/QTL for different root traits on the short arm of rye chromosome 1 (1RS) in bread wheat. This population consisted of induced homoeologous recombinants of 1RS with 1BS, each originating from a different crossover event and distinct from all other recombinants in the proportions of rye and wheat chromatin present. It provides a simple and powerful approach to detect even small QTL effects using fewer progeny. A promising empirical Bayes method was applied to estimate additive and epistatic effects for all possible marker pairs simultaneously in a single model. This method has an advantage for QTL analysis in minimizing the error variance and detecting interaction effects between loci with no main effect. A total of 15 QTL effects, 6 additive and 9 epistatic, were detected for different traits of root length and root weight in 1RS wheat. Epistatic interactions were further partitioned into inter-genomic (wheat and rye alleles) and intra-genomic (rye-rye or wheat-wheat alleles) interactions affecting various root traits. Four common regions were identified involving all the QTL for root traits. Two regions carried QTL for almost all the root traits and were responsible for all the epistatic interactions. Evidence for inter-genomic interactions is provided. Comparison of mean values supported the QTL detection.

A doubled haploid rye linkage map with a QTL affecting α-amylase activity

A rye doubled haploid (DH) mapping population (Amilo × Voima) segregating for pre-harvest sprouting (PHS) was generated through anther culture of F(1) plants. A linkage map was constructed using DHs, to our knowledge, for the first time in rye. The map was composed of 289 loci: amplified fragment length polymorphism (AFLP), microsatellite, random amplified polymorphic DNA (RAPD), retrotransposon-microsatellite amplified polymorphism (REMAP), inter-retrotransposon amplified polymorphism (IRAP), inter-simple sequence repeat (ISSR) and sequence-related amplified polymorphism (SRAP) markers, and extended altogether 732 cM (one locus in every 2.5 cM). All of the seven rye chromosomes and four unplaced groups were formed. Distorted segregation of markers (P ≤ 0.05) was detected on all chromosomes. One major quantitative trait locus (QTL) affecting α-amylase activity was found, which explained 16.1% of phenotypic variation. The QTL was localized on the long arm of chromosome 5R. Microsatellites SCM74, RMS1115, and SCM77, nearest to the QTL, can be used for marker-assisted selection as a part of a rye breeding program to decrease sprouting damage.

Comparative physical mapping between wheat chromosome arm 2BL and rice chromosome 4

Physical maps of chromosomes provide a framework for organizing and integrating diverse genetic information. DNA microarrays are a valuable technique for physical mapping and can also be used to facilitate the discovery of single feature polymorphisms (SFPs). Wheat chromosome arm 2BL was physically mapped using a Wheat Genome Array onto near-isogenic lines (NILs) with the aid of wheat-rice synteny and mapped wheat EST information. Using high variance probe set (HVP) analysis, 314 HVPs constituting genes present on 2BL were identified. The 314 HVPs were grouped into 3 categories: HVPs that match only rice chromosome 4 (298 HVPs), those that match only wheat ESTs mapped on 2BL (1), and those that match both rice chromosome 4 and wheat ESTs mapped on 2BL (15). All HVPs were converted into gene sets, which represented either unique rice gene models or mapped wheat ESTs that matched identified HVPs. Comparative physical maps were constructed for 16 wheat gene sets and 271 rice gene sets. Of the 271 rice gene sets, 257 were mapped to the 18-35 Mb regions on rice chromosome 4. Based on HVP analysis and sequence similarity between the gene models in the rice chromosomes and mapped wheat ESTs, the outermost rice gene model that limits the translocation breakpoint to orthologous regions was identified.

Dosage effect of the short arm of chromosome 1 of rye on root morphology and anatomy in bread wheat

The spontaneous translocation of the short arm of chromosome 1 of rye (1RS) in bread wheat is associated with higher root biomass and grain yield. Recent studies have confirmed the presence of QTL for different root morphological traits on the 1RS arm in bread wheat. This study was conducted to address two questions in wheat root genetics. First, does the presence of the 1RS arm in bread wheat affect its root anatomy? Second, how does root morphology and anatomy of bread wheat respond to different dosages of 1RS? Near-isogenic plants with a different number (0 to 4 dosages) of 1RS translocations were studied for root morphology and anatomy. The F(1) hybrid, with single doses of the 1RS and 1AS arms, showed heterosis for root and shoot biomass. In other genotypes, with 0, 2, or 4 doses of 1RS, root biomass was incremental with the increase in the dosage of 1RS in bread wheat. This study also provided evidence of the presence of gene(s) influencing root xylem vessel number, size, and distribution in bread wheat. It was found that root vasculature follows a specific developmental pattern along the length of the tap root and 1RS dosage tends to affect the transitions differentially in different positions. This study indicated that the inherent differences in root morphology and anatomy of different 1RS lines may be advantageous compared to normal bread wheat to survive under stress conditions.

Genomics of fungal disease resistance in tomato

Tomato (Solanum lycopersicum) is an important vegetable crop worldwide. Often times, its production is hindered by fungal diseases. Important fungal diseases limiting tomato production are late blight, caused by Phytophthora infestans, early blight, caused by Alternaria solanii, and septoria leaf spot, caused by Septoria lycopersici, fusarium wilt caused by Fusarium oxysporium fsp. oxysporium, and verticilium wilt caused by Verticilium dahlea. The Phytophthora infestans is the same fungus that caused the devastating loss of potato in Europe in 1845. A similar magnitude of crop loss in tomato has not occurred but Phytophthora infestans has caused the complete loss of tomato crops around the world on a small scale. Several attempts have been made through conventional breeding and the molecular biological approaches to understand the biology of host-pathogen interaction so that the disease can be managed and crop loss prevented. In this review, we present a comprehensive analysis of information produced by molecular genetic and genomic experiments on host-pathogen interactions of late blight, early blight, septoria leaf spot, verticilim wilt and fusarium wilt in tomato. Furthermore, approaches adopted to manage these diseases in tomato including genetic transformation are presented. Attempts made to link molecular markers with putative genes and their use in crop improvement are discussed.

Integrated genetic map and genetic analysis of a region associated with root traits on the short arm of rye chromosome 1 in bread wheat

A rye-wheat centric chromosome translocation 1RS.1BL has been widely used in wheat breeding programs around the world. Increased yield of translocation lines was probably a consequence of increased root biomass. In an effort to map loci-controlling root characteristics, homoeologous recombinants of 1RS with 1BS were used to generate a consensus genetic map comprised of 20 phenotypic and molecular markers, with an average spacing of 2.5 cM. Physically, all recombination events were located in the distal 40% of the arms. A total of 68 recombinants was used and recombination breakpoints were aligned and ordered over map intervals with all the markers, integrated together in a genetic map. This approach enabled dissection of genetic components of quantitative traits, such as root traits, present on 1S. To validate our hypothesis, phenotyping of 45-day-old wheat roots was performed in five lines including three recombinants representative of the entire short arm along with bread wheat parents 'Pavon 76' and Pavon 1RS.1BL. Individual root characteristics were ranked and the genotypic rank sums were subjected to Quade analysis to compare the overall rooting ability of the genotypes. It appears that the terminal 15% of the rye 1RS arm carries gene(s) for greater rooting ability in wheat.

Development and functional assessment of EST-derived 2RL-specific markers for 2BS.2RL translocations

ESTs-derived markers are useful for comparative genomic analysis and can also serve as phenotype-linked functional markers. Here, we report the development of EST-derived 2RL-specific markers and the evaluation of the possibility of functional assessment of markers tagging 2RL, which carries Hessian fly resistance genes (loci). To identify transcripts specific to 2RL, unigene sequences in combination with wheat progenitor genomes were used. Total 275 contigs mapped to the long arms of homoeologous group 2 chromosomes were downloaded. To obtain a cluster corresponding to each of the wheat 275 contigs, unigene sequences of wheat, rice, barley, and rye were pooled for cross-species clusters. Out of 275 clusters examined, it was possible to design 112 cross-species primer pairs for genome-specific amplifications. Out of 112 cross-species primer pairs, 45 primer pairs (40%) produced amplicons from at least one species (three wheat progenitors or rye). Among the 45 contigs, 73% were associated with one of known functions and 82% of the contigs associated with known functions were also associated with one of the GO categories. On the basis of the oligonucleotide sequence alignment of each of 45 genome-specific amplifications, 21 amplifications (47%) were suitable for designing RR genome-specific primers, which are specific to translocated rye chromatin 2RL. Six primer pairs (13%) successfully produced amplicons in the 2BS.2RL translocation lines and not in the non-2RLs. Functional assessment of one of the 2RL-specific markers, NSFT03P2_Contig4445, was performed on Hessian fly infested NILs. Under Hessian fly infestation, significantly high expression of a gene tagged by a 2RL-specific marker (NSFT03P2_Contig4445) was observed 1 day after infestation. EST-derived 2RL-specific marker development from this study provides a basis for the development of ESTs-derived markers for detecting wheat-rye translocations. In addition, these markers could be employed in elucidating functional analysis of genes on 2RL.