High levels of nucleotide diversity and fast decline of linkage disequilibrium in rye (Secale cereale L.) genes involved in frost response

Unveiling the distinctive traits of functional rye centromeres: minisatellites, retrotransposons, and R-loop formation

Centromeres play a vital role in cellular division by facilitating kinetochore assembly and spindle attachments. Despite their conserved functionality, centromeric DNA sequences exhibit rapid evolution, presenting diverse sizes and compositions across species. The functional significance of rye centromeric DNA sequences, particularly in centromere identity, remains unclear. In this study, we comprehensively characterized the sequence composition and organization of rye centromeres. Our findings revealed that these centromeres are primarily composed of long terminal repeat retrotransposons (LTR-RTs) and interspersed minisatellites. We systematically classified LTR-RTs into five categories, highlighting the prevalence of younger CRS1, CRS2, and CRS3 of CRSs (centromeric retrotransposons of Secale cereale) were primarily located in the core centromeres and exhibited a higher association with CENH3 nucleosomes. The minisatellites, mainly derived from retrotransposons, along with CRSs, played a pivotal role in establishing functional centromeres in rye. Additionally, we observed the formation of R-loops at specific regions of CRS1, CRS2, and CRS3, with both rye pericentromeres and centromeres exhibiting enrichment in R-loops. Notably, these R-loops selectively formed at binding regions of the CENH3 nucleosome in rye centromeres, suggesting a potential role in mediating the precise loading of CENH3 to centromeres and contributing to centromere specification. Our work provides insights into the DNA sequence composition, distribution, and potential function of R-loops in rye centromeres. This knowledge contributes valuable information to understanding the genetics and epigenetics of rye centromeres, offering implications for the development of synthetic centromeres in future plant modifications and beyond.

Improving Yield and Yield Stability in Winter Rye by Hybrid Breeding

Rye is the only cross-pollinating small-grain cereal. The unique reproduction biology results in an exceptional complexity concerning genetic improvement of rye by breeding. Rye is a close relative of wheat and has a strong adaptation potential that refers to its mating system, making this overlooked cereal readily adjustable to a changing environment. Rye breeding addresses the emerging challenges of food security associated with climate change. The systematic identification, management, and use of its valuable natural diversity became a feasible option in outbreeding rye only following the establishment of hybrid breeding late in the 20th century. In this article, we review the most recent technological advances to improve yield and yield stability in winter rye. Based on recently released reference genome sequences, SMART breeding approaches are described to counterbalance undesired linkage drag effects of major restorer genes on grain yield. We present the development of gibberellin-sensitive semidwarf hybrids as a novel plant breeding innovation based on an approach that is different from current methods of increasing productivity in rye and wheat. Breeding of new rye cultivars with improved performance and resilience is indispensable for a renaissance of this healthy minor cereal as a homogeneous commodity with cultural relevance in Europe that allows for comparatively smooth but substantial complementation of wheat with rye-based diets, supporting the necessary restoration of the balance between human action and nature.

Association mapping of autumn-seeded rye (Secale cereale L.) reveals genetic linkages between genes controlling winter hardiness and plant development

Winter field survival (WFS) in autumn-seeded winter cereals is a complex trait associated with low temperature tolerance (LTT), prostrate growth habit (PGH), and final leaf number (FLN). WFS and the three sub-traits were analyzed by a genome-wide association study of 96 rye (Secale cereal L.) genotypes of different origins and winter-hardiness levels. A total of 10,244 single nucleotide polymorphism (SNP) markers were identified by genotyping by sequencing and 259 marker-trait-associations (MTAs; p < 0.01) were revealed by association mapping. The ten most significant SNPs (p < 1.49e−04) associated with WFS corresponded to nine strong candidate genes: Inducer of CBF Expression 1 (ICE1), Cold-regulated 413-Plasma Membrane Protein 1 (COR413-PM1), Ice Recrystallization Inhibition Protein 1 (IRIP1), Jasmonate-resistant 1 (JAR1), BIPP2C1-like protein phosphatase, Chloroplast Unusual Positioning Protein-1 (CHUP1), FRIGIDA-like 4 (FRL4-like) protein, Chalcone Synthase 2 (CHS2), and Phenylalanine Ammonia-lyase 8 (PAL8). Seven of the candidate genes were also significant for one or several of the sub-traits supporting the hypothesis that WFS, LTT, FLN, and PGH are genetically interlinked. The winter-hardy rye genotypes generally carried additional allele variants for the strong candidate genes, which suggested allele diversity was a major contributor to cold acclimation efficiency and consistent high WFS under varying field conditions.

Identification of New QTLs for Dietary Fiber Content in Aegilops biuncialis

Grain dietary fiber content is an important health-promoting trait of bread wheat. A dominant dietary fiber component of wheat is the cell wall polysaccharide arabinoxylan and the goatgrass Aegilops biuncialis has high β-glucan content, which makes it an attractive gene source to develop wheat lines with modified fiber composition. In order to support introgression breeding, this work examined genetic variability in grain β-glucan, pentosan, and protein content in a collection of Ae. biuncialis. A large variation in grain protein and edible fiber content was revealed, reflecting the origin of Ae. biuncialis accessions from different eco-geographical habitats. Association analysis using DArTseq-derived SNPs identified 34 QTLs associated with β-glucan, pentosan, water-extractable pentosan, and protein content. Mapping the markers to draft chromosome assemblies of diploid progenitors of Ae. biuncialis underlined the role of genes on chromosomes 1M^b, 4M^b, and 5M^b in the formation of grain β-glucan content, while other QTLs on chromosome groups 3, 6, and 1 identified genes responsible for total- and water-extractable pentosan content. Functional annotation of the associated marker sequences identified fourteen genes, nine of which were identified in other monocots. The QTLs and genes identified in the present work are attractive targets for chromosome-mediated gene transfer to improve the health-promoting properties of wheat-derived foods.

Genomic prediction of preliminary yield trials in chickpea: Effect of functional annotation of SNPs and environment

Achieving yield potential in chickpea (Cicer arietinum L.) is limited by many constraints that include biotic and abiotic stresses. Combining next-generation sequencing technology with advanced statistical modeling has the potential to increase genetic gain efficiently. Whole genome resequencing data was obtained from 315 advanced chickpea breeding lines from the Australian chickpea breeding program resulting in more than 298,000 single nucleotide polymorphisms (SNPs) discovered. Analysis of population structure revealed a distinct group of breeding lines with many alleles that are absent from recently released Australian cultivars. Genome-wide association studies (GWAS) using these Australian breeding lines identified 20 SNPs significantly associated with grain yield in multiple field environments. A reduced level of nucleotide diversity and extended linkage disequilibrium suggested that some regions in these chickpea genomes may have been through selective breeding for yield or other traits. A large introgression segment that introduced from C. echinospermum for phytophthora root rot resistance was identified on chromosome 6, yet it also has unintended consequences of reducing yield due to linkage drag. We further investigated the effect of genotype by environment interaction on genomic prediction of yield. We found that the training set had better prediction accuracy when phenotyped under conditions relevant to the targeted environments. We also investigated the effect of SNP functional annotation on prediction accuracy using different subsets of SNPs based on their genomic locations: regulatory regions, exome, and alternative splice sites. Compared with the whole SNP dataset, a subset of SNPs did not significantly decrease prediction accuracy for grain yield despite consisting of a smaller number of SNPs.

Discovery of a novel powdery mildew (Blumeria graminis) resistance locus in rye (Secale cereale L.)

Powdery mildew is one of the most destructive diseases in the world, causing substantial grain yield losses and quality reduction in cereal crops. At present 23 powdery mildew resistance genes have been identified in rye, of which the majority are in wheat-rye translocation lines developed for wheat improvement. Here, we investigated the genetics underlying powdery mildew resistance in the Gülzow-type elite hybrid rye (Secale cereale L.) breeding germplasm. In total, 180 inbred breeding lines were genotyped using the state-of-the-art 600 K SNP array and phenotyped for infection type against three distinct field populations of B. graminis f. sp. secalis from Northern Germany (2013 and 2018) and Denmark (2020). We observed a moderate level of powdery mildew resistance in the non-restorer germplasm population, and by performing a genome-wide association study using 261,406 informative SNP markers, we identified a powdery mildew resistance locus, provisionally denoted PmNOS1, on the distal tip of chromosome arm 7RL. Using recent advances in rye genomic resources, we investigated whether nucleotide-binding leucine-rich repeat genes residing in the identified 17 Mbp block associated with PmNOS1 on recent reference genomes resembled known Pm genes.

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.

Gene-based mapping of trehalose biosynthetic pathway genes reveals association with source- and sink-related yield traits in a spring wheat panel

Trehalose 6‐phosphate (T6P) signalling regulates carbon use and allocation and is a target to improve crop yields. However, the specific contributions of trehalose phosphate synthase (TPS) and trehalose phosphate phosphatase (TPP) genes to source‐ and sink‐related traits remain largely unknown. We used enrichment capture sequencing on TPS and TPP genes to estimate and partition the genetic variation of yield‐related traits in a spring wheat (Triticum aestivum) breeding panel specifically built to capture the diversity across the 75,000 CIMMYT wheat cultivar collection. Twelve phenotypes were correlated to variation in TPS and TPP genes including plant height and biomass (source), spikelets per spike, spike growth and grain filling traits (sink) which showed indications of both positive and negative gene selection. Individual genes explained proportions of heritability for biomass and grain‐related traits. Three TPS1 homologues were particularly significant for trait variation. Epistatic interactions were found within and between the TPS and TPP gene families for both plant height and grain‐related traits. Gene‐based prediction improved predictive ability for grain weight when gene effects were combined with the whole‐genome markers. Our study has generated a wealth of information on natural variation of TPS and TPP genes related to yield potential which confirms the role for T6P in resource allocation and in affecting traits such as grain number and size confirming other studies which now opens up the possibility of harnessing natural genetic variation more widely to better understand the contribution of native genes to yield traits for incorporation into breeding programmes.

Quantitative trait loci and candidate genes associated with freezing tolerance of winter triticale (× Triticosecale Wittmack)

Freezing tolerance of triticale is a major trait contributing to its winter hardiness. The identification of genomic regions — quantitative trait loci (QTL) and molecular markers associated with freezing tolerance in winter hexaploid triticale — was the aim of this study. For that purpose, a new genetic linkage map was developed for the population of 92 doubled haploid lines derived from ‘Hewo’ × ‘Magnat’ F₁ hybrid. Those lines, together with parents were subjected to freezing tolerance test three times during two winter seasons. Plants were grown and cold-hardened under natural fall/winter conditions and then subjected to freezing in controlled conditions. Freezing tolerance was assessed as the plants recovery (REC), the electrolyte leakage (EL) from leaves and chlorophyll fluorescence parameters (JIP) after freezing. Three consistent QTL for several fluorescence parameters, electrolyte leakage, and the percentage of the survived plants were identified with composite interval mapping (CIM) and single marker analysis (SMA). The first locus Qfr.hm-7A.1 explained 9% of variation of both electrolyte leakage and plants recovery after freezing. Two QTL explaining up to 12% of variation in plants recovery and shared by selected chlorophyll fluorescence parameters were found on 4R and 5R chromosomes. Finally, main locus Qchl.hm-5A.1 was detected for chlorophyll fluorescence parameters that explained up to 19.6% of phenotypic variation. The co-located QTL on chromosomes 7A.1, 4R and 5R, clearly indicated physiological and genetic relationship of the plant survival after freezing with the ability to maintain optimal photochemical activity of the photosystem II and preservation of the cell membranes integrity. The genes located in silico within the identified QTL include those encoding BTR1-like protein, transmembrane helix proteins like potassium channel, and phosphoric ester hydrolase involved in response to osmotic stress as well as proteins involved in the regulation of the gene expression, chloroplast RNA processing, and pyrimidine salvage pathway. Additionally, our results confirm that the JIP test is a valuable tool to evaluate freezing tolerance of triticale under unstable winter environments.

Population genomics of Mediterranean oat (A. sativa) reveals high genetic diversity and three loci for heading date

Genomic analysis of Mediterranean oats reveals high genetic diversity and three loci for adaptation to this environment. This information together with phenotyping and passport data, gathered in an interactive map, will be a vital resource for oat genetic improvement. During the twentieth century, oat landraces have increasingly been replaced by modern cultivars, resulting in loss of genetic diversity. However, landraces have considerable potential to improve disease and abiotic stress tolerance and may outperform cultivars under low input systems. In this work, we assembled a panel of 669 oat landraces from Mediterranean rim and 40 cultivated oat varieties and performed the first large-scale population genetic analysis of both red and white oat types of Mediterranean origin. We created a public database associated with an interactive map to visualize information for each accession. The oat collection was genotyped with 17,288 single-nucleotide polymorphism (SNP) loci to evaluate population structure and linkage disequilibrium (LD); to perform a genome-wide association study (GWAs) for heading date, a key character closely correlated with performance in this drought-prone area. Population genetic analysis using both structure and PCA distinguished two main groups composed of the red and white oats, respectively. The white oat group was further divided into two subgroups. LD decay was slower within white lines in linkage groups Mrg01, 02, 04, 12, 13, 15, 23, 33, whereas it was slower within red lines in Mrg03, 05, 06, 11, 21, 24, and 28. Association analysis showed several significant markers associated with heading date on linkage group Mrg13 in white oats and on Mrg01 and Mrg08 in red oats.

Population genomics of Mediterranean oat (A. sativa) reveals high genetic diversity and three loci for heading date

KEY MESSAGE

Genomic analysis of Mediterranean oats reveals high genetic diversity and three loci for adaptation to this environment. This information together with phenotyping and passport data, gathered in an interactive map, will be a vital resource for oat genetic improvement. During the twentieth century, oat landraces have increasingly been replaced by modern cultivars, resulting in loss of genetic diversity. However, landraces have considerable potential to improve disease and abiotic stress tolerance and may outperform cultivars under low input systems. In this work, we assembled a panel of 669 oat landraces from Mediterranean rim and 40 cultivated oat varieties and performed the first large-scale population genetic analysis of both red and white oat types of Mediterranean origin. We created a public database associated with an interactive map to visualize information for each accession. The oat collection was genotyped with 17,288 single-nucleotide polymorphism (SNP) loci to evaluate population structure and linkage disequilibrium (LD); to perform a genome-wide association study (GWAs) for heading date, a key character closely correlated with performance in this drought-prone area. Population genetic analysis using both structure and PCA distinguished two main groups composed of the red and white oats, respectively. The white oat group was further divided into two subgroups. LD decay was slower within white lines in linkage groups Mrg01, 02, 04, 12, 13, 15, 23, 33, whereas it was slower within red lines in Mrg03, 05, 06, 11, 21, 24, and 28. Association analysis showed several significant markers associated with heading date on linkage group Mrg13 in white oats and on Mrg01 and Mrg08 in red oats.

Overexpression of RgPAL family genes involved in phenolic biosynthesis promotes the replanting disease development in Rehmannia glutinosa

Rehmannia glutinosa production is affected by the replanting disease, which involves autotoxic harm mediated by specific endogenous allelochemicals in root exudates. Many phenolics that act as allelochemical agents are mostly phenylpropanoid products of secondary metabolism in plants. Phenylalanine ammonia-lyase (PAL) is the first enzyme that catalyses the deamination of l-phenylalanine for entrance into the phenylpropanoid pathway. PAL family genes have been isolated and functionally characterized in many plant species. However, PAL family genes involved in phenolic biosynthesis remain largely uncharacterized in R. glutinosa. Here, we identified and characterized four PAL family genes (RgPAL2 to RgPAL5) in the species whose sequences exhibited highly conserved domains of PALs according to in silico analysis, implying their potential function in phenolic biosynthesis. Overexpression of RgPALs in R. glutinosa enhanced phenolic production, verifying that RgPAL family genes participate in phenolic biosynthesis pathways. Moreover, we found that the release of several allelopathic phenolics from the roots of RgPAL-overexpressing transgenic R. glutinosa increased, implying that the RgPALs positively promote their release. Importantly, under continuous monoculture stress, we found that the RgPAL transgenic plants exhibited more significant autotoxic harm than did non-transgenic (WT) plants by activating the phenolics/phenylpropanoid pathway, indicating that RgPAL family genes function as positive regulators of the replanting disease development in R. glutinosa. This study revealed that RgPAL family genes are involved in the biosynthesis and release of several phenolics and positively control the replanting disease development in R. glutinosa, laying a foundation for further clarification of the molecular mechanisms underlying the disease formation.

Deep sampling and pooled amplicon sequencing reveals hidden genic variation in heterogeneous rye accessions

BACKGROUND

Loss of genetic variation negatively impacts breeding efforts and food security. Genebanks house over 7 million accessions representing vast allelic diversity that is a resource for sustainable breeding. Discovery of DNA variations is an important step in the efficient use of these resources. While technologies have improved and costs dropped, it remains impractical to consider resequencing millions of accessions. Candidate genes are known for most agronomic traits, providing a list of high priority targets. Heterogeneity in seed stocks means that multiple samples from an accession need to be evaluated to recover available alleles. To address this we developed a pooled amplicon sequencing approach and applied it to the out-crossing cereal rye (Secale cereale L.).

RESULTS

Using the amplicon sequencing approach 95 rye accessions of different improvement status and worldwide origin, each represented by a pooled sample comprising DNA of 96 individual plants, were evaluated for sequence variation in six candidate genes with significant functions on biotic and abiotic stress resistance, and seed quality. Seventy-four predicted deleterious variants were identified using multiple algorithms. Rare variants were recovered including those found only in a low percentage of seed.

CONCLUSIONS

We conclude that this approach provides a rapid and flexible method for evaluating stock heterogeneity, probing allele diversity, and recovering previously hidden variation. A large extent of within-population heterogeneity revealed in the study provides an important point for consideration during rye germplasm conservation and utilization efforts.

Genetic structure of a germplasm for hybrid breeding in rye (Secale cereale L.)

Rye (Secale cereale L.) responds strongly to changes in heterozygosity with hybrids portraying strong heterosis effect on all developmental and yielding characteristics. In order to achieve the highest potential heterosis effect parental lines must originate from genetically distinct gene pools. Here we report the first comprehensive SNP-based population study of an elite germplasm using fertilization control system for hybrid breeding in rye that is genetically different to the predominating P-type. In total 376 inbred lines from Nordic Seed Germany GmbH were genotyped for 4419 polymorphic SNPs. The aim of this study was to confirm and quantify the genetic separation of parental populations, unveil their genetic characteristics and investigate underlying population structures. Through a palette of complimenting analysis, we confirmed a strong genetic differentiation (FST = 0.332) of parental populations validating the germplasms suitability for hybrid breeding. These were, furthermore, found to diverge considerably in several features with the maternal population portraying a strong population structure characterized by a narrow genetic profile, small effective population size and high genome-wise linkage disequilibrium. We propose that the employed male-sterility system putatively constitutes a population determining parameter by influencing the rate of introducing novel genetic variation to the parental populations. Functional analysis of linkage blocks led to identification of a conserved segment on the distal 4RL chromosomal region annotated to the Rfp3 male-fertility restoration 'Pampa' type gene. Findings of our study emphasized the immediate value of comprehensive population studies on elite breeding germplasms as a pre-requisite for application of genomic-based breeding techniques, introgression of novel material and to support breeder decision-making.

Patterns of nucleotide diversity for different domains of centromeric histone H3 (CENH3) gene in Secale L

Rye (Secale) is among staple cereals along with other members of the Triticeae tribe: wheat and barley. The genus Secale includes perennial and annual, cross-pollinating and self-pollinating species, and they can be donors of valuable genes in wheat and rye breeding programs. Studies of the structure of the gene for centromeric histone H3 (CENH3), essential for centromere functions, are relevant to the breeding of agronomically important crops. We have investigated the nucleotide diversity of sequences of two variants of the rye CENH3 gene inside the N-terminal tail (NTT) and the conservative HFD (histone fold domain) domain in the genus Secale. The mean values of nucleotide diversity in the NTT and HFD of wild cross- and self-pollinating taxa are close in αCENH3: πtot = 0.0176–0.0090 and 0.0136–0. 0052, respectively. In the case of βCENH3, the mean values for NTT (πtot = 0.0168–0.0062) are lower than for HFD (πtot = 0.0259–0.084). The estimates of nucleotide and haplotype diversity per site for the CENH3 domains are considerably lower in taxa with narrow geographic ranges: S. cereale subsp. dighoricum and S. strictum subsp. kuprijanovii. Commercial breeding reduces the nucleotide sequence variability in αCENH3 and βCENH3. Cultivated rye varieties have π values within 0.0122–0.0014. The nucleotide and haplotype diversity values in αCENH3 and βCENH3 are close in S. sylvestre, which is believed to be the oldest rye species. The results of this study prove that the frequency of single nucleotide polymorphisms and nucleotide diversity of sequences in genes for CENH3 in Secale species are influenced by numerous factors, including reproduction habits, the geographic isolation of taxa, breeding, and the evolutionary age of species.

Volatile composition and classification of Lilium flower aroma types and identification, polymorphisms, and alternative splicing of their monoterpene synthase genes

Lily is a well-known ornamental plant with a diversity of fragrant types. Basic information on lily floral scent compounds has been obtained for only a few accessions, and little is known about Lilium aroma types, the terpene synthase genes that may play roles in the production of key volatiles, or the range of monoterpenes that these genes produce. In this study, 41 cultivars were analyzed for volatile emissions, and a total of 46 individual volatile compounds were identified, 16 for the first time in lilies. Lily accessions were classified into six groups according to the composition of major scent components: faint-scented, cool, fruity, musky, fruity-honey, and lily. Monoterpenes were one of the main groups of volatiles identified, and attention was focused on terpene synthase (TPS) genes, which encode enzymes that catalyze the last steps in monoterpene synthesis. Thirty-two candidate monoterpene synthase cDNAs were obtained from 66 lily cultivars, and 64 SNPs were identified. Two InDels were also shown to result from variable splicing, and sequence analysis suggested that different transcripts arose from the same gene. All identified nucleotide substitution sites were highly correlated with the amounts of myrcene emitted, and InDel site 230 was highly correlated with the emission of all major monoterpenoid components, especially (E)-β-ocimene. Heterologous expression of five cDNAs cloned from faint-scented and strong-scented lilies showed that their corresponding enzymes could convert geranyl diphosphate to (E)-β-ocimene, α-pinene, and limonene. The findings from this study provide a major resource for the assessment of lily scent volatiles and will be helpful in breeding of improved volatile components.

Association genetics studies on frost tolerance in wheat (Triticum aestivum L.) reveal new highly conserved amino acid substitutions in CBF-A3, CBF-A15, VRN3 and PPD1 genes

BACKGROUND

Understanding the genetic basis of frost tolerance (FT) in wheat (Triticum aestivum L.) is essential for preventing yield losses caused by frost due to cellular damage, dehydration and reduced metabolism. FT is a complex trait regulated by a number of genes and several gene families. Availability of the wheat genomic sequence opens new opportunities for exploring candidate genes diversity for FT. Therefore, the objectives of this study were to identity SNPs and insertion-deletion (indels) in genes known to be involved in frost tolerance and to perform association genetics analysis of respective SNPs and indels on FT.

RESULTS

Here we report on the sequence analysis of 19 candidate genes for FT in wheat assembled using the Chinese Spring IWGSC RefSeq v1.0. Out of these, the tandem duplicated C-repeat binding factors (CBF), i.e. CBF-A3, CBF-A5, CBF-A10, CBF-A13, CBF-A14, CBF-A15, CBF-A18, the vernalisation response gene VRN-A1, VRN-B3, the photoperiod response genes PPD-B1 and PPD-D1 revealed association to FT in 235 wheat cultivars. Within six genes (CBF-A3, CBF-A15, VRN-A1, VRN-B3, PPD-B1 and PPD-D1) amino acid (AA) substitutions in important protein domains were identified. The amino acid substitution effect in VRN-A1 on FT was confirmed and new AA substitutions in CBF-A3, CBF-A15, VRN-B3, PPD-B1 and PPD-D1 located at highly conserved sites were detected. Since these results rely on phenotypic data obtained at five locations in 2 years, detection of significant associations of FT to AA changes in CBF-A3, CBF-A15, VRN-A1, VRN-B3, PPD-B1 and PPD-D1 may be exploited in marker assisted breeding for frost tolerance in winter wheat.

CONCLUSIONS

A set of 65 primer pairs for the genes mentioned above from a previous study was BLASTed against the IWGSC RefSeq resulting in the identification of 39 primer combinations covering the full length of 19 genes. This work demonstrates the usefulness of the IWGSC RefSeq in specific primer development for highly conserved gene families in hexaploid wheat and, that a candidate gene association genetics approach based on the sequence data is an efficient tool to identify new alleles of genes important for the response to abiotic stress in wheat.

Investigating Drought Tolerance in Chickpea Using Genome-Wide Association Mapping and Genomic Selection Based on Whole-Genome Resequencing Data

Drought tolerance is a complex trait that involves numerous genes. Identifying key causal genes or linked molecular markers can facilitate the fast development of drought tolerant varieties. Using a whole-genome resequencing approach, we sequenced 132 chickpea varieties and advanced breeding lines and found more than 144,000 single nucleotide polymorphisms (SNPs). We measured 13 yield and yield-related traits in three drought-prone environments of Western Australia. The genotypic effects were significant for all traits, and many traits showed highly significant correlations, ranging from 0.83 between grain yield and biomass to -0.67 between seed weight and seed emergence rate. To identify candidate genes, the SNP and trait data were incorporated into the SUPER genome-wide association study (GWAS) model, a modified version of the linear mixed model. We found that several SNPs from auxin-related genes, including auxin efflux carrier protein (PIN3), p-glycoprotein, and nodulin MtN21/EamA-like transporter, were significantly associated with yield and yield-related traits under drought-prone environments. We identified four genetic regions containing SNPs significantly associated with several different traits, which was an indication of pleiotropic effects. We also investigated the possibility of incorporating the GWAS results into a genomic selection (GS) model, which is another approach to deal with complex traits. Compared to using all SNPs, application of the GS model using subsets of SNPs significantly associated with the traits under investigation increased the prediction accuracies of three yield and yield-related traits by more than twofold. This has important implication for implementing GS in plant breeding programs.

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.

Towards a whole-genome sequence for rye (Secale cereale L.)

We report on a whole-genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole-genome shotgun sequencing of the 7.9-Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low-copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high-density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high-density genetic map and establishing a synteny-based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome-wide map of selection signals differentiating the divergent gene pools. This rye whole-genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome-based breeding in rye.

Association Mapping in Turkish Olive Cultivars Revealed Significant Markers Related to Some Important Agronomic Traits

Olive (Olea europaea L.) is one of the most important fruit trees especially in the Mediterranean countries due to high consumption of table olive and olive oil. In olive breeding, the phenotypic traits associated to fruit are the key factors that determine productivity. Association mapping has been used in some tree species and a lot of crop plant species, and here, we perform an initial effort to detect marker-trait associations in olive tree. In the current study, a total of 96 olive genotypes, including both oil and table olive genotypes from Turkish Olive GenBank Resources, were used to examine marker-trait associations. For olive genotyping, SNP, AFLP, and SSR marker data were selected from previously published study and association analysis was performed between these markers and 5 yield-related traits. Three different approaches were used to check for false-positive results in association tests, and association results obtained from these models were compared. Using the model utilizing both population structure and relative kinship, eleven associations were significant with FDR ≤ 0.05. The largest number of significant associations was detected for fruit weight and stone weight. Our results suggested that association mapping could be an effective approach for identifying marker-trait associations in olive genotypes, without the development of mapping populations. This study shows for the first time the use of association mapping for identifying molecular markers linked to important traits in olive tree.

Geographical distribution of genetic diversity in Secale landrace and wild accessions

BACKGROUND

Rye, Secale cereale L., has historically been a crop of major importance and is still a key cereal in many parts of Europe. Single populations of cultivated rye have been shown to capture a large proportion of the genetic diversity present in the species, but the distribution of genetic diversity in subspecies and across geographical areas is largely unknown. Here we explore the structure of genetic diversity in landrace rye and relate it to that of wild and feral relatives.

RESULTS

A total of 567 SNPs were analysed in 434 individuals from 76 accessions of wild, feral and cultivated rye. Genetic diversity was highest in cultivated rye, slightly lower in feral rye taxa and significantly lower in the wild S. strictum Presl. and S. africanum Stapf. Evaluation of effects from ascertainment bias suggests underestimation of diversity primarily in S. strictum and S. africanum. Levels of ascertainment bias, STRUCTURE and principal component analyses all supported the proposed classification of S. africanum and S. strictum as a separate species from S. cereale. S. afghanicum (Vav.) Roshev, S. ancestrale Zhuk., S. dighoricum (Vav.) Roshev, S. segetale (Zhuk.) Roshev and S. vavilovii Grossh. seemed, in contrast, to share the same gene pool as S. cereale and their genetic clustering was more dependent on geographical origin than taxonomic classification. S. vavilovii was found to be the most likely wild ancestor of cultivated rye. Among cultivated rye landraces from Europe, Asia and North Africa five geographically discrete genetic clusters were identified. These had only limited overlap with major agro-climatic zones. Slash-and-burn rye from the Finnmark area in Scandinavia formed a distinct cluster with little similarity to other landrace ryes. Regional studies of Northern and South-West Europe demonstrate different genetic distribution patterns as a result of varying cultivation intensity.

CONCLUSIONS

With the exception of S. strictum and S. africanum different rye taxa share the majority of the genetic variation. Due to the vast sharing of genetic diversity within the S. cereale clade, ascertainment bias seems to be a lesser problem in rye than in predominantly selfing species. By exploiting within accession diversity geographic structure can be shown on a much finer scale than previously reported.

A horizontally transferred nuclear gene is associated with microhabitat variation in a natural plant population

Horizontal gene transfer involves the non-sexual interspecific transmission of genetic material. Even if they are initially functional, horizontally transferred genes are expected to deteriorate into non-expressed pseudogenes, unless they become adaptively relevant in the recipient organism. However, little is known about the distributions of natural transgenes within wild species or the adaptive significance of natural transgenes within wild populations. Here, we examine the distribution of a natural plant-to-plant nuclear transgene in relation to environmental variation within a wild population. Festuca ovina is polymorphic for an extra (second) expressed copy of the nuclear gene (PgiC) encoding cytosolic phosphoglucose isomerase, with the extra PgiC locus having been acquired horizontally from the distantly related grass genus Poa. We investigated variation at PgiC in samples of F. ovina from a fine-scale, repeating patchwork of grassland microhabitats, replicated within spatially separated sites. Even after accounting for spatial effects, the distributions of F. ovina individuals carrying the additional PgiC locus, and one of the enzyme products encoded by the locus, are significantly associated with fine-scale habitat variation. Our results suggest that the PgiC transgene contributes, together with the unlinked 'native' PgiC locus, to local adaptation to a fine-scale mosaic of edaphic and biotic grassland microhabitats.

An Efficient Approach for the Development of Locus Specific Primers in Bread Wheat (Triticum aestivum L.) and Its Application to Re-Sequencing of Genes Involved in Frost Tolerance

Recent declines in costs accelerated sequencing of many species with large genomes, including hexaploid wheat (Triticum aestivum L.). Although the draft sequence of bread wheat is known, it is still one of the major challenges to developlocus specific primers suitable to be used in marker assisted selection procedures, due to the high homology of the three genomes. In this study we describe an efficient approach for the development of locus specific primers comprising four steps, i.e. (i) identification of genomic and coding sequences (CDS) of candidate genes, (ii) intron- and exon-structure reconstruction, (iii) identification of wheat A, B and D sub-genome sequences and primer development based on sequence differences between the three sub-genomes, and (iv); testing of primers for functionality, correct size and localisation. This approach was applied to single, low and high copy genes involved in frost tolerance in wheat. In summary for 27 of these genes for which sequences were derived from Triticum aestivum, Triticum monococcum and Hordeum vulgare, a set of 119 primer pairs was developed and after testing on Nulli-tetrasomic (NT) lines, a set of 65 primer pairs (54.6%), corresponding to 19 candidate genes, turned out to be specific. Out of these a set of 35 fragments was selected for validation via Sanger's amplicon re-sequencing. All fragments, with the exception of one, could be assigned to the original reference sequence. The approach presented here showed a much higher specificity in primer development in comparison to techniques used so far in bread wheat and can be applied to other polyploid species with a known draft sequence.

Employing genome-wide SNP discovery and genotyping strategy to extrapolate the natural allelic diversity and domestication patterns in chickpea

The genome-wide discovery and high-throughput genotyping of SNPs in chickpea natural germplasm lines is indispensable to extrapolate their natural allelic diversity, domestication, and linkage disequilibrium (LD) patterns leading to the genetic enhancement of this vital legume crop. We discovered 44,844 high-quality SNPs by sequencing of 93 diverse cultivated desi, kabuli, and wild chickpea accessions using reference genome- and de novo-based GBS (genotyping-by-sequencing) assays that were physically mapped across eight chromosomes of desi and kabuli. Of these, 22,542 SNPs were structurally annotated in different coding and non-coding sequence components of genes. Genes with 3296 non-synonymous and 269 regulatory SNPs could functionally differentiate accessions based on their contrasting agronomic traits. A high experimental validation success rate (92%) and reproducibility (100%) along with strong sensitivity (93-96%) and specificity (99%) of GBS-based SNPs was observed. This infers the robustness of GBS as a high-throughput assay for rapid large-scale mining and genotyping of genome-wide SNPs in chickpea with sub-optimal use of resources. With 23,798 genome-wide SNPs, a relatively high intra-specific polymorphic potential (49.5%) and broader molecular diversity (13-89%)/functional allelic diversity (18-77%) was apparent among 93 chickpea accessions, suggesting their tremendous applicability in rapid selection of desirable diverse accessions/inter-specific hybrids in chickpea crossbred varietal improvement program. The genome-wide SNPs revealed complex admixed domestication pattern, extensive LD estimates (0.54-0.68) and extended LD decay (400-500 kb) in a structured population inclusive of 93 accessions. These findings reflect the utility of our identified SNPs for subsequent genome-wide association study (GWAS) and selective sweep-based domestication trait dissection analysis to identify potential genomic loci (gene-associated targets) specifically regulating important complex quantitative agronomic traits in chickpea. The numerous informative genome-wide SNPs, natural allelic diversity-led domestication pattern, and LD-based information generated in our study have got multidimensional applicability with respect to chickpea genomics-assisted breeding.

Development of genetic markers in Eucalyptus species by target enrichment and exome sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.

Assessment of genetic variation within a global collection of lentil (Lens culinaris Medik.) cultivars and landraces using SNP markers

Background

Lentil is a self-pollinated annual diploid (2n = 2× = 14) crop with a restricted history of genetic improvement through breeding, particularly when compared to cereal crops. This limited breeding has probably contributed to the narrow genetic base of local cultivars, and a corresponding potential to continue yield increases and stability. Therefore, knowledge of genetic variation and relationships between populations is important for understanding of available genetic variability and its potential for use in breeding programs. Single nucleotide polymorphism (SNP) markers provide a method for rapid automated genotyping and subsequent data analysis over large numbers of samples, allowing assessment of genetic relationships between genotypes.

Results

In order to investigate levels of genetic diversity within lentil germplasm, 505 cultivars and landraces were genotyped with 384 genome-wide distributed SNP markers, of which 266 (69.2%) obtained successful amplification and detected polymorphisms. Gene diversity and PIC values varied between 0.108-0.5 and 0.102-0.375, with averages of 0.419 and 0.328, respectively. On the basis of clarity and interest to lentil breeders, the genetic structure of the germplasm collection was analysed separately for cultivars and landraces. A neighbour-joining (NJ) dendrogram was constructed for commercial cultivars, in which lentil cultivars were sorted into three major groups (G-I, G-II and G-III). These results were further supported by principal coordinate analysis (PCoA) and STRUCTURE, from which three clear clusters were defined based on differences in geographical location. In the case of landraces, a weak correlation between geographical origin and genetic relationships was observed. The landraces from the Mediterranean region, predominantly Greece and Turkey, revealed very high levels of genetic diversity.

Conclusions

Lentil cultivars revealed clear clustering based on geographical origin, but much more limited correlation between geographic origin and genetic diversity was observed for landraces. These results suggest that selection of divergent parental genotypes for breeding should be made actively on the basis of systematic assessment of genetic distance between genotypes, rather than passively based on geographical distance.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0150-3) contains supplementary material, which is available to authorized users.

Genome-wide characterization of genetic diversity and population structure in Secale

BACKGROUND

Numerous rye accessions are stored in ex situ genebanks worldwide. Little is known about the extent of genetic diversity contained in any of them and its relation to contemporary varieties, since to date rye genetic diversity studies had a very limited scope, analyzing few loci and/ or few accessions. Development of high throughput genotyping methods for rye opened the possibility for genome wide characterizations of large accessions sets. In this study we used 1054 Diversity Array Technology (DArT) markers with defined chromosomal location to characterize genetic diversity and population structure in a collection of 379 rye accessions including wild species, landraces, cultivated materials, historical and contemporary rye varieties.

RESULTS

Average genetic similarity (GS) coefficients and average polymorphic information content (PIC) values varied among chromosomes. Comparison of chromosome specific average GS within and between germplasm sub-groups indicated regions of chromosomes 1R and 4R as being targeted by selection in current breeding programs. Bayesian clustering, principal coordinate analysis and Neighbor Joining clustering demonstrated that source and improvement status contributed significantly to the structure observed in the analyzed set of Secale germplasm. We revealed a relatively limited diversity in improved rye accessions, both historical and contemporary, as well as lack of correlation between clustering of improved accessions and geographic origin, suggesting common genetic background of rye accessions from diverse geographic regions and extensive germplasm exchange. Moreover, contemporary varieties were distinct from the remaining accessions.

CONCLUSIONS

Our results point to an influence of reproduction methods on the observed diversity patterns and indicate potential of ex situ collections for broadening the genetic diversity in rye breeding programs. Obtained data show that DArT markers provide a realistic picture of the genetic diversity and population structure present in the collection of 379 rye accessions and are an effective platform for rye germplasm characterization and association mapping studies.

Association analysis of single nucleotide polymorphisms in candidate genes with root traits in maize (Zea mays L.) seedlings

Several genes involved in maize root development have been isolated. Identification of SNPs associated with root traits would enable the selection of maize lines with better root architecture that might help to improve N uptake, and consequently plant growth particularly under N deficient conditions. In the present study, an association study (AS) panel consisting of 74 maize inbred lines was screened for seedling root traits in 6, 10, and 14-day-old seedlings. Allele re-sequencing of candidate root genes Rtcl, Rth3, Rum1, and Rul1 was also carried out in the same AS panel lines. All four candidate genes displayed different levels of nucleotide diversity, haplotype diversity and linkage disequilibrium. Gene based association analyses were carried out between individual polymorphisms in candidate genes, and root traits measured in 6, 10, and 14-day-old maize seedlings. Association analyses revealed several polymorphisms within the Rtcl, Rth3, Rum1, and Rul1 genes associated with seedling root traits. Several nucleotide polymorphisms in Rtcl, Rth3, Rum1, and Rul1 were significantly (P<0.05) associated with seedling root traits in maize suggesting that all four tested genes are involved in the maize root development. Thus considerable allelic variation present in these root genes can be exploited for improving maize root characteristics.

Expression and molecular evolution of two DREB1 genes in black poplar (Populus nigra)

Environmental stresses such as low temperature, drought, and high salinity significantly affect plant growth and yield. As selective forces, these adverse factors play essential roles in shaping phenotypic variation in plant populations. Black poplar (Populus nigra) is an economically and ecologically important forest tree species with widely distributed populations and is thus suitable for experiments detecting evolutionary footprints left by stress. Here, we performed expression and evolutionary analysis of two duplicated DREB A1-subgroup (DREB1) genes, PnDREB68 and PnDREB69, encoding transcription factors that are involved in stress responses. The two genes showed partially overlapping but distinct expression patterns in response to stresses. These genes were strongly and rapidly induced by cold stress in leaves, stems, and roots. In leaf tissue, dehydration stress induced the expression of PnDREB68 but not PnDREB69. PnDREB69 displayed more rapid responses and longer expression durations than PnDREB68 under salt and ABA stress, respectively. Based on single nucleotide polymorphism (SNP) analysis, we found significant population genetic differentiation, with a greater F_ST value (0.09189) for PnDREB69 than for PnDREB68 (0.07743). Nucleotide diversity analysis revealed a two-fold higher π_T for PnDREB68 than for PnDREB69 (0.00563 vs. 0.00243), reflecting strong purifying selection acting on the former. The results suggest that positive selection acted on PnDREB69, as evidenced by neutral testing using Tajima’s D statistic. The distinct selective forces to which each of the genes was subjected may be associated with expression divergence. Linkage disequilibrium (LD) was low for the sequenced region, with a higher level for PnDREB68 than for PnDREB69. Additionally, analysis of the relationship among carbon isotope ratios, SNP classes and gene expression, together with motif and domain analysis, suggested that 14 polymorphisms within the two genes may be candidates for an association study of important traits such as water use efficiency/drought tolerance in black poplar.

A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the 'true citrus fruit trees' group (Citrinae, Rutaceae) and the origin of cultivated species

BACKGROUND AND AIMS

Despite differences in morphology, the genera representing 'true citrus fruit trees' are sexually compatible, and their phylogenetic relationships remain unclear. Most of the important commercial 'species' of Citrus are believed to be of interspecific origin. By studying polymorphisms of 27 nuclear genes, the average molecular differentiation between species was estimated and some phylogenetic relationships between 'true citrus fruit trees' were clarified.

METHODS

Sanger sequencing of PCR-amplified fragments from 18 genes involved in metabolite biosynthesis pathways and nine putative genes for salt tolerance was performed for 45 genotypes of Citrus and relatives of Citrus to mine single nucleotide polymorphisms (SNPs) and indel polymorphisms. Fifty nuclear simple sequence repeats (SSRs) were also analysed.

KEY RESULTS

A total of 16 238 kb of DNA was sequenced for each genotype, and 1097 single nucleotide polymorphisms (SNPs) and 50 indels were identified. These polymorphisms were more valuable than SSRs for inter-taxon differentiation. Nuclear phylogenetic analysis revealed that Citrus reticulata and Fortunella form a cluster that is differentiated from the clade that includes three other basic taxa of cultivated citrus (C. maxima, C. medica and C. micrantha). These results confirm the taxonomic subdivision between the subgenera Metacitrus and Archicitrus. A few genes displayed positive selection patterns within or between species, but most of them displayed neutral patterns. The phylogenetic inheritance patterns of the analysed genes were inferred for commercial Citrus spp.

CONCLUSIONS

Numerous molecular polymorphisms (SNPs and indels), which are potentially useful for the analysis of interspecific genetic structures, have been identified. The nuclear phylogenetic network for Citrus and its sexually compatible relatives was consistent with the geographical origins of these genera. The positive selection observed for a few genes will help further works to analyse the molecular basis of the variability of the associated traits. This study presents new insights into the origin of C. sinensis.

Nucleotide diversity and linkage disequilibrium in antioxidant genes of Brachypodium distachyon

Brachypodium distachyon (Brachypodium) is a powerful model system for studying cereal, bioenergy, forage, and turf grasses. Nucleotide diversity (π) and linkage disequilibrium (LD) in candidate genes involved in the antioxidative pathways in this species are not known. The average π for CAT encoding catalase, GPX encoding glutathione peroxidase, DHAR encoding dehydroascorbate reductase, MDHAR encoding monodehydroascorbate reductase, and APX ecoding ascorbate peroxidase was 0.0027 among 19 accessions contrasting for drought tolerance. The highest value of π was found in APX (0.0046) and the lowest π was in MDHAR (0.0006). The average single nucleotide polymorphism (SNP) frequency across these five genes was one SNP per 131 bp between two randomly sampled sequences for the five genes in the sequence length ranging from 1,447 bp to 1,701 bp. The LD decay was slow and extended to a distance of more than 1.2kb for all genes. The neighbor-joining tree analyses of DHAR, MDHAR, and CAT generally separated accessions differing in drought tolerance. The results indicate a putative role of these candidate genes in increasing general fitness of Brachypodium.

High-throughput and parallel SNP discovery in selected candidate genes in Eucalyptus camaldulensis using Illumina NGS platform

Next generation sequencing (NGS) technologies have revolutionized the pace and scale of genomics- and transcriptomics-based SNP discovery across different plant and animal species. Herein, 72-base paired-end Illumina sequencing was employed for high-throughput, parallel and large-scale SNP discovery in 41 growth-related candidate genes in Eucalyptus camaldulensis. Approximately 100 kb of genome from 96 individuals was amplified and sequenced using a hierarchical DNA/PCR pooling strategy and assembled over corresponding E. grandis reference. A total of 1191 SNPs (minimum 5% other allele frequency) were identified with an average frequency of 1 SNP/83.9 bp, whereas in exons and introns, it was 1 SNP/108.4 bp and 1 SNP/65.6 bp, respectively. A total of 75 insertions and 89 deletions were detected of which approximately 15% were exonic. Transitions (Tr) were in excess than transversions (Tv) (Tr/Tv: 1.89), but exceeded in exons (Tr/Tv: 2.73). In exons, synonymous SNPs (Ka) prevailed over the non-synonymous SNPs (Ks; average Ka/Ks ratio: 0.72, range: 0-3.00 across genes). Many of the exonic SNPs/indels had potential to change amino acid sequence of respective genes. Transcription factors appeared more conserved, whereas enzyme coding genes appeared under relaxed control. Further, 541 SNPs were classified into 196 'equal frequency' (EF) blocks with almost similar minor allele frequencies to facilitate selection of one tag-SNP/EF-block. There were 241 (approximately 20%) 'zero-SNP' blocks with absence of SNPs in surrounding ±60 bp windows. The data thus indicated enormous extant and unexplored diversity in E. camaldulensis in the studied genes with potential applications for marker-trait associations.

Genome-wide evaluation of genetic diversity and linkage disequilibrium in winter and spring triticale (x Triticosecale Wittmack)

Background

Recent advances in genotyping with high-density markers nowadays enable genome-wide genomic analyses in crops. A detailed characterisation of the population structure and linkage disequilibrium (LD) is essential for the application of genomic approaches and consequently for knowledge-based breeding. In this study we used the triticale-specific DArT array to analyze population structure, genetic diversity, and LD in a worldwide set of 161 winter and spring triticale lines.

Results

The principal coordinate analysis revealed that the first principal coordinate divides the triticale population into two clusters according to their growth habit. The density distributions of the first ten principal coordinates revealed that several show a distribution indicative of population structure. In addition, we observed relatedness within growth habits which was higher among the spring types than among the winter types. The genome-wide analysis of polymorphic information content (PIC) showed that the PIC is variable among and along chromosomes and that especially the R genome of spring types possesses a reduced genetic diversity. We also found that several chromosomes showed regions of high genetic distance between the two growth habits, indicative of divergent selection. Regarding linkage disequilibrium, the A and B genomes showed a similar LD of 0.24 for closely linked markers and a decay within approximately 12 cM. LD in the R genome was lower with 0.19 and decayed within a shorter map distance of approximately 5 cM. The extent of LD was generally higher for the spring types compared to the winter types. In addition, we observed strong variability of LD along the chromosomes.

Conclusions

Our results confirm winter and spring growth habit are the major contributors to population structure in triticale, and a family structure exists in both growth types. The specific patterns of genetic diversity observed within these types, such as the low diversity on some rye chromosomes of spring habits, provide a basis for targeted broadening of the available breeding germplasm. In addition, the genome-wide analysis of the extent and the pattern of LD will assist scientists and breeders alike in the implementation and the interpretation of association mapping in triticale.

Species discrimination, population structure and linkage disequilibrium in Eucalyptus camaldulensis and Eucalyptus tereticornis using SSR markers

Eucalyptus camaldulensis and E. tereticornis are closely related species commonly cultivated for pulp wood in many tropical countries including India. Understanding the genetic structure and linkage disequilibrium (LD) existing in these species is essential for the improvement of industrially important traits. Our goal was to evaluate the use of simple sequence repeat (SSR) loci for species discrimination, population structure and LD analysis in these species. Investigations were carried out with the most common alleles in 93 accessions belonging to these two species using 62 SSR markers through cross amplification. The polymorphic information content (PIC) ranged from 0.44 to 0.93 and 0.36 to 0.93 in E. camaldulensis and E. tereticornis respectively. A clear delineation between the two species was evident based on the analysis of population structure and species-specific alleles. Significant genotypic LD was found in E. camaldulensis, wherein out of 135 significant pairs, 17 pairs showed r²≥0.1. Similarly, in E. tereticornis, out of 136 significant pairs, 18 pairs showed r²≥0.1. The extent of LD decayed rapidly showing the significance of association analyses in eucalypts with higher resolution markers. The availability of whole genome sequence for E. grandis and the synteny and co-linearity in the genome of eucalypts, will allow genome-wide genotyping using microsatellites or single nucleotide polymorphims.

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.

Association analysis of frost tolerance in rye using candidate genes and phenotypic data from controlled, semi-controlled, and field phenotyping platforms

BACKGROUND

Frost is an important abiotic stress that limits cereal production in the temperate zone. As the most frost tolerant small grain cereal, rye (Secale cereale L.) is an ideal cereal model for investigating the genetic basis of frost tolerance (FT), a complex trait with polygenic inheritance. Using 201 genotypes from five Eastern and Middle European winter rye populations, this study reports a multi-platform candidate gene-based association analysis in rye using 161 single nucleotide polymorphisms (SNPs) and nine insertion-deletion (Indel) polymorphisms previously identified from twelve candidate genes with a putative role in the frost responsive network.

RESULTS

Phenotypic data analyses of FT in three different phenotyping platforms, controlled, semi-controlled and field, revealed significant genetic variations in the plant material under study. Statistically significant (P < 0.05) associations between FT and SNPs/haplotypes of candidate genes were identified. Two SNPs in ScCbf15 and one in ScCbf12, all leading to amino acid exchanges, were significantly associated with FT over all three phenotyping platforms. Distribution of SNP effect sizes expressed as percentage of the genetic variance explained by individual SNPs was highly skewed towards zero with a few SNPs obtaining large effects. Two-way epistasis was found between 14 pairs of candidate genes. Relatively low to medium empirical correlations of SNP-FT associations were observed across the three platforms underlining the need for multi-level experimentation for dissecting complex associations between genotypes and FT in rye.

CONCLUSIONS

Candidate gene based-association studies are a powerful tool for investigating the genetic basis of FT in rye. Results of this study support the findings of bi-parental linkage mapping and expression studies that the Cbf gene family plays an essential role in FT.