Rapid reorganization of resistance gene homologues in cereal genomes

Resistance strategies for defense against Albugo candida causing white rust disease

Albugo candida, the causal organism of white rust, is an oomycete obligate pathogen infecting crops of Brassicaceae family occurred on aerial part, including vegetable and oilseed crops at all growth stages. The disease expression is characterized by local infection appearing on the abaxial region developing white or creamy yellow blister (sori) on leaves and systemic infections cause hypertrophy and hyperplasia leading to stag-head of reproductive organ. To overcome this problem, several disease management strategies like fungicide treatments were used in the field and disease-resistant varieties have also been developed using conventional and molecular breeding. Due to high variability among A. candida isolates, there is no single approach available to understand the diverse spectrum of disease symptoms. In absence of resistance sources against pathogen, repetitive cultivation of genetically-similar varieties locally tends to attract oomycete pathogen causing heavy yield losses. In the present review, a deep insight into the underlying role of the non-host resistance (NHR) defence mechanism available in plants, and the strategies to exploit available gene pools from plant species that are non-host to A. candida could serve as novel sources of resistance. This work summaries the current knowledge pertaining to the resistance sources available in non-host germ plasm, the understanding of defence mechanisms and the advance strategies covers molecular, biochemical and nature-based solutions in protecting Brassica crops from white rust disease.

Novel Fusarium wilt resistance genes uncovered in natural and cultivated strawberry populations are found on three non-homoeologous chromosomes

Several Fusarium wilt resistance genes were discovered, genetically and physically mapped, and rapidly deployed via marker-assisted selection to develop cultivars resistant to Fusarium oxysporum f. sp. fragariae, a devastating soil-borne pathogen of strawberry. Fusarium wilt, a soilborne disease caused by Fusarium oxysporum f. sp. fragariae, poses a significant threat to strawberry (Fragaria [Formula: see text] ananassa) production in many parts of the world. This pathogen causes wilting, collapse, and death in susceptible genotypes. We previously identified a dominant gene (FW1) on chromosome 2B that confers resistance to race 1 of the pathogen, and hypothesized that gene-for-gene resistance to Fusarium wilt was widespread in strawberry. To explore this, a genetically diverse collection of heirloom and modern cultivars and octoploid ecotypes were screened for resistance to Fusarium wilt races 1 and 2. Here, we show that resistance to both races is widespread in natural and domesticated populations and that resistance to race 1 is conferred by partially to completely dominant alleles among loci (FW1, FW2, FW3, FW4, and FW5) found on three non-homoeologous chromosomes (1A, 2B, and 6B). The underlying genes have not yet been cloned and functionally characterized; however, plausible candidates were identified that encode pattern recognition receptors or other proteins known to confer gene-for-gene resistance in plants. High-throughput genotyping assays for SNPs in linkage disequilibrium with FW1-FW5 were developed to facilitate marker-assisted selection and accelerate the development of race 1 resistant cultivars. This study laid the foundation for identifying the genes encoded by FW1-FW5, in addition to exploring the genetics of resistance to race 2 and other races of the pathogen, as a precaution to averting a Fusarium wilt pandemic.

A high-density linkage map of finger millet provides QTL for blast resistance and other agronomic traits

Finger millet [Eleusine coracana (L.) Gaertn.] is a critical subsistence crop in eastern Africa and southern Asia but has few genomic resources and modern breeding programs. To aid in the understanding of finger millet genomic organization and genes underlying disease resistance and agronomically important traits, we generated a F_2:3 population from a cross between E. coracana (L.) Gaertn. subsp. coracana accession ACC 100007 and E. coracana (L.) Gaertn. subsp. africana , accession GBK 030647. Phenotypic data on morphology, yield, and blast (Magnaporthe oryzae) resistance traits were taken on a subset of the F_2:3 population in a Kenyan field trial. The F_2:3 population was genotyped via genotyping-by-sequencing (GBS) and the UGbS-Flex pipeline was used for sequence alignment, nucleotide polymorphism calling, and genetic map construction. An 18-linkage-group genetic map consisting of 5,422 markers was generated that enabled comparative genomic analyses with rice (Oryza sativa L.), foxtail millet [Setaria italica (L.) P. Beauv.], and sorghum [Sorghum bicolor (L.) Moench]. Notably, we identified conserved acrocentric homoeologous chromosomes (4A and 4B in finger millet) across all species. Significant quantitative trait loci (QTL) were discovered for flowering date, plant height, panicle number, and blast incidence and severity. Sixteen putative candidate genes that may underlie trait variation were identified. Seven LEUCINE-RICH REPEAT-CONTAINING PROTEIN genes, with homology to nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance proteins, were found on three chromosomes under blast resistance QTL. This high-marker-density genetic map provides an important tool for plant breeding programs and identifies genomic regions and genes of critical interest for agronomic traits and blast resistance.

Plant pan-genomes are the new reference

Recent years have seen a surge in plant genome sequencing projects and the comparison of multiple related individuals. The high degree of genomic variation observed led to the realization that single reference genomes do not represent the diversity within a species, and led to the expansion of the pan-genome concept. Pan-genomes represent the genomic diversity of a species and includes core genes, found in all individuals, as well as variable genes, which are absent in some individuals. Variable gene annotations often show similarities across plant species, with genes for biotic and abiotic stress commonly enriched within variable gene groups. Here we review the growth of pan-genomics in plants, explore the origins of gene presence and absence variation, and show how pan-genomes can support plant breeding and evolution studies.

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

•

Species-wide NLR diversity is high but not unlimited

•

A large fraction of NLR diversity is recovered with 40–50 accessions

•

Presence/absence variation in NLRs is widespread, resulting in a mosaic population

•

A high diversity of NLR-integrated domains favor known virulence targets

A Species-Wide Inventory of NLR Genes and Alleles in Arabidopsis thaliana

Infectious disease is both a major force of selection in nature and a prime cause of yield loss in agriculture. In plants, disease resistance is often conferred by nucleotide-binding leucine-rich repeat (NLR) proteins, intracellular immune receptors that recognize pathogen proteins and their effects on the host. Consistent with extensive balancing and positive selection, NLRs are encoded by one of the most variable gene families in plants, but the true extent of intraspecific NLR diversity has been unclear. Here, we define a nearly complete species-wide pan-NLRome in Arabidopsis thaliana based on sequence enrichment and long-read sequencing. The pan-NLRome largely saturates with approximately 40 well-chosen wild strains, with half of the pan-NLRome being present in most accessions. We chart NLR architectural diversity, identify new architectures, and quantify selective forces that act on specific NLRs and NLR domains. Our study provides a blueprint for defining pan-NLRomes.

BvcZR3 and BvHs1pro-1 Genes Pyramiding Enhanced Beet Cyst Nematode (Heterodera schachtii Schm.) Resistance in Oilseed Rape (Brassica napus L.)

Beet cyst nematode (Heterodera schachtii Schm.) is one of the most damaging pests in sugar beet growing areas around the world. The Hs1^pro-1 and cZR3 genes confer resistance to the beet cyst nematode, and both were cloned from sugar beet translocation line (A906001). The translocation line carried the locus from B. procumbens chromosome 1 including Hs1^pro-1 gene and resistance gene analogs (RGA), which confer resistance to Heterodera schachtii. In this research, BvHs1^pro-1 and BvcZR3 genes were transferred into oilseed rape to obtain different transgenic lines by A. tumefaciens mediated transformation method. The cZR3Hs1^pro-1 gene was pyramided into the same plants by crossing homozygous cZR3 and Hs1^pro-1 plants to identify the function and interaction of cZR3 and Hs1^pro-1 genes. In vitro and in vivo cyst nematode resistance tests showed that cZR3 and Hs1^pro-1 plants could be infested by beet cyst nematode (BCN) juveniles, however a large fraction of penetrated nematode juveniles was not able to develop normally and stagnated in roots of transgenic plants, consequently resulting in a significant reduction in the number of developed nematode females. A higher efficiency in inhibition of nematode females was observed in plants expressing pyramiding genes than in those only expressing a single gene. Molecular analysis demonstrated that BvHs1^pro-1 and BvcZR3 gene expressions in oilseed rape constitutively activated transcription of plant-defense related genes such as NPR1 (non-expresser of PR1), SGT1b (enhanced disease resistance 1) and RAR1 (suppressor of the G2 allele of skp1). Transcript of NPR1 gene in transgenic cZR3 and Hs1^pro-1 plants were slightly up-regulated, while its expression was considerably enhanced in cZR3Hs1^pro-1 gene pyramiding plants. The expression of EDS1 gene did not change significantly among transgenic cZR3, Hs1^pro-1 and cZR3Hs1^pro-1 gene pyramiding plants and wild type. The expression of SGT1b gene was slightly up-regulated in transgenic cZR3 and Hs1^pro-1 plants compared with the wild type, however, its expression was not changed in cZR3Hs1^pro-1 gene pyramiding plant and had no interaction effect. RAR1 gene expression was significantly up-regulated in transgenic cZR3 and cZR3Hs1^pro-1 genes pyramiding plants, but almost no expression was found in Hs1^pro-1 transgenic plants. These results show that nematode resistance genes from sugar beet were functional in oilseed rape and conferred BCN resistance by activation of a CC-NBS-LRR R gene mediated resistance response. The gene pyramiding had enhanced resistance, thus offering a novel approach for the BCN control by preventing the propagation of BCN in oilseed rape. The transgenic oilseed rape could be used as a trap crop to offer an alternative method for beet cyst nematode control.

Prospects of pan-genomics in barley

The concept of a pan-genome refers to intraspecific diversity in genome content and structure, encompassing both genes and intergenic space. Pan-genomic studies employ a combination of de novo sequence assembly and reference-based alignment to discover and genotype structural variants. The large size and complex structure of Triticeae genomes were for a long time an obstacle for genomic research in barley and its relatives. Now that a reference genome is available, computational pipelines for high-quality sequence assembly are in place, and sequence costs continue to drop, investigations into the structural diversity of the barley genome seem within reach. Here, we review the recent progress on pan-genomics in the model grass Brachypodium distachyon, and the cereal crops rice and maize, and devise a multi-tiered strategy for a pan-genome project in barley. Our design involves: (1) the construction of high-quality de novo sequence assemblies for a small core set of representative genotypes, (2) short-read sequencing of a large diversity panel of genebank accessions to medium coverage and (3) the use of complementary methods such as chromosome-conformation capture sequencing and k-mer-based association genetics. The in silico representation of the barley pan-genome may inform about the mechanisms of structural genome evolution in the Triticeae and supplement quantitative genetics models of crop performance for better accuracy and predictive ability.

Detecting Large Chromosomal Modifications Using Short Read Data From Genotyping-by-Sequencing

Markers linked to agronomic traits are of the prerequisite for molecular breeding. Genotyping-by-sequencing (GBS) data enables to detect small polymorphisms including single nucleotide polymorphisms (SNPs) and short insertions or deletions (InDels) that can be used, for instance, for marker-assisted selection, population genetics, and genome-wide association studies (GWAS). Here, we aim at detecting large chromosomal modifications in barley and wheat based on GBS data. These modifications could be duplications, deletions, substitutions including introgressions as well as alterations of DNA methylation. We demonstrate that GBS coverage analysis is capable to detect Hordeum vulgare/Hordeum bulbosum introgression lines. Furthermore, we identify large chromosomal modifications in barley and wheat collections. Hence, large chromosomal modifications, including introgressions and copy number variations (CNV), can be detected easily and can be used as markers in research and breeding without additional wet-lab experiments.

Dominant integration locus drives continuous diversification of plant immune receptors with exogenous domain fusions

BACKGROUND

The plant immune system is innate and encoded in the germline. Using it efficiently, plants are capable of recognizing a diverse range of rapidly evolving pathogens. A recently described phenomenon shows that plant immune receptors are able to recognize pathogen effectors through the acquisition of exogenous protein domains from other plant genes.

RESULTS

We show that plant immune receptors with integrated domains are distributed unevenly across their phylogeny in grasses. Using phylogenetic analysis, we uncover a major integration clade, whose members underwent repeated independent integration events producing diverse fusions. This clade is ancestral in grasses with members often found on syntenic chromosomes. Analyses of these fusion events reveals that homologous receptors can be fused to diverse domains. Furthermore, we discover a 43 amino acid long motif associated with this dominant integration clade which is located immediately upstream of the fusion site. Sequence analysis reveals that DNA transposition and/or ectopic recombination are the most likely mechanisms of formation for nucleotide binding leucine rich repeat proteins with integrated domains.

CONCLUSIONS

The identification of this subclass of plant immune receptors that is naturally adapted to new domain integration will inform biotechnological approaches for generating synthetic receptors with novel pathogen "baits."

Copy number variation and disease resistance in plants

Plant genome diversity varies from single nucleotide polymorphisms to large-scale deletions, insertions, duplications, or re-arrangements. These re-arrangements of sequences resulting from duplication, gains or losses of DNA segments are termed copy number variations (CNVs). During the last decade, numerous studies have emphasized the importance of CNVs as a factor affecting human phenotype; in particular, CNVs have been associated with risks for several severe diseases. In plants, the exploration of the extent and role of CNVs in resistance against pathogens and pests is just beginning. Since CNVs are likely to be associated with disease resistance in plants, an understanding of the distribution of CNVs could assist in the identification of novel plant disease-resistance genes. In this paper, we review existing information about CNVs; their importance, role and function, as well as their association with disease resistance in plants.

NLR diversity, helpers and integrated domains: making sense of the NLR IDentity

Plant innate immunity relies on genetically predetermined repertoires of immune receptors to detect pathogens and trigger an effective immune response. A large proportion of these receptors are from the Nucletoide Binding Leucine Rich Repeat (NLR) gene family. As plants live longer than most pathogens, maintaining diversity of NLRs and deploying efficient 'pathogen traps' is necessary to withstand the evolutionary battle. In this review, we summarize the sources of diversity in NLR plant immune receptors giving an overview of genomic, regulatory as well as functional studies, including the latest concepts of NLR helpers and NLRs with integrated domains.

Isolation, characterization and genetic diversity of NBS-LRR class disease-resistant gene analogs in multiple virus resistant line of chilli (Capsicum annuum L.)

Viruses are serious threat to chilli crop production worldwide. Resistance screening against several viruses resulted in identifying a multiple virus resistant genotype 'IHR 2451'. Degenerate primers based on the conserved regions between P-Loop and GLPL of Resistance genes (R-genes) were used to amplify nucleotide binding sites (NBS)-encoding regions from genotype 'IHR 2451'. Alignment of deduced amino acid sequences and phylogenetic analyses of isolated sequences distinguished into two groups representing toll interleukin-1 receptor (TIR) and non-TIR, and different families within the group confirming the hypotheses that dicots have both the types of NBS-LRR genes. The alignment of deduced amino acid sequences revealed conservation of subdomains P-loop, RNBS-A, kinase2, RNBS-B, and GLPL. The distinctive five RGAs showing specific conserved motifs were subjected to BLASTp and indicated high homology at deduced amino acid level with R genes identified such as Pvr9 gene for potyvirus resistance, putative late blight resistance protein homolog R1B-23 and other disease resistance genes suggesting high correlation with resistance to different pathogens. These pepper RGAs could be regarded as candidate sequences of resistant genes for marker development.

ARGOS8 variants generated by CRISPR-Cas9 improve maize grain yield under field drought stress conditions

Maize ARGOS8 is a negative regulator of ethylene responses. A previous study has shown that transgenic plants constitutively overexpressing ARGOS8 have reduced ethylene sensitivity and improved grain yield under drought stress conditions. To explore the targeted use of ARGOS8 native expression variation in drought-tolerant breeding, a diverse set of over 400 maize inbreds was examined for ARGOS8 mRNA expression, but the expression levels in all lines were less than that created in the original ARGOS8 transgenic events. We then employed a CRISPR-Cas-enabled advanced breeding technology to generate novel variants of ARGOS8. The native maize GOS2 promoter, which confers a moderate level of constitutive expression, was inserted into the 5'-untranslated region of the native ARGOS8 gene or was used to replace the native promoter of ARGOS8. Precise genomic DNA modification at the ARGOS8 locus was verified by PCR and sequencing. The ARGOS8 variants had elevated levels of ARGOS8 transcripts relative to the native allele and these transcripts were detectable in all the tissues tested, which was the expected results using the GOS2 promoter. A field study showed that compared to the WT, the ARGOS8 variants increased grain yield by five bushels per acre under flowering stress conditions and had no yield loss under well-watered conditions. These results demonstrate the utility of the CRISPR-Cas9 system in generating novel allelic variation for breeding drought-tolerant crops.

Quantitative Trait Loci and Resistance Gene Analogs Associated with Maize White Spot Resistance

Maize white spot (MWS), caused by the bacterium Pantoea ananatis, is one of the most important maize foliar diseases in tropical and subtropical regions, causing significant yield losses. Despite its economic importance, genetic studies of MWS are scarce. The aim of this study was to map quantitative trait loci (QTL) associated with MWS resistance and to identify resistance gene analogs (RGA) underlying these QTL. QTL mapping was performed in a tropical maize F_2:3 population, which was genotyped with simple-sequence repeat and RGA-tagged markers and phenotyped for the response to MWS in two Brazilian southeastern locations. Nine QTL explained approximately 70% of the phenotypic variance for MWS resistance at each location, with two of them consistently detected in both environments. Data mining using 112 resistance genes cloned from different plant species revealed 1,697 RGA distributed in clusters within the maize genome. The RGA Pto19, Pto20, Pto99, and Xa26.151.4 were genetically mapped within MWS resistance QTL on chromosomes 4 and 8 and were preferentially expressed in the resistant parental line at locations where their respective QTL occurred. The consistency of QTL mapping, in silico prediction, and gene expression analyses revealed RGA and genomic regions suitable for marker-assisted selection to improve MWS resistance.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

Fine Mapping of Two Wheat Powdery Mildew Resistance Genes Located at the Cluster

Powdery mildew caused by (DC.) f. sp. () is a globally devastating foliar disease of wheat ( L.). More than a dozen genes against this disease, identified from wheat germplasms of different ploidy levels, have been mapped to the region surrounding the locus on the long arm of chromosome 7A, which forms a resistance ()-gene cluster. and from einkorn wheat ( L.) were two of the genes belonging to this cluster. This study was initiated to fine map these two genes toward map-based cloning. Comparative genomics study showed that macrocolinearity exists between L. chromosome 1 (Bd1) and the - region, which allowed us to develop markers based on the wheat sequences orthologous to genes contained in the Bd1 region. With these and other newly developed and published markers, high-resolution maps were constructed for both and using large F populations. Moreover, a physical map of was constructed through chromosome walking with bacterial artificial chromosome (BAC) clones and comparative mapping. Eventually, and were restricted to a 0.12- and 0.86-cM interval, respectively. Based on the closely linked common markers, , , and (another powdery mildew resistance gene in the cluster) were not allelic to one another. Severe recombination suppression and disruption of synteny were noted in the region encompassing . These results provided useful information for map-based cloning of the genes in the cluster and interpretation of their evolution.

Molecular phylogeny and dynamic evolution of disease resistance genes in the legume family

Background

Legumes are the second-most important crop family in agriculture for its economic and nutritional values. Disease resistance (R-) genes play an important role in responding to pathogen infections in plants. To further increase the yield of legume crops, we need a comprehensive understanding of the evolution of R-genes in the legume family.

Results

In this study, we developed a robust pipeline and identified a total of 4,217 R-genes in the genomes of seven sequenced legume species. A dramatic diversity of R-genes with structural variances indicated a rapid birth-and-death rate during the R-gene evolution in legumes. The number of R-genes transiently expanded and then quickly contracted after whole-genome duplications, which meant that R-genes were sensitive to subsequent diploidization. R proteins with the Coiled-coil (CC) domain are more conserved than others in legumes. Meanwhile, other types of legume R proteins with only one or two typical domains were subjected to higher rates of loss during evolution. Although R-genes evolved quickly in legumes, they tended to undergo purifying selection instead of positive selection during evolution. In addition, domestication events in some legume species preferentially selected for the genes directly involved in the plant-pathogen interaction pathway while suppressing those R-genes with low occurrence rates.

Conclusions

Our results provide insights into the dynamic evolution of R-genes in the legume family, which will be valuable for facilitating genetic improvements in the disease resistance of legume cultivars.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2736-9) contains supplementary material, which is available to authorized users.

Morphological features of different polyploids for adaptation and molecular characterization of CC-NBS-LRR and LEA gene families in Agave L

Polyploidy has been widely described in many Agave L. species, but its influence on environmental response to stress is still unknown. With the objective of knowing the morphological adaptations and regulation responses of genes related to biotic (LEA) and abiotic (NBS-LRR) stress in species of Agave with different levels of ploidy, and how these factors contribute to major response of Agave against environmental stresses, we analyzed 16 morphological trials on five accessions of three species (Agave tequilana Weber, Agave angustifolia Haw. and Agave fourcroydes Lem.) with different ploidy levels (2n=2x=60 2n=3x=90, 2n=5x=150, 2n=6x=180) and evaluated the expression of NBS-LRR and LEA genes regulated by biotic and abiotic stress. It was possible to associate some morphological traits (spines, nuclei, and stomata) to ploidy level. The genetic characterization of stress-related genes NBS-LRR induced by pathogenic infection and LEA by heat or saline stresses indicated that amino acid sequence analysis in these genes showed more substitutions in higher ploidy level accessions of A. fourcroydes Lem. 'Sac Ki' (2n=5x=150) and A. angustifolia Haw. 'Chelem Ki' (2n=6x=180), and a higher LEA and NBS-LRR representativeness when compared to their diploid and triploid counterparts. In all studied Agave accessions expression of LEA and NBS-LRR genes was induced by saline or heat stresses or by infection with Erwinia carotovora, respectively. The transcriptional activation was also higher in A. angustifolia Haw. 'Chelem Ki' (2n=6x=180) and A. fourcroydes 'Sac Ki' (2n=5x=150) than in their diploid and triploid counterparts, which suggests higher adaptation to stress. Finally, the diploid accession A. tequilana Weber 'Azul' showed a differentiated genetic profile relative to other Agave accessions. The differences include similar or higher genetic representativeness and transcript accumulation of LEA and NBS-LRR genes than in polyploid (2n=5x=150 and 2n=6x=180) Agave accessions, thus suggesting a differentiated selection pressure for overcoming the lower ploidy level of the diploid A. tequilana Weber 'Azul'.

Assessing the Barley Genome Zipper and Genomic Resources for Breeding Purposes

The aim of this study was to estimate the accuracy and convergence of newly developed barley (Hordeum vulgare L.) genomic resources, primarily genome zipper (GZ) and population sequencing (POPSEQ), at the genome-wide level and to assess their usefulness in applied barley breeding by analyzing seven known loci. Comparison of barley GZ and POPSEQ maps to a newly developed consensus genetic map constructed with data from 13 individual linkage maps yielded an accuracy of 97.8% (GZ) and 99.3% (POPSEQ), respectively, regarding the chromosome assignment. The percentage of agreement in marker position indicates that on average only 3.7% GZ and 0.7% POPSEQ positions are not in accordance with their centimorgan coordinates in the consensus map. The fine-scale comparison involved seven genetic regions on chromosomes 1H, 2H, 4H, 6H, and 7H, harboring major genes and quantitative trait loci (QTL) for disease resistance. In total, 179 GZ loci were analyzed and 64 polymorphic markers were developed. Entirely, 89.1% of these were allocated within the targeted intervals and 84.2% followed the predicted order. Forty-four markers showed a match to a POPSEQ-anchored contig, the percentage of collinearity being 93.2%, on average. Forty-four markers allowed the identification of twenty-five fingerprinted contigs (FPCs) and a more clear delimitation of the physical regions containing the traits of interest. Our results demonstrate that an increase in marker density of barley maps by using new genomic data significantly improves the accuracy of GZ. In addition, the combination of different barley genomic resources can be considered as a powerful tool to accelerate barley breeding.

Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B

Background

Bread wheat is not only an important crop, but its large (17 Gb), highly repetitive, and hexaploid genome makes it a good model to study the organization and evolution of complex genomes. Recently, we produced a high quality reference sequence of wheat chromosome 3B (774 Mb), which provides an excellent opportunity to study the evolutionary dynamics of a large and polyploid genome, specifically the impact of single gene duplications.

Results

We find that 27 % of the 3B predicted genes are non-syntenic with the orthologous chromosomes of Brachypodium distachyon, Oryza sativa, and Sorghum bicolor, whereas, by applying the same criteria, non-syntenic genes represent on average only 10 % of the predicted genes in these three model grasses. These non-syntenic genes on 3B have high sequence similarity to at least one other gene in the wheat genome, indicating that hexaploid wheat has undergone massive small-scale interchromosomal gene duplications compared to other grasses. Insertions of non-syntenic genes occurred at a similar rate along the chromosome, but these genes tend to be retained at a higher frequency in the distal, recombinogenic regions. The ratio of non-synonymous to synonymous substitution rates showed a more relaxed selection pressure for non-syntenic genes compared to syntenic genes, and gene ontology analysis indicated that non-syntenic genes may be enriched in functions involved in disease resistance.

Conclusion

Our results highlight the major impact of single gene duplications on the wheat gene complement and confirm the accelerated evolution of the Triticeae lineage among grasses.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0754-6) contains supplementary material, which is available to authorized users.

Identification of a novel major locus for gray leaf spot resistance in Italian ryegrass (Lolium multiflorum Lam.)

BACKGROUND

Gray leaf spot (GLS), caused by Magnaporthe oryzae (anamorph Pyricularia oryzae), in ryegrasses is a very serious problem. Heavily infected small seedlings die within a matter of days, and stands of the grasses are seriously damaged by the disease. Thus, the development of GLS-resistant cultivars has become a concern in ryegrass breeding.

RESULTS

Phenotypic segregations in a single cross-derived F1 population of Italian ryegrass (Lolium multiflorum Lam.) indicated that the GLS resistance in the population was possibly controlled by one or two dominant genes with 66.5-77.9% of broad-sense heritability. In bulked segregant analyses, two simple sequence repeat (SSR) markers, which have so far been reported to locate on linkage group (LG) 3 of Italian ryegrass, showed specific signals in the resistant parent and resistant bulk, indicating that the resistance gene locus was possibly in the LG 3. We thus constructed a genetic linkage map of the LG 3 covering 133.6 centimorgan with other SSR markers of the LG 3 of Italian ryegrass and grass anchor probes that have previously been assigned to LG 3 of ryegrasses, and with rice expressed sequence tag (EST)-derived markers selected from a rice EST map of chromosome (Chr) 1 since LG 3 of ryegrasses are syntenic to rice Chr 1. Quantitative trait locus (QTL) analysis with the genetic linkage map and phenotypic data of the F1 population detected a major locus for GLS resistance. Proportions of phenotypic variance explained by the QTL at the highest logarithm of odds scores were 61.0-69.5%.

CONCLUSIONS

A resistance locus was confirmed as novel for GLS resistance, because its genetic position was different from other known loci for GLS resistance. Broad-sense heritability and the proportion of phenotypic variance explained by the QTL were similar, suggesting that most of the genetic factors for the resistance phenotype against GLS in the F1 population can be explained by a function of the single resistance locus. We designated the putative gene for the novel resistance locus as LmPi2. LmPi2 will be useful for future development of GLS-resistant cultivars in combination with other resistance genes.

Whole-genome profiling and shotgun sequencing delivers an anchored, gene-decorated, physical map assembly of bread wheat chromosome 6A

Bread wheat (Triticum aestivum L.) is the most important staple food crop for 35% of the world's population. International efforts are underway to facilitate an increase in wheat production, of which the International Wheat Genome Sequencing Consortium (IWGSC) plays an important role. As part of this effort, we have developed a sequence-based physical map of wheat chromosome 6A using whole-genome profiling (WGP™). The bacterial artificial chromosome (BAC) contig assembly tools fingerprinted contig (fpc) and linear topological contig (ltc) were used and their contig assemblies were compared. A detailed investigation of the contigs structure revealed that ltc created a highly robust assembly compared with those formed by fpc. The ltc assemblies contained 1217 contigs for the short arm and 1113 contigs for the long arm, with an L50 of 1 Mb. To facilitate in silico anchoring, WGP™ tags underlying BAC contigs were extended by wheat and wheat progenitor genome sequence information. Sequence data were used for in silico anchoring against genetic markers with known sequences, of which almost 79% of the physical map could be anchored. Moreover, the assigned sequence information led to the 'decoration' of the respective physical map with 3359 anchored genes. Thus, this robust and genetically anchored physical map will serve as a framework for the sequencing of wheat chromosome 6A, and is of immediate use for map-based isolation of agronomically important genes/quantitative trait loci located on this chromosome.

Characterization of novel wheat NBS domain-containing sequences and their utilization, in silico, for genome-scale R-gene mining

In crop improvement, the isolation, cloning and transfer of disease resistance genes (R-genes) is an ultimate goal usually starting from tentative R-gene analogs (RGAs) that are identified on the basis of their structure. For bread wheat, recent advances in genome sequencing are supporting the efforts of wheat geneticists worldwide. Among wheat R-genes, nucleotide-binding site (NBS)-encoding ones represent a major class. In this study, we have used a polymerase chain reaction-based approach to amplify and clone NBS-type RGAs from a bread wheat cultivar, 'Salambo 80.' Four novel complete ORF sequences showing similarities to previously reported R-genes/RGAs were used for in silico analyses. In a first step, where analyses were focused on the NBS domain, these sequences were phylogenetically assigned to two distinct groups: a first group close to leaf rust Lr21 resistance proteins; and a second one similar to cyst nematode resistance proteins. In a second step, sequences were used as initial seeds to walk up and downstream the NBS domain. This procedure enabled identifying 8 loci ranging in size between 2,115 and 7,653 bp. Ab initio gene prediction identified 8 gene models, among which two had complete ORFs. While GenBank survey confirmed the belonging of sequences to two groups, subsequent characterization using IWGSC genomic and proteomic data showed that the 8 gene models, reported in this study, were unique and their loci matched scaffolds on chromosome arms 1AS, 1BS, 4BS and 1DS. The gene model located on 1DS is a pseudo-Lr21 that was shown to have an NBS-LRR domain structure, while the potential association of the RGAs, here reported, is discussed. This study has produced novel R-gene-like loci and models in the wheat genome and provides the first steps toward further elucidation of their role in wheat disease resistance.

Paleo-evolutionary plasticity of plant disease resistance genes

Background

The recent access to a large set of genome sequences, combined with a robust evolutionary scenario of modern monocot (i.e. grasses) and eudicot (i.e. rosids) species from their founder ancestors, offered the opportunity to gain insights into disease resistance genes (R-genes) evolutionary plasticity.

Results

We unravel in the current article (i) a R-genes repertoire consisting in 7883 for monocots and 15758 for eudicots, (ii) a contrasted R-genes conservation with 23.8% for monocots and 6.6% for dicots, (iii) a minimal ancestral founder pool of 384 R-genes for the monocots and 150 R-genes for the eudicots, (iv) a general pattern of organization in clusters accounting for more than 60% of mapped R-genes, (v) a biased deletion of ancestral duplicated R-genes between paralogous blocks possibly compensated by clusterization, (vi) a bias in R-genes clusterization where Leucine-Rich Repeats act as a ‘glue’ for domain association, (vii) a R-genes/miRNAs interome enriched toward duplicated R-genes.

Conclusions

Together, our data may suggest that R-genes family plasticity operated during plant evolution (i) at the structural level through massive duplicates loss counterbalanced by massive clusterization following polyploidization; as well as at (ii) the regulation level through microRNA/R-gene interactions acting as a possible source of functional diploidization of structurally retained R-genes duplicates. Such evolutionary shuffling events leaded to CNVs (i.e. Copy Number Variation) and PAVs (i.e. Presence Absence Variation) between related species operating in the decay of R-genes colinearity between plant species.

The physical map of wheat chromosome 1BS provides insights into its gene space organization and evolution

Background

The wheat genome sequence is an essential tool for advanced genomic research and improvements. The generation of a high-quality wheat genome sequence is challenging due to its complex 17 Gb polyploid genome. To overcome these difficulties, sequencing through the construction of BAC-based physical maps of individual chromosomes is employed by the wheat genomics community. Here, we present the construction of the first comprehensive physical map of chromosome 1BS, and illustrate its unique gene space organization and evolution.

Results

Fingerprinted BAC clones were assembled into 57 long scaffolds, anchored and ordered with 2,438 markers, covering 83% of chromosome 1BS. The BAC-based chromosome 1BS physical map and gene order of the orthologous regions of model grass species were consistent, providing strong support for the reliability of the chromosome 1BS assembly. The gene space for chromosome 1BS spans the entire length of the chromosome arm, with 76% of the genes organized in small gene islands, accompanied by a two-fold increase in gene density from the centromere to the telomere.

Conclusions

This study provides new evidence on common and chromosome-specific features in the organization and evolution of the wheat genome, including a non-uniform distribution of gene density along the centromere-telomere axis, abundance of non-syntenic genes, the degree of colinearity with other grass genomes and a non-uniform size expansion along the centromere-telomere axis compared with other model cereal genomes. The high-quality physical map constructed in this study provides a solid basis for the assembly of a reference sequence of chromosome 1BS and for breeding applications.

Identification and phylogenetic analysis of a CC-NBS-LRR encoding gene assigned on chromosome 7B of wheat

Hexaploid wheat displays limited genetic variation. As a direct A and B genome donor of hexaploid wheat, tetraploid wheat represents an important gene pool for cultivated bread wheat. Many disease resistant genes express conserved domains of the nucleotide-binding site and leucine-rich repeats (NBS-LRR). In this study, we isolated a CC-NBS-LRR gene locating on chromosome 7B from durum wheat variety Italy 363, and designated it TdRGA-7Ba. Its open reading frame was 4014 bp, encoding a 1337 amino acid protein with a complete NBS domain and 18 LRR repeats, sharing 44.7% identity with the PM3B protein. TdRGA-7Ba expression was continuously seen at low levels and was highest in leaves. TdRGA-7Ba has another allele TdRGA-7Bb with a 4 bp deletion at position +1892 in other cultivars of tetraploid wheat. In Ae. speltoides, as a B genome progenitor, both TdRGA-7Ba and TdRGA-7Bb were detected. In all six species of hexaploid wheats (AABBDD), only TdRGA-7Bb existed. Phylogenic analysis showed that all TdRGA-7Bb type genes were grouped in one sub-branch. We speculate that TdRGA-7Bb was derived from a TdRGA-7Ba mutation, and it happened in Ae. speltoides. Both types of TdRGA-7B participated in tetraploid wheat formation. However, only the TdRGA-7Bb was retained in hexaploid wheat.

Isolation and characterization of nucleotide-binding site resistance gene homologues in common bean (Phaseolus vulgaris)

Common bean production is constrained by many fungal, viral, and bacterial pathogens. Thus, the identification of resistance (R) genes is an important focal point of common bean research. The main goal of our study was to identify resistance gene homologues (RGH) in the crop, using degenerate primers designed from conserved sequences in the nucleotide-binding site (NBS) domains of R-genes from the model legume Medicago truncatula. Total DNA of the Andean common bean genotype G19833 was used for amplification of over 500 primer combinations. Sequencing of amplicons showed that 403 cloned fragments had uninterrupted open reading frames and were considered representative of functional RGH genes. The sequences were grouped at two levels of nucleotide identity (90 and 80%) and representative sequences of each group were used for phylogenetic analyses. The RGH sequence diversity of common bean was divided into TIR and non-TIR families, each with different clusters. The TIR sequences grouped into 14 clades while non-TIR sequences grouped into seven clades. Pairwise comparisons showed purifying selection, although some sequences may have been the result of diversifying selection. Knowledge about RGH genes in common bean can allow the design of molecular markers for pyramiding of resistance genes against various pathogens.

Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family

The proper use of resistance genes (R genes) requires a comprehensive understanding of their genomics and evolution. We analyzed genes encoding nucleotide-binding sites and leucine-rich repeats in the genomes of rice (Oryza sativa), maize (Zea mays), sorghum (Sorghum bicolor), and Brachypodium distachyon. Frequent deletions and translocations of R genes generated prevalent presence/absence polymorphism between different accessions/species. The deletions were caused by unequal crossover, homologous repair, nonhomologous repair, or other unknown mechanisms. R gene loci identified from different genomes were mapped onto the chromosomes of rice cv Nipponbare using comparative genomics, resulting in an integrated map of 495 R loci. Sequence analysis of R genes from the partially sequenced genomes of an African rice cultivar and 10 wild accessions suggested that there are many additional R gene lineages in the AA genome of Oryza. The R genes with chimeric structures (termed type I R genes) are diverse in different rice accessions but only account for 5.8% of all R genes in the Nipponbare genome. In contrast, the vast majority of R genes in the rice genome are type II R genes, which are highly conserved in different accessions. Surprisingly, pseudogene-causing mutations in some type II lineages are often conserved, indicating that their conservations were not due to their functions. Functional R genes cloned from rice so far have more type II R genes than type I R genes, but type I R genes are predicted to contribute considerable diversity in wild species. Type I R genes tend to reduce the microsynteny of their flanking regions significantly more than type II R genes, and their flanking regions have slightly but significantly lower G/C content than those of type II R genes.

Deconstruction of the (paleo)polyploid grapevine genome based on the analysis of transposition events involving NBS resistance genes

Plants have followed a reticulate type of evolution and taxa have frequently merged via allopolyploidization. A polyploid structure of sequenced genomes has often been proposed, but the chromosomes belonging to putative component genomes are difficult to identify. The 19 grapevine chromosomes are evolutionary stable structures: their homologous triplets have strongly conserved gene order, interrupted by rare translocations. The aim of this study is to examine how the grapevine nucleotide-binding site (NBS)-encoding resistance (NBS-R) genes have evolved in the genomic context and to understand mechanisms for the genome evolution. We show that, in grapevine, i) helitrons have significantly contributed to transposition of NBS-R genes, and ii) NBS-R gene cluster similarity indicates the existence of two groups of chromosomes (named as Va and Vc) that may have evolved independently. Chromosome triplets consist of two Va and one Vc chromosomes, as expected from the tetraploid and diploid conditions of the two component genomes. The hexaploid state could have been derived from either allopolyploidy or the separation of the Va and Vc component genomes in the same nucleus before fusion, as known for Rosaceae species. Time estimation indicates that grapevine component genomes may have fused about 60 mya, having had at least 40–60 mya to evolve independently. Chromosome number variation in the Vitaceae and related families, and the gap between the time of eudicot radiation and the age of Vitaceae fossils, are accounted for by our hypothesis.

High-density genetic and physical bin mapping of wheat chromosome 1D reveals that the powdery mildew resistance gene Pm24 is located in a highly recombinogenic region

Genetic maps of wheat chromosome 1D consisting of 57 microsatellite marker loci were constructed using Chinese Spring (CS) × Chiyacao F(2) and the International Triticeae Mapping Initiative (ITMI) recombinant inbred lines (RILs) mapping populations. Marker order was consistent, but genetic distances of neighboring markers were different in two populations. Physical bin map of 57 microsatellite marker loci was generated by means of 10 CS 1D deletion lines. The physical bin mapping indicated that microsatellite marker loci were not randomly distributed on chromosome 1D. Nineteen of the 24 (79.2%) microsatellite markers were mapped in the distal 30% genomic region of 1DS, whereas 25 of the 33 (75.8%) markers were assigned to the distal 59% region of 1DL. The powdery mildew resistance gene Pm24, originating from the Chinese wheat landrace Chiyacao, was previously mapped in the vicinity of the centromere on the short arm of chromosome 1D. A high density genetic map of chromosome 1D was constructed, consisting of 36 markers and Pm24, with a total map length of 292.7 cM. Twelve marker loci were found to be closely linked to Pm24. Pm24 was flanked by Xgwm789 (Xgwm603) and Xbarc229 with genetic distances of 2.4 and 3.6 cM, respectively, whereas a microsatellite marker Xgwm1291 co-segregated with Pm24. The microsatellite marker Xgwm1291 was assigned to the bin 1DS5-0.70-1.00 of the chromosome arm 1DS. It could be concluded that Pm24 is located in the '1S0.8 gene-rich region', a highly recombinogenic region of wheat. The results presented here would provide a start point for the map-based cloning of Pm24.

Relative evolutionary rates of NBS-encoding genes revealed by soybean segmental duplication

It is well known that nucleotide binding site (NBS)-encoding genes are duplicate-rich and fast-evolving genes. However, there is little information on the relative importance of tandem and segmental NBS duplicates and their exact evolutionary rates. The two rounds of large-scale duplication that have occurred in soybean provide a unique opportunity to investigate these issues. Comparison of NBS and non-NBS genes on segments of syntenic homoeologs shows that NBS-encoding genes evolve at least 1.5-fold faster (~1.5-fold higher synonymous and approximately 2.3-fold higher nonsynonymous substitution rates) and lose their genes approximately twofold faster than the flanking non-NBS genes. Compared with segmental duplicates, tandem NBS duplicates are more abundant in soybean, suggesting that tandem duplication is the major driving force in the expansion of NBS genes. Notably, significant sequence exchanges along with significantly positive selection were detected in most tandem-duplicated NBS gene families. The results suggest that the rapid evolution of NBS genes may be due to the combined effects of diversifying selection and frequent sequence exchanges. Interestingly, TIR-NBS-LRR genes (TNLs) have a higher nucleotide substitution rate than non-TNLs, indicating that these types of NBS genes may have a rather different evolutionary pattern. It is important to determine the exact relative evolutionary rates of TNL, non-TNL, and non-NBS genes in order to understand how fast the host plant can adjust its response to rapidly evolving pathogens in a coevolutionary context.

Identification and characterization of potential NBS-encoding resistance genes and induction kinetics of a putative candidate gene associated with downy mildew resistance in Cucumis

BACKGROUND

Due to the variation and mutation of the races of Pseudoperonospora cubensis, downy mildew has in recent years become the most devastating leaf disease of cucumber worldwide. Novel resistance to downy mildew has been identified in the wild Cucumis species, C. hystrix Chakr. After the successful hybridization between C. hystrix and cultivated cucumber (C. sativus L.), an introgression line (IL5211S) was identified as highly resistant to downy mildew. Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes are the largest class of disease resistance genes cloned from plant with highly conserved domains, which can be used to facilitate the isolation of candidate genes associated with downy mildew resistance in IL5211S.

RESULTS

Degenerate primers that were designed based on the conserved motifs in the NBS domain of resistance (R) proteins were used to isolate NBS-type sequences from IL5211S. A total of 28 sequences were identified and named as cucumber (C. sativus = CS) resistance gene analogs as CSRGAs. Polygenetic analyses separated these sequences into four different classes. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that these CSRGAs expressed at different levels in leaves, roots, and stems. In addition, introgression from C. hystrix induced expression of the partial CSRGAs in cultivated cucumber, especially CSRGA23, increased four-fold when compared to the backcross parent CC3. Furthermore, the expression of CSRGA23 under P. cubensis infection and abiotic stresses was also analyzed at different time points. Results showed that the P. cubensis treatment and four tested abiotic stimuli, MeJA, SA, ABA, and H2O2, triggered a significant induction of CSRGA23 within 72 h of inoculation. The results indicate that CSRGA23 may play a critical role in protecting cucumber against P. cubensis through a signaling the pathway triggered by these molecules.

CONCLUSIONS

Four classes of NBS-type RGAs were successfully isolated from IL5211S, and the possible involvement of CSRGA23 in the active defense response to P. cubensis was demonstrated. These results will contribute to develop analog-based markers related to downy mildew resistance gene and elucidate the molecular mechanisms causing resistance in IL5211S in the future.

Combined meta-genomics analyses unravel candidate genes for the grain dietary fiber content in bread wheat (Triticum aestivum L.)

Grain dietary fiber content in wheat not only affects its end use and technological properties including milling, baking and animal feed but is also of great importance for health benefits. In this study, integration of association genetics (seven detected loci on chromosomes 1B, 3A, 3D, 5B, 6B, 7A, 7B) and meta-QTL (three consensus QTL on chromosomes 1B, 3D and 6B) analyses allowed the identification of seven chromosomal regions underlying grain dietary fiber content in bread wheat. Based either on a diversity panel or on bi-parental populations, we clearly demonstrate that this trait is mainly driven by a major locus located on chromosome 1B associated with a log of p value >13 and a LOD score >8, respectively. In parallel, we identified 73 genes differentially expressed during the grain development and between genotypes with contrasting grain fiber contents. Integration of quantitative genetics and transcriptomic data allowed us to propose a short list of candidate genes that are conserved in the rice, sorghum and Brachypodium chromosome regions orthologous to the seven wheat grain fiber content QTL and that can be considered as major candidate genes for future improvement of the grain dietary fiber content in bread wheat breeding programs.

Identification and mapping of PmG16, a powdery mildew resistance gene derived from wild emmer wheat

The gene-pool of wild emmer wheat, Triticum turgidum ssp. dicoccoides, harbors a rich allelic repertoire for disease resistance. In the current study, we made use of tetraploid wheat mapping populations derived from a cross between durum wheat (cv. Langdon) and wild emmer (accession G18-16) to identify and map a new powdery mildew resistance gene derived from wild emmer wheat. Initially, the two parental lines were screened with a collection of 42 isolates of Blumeria graminis f. sp. tritici (Bgt) from Israel and 5 isolates from Switzerland. While G18-16 was resistant to 34 isolates, Langdon was resistant only to 5 isolates and susceptible to 42 isolates. Isolate Bgt#15 was selected to differentiate between the disease reactions of the two genotypes. Segregation ratio of F(2-3) and recombinant inbreed line (F(7)) populations to inoculation with isolate Bgt#15 indicated the role of a single dominant gene in conferring resistance to Bgt#15. This gene, temporarily designated PmG16, was located on the distal region of chromosome arm 7AL. Genetic map of PmG16 region was assembled with 32 simple sequence repeat (SSR), sequence tag site (STS), Diversity array technology (DArT) and cleaved amplified polymorphic sequence (CAPS) markers and assigned to the 7AL physical bin map (7AL-16). Using four DNA markers we established colinearity between the genomic region spanning the PmG16 locus within the distal region of chromosome arm 7AL and the genomic regions on rice chromosome 6 and Brachypodium Bd1. A comparative analysis was carried out between PmG16 and other known Pm genes located on chromosome arm 7AL. The identified PmG16 may facilitate the use of wild alleles for improvement of powdery mildew resistance in elite wheat cultivars via marker-assisted selection.

Comparative genomic sequence analysis of strawberry and other rosids reveals significant microsynteny

BACKGROUND

Fragaria belongs to the Rosaceae, an economically important family that includes a number of important fruit producing genera such as Malus and Prunus. Using genomic sequences from 50 Fragaria fosmids, we have examined the microsynteny between Fragaria and other plant models.

RESULTS

In more than half of the strawberry fosmids, we found syntenic regions that are conserved in Populus, Vitis, Medicago and/or Arabidopsis with Populus containing the greatest number of syntenic regions with Fragaria. The longest syntenic region was between LG VIII of the poplar genome and the strawberry fosmid 72E18, where seven out of twelve predicted genes were collinear. We also observed an unexpectedly high level of conserved synteny between Fragaria (rosid I) and Vitis (basal rosid). One of the strawberry fosmids, 34E24, contained a cluster of R gene analogs (RGAs) with NBS and LRR domains. We detected clusters of RGAs with high sequence similarity to those in 34E24 in all the genomes compared. In the phylogenetic tree we have generated, all the NBS-LRR genes grouped together with Arabidopsis CNL-A type NBS-LRR genes. The Fragaria RGA grouped together with those of Vitis and Populus in the phylogenetic tree.

CONCLUSIONS

Our analysis shows considerable microsynteny between Fragaria and other plant genomes such as Populus, Medicago, Vitis, and Arabidopsis to a lesser degree. We also detected a cluster of NBS-LRR type genes that are conserved in all the genomes compared.

Numbers of genes in the NBS and RLK families vary by more than four-fold within a plant species and are regulated by multiple factors

Many genes exist in the form of families; however, little is known about their size variation, evolution and biology. Here, we present the size variation and evolution of the nucleotide-binding site (NBS)-encoding gene family and receptor-like kinase (RLK) gene family in Oryza, Glycine and Gossypium. The sizes of both families vary by numeral fold, not only among species, surprisingly, also within a species. The size variations of the gene families are shown to correlate with each other, indicating their interactions, and driven by natural selection, artificial selection and genome size variation, but likely not by polyploidization. The numbers of genes in the families in a polyploid species are similar to those of one of its diploid donors, suggesting that polyploidization plays little roles in the expansion of the gene families and that organisms tend not to maintain their ‘surplus’ genes in the course of evolution. Furthermore, it is found that the size variations of both gene families are associated with organisms’ phylogeny, suggesting their roles in speciation and evolution. Since both selection and speciation act on organism’s morphological, physiological and biological variation, our results indicate that the variation of gene family size provides a source of genetic variation and evolution.

Phylogenetic analyses of peanut resistance gene candidates and screening of different genotypes for polymorphic markers

The nucleotide-binding-site-leucine-rich-repeat (NBS-LRR)-encoding gene family has attracted much research interest because approximately 75% of the plant disease resistance genes that have been cloned to date are from this gene family. Here, we describe a collection of peanut NBS-LRR resistance gene candidates (RGCs) isolated from peanut (Arachis) species by mining Gene Bank data base. NBS-LRR sequences assembled into TIR-NBS-LRR (75.4%) and non-TIR-NBS-LRR (24.6%) subfamilies. Total of 20 distinct clades were identified and showed a high level of sequence divergence within TIR-NBS and non-TIR-NBS subfamilies. Thirty-four primer pairs were designed from these RGC sequences and used for screening different genotypes belonging to wild and cultivated peanuts. Therefore, peanut RGC identified in this study will provide useful tools for developing DNA markers and cloning the genes for resistance to different pathogens in peanut.

Molecular genetic description of the cryptic wheat-Aegilops geniculata introgression carrying rust resistance genes Lr57 and Yr40 using wheat ESTs and synteny with rice

The cryptic wheat-alien translocation T5DL.5DS-5MgS(0.95), with leaf rust and stripe rust resistance genes Lr57 and Yr40 transferred from Aegilops geniculata (UgMg) into common wheat, was further analyzed. Molecular genetic analysis using physically mapped ESTs showed that the alien segment in T5DL.5DS-5MgS(0.95) represented only a fraction of the wheat deletion bin 5DS2-0.78-1.00 and was less than 3.3 cM in length in the diploid wheat genetic map. Comparative genomic analysis indicated a high level of colinearity between the distal region of the long arm of chromosome 12 of rice and the genomic region spanning the Lr57 and Yr40 genes in wheat. The alien segment with genes Lr57 and Yr40 corresponds to fewer than four overlapping BAC or PAC clones of the syntenic rice chromosome arm 12L. The wheat-alien translocation breakpoint in T5DL.5DS-5MgS(0.95) was further localized to a single BAC clone of the syntenic rice genomic sequence. The small size of the terminal wheat-alien translocation, as established precisely with respect to Chinese Spring deletion bins and the syntenic rice genomic sequence, further confirmed the escaping nature of cryptic wheat-alien translocations in introgressive breeding. The molecular genetic resources and information developed in the present study will facilitate further fine-scale physical mapping and map-based cloning of the Lr57 and Yr40 genes.

A nomadic subtelomeric disease resistance gene cluster in common bean

The B4 resistance (R) gene cluster is one of the largest clusters known in common bean (Phaseolus vulgaris [Pv]). It is located in a peculiar genomic environment in the subtelomeric region of the short arm of chromosome 4, adjacent to two heterochromatic blocks (knobs). We sequenced 650 kb spanning this locus and annotated 97 genes, 26 of which correspond to Coiled-Coil-Nucleotide-Binding-Site-Leucine-Rich-Repeat (CNL). Conserved microsynteny was observed between the Pv B4 locus and corresponding regions of Medicago truncatula and Lotus japonicus in chromosomes Mt6 and Lj2, respectively. The notable exception was the CNL sequences, which were completely absent in these regions. The origin of the Pv B4-CNL sequences was investigated through phylogenetic analysis, which reveals that, in the Pv genome, paralogous CNL genes are shared among nonhomologous chromosomes (4 and 11). Together, our results suggest that Pv B4-CNL was derived from CNL sequences from another cluster, the Co-2 cluster, through an ectopic recombination event. Integration of the soybean (Glycine max) genome data enables us to date more precisely this event and also to infer that a single CNL moved from the Co-2 to the B4 cluster. Moreover, we identified a new 528-bp satellite repeat, referred to as khipu, specific to the Phaseolus genus, present both between B4-CNL sequences and in the two knobs identified at the B4 R gene cluster. The khipu repeat is present on most chromosomal termini, indicating the existence of frequent ectopic recombination events in Pv subtelomeric regions. Our results highlight the importance of ectopic recombination in R gene evolution.

Fine mapping, physical mapping and development of diagnostic markers for the Rrs2 scald resistance gene in barley

The Rrs2 gene confers resistance to the fungal pathogen Rhynchosporium secalis which causes leaf scald, a major barley disease. The Rrs2 gene was fine mapped to an interval of 0.08 cM between markers 693M6_6 and P1D23R on the distal end of barley chromosome 7HS using an Atlas (resistant) x Steffi (susceptible) mapping population of 9,179 F(2)-plants. The establishment of a physical map of the Rrs2 locus led to the discovery that Rrs2 is located in an area of suppressed recombination within this mapping population. The analysis of 58 barley genotypes revealed a large linkage block at the Rrs2 locus extending over several hundred kb which is present only in Rrs2 carrying cultivars. Due to the lack of recombination in the mapping population and the presence of a Rrs2-specific linkage block, we assume a local chromosomal rearrangement (alien introgression or inversion) in Rrs2 carrying varieties. The variety analysis led to the discovery of eight SNPs which were diagnostic for the Rrs2 phenotype. Based on these SNPs diagnostic molecular markers (CAPS and pyrosequencing markers) were developed which are highly useful for marker-assisted selection in resistance gene pyramiding programmes for Rhynchosporium secalis resistance in barley.

Physical mapping, expression analysis and polymorphism survey of resistance gene analogues on chromosome 11 of rice

Rice is the first cereal genome with a finished sequence and a model crop that has important syntenic relationships with other cereal species. The objectives of our study were to identify resistance gene analogue (RGA) sequences from chromosome 11 of rice, understand their expression in other cereals and dicots by in silico analysis, determine their presence on other rice chromosomes, and evaluate the extent of polymorphism and actual expression in a set of rice genotypes. A total of 195 RGAs were predicted and physically localised. Of these, 91.79% expressed in rice, and 51.28% expressed in wheat, which was the highest among other cereals. Among monocots, sugarcane showed the highest (78.92%) expression, while among dicots, RGAs were maximally expressed in Arabidopsis (11.79%). Interestingly, two of the chromosome 11-specific RGAs were found to be expressing in all the organisms studied. Eighty RGAs of chromosome 11 had significant homology with chromosome 12, which was the maximum among all the rice chromosomes. Thirty-one per cent of the RGAs used in polymerase chain reaction (PCR) amplification showed polymorphism in a set of rice genotypes. Actual gene expression analysis revealed post-inoculation induction of one RGA in the rice line IRBB-4 carrying the bacterial blight resistance gene Xa-4. Our results have implications for the development of sequence-based markers and functional validation of specific RGAs in rice.

Genomics in cereals: from genome-wide conserved orthologous set (COS) sequences to candidate genes for trait dissection

Recent updates in comparative genomics among cereals have provided the opportunity to identify conserved orthologous set (COS) DNA sequences for cross-genome map-based cloning of candidate genes underpinning quantitative traits. New tools are described that are applicable to any cereal genome of interest, namely, alignment criterion for orthologous couples identification, as well as the Intron Spanning Marker software to automatically select intron-spanning primer pairs. In order to test the software, it was applied to the bread wheat genome, and 695 COS markers were assigned to 1,535 wheat loci (on average one marker/2.6 cM) based on 827 robust rice-wheat orthologs. Furthermore, 31 of the 695 COS markers were selected to fine map a pentosan viscosity quantitative trait loci (QTL) on wheat chromosome 7A. Among the 31 COS markers, 14 (45%) were polymorphic between the parental lines and 12 were mapped within the QTL confidence interval with one marker every 0.6 cM defining candidate genes among the rice orthologous region.

Mapping of STS markers obtained from oat resistance gene analog sequences

Two previously isolated resistance gene analogs (RGAs) of oat have been located as RFLPs in the reference map of Avena byzantina ‘Kanota’ × Avena sativa ‘Ogle’ in regions either homologous or homoeologous to loci for resistance to Puccinia coronata , the causal agent of crown rust. In this study, the RGAs were mapped in two recombinant inbred line (RIL) populations that segregate for crown rust resistance: the diploid Avena strigosa × Avena wiestii RIL population (Asw), which has been used for mapping the complex locus PcA, and the hexaploid MN841801-1 × Noble-2 RIL population (MN), in which QTLs have been located. To obtain single-locus markers, RGAs were converted to sequence tagged site (STS) markers using a procedure involving extension of the original RGA sequence lengths by PCR genome walking, amplification and cloning of the parental fragments, and identification of single nucleotide polymorphisms. The procedure successfully obtained STSs from different members of the L7M2 family of sequences, the initial NBS of which have nucleotide similarities of >83%. However, for RGA III2.18, the parental lines were not polymorphic for the STSs assayed. A sequence characterized amplified region (SCAR) marker with features of an RGA had been previously identified for gene Pc94. This marker was also mapped in the above RIL populations. Markers based on RGA L7M2 co-localized with markers defining the QTL Prq1a in linkage group MN3, and were located 15.2 cM from PcA in linkage group AswAC. The SCAR marker for Pc94 was also located in the QTL Prq1a but at 39.5 cM from PcA in AswAC, indicating that the NBS-LRR sequence represented by this marker is not related to PcA. L7M2 was also excluded as a member of the PcA cluster, although it could be an appropriate marker for the Prq1a cluster if chromosome rearrangements are postulated.

The fractionated orthology of Bs2 and Rx/Gpa2 supports shared synteny of disease resistance in the Solanaceae

Comparative genomics provides a powerful tool for the identification of genes that encode traits shared between crop plants and model organisms. Pathogen resistance conferred by plant R genes of the nucleotide-binding-leucine-rich-repeat (NB-LRR) class is one such trait with great agricultural importance that occupies a critical position in understanding fundamental processes of pathogen detection and coevolution. The proposed rapid rearrangement of R genes in genome evolution would make comparative approaches tenuous. Here, we test the hypothesis that orthology is predictive of R-gene genomic location in the Solanaceae using the pepper R gene Bs2. Homologs of Bs2 were compared in terms of sequence and gene and protein architecture. Comparative mapping demonstrated that Bs2 shared macrosynteny with R genes that best fit criteria determined to be its orthologs. Analysis of the genomic sequence encompassing solanaceous R genes revealed the magnitude of transposon insertions and local duplications that resulted in the expansion of the Bs2 intron to 27 kb and the frequently detected duplications of the 5'-end of R genes. However, these duplications did not impact protein expression or function in transient assays. Taken together, our results support a conservation of synteny for NB-LRR genes and further show that their distribution in the genome has been consistent with global rearrangements.

Molecular characterisation and genetic mapping of candidate genes for qualitative disease resistance in perennial ryegrass (Lolium perenne L.)

BACKGROUND

Qualitative pathogen resistance in both dicotyledenous and monocotyledonous plants has been attributed to the action of resistance (R) genes, including those encoding nucleotide binding site--leucine rich repeat (NBS-LRR) proteins and receptor-like kinase enzymes. This study describes the large-scale isolation and characterisation of candidate R genes from perennial ryegrass. The analysis was based on the availability of an expressed sequence tag (EST) resource and a functionally-integrated bioinformatics database.

RESULTS

Amplification of R gene sequences was performed using template EST data and information from orthologous candidate using a degenerate consensus PCR approach. A total of 102 unique partial R genes were cloned, sequenced and functionally annotated. Analysis of motif structure and R gene phylogeny demonstrated that Lolium R genes cluster with putative ortholoci, and evolved from common ancestral origins. Single nucleotide polymorphisms (SNPs) predicted through resequencing of amplicons from the parental genotypes of a genetic mapping family were validated, and 26 distinct R gene loci were assigned to multiple genetic maps. Clusters of largely non-related NBS-LRR genes were located at multiple distinct genomic locations and were commonly found in close proximity to previously mapped defence response (DR) genes. A comparative genomics analysis revealed the co-location of several candidate R genes with disease resistance quantitative trait loci (QTLs).

CONCLUSION

This study is the most comprehensive analysis to date of qualitative disease resistance candidate genes in perennial ryegrass. SNPs identified within candidate genes provide a valuable resource for mapping in various ryegrass pair cross-derived populations and further germplasm analysis using association genetics. In parallel with the use of specific pathogen virulence races, such resources provide the means to identify gene-for-gene mechanisms for multiple host pathogen-interactions and ultimately to obtain durable field-based resistance.

Integration of new CAPS and dCAPS-RGA markers into a composite chickpea genetic map and their association with disease resistance

A composite linkage map was constructed based on two interspecific recombinant inbred line populations derived from crosses between Cicer arietinum (ILC72 and ICCL81001) and Cicer reticulatum (Cr5-10 or Cr5-9). These mapping populations segregate for resistance to ascochyta blight (caused by Ascochyta rabiei), fusarium wilt (caused by Fusarium oxysporum f. sp. ciceris) and rust (caused by Uromyces ciceris-arietini). The presence of single nucleotide polymorphisms in ten resistance gene analogs (RGAs) previously isolated and characterized was exploited. Six out of the ten RGAs were novel sequences. In addition, classes RGA05, RGA06, RGA07, RGA08, RGA09 and RGA10 were considerate putatively functional since they matched with several legume expressed sequences tags (ESTs) obtained under infection conditions. Seven RGA PCR-based markers (5 CAPS and 2 dCAPS) were developed and successfully genotyped in the two progenies. Six of them have been mapped in different linkage groups where major quantitative trait loci conferring resistance to ascochyta blight and fusarium wilt have been reported. Genomic locations of RGAs were compared with those of known Cicer R-genes and previously mapped RGAs. Association was detected between RGA05 and genes controlling resistance to fusarium wilt caused by races 0 and 5.

Dynamic evolution of oryza genomes is revealed by comparative genomic analysis of a genus-wide vertical data set

Oryza (23 species; 10 genome types) contains the world's most important food crop - rice. Although the rice genome serves as an essential tool for biological research, little is known about the evolution of the other Oryza genome types. They contain a historical record of genomic changes that led to diversification of this genus around the world as well as an untapped reservoir of agriculturally important traits. To investigate the evolution of the collective Oryza genome, we sequenced and compared nine orthologous genomic regions encompassing the Adh1-Adh2 genes (from six diploid genome types) with the rice reference sequence. Our analysis revealed the architectural complexities and dynamic evolution of this region that have occurred over the past approximately 15 million years. Of the 46 intact genes and four pseudogenes in the japonica genome, 38 (76%) fell into eight multigene families. Analysis of the evolutionary history of each family revealed independent and lineage-specific gain and loss of gene family members as frequent causes of synteny disruption. Transposable elements were shown to mediate massive replacement of intergenic space (>95%), gene disruption, and gene/gene fragment movement. Three cases of long-range structural variation (inversions/deletions) spanning several hundred kilobases were identified that contributed significantly to genome diversification.