Genes and pathways revealed by whole transcriptome analysis of milk derived bovine mammary epithelial cells after Escherichia coli challenge.

Mastitis, inflammation of the mammary gland, is the costliest disease in dairy cattle and a major animal welfare concern. Mastitis is usually caused by bacteria, of which staphylococci, streptococci and Escherichia coli are most frequently isolated from bovine mastitis. Bacteria activate the mammary immune system in variable ways, thereby influencing the severity of the disease. Escherichia coli is a common cause of mastitis in cattle causing both subclinical and clinical mastitis. Understanding of the molecular mechanisms that activate and regulate the host response would be central to effective prevention of mastitis and breeding of cows more resistant to mastitis. We used primary bovine mammary epithelial cell cultures extracted noninvasively from bovine milk samples to monitor the cellular responses to Escherichia coli challenge. Differences in gene expression between control and challenged cells were studied by total RNA-sequencing at two time points post-challenge. In total, 150 and 440 (Padj < 0.05) differentially expressed genes were identified at 3 h and 24 h post-challenge, respectively. The differentially expressed genes were mostly upregulated at 3 h (141/150) and 24 h (424/440) post-challenge. Our results are in line with known effects of E. coli infection, with a strong early inflammatory response mediated by pathogen receptor families. Among the most significantly enriched early KEGG pathways were the TNF signalling pathway, the cytokine-cytokine receptor interaction, and the NF-kappa B signalling pathway. At 24 h post-challenge, most significantly enriched were the Influenza A, the NOD-like receptor signalling, and the IL-17 signaling pathway.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Biochemical and structural analyses suggest that plasminogen activators coevolved with their cognate protein substrates and inhibitors.

Protein sequences of members of the plasminogen activation system are present throughout the entire vertebrate phylum. This important and well-described proteolytic cascade is governed by numerous protease-substrate and protease-inhibitor interactions whose conservation is crucial to maintaining unchanged protein function throughout evolution. The pressure to preserve protein-protein interactions may lead to either co-conservation or covariation of binding interfaces. Here, we combined covariation analysis and structure-based prediction to analyze the binding interfaces of urokinase (uPA):plasminogen activator inhibitor-1 (PAI-1) and uPA:plasminogen complexes. We detected correlated variation between the S3-pocket-lining residues of uPA and the P3 residue of both PAI-1 and plasminogen. These residues are known to form numerous polar interactions in the human uPA:PAI-1 Michaelis complex. To test the effect of mutations that correlate with each other and have occurred during mammalian diversification on protein-protein interactions, we produced uPA, PAI-1, and plasminogen from human and zebrafish to represent mammalian and nonmammalian orthologs. Using single amino acid point substitutions in these proteins, we found that the binding interfaces of uPA:plasminogen and uPA:PAI-1 may have coevolved to maintain tight interactions. Moreover, we conclude that although the interaction areas between protease-substrate and protease-inhibitor are shared, the two interactions are mechanistically different. Compared with a protease cleaving its natural substrate, the interaction between a protease and its inhibitor is more complex and involves a more fine-tuned mechanism. Understanding the effects of evolution on specific protein interactions may help further pharmacological interventions of the plasminogen activation system and other proteolytic systems.

Origin and diversification of the plasminogen activation system among chordates.

BACKGROUND: The plasminogen (PLG) activation system is composed by a series of serine proteases, inhibitors and several binding proteins, which together control the temporal and spatial generation of the active serine protease plasmin. As this proteolytic system plays a central role in human physiology and pathophysiology it has been extensively studied in mammals. The serine proteases of this system are believed to originate from an ancestral gene by gene duplications followed by domain gains and deletions. However, the identification of ancestral forms in primitive chordates supporting these theories remains elusive. In addition, evolutionary studies of the non-proteolytic members of this system are scarce. RESULTS: Our phylogenetic analyses place lamprey PLG at the root of the vertebrate PLG-group, while lamprey PLG-related growth factors represent the ancestral forms of the jawed-vertebrate orthologues. Furthermore, we find that the earliest putative orthologue of the PLG activator group is the hyaluronan binding protein 2 (HABP2) gene found in lampreys. The prime plasminogen activators (tissue- and urokinase-type plasminogen activator, tPA and uPA) first occur in cartilaginous fish and phylogenetic analyses confirm that all orthologues identified compose monophyletic groups to their mammalian counterparts. Cartilaginous fishes exhibit the most ancient vitronectin of all vertebrates, while plasminogen activator inhibitor 1 (PAI-1) appears for the first time in cartilaginous fishes and is conserved in the rest of jawed vertebrate clades. PAI-2 appears for the first time in the common ancestor of reptiles and mammals, and represents the latest appearing plasminogen activator inhibitor. Finally, we noted that the urokinase-type plasminogen activator receptor (uPAR)-and three-LU domain containing genes in general-occurred later in evolution and was first detectable after coelacanths. CONCLUSIONS: This study identifies several primitive orthologues of the mammalian plasminogen activation system. These ancestral forms provide clues to the origin and diversification of this enzyme system. Further, the discovery of several members-hitherto unknown in mammals-provide new perspectives on the evolution of this important enzyme system.

MicroRNA-guided prioritization of genome-wide association signals reveals the importance of microRNA-target gene networks for complex traits in cattle.

MicroRNAs (miRNA) are key modulators of gene expression and so act as putative fine-tuners of complex phenotypes. Here, we hypothesized that causal variants of complex traits are enriched in miRNAs and miRNA-target networks. First, we conducted a genome-wide association study (GWAS) for seven functional and milk production traits using imputed sequence variants (13~15 million) and >10,000 animals from three dairy cattle breeds, i.e., Holstein (HOL), Nordic red cattle (RDC) and Jersey (JER). Second, we analyzed for enrichments of association signals in miRNAs and their miRNA-target networks. Our results demonstrated that genomic regions harboring miRNA genes were significantly (P < 0.05) enriched with GWAS signals for milk production traits and mastitis, and that enrichments within miRNA-target gene networks were significantly higher than in random gene-sets for the majority of traits. Furthermore, most between-trait and across-breed correlations of enrichments with miRNA-target networks were significantly greater than with random gene-sets, suggesting pleiotropic effects of miRNAs. Intriguingly, genes that were differentially expressed in response to mammary gland infections were significantly enriched in the miRNA-target networks associated with mastitis. All these findings were consistent across three breeds. Collectively, our observations demonstrate the importance of miRNAs and their targets for the expression of complex traits.

Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals.

Stature is affected by many polymorphisms of small effect in humans 1 . In contrast, variation in dogs, even within breeds, has been suggested to be largely due to variants in a small number of genes2,3. Here we use data from cattle to compare the genetic architecture of stature to those in humans and dogs. We conducted a meta-analysis for stature using 58,265 cattle from 17 populations with 25.4 million imputed whole-genome sequence variants. Results showed that the genetic architecture of stature in cattle is similar to that in humans, as the lead variants in 163 significantly associated genomic regions (P < 5 × 10-8) explained at most 13.8% of the phenotypic variance. Most of these variants were noncoding, including variants that were also expression quantitative trait loci (eQTLs) and in ChIP-seq peaks. There was significant overlap in loci for stature with humans and dogs, suggesting that a set of common genes regulates body size in mammals.

The first draft reference genome of the American mink (Neovison vison).

The American mink (Neovison vison) is a semiaquatic species of mustelid native to North America. It's an important animal for the fur industry. Many efforts have been made to locate genes influencing fur quality and color, but this search has been impeded by the lack of a reference genome. Here we present the first draft genome of mink. In our study, two mink individuals were sequenced by Illumina sequencing with 797 Gb sequence generated. Assembly yielded 7,175 scaffolds with an N50 of 6.3 Mb and length of 2.4 Gb including gaps. Repeat sequences constitute around 31% of the genome, which is lower than for dog and cat genomes. The alignments of mink, ferret and dog genomes help to illustrate the chromosomes rearrangement. Gene annotation identified 21,053 protein-coding sequences present in mink genome. The reference genome's structure is consistent with the microsatellite-based genetic map. Mapping of well-studied genes known to be involved in coat quality and coat color, and previously located fur quality QTL provide new knowledge about putative candidate genes for fur traits. The draft genome shows great potential to facilitate genomic research towards improved breeding for high fur quality animals and strengthen our understanding on evolution of Carnivora.

Latency transition of plasminogen activator inhibitor type 1 is evolutionarily conserved.

Plasminogen activator inhibitor type 1 (PAI-1) is a central regulator of fibrinolysis and tissue remodelling. PAI-1 belongs to the serpin superfamily and unlike other inhibitory serpins undergoes a spontaneous inactivation process under physiological conditions, termed latency transition. During latency transition the solvent exposed reactive centre loop is inserted into the central ß-sheet A of the molecule, and is no longer accessible to reaction with the protease. More than three decades of research on mammalian PAI-1 has not been able to clarify the evolutionary advantage and physiological relevance of latency transition. In order to study the origin of PAI-1 latency transition, we produced PAI-1 from Spiny dogfish shark (Squalus acanthias) and African lungfish (Protopterus sp.), which represent central species in the evolution of vertebrates. Although human PAI-1 and the non-mammalian PAI-1 variants share only approximately 50 % sequence identity, our results showed that all tested PAI-1 variants undergo latency transition with a similar rate. Since the functional stability of PAI-1 can be greatly increased by substitution of few amino acid residues, we conclude that the ability to undergo latency transition must have been a specific selection criterion for the evolution of PAI-1. It appears that all PAI-1 molecules must harbour latency transition to fulfil their physiological function, stressing the importance to further pursue a complete understanding of this molecular phenomenon with possible implication to pharmacological intervention. Our results provide the next step in understanding how the complete role of this important protease inhibitor evolved along with the fibrinolytic system.

Multi-tissue RNA-seq and transcriptome characterisation of the spiny dogfish shark (Squalus acanthias) provides a molecular tool for biological research and reveals new genes involved in osmoregulation.

The spiny dogfish shark (Squalus acanthias) is one of the most commonly used cartilaginous fishes in biological research, especially in the fields of nitrogen metabolism, ion transporters and osmoregulation. Nonetheless, transcriptomic data for this organism is scarce. In the present study, a multi-tissue RNA-seq experiment and de novo transcriptome assembly was performed in four different spiny dogfish tissues (brain, liver, kidney and ovary), providing an annotated sequence resource. The characterization of the transcriptome greatly increases the scarce sequence information for shark species. Reads were assembled with the Trinity de novo assembler both within each tissue and across all tissues combined resulting in 362,690 transcripts in the combined assembly which represent 289,515 Trinity genes. BUSCO analysis determined a level of 87% completeness for the combined transcriptome. In total, 123,110 proteins were predicted of which 78,679 and 83,164 had significant hits against the SwissProt and Uniref90 protein databases, respectively. Additionally, 61,215 proteins aligned to known protein domains, 7,208 carried a signal peptide and 15,971 possessed at least one transmembrane region. Based on the annotation, 81,582 transcripts were assigned to gene ontology terms and 42,078 belong to known clusters of orthologous groups (eggNOG). To demonstrate the value of our molecular resource, we show that the improved transcriptome data enhances the current possibilities of osmoregulation research in spiny dogfish by utilizing the novel gene and protein annotations to investigate a set of genes involved in urea synthesis and urea, ammonia and water transport, all of them crucial in osmoregulation. We describe the presence of different gene copies and isoforms of key enzymes involved in this process, including arginases and transporters of urea and ammonia, for which sequence information is currently absent in the databases for this model species. The transcriptome assemblies and the derived annotations generated in this study will support the ongoing research for this particular animal model and provides a new molecular tool to assist biological research in cartilaginous fishes.

Deep sequencing of Danish Holstein dairy cattle for variant detection and insight into potential loss-of-function variants in protein coding genes.

BACKGROUND: Over the last few years, continuous development of high-throughput sequencing platforms and sequence analysis tools has facilitated reliable identification and characterization of genetic variants in many cattle breeds. Deep sequencing of entire genomes within a cattle breed that has not been thoroughly investigated would be imagined to discover functional variants that are underlying phenotypic differences. Here, we sequenced to a high coverage the Danish Holstein cattle breed to detect and characterize single nucleotide polymorphisms (SNPs), insertion/deletions (Indels), and loss-of-function (LoF) variants in protein-coding genes in order to provide a comprehensive resource for subsequent detection of causal variants for recessive traits. RESULTS: We sequenced four genetically unrelated Danish Holstein cows with a mean coverage of 27X using an Illumina Hiseq 2000. Multi-sample SNP calling identified 10,796,794 SNPs and 1,295,036 indels whereof 482,835 (4.5 %) SNPs and 231,359 (17.9 %) indels were novel. A comparison between sequencing-derived SNPs and genotyping from the BovineHD BeadChip revealed a concordance rate of 99.6-99.8 % for homozygous SNPs and 93.3-96.5 % for heterozygous SNPs. Annotation of the SNPs discovered 74,886 SNPs and 1937 indels affecting coding sequences with 2145 being LoF mutations. The frequency of LoF variants differed greatly across the genome, a hot spot with a strikingly high density was observed in a 6 Mb region on BTA18. LoF affected genes were enriched for functional categories related to olfactory reception and underrepresented for genes related to key cellular constituents and cellular and biological process regulation. Filtering using sequence derived genotype data for 288 Holstein animals from the 1000 bull genomes project removing variants containing homozygous individuals retained 345 of the LoF variants as putatively deleterious. A substantial number of the putative deleterious LoF variants had a minor allele frequency >0.05 in the 1000 bull genomes data set. CONCLUSIONS: Deep sequencing of Danish Holstein genomes enabled us to identify 12.1 million variants. An investigation into LoF variants discovered a set of variants predicted to disrupt protein-coding genes. This catalog of variants will be a resource for future studies to understand variation underlying important phenotypes, particularly recessively inherited lethal phenotypes.

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.

Draft Genome Sequence of Clostridium tyrobutyricum Strain DIVETGP, Isolated from Cow's Milk for Grana Padano Production.

We announce the draft genome sequence of Clostridium tyrobutyricum strain DIVETGP. This strain was isolated from cow's milk used for Grana Padano cheese production. The genome was obtained using Illumina HiSeq technology and comprises 45 contigs for 3,018,999 bp, with a G+C content of 30.8%.

A 660-Kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in Nordic Red cattle: additional evidence for the common occurrence of balancing selection in livestock.

In dairy cattle, the widespread use of artificial insemination has resulted in increased selection intensity, which has led to spectacular increase in productivity. However, cow fertility has concomitantly severely declined. It is generally assumed that this reduction is primarily due to the negative energy balance of high-producing cows at the peak of lactation. We herein describe the fine-mapping of a major fertility QTL in Nordic Red cattle, and identify a 660-kb deletion encompassing four genes as the causative variant. We show that the deletion is a recessive embryonically lethal mutation. This probably results from the loss of RNASEH2B, which is known to cause embryonic death in mice. Despite its dramatic effect on fertility, 13%, 23% and 32% of the animals carry the deletion in Danish, Swedish and Finnish Red Cattle, respectively. To explain this, we searched for favorable effects on other traits and found that the deletion has strong positive effects on milk yield. This study demonstrates that embryonic lethal mutations account for a non-negligible fraction of the decline in fertility of domestic cattle, and that associated positive effects on milk yield may account for part of the negative genetic correlation. Our study adds to the evidence that structural variants contribute to animal phenotypic variation, and that balancing selection might be more common in livestock species than previously appreciated.

Integration of mate pair sequences to improve shotgun assemblies of flow-sorted chromosome arms of hexaploid wheat.

BACKGROUND: The assembly of the bread wheat genome sequence is challenging due to allohexaploidy and extreme repeat content (>80%). Isolation of single chromosome arms by flow sorting can be used to overcome the polyploidy problem, but the repeat content cause extreme assembly fragmentation even at a single chromosome level. Long jump paired sequencing data (mate pairs) can help reduce assembly fragmentation by joining multiple contigs into single scaffolds. The aim of this work was to assess how mate pair data generated from multiple displacement amplified DNA of flow-sorted chromosomes affect assembly fragmentation of shotgun assemblies of the wheat chromosomes. RESULTS: Three mate pair (MP) libraries (2 Kb, 3 Kb, and 5 Kb) were sequenced to a total coverage of 89x and 64x for the short and long arm of chromosome 7B, respectively. Scaffolding using SSPACE improved the 7B assembly contiguity and decreased gene space fragmentation, but the degree of improvement was greatly affected by scaffolding stringency applied. At the lowest stringency the assembly N50 increased by ~7 fold, while at the highest stringency N50 was only increased by ~1.5 fold. Furthermore, a strong positive correlation between estimated scaffold reliability and scaffold assembly stringency was observed. A 7BS scaffold assembly with reduced MP coverage proved that assembly contiguity was affected only to a small degree down to ~50% of the original coverage. CONCLUSION: The effect of MP data integration into pair end shotgun assemblies of wheat chromosome was moderate; possibly due to poor contig assembly contiguity, the extreme repeat content of wheat, and the use of amplified chromosomal DNA for MP library construction.

Genome wide allele frequency fingerprints (GWAFFs) of populations via genotyping by sequencing.

Genotyping-by-Sequencing (GBS) is an excellent tool for characterising genetic variation between plant genomes. To date, its use has been reported only for genotyping of single individuals. However, there are many applications where resolving allele frequencies within populations on a genome-wide scale would be very powerful, examples include the breeding of outbreeding species, varietal protection in outbreeding species, monitoring changes in population allele frequencies. This motivated us to test the potential to use GBS to evaluate allele frequencies within populations. Perennial ryegrass is an outbreeding species, and breeding programs are based upon selection on populations. We tested two restriction enzymes for their efficiency in complexity reduction of the perennial ryegrass genome. The resulting profiles have been termed Genome Wide Allele Frequency Fingerprints (GWAFFs), and we have shown how these fingerprints can be used to distinguish between plant populations. Even at current costs and throughput, using sequencing to directly evaluate populations on a genome-wide scale is viable. GWAFFs should find many applications, from varietal development in outbreeding species right through to playing a role in protecting plant breeders' rights.

The genome and transcriptome of perennial ryegrass mitochondria.

BACKGROUND: Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass. RESULTS: Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes - nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 - all showed high expression levels both in absolute numbers and after normalization. CONCLUSIONS: The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass.

Spatio-temporal regulation of ADAR editing during development in porcine neural tissues.

Editing by ADAR enzymes is essential for mammalian life. Still, knowledge of the spatio-temporal editing patterns in mammals is limited. By use of 454 amplicon sequencing we examined the editing status of 12 regionally extracted mRNAs from porcine developing brain encompassing a total of 64 putative ADAR editing sites. In total 24 brain tissues, dissected from up to five regions from embryonic gestation day 23, 42, 60, 80, 100 and 115, were examined for editing.   Generally, editing increased during embryonic development concomitantly with an increase in ADAR2 mRNA level. Notably, the Gria2 (GluR-B) Q/R site, reported to be ~100% edited in previous studies, is only 54% edited at embryonic day 23. Transcripts with multiple editing sites in close proximity to each other exhibit coupled editing and an extraordinary incidence of long-range coupling of editing events more than 32 kb apart is observed for the kainate glutamate receptor 2 transcript, Grik2. Our study reveals complex spatio-temporal ADAR editing patterns of coordinated editing events that may play important roles in the development of the mammalian brain.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus).

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

Illegal, Unreported and Unregulated fishing has had a major role in the overexploitation of global fish populations. In response, international regulations have been imposed and many fisheries have been 'eco-certified' by consumer organizations, but methods for independent control of catch certificates and eco-labels are urgently needed. Here we show that, by using gene-associated single nucleotide polymorphisms, individual marine fish can be assigned back to population of origin with unprecedented high levels of precision. By applying high differentiation single nucleotide polymorphism assays, in four commercial marine fish, on a pan-European scale, we find 93-100% of individuals could be correctly assigned to origin in policy-driven case studies. We show how case-targeted single nucleotide polymorphism assays can be created and forensically validated, using a centrally maintained and publicly available database. Our results demonstrate how application of gene-associated markers will likely revolutionize origin assignment and become highly valuable tools for fighting illegal fishing and mislabelling worldwide.

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.