Modeling Flexible Protein Structure With AlphaFold2 and Crosslinking Mass Spectrometry.

We propose a pipeline that combines AlphaFold2 (AF2) and crosslinking mass spectrometry (XL-MS) to model the structure of proteins with multiple conformations. The pipeline consists of two main steps: ensemble generation using AF2 and conformer selection using XL-MS data. For conformer selection, we developed two scores-the monolink probability score (MP) and the crosslink probability score (XLP)-both of which are based on residue depth from the protein surface. We benchmarked MP and XLP on a large dataset of decoy protein structures and showed that our scores outperform previously developed scores. We then tested our methodology on three proteins having an open and closed conformation in the Protein Data Bank: Complement component 3 (C3), luciferase, and glutamine-binding periplasmic protein, first generating ensembles using AF2, which were then screened for the open and closed conformations using experimental XL-MS data. In five out of six cases, the most accurate model within the AF2 ensembles-or a conformation within 1 Å of this model-was identified using crosslinks, as assessed through the XLP score. In the remaining case, only the monolinks (assessed through the MP score) successfully identified the open conformation of glutamine-binding periplasmic protein, and these results were further improved by including the "occupancy" of the monolinks. This serves as a compelling proof-of-concept for the effectiveness of monolinks. In contrast, the AF2 assessment score was only able to identify the most accurate conformation in two out of six cases. Our results highlight the complementarity of AF2 with experimental methods like XL-MS, with the MP and XLP scores providing reliable metrics to assess the quality of the predicted models. The MP and XLP scoring functions mentioned above are available at https://gitlab.com/topf-lab/xlms-tools.

The Systems Architecture of Molecular Memory in Poplar after Abiotic Stress.

Throughout the temperate zones, plants face combined drought and heat spells in increasing frequency and intensity. Here, we compared periodic (intermittent, i.e., high-frequency) versus chronic (continuous, i.e., high-intensity) drought-heat stress scenarios in gray poplar (Populus× canescens) plants for phenotypic and transcriptomic effects during stress and after recovery. Photosynthetic productivity after stress recovery exceeded the performance of poplar trees without stress experience. We analyzed the molecular basis of this stress-related memory phenotype and investigated gene expression responses across five major tree compartments including organs and wood tissues. For each of these tissue samples, transcriptomic changes induced by the two stress scenarios were highly similar during the stress phase but strikingly divergent after recovery. Characteristic molecular response patterns were found across tissues but involved different genes in each tissue. Only a small fraction of genes showed similar stress and recovery expression profiles across all tissues, including type 2C protein phosphatases, the LATE EMBRYOGENESIS ABUNDANT PROTEIN4-5 genes, and homologs of the Arabidopsis (Arabidopsis thaliana) transcription factor HOMEOBOX7. Analysis of the predicted transcription factor regulatory networks for these genes suggested that a complex interplay of common and tissue-specific components contributes to the coordination of post-recovery responses to stress in woody plants.

SSB1/SSB2 Proteins Safeguard B Cell Development by Protecting the Genomes of B Cell Precursors.

Induction of programmed DNA damage and its recognition and repair are fundamental for B cell development. The ssDNA-binding protein SSB1 has been described in human cells as essential for the recognition and repair of DNA damage. To study its relevance for B cells, we recently developed Ssb1 -/- and conditional Ssb1 -/- mice. Although SSB1 loss did not affect B cell development, Ssb1 -/- cells exhibited compensatory expression of its homolog SSB2. We have now generated Ssb2 -/- mice and show in this study that SSB2 is also dispensable for B cell development and DNA damage response activation. In contrast to the single loss of Ssb1 or Ssb2, however, combined SSB1/2 deficiency caused a defect in early B cell development. We relate this to the sensitivity of B cell precursors as mature B cells largely tolerated their loss. Toxicity of combined genetic SSB1/2 loss can be rescued by ectopic expression of either SSB1 or SSB2, mimicked by expression of SSB1 ssDNA-binding mutants, and attenuated by BCL2-mediated suppression of apoptosis. SSB1/2 loss in B cell precursors further caused increased exposure of ssDNA associated with disruption of genome fragile sites, inefficient cell cycle progression, and increased DNA damage if apoptosis is suppressed. As such, our results establish SSB1/2 as safeguards of B cell development and unveil their differential requirement in immature and mature B lymphocytes.

Sensitivity to PI3K and AKT inhibitors is mediated by divergent molecular mechanisms in subtypes of DLBCL.

Activated B-cell-like (ABC) and germinal center B-cell-like diffuse large B-cell lymphoma (DLBCL) represent the 2 major molecular DLBCL subtypes. They are characterized by differences in clinical course and by divergent addiction to oncogenic pathways. To determine activity of novel compounds in these 2 subtypes, we conducted an unbiased pharmacologic in vitro screen. The phosphatidylinositol-3-kinase (PI3K) α/δ (PI3Kα/δ) inhibitor AZD8835 showed marked potency in ABC DLBCL models, whereas the protein kinase B (AKT) inhibitor AZD5363 induced apoptosis in PTEN-deficient DLBCLs irrespective of their molecular subtype. These in vitro results were confirmed in various cell line xenograft and patient-derived xenograft mouse models in vivo. Treatment with AZD8835 induced inhibition of nuclear factor κB signaling, prompting us to combine AZD8835 with the Bruton's tyrosine kinase inhibitor ibrutinib. This combination was synergistic and effective both in vitro and in vivo. In contrast, the AKT inhibitor AZD5363 was effective in PTEN-deficient DLBCLs through downregulation of the oncogenic transcription factor MYC. Collectively, our data suggest that patients should be stratified according to their oncogenic dependencies when treated with PI3K and AKT inhibitors.

Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation.

About 8,000 years ago in the Fertile Crescent, a spontaneous hybridization of the wild diploid grass Aegilops tauschii (2n = 14; DD) with the cultivated tetraploid wheat Triticum turgidum (2n = 4x = 28; AABB) resulted in hexaploid wheat (T. aestivum; 2n = 6x = 42; AABBDD). Wheat has since become a primary staple crop worldwide as a result of its enhanced adaptability to a wide range of climates and improved grain quality for the production of baker's flour. Here we describe sequencing the Ae. tauschii genome and obtaining a roughly 90-fold depth of short reads from libraries with various insert sizes, to gain a better understanding of this genetically complex plant. The assembled scaffolds represented 83.4% of the genome, of which 65.9% comprised transposable elements. We generated comprehensive RNA-Seq data and used it to identify 43,150 protein-coding genes, of which 30,697 (71.1%) were uniquely anchored to chromosomes with an integrated high-density genetic map. Whole-genome analysis revealed gene family expansion in Ae. tauschii of agronomically relevant gene families that were associated with disease resistance, abiotic stress tolerance and grain quality. This draft genome sequence provides insight into the environmental adaptation of bread wheat and can aid in defining the large and complicated genomes of wheat species.

Analysis of the bread wheat genome using whole-genome shotgun sequencing.

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

Modulation of Ambient Temperature-Dependent Flowering in Arabidopsis thaliana by Natural Variation of FLOWERING LOCUS M.

Plants integrate seasonal cues such as temperature and day length to optimally adjust their flowering time to the environment. Compared to the control of flowering before and after winter by the vernalization and day length pathways, mechanisms that delay or promote flowering during a transient cool or warm period, especially during spring, are less well understood. Due to global warming, understanding this ambient temperature pathway has gained increasing importance. In Arabidopsis thaliana, FLOWERING LOCUS M (FLM) is a critical flowering regulator of the ambient temperature pathway. FLM is alternatively spliced in a temperature-dependent manner and the two predominant splice variants, FLM-ß and FLM-δ, can repress and activate flowering in the genetic background of the A. thaliana reference accession Columbia-0. The relevance of this regulatory mechanism for the environmental adaptation across the entire range of the species is, however, unknown. Here, we identify insertion polymorphisms in the first intron of FLM as causative for accelerated flowering in many natural A. thaliana accessions, especially in cool (15°C) temperatures. We present evidence for a potential adaptive role of this structural variation and link it specifically to changes in the abundance of FLM-ß. Our results may allow predicting flowering in response to ambient temperatures in the Brassicaceae.

A synteny-based draft genome sequence of the forage grass Lolium perenne.

Here we report the draft genome sequence of perennial ryegrass (Lolium perenne), an economically important forage and turf grass species that is widely cultivated in temperate regions worldwide. It is classified along with wheat, barley, oats and Brachypodium distachyon in the Pooideae sub-family of the grass family (Poaceae). Transcriptome data was used to identify 28,455 gene models, and we utilized macro-co-linearity between perennial ryegrass and barley, and synteny within the grass family, to establish a synteny-based linear gene order. The gametophytic self-incompatibility mechanism enables the pistil of a plant to reject self-pollen and therefore promote out-crossing. We have used the sequence assembly to characterize transcriptional changes in the stigma during pollination with both compatible and incompatible pollen. Characterization of the pollen transcriptome identified homologs to pollen allergens from a range of species, many of which were expressed to very high levels in mature pollen grains, and are potentially involved in the self-incompatibility mechanism. The genome sequence provides a valuable resource for future breeding efforts based on genomic prediction, and will accelerate the development of new varieties for more productive grasslands.

PTEN loss defines a PI3K/AKT pathway-dependent germinal center subtype of diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) represents a heterogeneous diagnostic category with distinct molecular subtypes that can be defined by gene expression profiling. However, even within these defined subtypes, heterogeneity prevails. To further elucidate the pathogenesis of these entities, we determined the expression of the tumor suppressor phosphatase and tensin homolog (PTEN) in 248 primary DLBCL patient samples. These analyses revealed that loss of PTEN was detectable in 55% of germinal center B-cell-like (GCB) DLBCLs, whereas this abnormality was found in only 14% of non-GCB DLBCL patient samples. In GCB DLBCL, the PTEN status was inversely correlated with activation of the oncogenic PI3K/protein kinase B (AKT) pathway in both DLBCL cell lines and primary patient samples. Reexpression of PTEN induced cytotoxicity in PTEN-deficient GCB DLBCL cell line models by inhibiting PI3K/AKT signaling, indicating an addiction to this pathway in this subset of GCB DLBCLs. PI3K/AKT inhibition induced down-regulation of the transcription factor MYC. Reexpression of MYC rescued GCB DLBCL cells from PTEN-induced toxicity, identifying a regulatory mechanism of MYC expression in DLBCL. Finally, pharmacologic PI3K inhibition resulted in toxicity selectively in PTEN-deficient GCB DLBCL lines. Collectively, our results indicate that PTEN loss defines a PI3K/AKT-dependent GCB DLBCL subtype that is addicted to PI3K and MYC signaling and suggest that pharmacologic inhibition of PI3K might represent a promising therapeutic approach in these lymphomas.

IκB-ζ controls the constitutive NF-κB target gene network and survival of ABC DLBCL.

Constitutive activation of the nuclear factor-κ B (NF-κB) pathway is a hallmark of the activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL). Recurrent mutations of NF-κB regulators that cause constitutive activity of this oncogenic pathway have been identified. However, it remains unclear how specific target genes are regulated. We identified the atypical nuclear IκB protein IκB-ζ to be upregulated in ABC compared with germinal center B-cell-like (GCB) DLBCL primary patient samples. Knockdown of IκB-ζ by RNA interference was toxic to ABC but not to GCB DLBCL cell lines. Gene expression profiling after IκB-ζ knockdown demonstrated a significant downregulation of a large number of known NF-κB target genes, indicating an essential role of IκB-ζ in regulating a specific set of NF-κB target genes. To further investigate how IκB-ζ mediates NF-κB activity, we performed immunoprecipitations and detected a physical interaction of IκB-ζ with both p50 and p52 NF-κB subunits, indicating that IκB-ζ interacts with components of both the canonical and the noncanonical NF-κB pathway in ABC DLBCL. Collectively, our data demonstrate that IκB-ζ is essential for nuclear NF-κB activity in ABC DLBCL, and thus might represent a promising molecular target for future therapies.

MIPS PlantsDB: a database framework for comparative plant genome research.

The rapidly increasing amount of plant genome (sequence) data enables powerful comparative analyses and integrative approaches and also requires structured and comprehensive information resources. Databases are needed for both model and crop plant organisms and both intuitive search/browse views and comparative genomics tools should communicate the data to researchers and help them interpret it. MIPS PlantsDB (http://mips.helmholtz-muenchen.de/plant/genomes.jsp) was initially described in NAR in 2007 [Spannagl,M., Noubibou,O., Haase,D., Yang,L., Gundlach,H., Hindemitt, T., Klee,K., Haberer,G., Schoof,H. and Mayer,K.F. (2007) MIPSPlantsDB-plant database resource for integrative and comparative plant genome research. Nucleic Acids Res., 35, D834-D840] and was set up from the start to provide data and information resources for individual plant species as well as a framework for integrative and comparative plant genome research. PlantsDB comprises database instances for tomato, Medicago, Arabidopsis, Brachypodium, Sorghum, maize, rice, barley and wheat. Building up on that, state-of-the-art comparative genomics tools such as CrowsNest are integrated to visualize and investigate syntenic relationships between monocot genomes. Results from novel genome analysis strategies targeting the complex and repetitive genomes of triticeae species (wheat and barley) are provided and cross-linked with model species. The MIPS Repeat Element Database (mips-REdat) and Catalog (mips-REcat) as well as tight connections to other databases, e.g. via web services, are further important components of PlantsDB.

Barley whole exome capture: a tool for genomic research in the genus Hordeum and beyond.

Advanced resources for genome-assisted research in barley (Hordeum vulgare) including a whole-genome shotgun assembly and an integrated physical map have recently become available. These have made possible studies that aim to assess genetic diversity or to isolate single genes by whole-genome resequencing and in silico variant detection. However such an approach remains expensive given the 5 Gb size of the barley genome. Targeted sequencing of the mRNA-coding exome reduces barley genomic complexity more than 50-fold, thus dramatically reducing this heavy sequencing and analysis load. We have developed and employed an in-solution hybridization-based sequence capture platform to selectively enrich for a 61.6 megabase coding sequence target that includes predicted genes from the genome assembly of the cultivar Morex as well as publicly available full-length cDNAs and de novo assembled RNA-Seq consensus sequence contigs. The platform provides a highly specific capture with substantial and reproducible enrichment of targeted exons, both for cultivated barley and related species. We show that this exome capture platform provides a clear path towards a broader and deeper understanding of the natural variation residing in the mRNA-coding part of the barley genome and will thus constitute a valuable resource for applications such as mapping-by-sequencing and genetic diversity analyzes.

Ragweed (Ambrosia artemisiifolia) pollen allergenicity: SuperSAGE transcriptomic analysis upon elevated CO2 and drought stress.

BACKGROUND: Pollen of common ragweed (Ambrosia artemisiifolia) is a main cause of allergic diseases in Northern America. The weed has recently become spreading as a neophyte in Europe, while climate change may also affect the growth of the plant and additionally may also influence pollen allergenicity. To gain better insight in the molecular mechanisms in the development of ragweed pollen and its allergenic proteins under global change scenarios, we generated SuperSAGE libraries to identify differentially expressed transcripts. RESULTS: Ragweed plants were grown in a greenhouse under 380 ppm CO2 and under elevated level of CO2 (700 ppm). In addition, drought experiments under both CO2 concentrations were performed. The pollen viability was not altered under elevated CO2, whereas drought stress decreased its viability. Increased levels of individual flavonoid metabolites were found under elevated CO2 and/or drought. Total RNA was isolated from ragweed pollen, exposed to the four mentioned scenarios and four SuperSAGE libraries were constructed. The library dataset included 236,942 unique sequences, showing overlapping as well as clear differently expressed sequence tags (ESTs). The analysis targeted ESTs known in Ambrosia, as well as in pollen of other plants. Among the identified ESTs, those encoding allergenic ragweed proteins (Amb a) increased under elevated CO2 and drought stress. In addition, ESTs encoding allergenic proteins in other plants were also identified. CONCLUSIONS: The analysis of changes in the transcriptome of ragweed pollen upon CO2 and drought stress using SuperSAGE indicates that under global change scenarios the pollen transcriptome was altered, and impacts the allergenic potential of ragweed pollen.

chromoWIZ: a web tool to query and visualize chromosome-anchored genes from cereal and model genomes.

BACKGROUND: Over the last years reference genome sequences of several economically and scientifically important cereals and model plants became available. Despite the agricultural significance of these crops only a small number of tools exist that allow users to inspect and visualize the genomic position of genes of interest in an interactive manner. DESCRIPTION: We present chromoWIZ, a web tool that allows visualizing the genomic positions of relevant genes and comparing these data between different plant genomes. Genes can be queried using gene identifiers, functional annotations, or sequence homology in four grass species (Triticum aestivum, Hordeum vulgare, Brachypodium distachyon, Oryza sativa). The distribution of the anchored genes is visualized along the chromosomes by using heat maps. Custom gene expression measurements, differential expression information, and gene-to-group mappings can be uploaded and can be used for further filtering. CONCLUSIONS: This tool is mainly designed for breeders and plant researchers, who are interested in the location and the distribution of candidate genes as well as in the syntenic relationships between different grass species. chromoWIZ is freely available and online accessible at http://mips.helmholtz-muenchen.de/plant/chromoWIZ/index.jsp.

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

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

Critical role of PI3K signaling for NF-kappaB-dependent survival in a subset of activated B-cell-like diffuse large B-cell lymphoma cells.

The activated B-cell-like (ABC) subtype of diffuse large B-cell lymphoma (DLBCL) represents a very aggressive human lymphoma entity. Constitutive NF-κB activation caused by chronic active B-cell receptor (BCR) signaling is common feature of many ABC DLBCL cells; however, the pathways linking BCR signaling to the NF-κB prosurvival network are largely unknown. Here we report that constitutive activity of PI3K and the downstream kinase PDK1 are essential for the viability of two ABC DLBCL cell lines that carry mutations in the BCR proximal signaling adaptor CD79B. In these cells, PI3K inhibition reduces NF-κB activity and decreases the expression of NF-κB target genes. Furthermore, PI3K and PDK1 are required for maintaining MALT1 protease activity, which promotes survival of the affected ABC DLBCL cells. These results demonstrate a critical function of PI3K-PDK1 signaling upstream of MALT1 protease and NF-κB in distinct ABC DLBCL cells and provide a rationale for the pharmacologic use of PI3K inhibitors in DLBCL therapy.

Overcoming osimertinib resistance with AKT inhibition in EGFRm-driven Non-Small-Cell-Lung-Cancer with PIK3CA/PTEN alterations.

PURPOSE: Osimertinib is an epidermal growth factor receptor tyrosine kinase inhibitor (EGFR-TKI) indicated for the treatment of EGFR mutated (EGFRm)-driven lung adenocarcinomas. Osimertinib significantly improves progression-free survival in first-line treated patients with EGFRm advanced NSCLC. Despite the durable disease control, the majority of patients receiving osimertinib eventually develop disease progression. EXPERIMENTAL DESIGN: ctDNA profiling analysis on-progression plasma samples from patients treated with osimertinib in both first (Phase 3, FLAURA trial) and second-line trials (Phase 3, AURA3 trial) revealed a high prevalence of PIK3CA/AKT/PTEN alterations. In vitro and in vivo evidence using CRISPR engineered NSCLC cell lines and PXD models support a functional role for PIK3CA and PTEN mutations in the development of osimertinib resistance. RESULTS: These alterations are functionally relevant as EGFRm NSCLC cells with engineered PIK3CA/AKT/PTEN alterations develop resistance to osimertinib and can be re-sensitized by treatment with the combination of osimertinib and the AKT inhibitor capivasertib. Moreover, xenograft and PDX in vivo models with PIK3CA/AKT/PTEN alterations display limited sensitivity to osimertinib relative to models without alteration, and in these double mutant models capivasertib and osimertinib combination elicits an improved anti-tumor effect versus osimertinib alone. CONCLUSIONS: Together, this approach offers a potential treatment strategy for patients with EGFRm-driven NSCLC that have a sub-optimal response, or develop resistance, to osimertinib through PIK3CA/AKT/PTEN alterations.

Genome-wide CRISPR screens identify the YAP/TEAD axis as a driver of persister cells in EGFR mutant lung cancer.

Most lung cancer patients with metastatic cancer eventually relapse with drug-resistant disease following treatment and EGFR mutant lung cancer is no exception. Genome-wide CRISPR screens, to either knock out or overexpress all protein-coding genes in cancer cell lines, revealed the landscape of pathways that cause resistance to the EGFR inhibitors osimertinib or gefitinib in EGFR mutant lung cancer. Among the most recurrent resistance genes were those that regulate the Hippo pathway. Following osimertinib treatment a subpopulation of cancer cells are able to survive and over time develop stable resistance. These 'persister' cells can exploit non-genetic (transcriptional) programs that enable cancer cells to survive drug treatment. Using genetic and pharmacologic tools we identified Hippo signalling as an important non-genetic mechanism of cell survival following osimertinib treatment. Further, we show that combinatorial targeting of the Hippo pathway and EGFR is highly effective in EGFR mutant lung cancer cells and patient-derived organoids, suggesting a new therapeutic strategy for EGFR mutant lung cancer patients.

Next-generation sequencing and syntenic integration of flow-sorted arms of wheat chromosome 4A exposes the chromosome structure and gene content.

Wheat is the third most important crop for human nutrition in the world. The availability of high-resolution genetic and physical maps and ultimately a complete genome sequence holds great promise for breeding improved varieties to cope with increasing food demand under the conditions of changing global climate. However, the large size of the bread wheat (Triticum aestivum) genome (approximately 17 Gb/1C) and the triplication of genic sequence resulting from its hexaploid status have impeded genome sequencing of this important crop species. Here we describe the use of mitotic chromosome flow sorting to separately purify and then shotgun-sequence a pair of telocentric chromosomes that together form chromosome 4A (856 Mb/1C) of wheat. The isolation of this much reduced template and the consequent avoidance of the problem of sequence duplication, in conjunction with synteny-based comparisons with other grass genomes, have facilitated construction of an ordered gene map of chromosome 4A, embracing ≥85% of its total gene content, and have enabled precise localization of the various translocation and inversion breakpoints on chromosome 4A that differentiate it from its progenitor chromosome in the A genome diploid donor. The gene map of chromosome 4A, together with the emerging sequences of homoeologous wheat chromosome groups 4, 5 and 7, represent unique resources that will allow us to obtain new insights into the evolutionary dynamics between homoeologous chromosomes and syntenic chromosomal regions.

A transcriptome map of perennial ryegrass (Lolium perenne L.).

BACKGROUND: Single nucleotide polymorphisms (SNPs) are increasingly becoming the DNA marker system of choice due to their prevalence in the genome and their ability to be used in highly multiplexed genotyping assays. Although needed in high numbers for genome-wide marker profiles and genomics-assisted breeding, a surprisingly low number of validated SNPs are currently available for perennial ryegrass. RESULTS: A perennial ryegrass unigene set representing 9,399 genes was used as a reference for the assembly of 802,156 high quality reads generated by 454 transcriptome sequencing and for in silico SNP discovery. Out of more than 15,433 SNPs in 1,778 unigenes fulfilling highly stringent assembly and detection parameters, a total of 768 SNP markers were selected for GoldenGate genotyping in 184 individuals of the perennial ryegrass mapping population VrnA, a population being previously evaluated for important agronomic traits. A total of 592 (77%) of the SNPs tested were successfully called with a cluster separation above 0.9. Of these, 509 (86%) genic SNP markers segregated in the VrnA mapping population, out of which 495 were assigned to map positions. The genetic linkage map presented here comprises a total of 838 DNA markers (767 gene-derived markers) and spans 750 centi Mogan (cM) with an average marker interval distance of less than 0.9 cM. Moreover, it locates 732 expressed genes involved in a broad range of molecular functions of different biological processes in the perennial ryegrass genome. CONCLUSIONS: Here, we present an efficient approach of using next generation sequencing (NGS) data for SNP discovery, and the successful design of a 768-plex Illumina GoldenGate genotyping assay in a complex genome. The ryegrass SNPs along with the corresponding transcribed sequences represent a milestone in the establishment of genetic and genomics resources available for this species and constitute a further step towards molecular breeding strategies. Moreover, the high density genetic linkage map predominantly based on gene-associated DNA markers provides an important tool for the assignment of candidate genes to quantitative trait loci (QTL), functional genomics and the integration of genetic and physical maps in perennial ryegrass, one of the most important temperate grassland species.