Search

During myogenesis the satellite cells undergo a series of highly ordered events that include division and proliferation of mononucleated myoblasts, irreversible cell cycle withdrawal, elongation of myoblasts and finally cell-cell fusion, leading to the formation of multinucleated myotubes. These cellular processes are strictly regulated in a spatiotemporally manner through a coordinated patterns of gene expression governed by specific transcriptional activators and networks of signal transduction. Here we analyzed the genome-wide changes in gene transcription occurring during growth and differentiation of porcine primary myoblast cells. Analysis for differential expression revealed 1656 gene transcripts with significant changes in expression levels during myogenesis. These analyses expand our knowledge of the temporal expression patterns and provide an insight into known and unknown genes and pathways that are transcriptionally regulated and, therefore, likely to be critically implicated in myogenesis. Intriguingly, a remarkable number of the differentially expressed genes encode proteins engaged in histone modification and remodeling of nucleosome structure, emphasizing that extensive changes in chromatin structure are required for regulated gene transcription in the transition from proliferating myoblasts to terminally differentiated myotubes.

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

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.

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. Electronic supplementary material The online version of this article (10.1186/s12862-019-1353-z) contains supplementary material, which is available to authorized users.

Rennet-induced milk coagulation is an important trait for cheese production. Recent studies have reported an alarming frequency of cows producing poorly coagulating milk unsuitable for cheese production. Several genetic factors are known to affect milk coagulation, including variation in the major milk proteins; however, recent association studies indicate genetic effects from other genomic regions as well. The aim of this study was to detect genetic variation affecting milk coagulation properties, measured as curd-firming rate (CFR) and milk pH. This was achieved by examining allele frequency differences between pooled whole-genome sequences of phenotypically extreme samples (pool-seq).. Curd-firming rate and raw milk pH were measured for 415 Danish Holstein cows, and each animal was sequenced at low coverage. Pools were created containing whole genome sequence reads from samples with "extreme" values (high or low) for both phenotypic traits. A total of 6,992,186 and 5,295,501 SNP were assessed in relation to CFR and milk pH, respectively. Allele frequency differences were calculated between pools and 32 significantly different SNP were detected, 1 for milk pH and 31 for CFR, of which 19 are located on chromosome 6. A total of 9 significant SNP, which were selected based on the possible function of proximal candidate genes, were genotyped in the entire sample set ( = 415) to test for an association. The most significant SNP was located proximal to , explaining 33% of the phenotypic variance. , coding for κ-casein, is the most studied in relation to milk coagulation due to its position on the surface of the casein micelles and the direct involvement in milk coagulation. Three additional SNP located on chromosome 6 showed significant associations explaining 7, 3.6, and 1.3% of the phenotypic variance of CFR. The significant SNP on chromosome 6 were shown to be in linkage disequilibrium with the SNP peaking proximal to ; however, after accounting for the genotype of the peak SNP within this QTL, significant effects (-value < 0.1) could still be detected for 2 of the SNP accounting for 2 and 1% of the phenotypic variance. These 2 interesting SNP were located within introns or proximal to the candidate genes-solute carrier family 4 (sodium bicarbonate cotransporter), member 4 () and LIM and calponin homology domains 1 (), respectively-making them interesting targets for further analysis.

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

peer-reviewed H.D.D., A.J.C., P.J.B. and B.J.H. would like to acknowledge the Dairy Futures Cooperative Research Centre for funding. H.P. and R.F. acknowledge funding from the German Federal Ministry of Education and Research (BMBF) within the AgroClustEr ‘Synbreed—Synergistic Plant and Animal Breeding’ (grant 0315527B). H.P., R.F., R.E. and K.-U.G. acknowledge the Arbeitsgemeinschaft Süddeutscher Rinderzüchter, the Arbeitsgemeinschaft Österreichischer Fleckviehzüchter and ZuchtData EDV Dienstleistungen for providing genotype data. A. Bagnato acknowledges the European Union (EU) Collaborative Project LowInputBreeds (grant agreement 222623) for providing Brown Swiss genotypes. Braunvieh Schweiz is acknowledged for providing Brown Swiss phenotypes. H.P. and R.F. acknowledge the German Holstein Association (DHV) and the Confederación de Asociaciones de Frisona Española (CONCAFE) for sharing genotype data. H.P. was financially supported by a postdoctoral fellowship from the Deutsche Forschungsgemeinschaft (DFG) (grant PA 2789/1-1). D.B. and D.C.P. acknowledge funding from the Research Stimulus Fund (11/S/112) and Science Foundation Ireland (14/IA/2576). M.S. and F.S.S. acknowledge the Canadian Dairy Network (CDN) for providing the Holstein genotypes. P.S. acknowledges funding from the Genome Canada project entitled ‘Whole Genome Selection through Genome Wide Imputation in Beef Cattle’ and acknowledges WestGrid and Compute/Calcul Canada for providing computing resources. J.F.T. was supported by the National Institute of Food and Agriculture, US Department of Agriculture, under awards 2013-68004-20364 and 2015-67015-23183. A. Bagnato, F.P., M.D. and J.W. acknowledge EU Collaborative Project Quantomics (grant 516 agreement 222664) for providing Brown Swiss and Finnish Ayrshire sequences and genotypes. A.C.B. and R.F.V. acknowledge funding from the public–private partnership ‘Breed4Food’ (code BO-22.04-011- 001-ASG-LR) and EU FP7 IRSES SEQSEL (grant 317697). A.C.B. and R.F.V. acknowledge CRV (Arnhem, the Netherlands) for providing data on Dutch and New Zealand Holstein and Jersey bulls. Stature is affected by many polymorphisms of small effect in humans1. 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.

Here we report the draft genome sequence of perennial ryegrass (Lolium perenne), an economically important forage and turf grass species 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 utilize 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 (SI) mechanism enables the pistil of a plant to reject self-pollen and therefore promote outcrossing. We have used the sequence assembly to characterise transcriptional changes in the stigma during pollination with both compatible and incompatible pollen. Characterisation 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 potentially involved in the SI mechanism. The genome sequence provides a valuable resource for future breeding efforts based on genomic prediction, and will accelerate the development of varieties for more productive grasslands. 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.

SummaryPlasminogen activator inhibitor type 1 (PAI-1) is a central regulator of fibrinolysis and tissue remodelling. PAI-1 belongs to the serpin super-family 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.Supplementary Material to this article is available online at www.thrombosis-online.com.

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.

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.

Symbiotic associations between bacteria and animal hosts shape the ecology and evolution of both partners (Dubilier et al. 2008; Moran et al. 2008; Moya et al. 2008). Transition of free-living bacteria to obligate endosymbionts profoundly affects bacterial genome evolution (Moran et al. 2008; Moya et al. 2008; Bright and Bulgheresi 2010; Toft and Andersson 2010). Periodic population bottlenecks during vertical transmission magnify the effects of genetic drift. A life cycle spent in the stable and protected host environment also reduces the scope for genetic mixing. These factors lessen the efficacy of natural selection and explain genomic erosion, that is, accelerated substitution rates, biased nucleotide base composition, pseudogenization, and gene loss in symbiont genomes. Genome erosion is best documented for obligate intracellular endosymbionts of insects (Moran et al. 2008, 2009; Toft and Andersson 2010) and has become the paradigm for genome evolution of such heritable symbionts (Andersson and Kurland 1998). In contrast, little (and conflicting) information is available for extracellular bacterial endosymbionts (Dubilier et al. 2008); while those of stinkbugs show vertical transmission, cospeciation, and reductive genome evolution in line with intracellular insect endosymbionts (Hosokawa et al. 2006; Kikuchi et al. 2009), the genomes of gutless oligochaete symbionts are similar in size and G + C content to their free-living relatives; however, their mode of transmission has not been fully resolved (Woyke et al. 2006). The genome of the extracellular earthworm endosymbiont Verminephrobacter eiseniae (Pinel et al. 2008) shows no signs of genome reduction either (Pinel 2009). Verminephrobacter (Betaproteobacteria) are species specific, extracellular endosymbionts of most lumbricid earthworm species (Lund et al. 2010). They colonize the earthworm's excretion organs, nephridia (Schramm et al. 2003), and are transmitted vertically: symbionts are deposited along with eggs and sperm in the egg capsules (Davidson and Stahl 2006) and during embryogenesis, selectively recruited by the earthworm (Davidson and Stahl 2008). Fitness experiments have shown that the symbionts increase the reproductive success of their host (Lund, Holmstrup, et al. 2010); the earthworm–Verminephrobacter symbiosis is therefore considered mutually beneficial. Although host and symbiont can be separated under laboratory conditions (Lund, Holmstrup, et al. 2010), the host-associated lifestyle of Verminephrobacter is evolutionary stable and as old as that of many heritable intracellular endosymbionts of insects (Moran et al. 2008). Strict host fidelity and a pattern of cospeciation indicate an origin of the Verminephrobacter–earthworm symbiosis in the most recent common ancestor of lumbricid earthworms, 60–130 Myr (Lund et al. 2010). Here, we obtain a draft genome from a second Verminephrobacter species (Lund et al. 2012). We contrast the strength of purifying selection and the level of positive selection in the two Verminephrobacter genomes with those found in two free-living bacteria from the sister genus Acidovorax to study how the ancient transition to a symbiotic lifestyle affected genome evolution in these extracellular heritable endosymbionts.

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model., The research leading to these results has received funding from the European Research Council under the European Community’s Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. 249894 (SelSweep)

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. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2249-y) contains supplementary material, which is available to authorized users.

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%.

The extensive use of a limited number of elite bulls in cattle breeding can lead to rapid spread of recessively inherited disorders. A recent example is the globally distributed syndrome Complex Vertebral Malformation (CVM), which is characterized by misshapen and fused vertebrae around the cervico-thoracic junction. Here, we show that CVM is caused by a mutation in the Golgi-resident nucleotide-sugar transporter encoded by SLC35A3. Thus, the disease showed complete cosegregation with the mutation in a homozygous state, and proteome patterns indicated abnormal protein glycosylation in tissues of affected animals. In addition, a yeast mutant that is deficient in the transport of UDP-N-acetylglucosamine into its Golgi lumen can be rescued by the wild-type SLC35A3 gene, but not by the mutated gene. These results provide the first demonstration of a genetic disorder associated with a defective SLC35A3 gene, and reveal a new mechanism for malformation of the vertebral column caused by abnormal nucleotide-sugar transport into the Golgi apparatus.

Next-generation sequencing (NGS) technologies now offer the opportunity to move from population genetics to population genomics, even in non-model species. In the present study, two NGS platforms, the Roche 454 FLX and the Illumina Genome Analyzer were used to discover Single Nucleotide Polymorphisms (SNPs) in the muscle transcriptome of the European hake (Merluccius merluccius), which represents one of the most important demersal fisheries in Europe. Transcriptome sequencing and de novo assembly into unique contigs, in silico SNP detection, and validation by high-throughput genotyping yielded a set of 395 polymorphic SNP loci. Of these, 381 nuclear SNPs were used to scan genomic variation in 19 Atlantic and Mediterranean populations of M. merluccius. As all SNPs reside in the functional (transcribed) portion of the hake genome, it is expected that a subset of sequence variants might be located in genomic regions potentially involved in local adaptation, conveying greater resolving power for the detection of spatial structuring on a large and/or local scale. Combining genome scans and correlation analysis between SNP allele frequencies and environmental factors, several candidate loci that are potentially under divergent selection were identified, mainly genes involved in muscles structural organization/function and protein synthesis regulation, or implicated in energy metabolism. Such “outlier” loci showed high performance in resolving population structure, unveiling distinct genetic clusters at multiple spatial scales, which possibly correspond to locally adapted populations. Following a precautionary approach, these results suggest the need to revise the current management strategies of European hake stocks, particularly within basins.

Identification of genetic variation in cattle breeds using next-generation sequencing technology has focused on the modern production cattle breeds. We focused on one of the oldest indigenous breeds, the Icelandic cattle breed.Sequencing of two individuals enabled identification of more than 8 million SNPs and more than one million short indels. Annotation of the genetic variants identified a substantial number of functional SNPs and variants. The number of genetic variants identified in the Icelandic cattle breed is on the same level as previously seen in other studies on Holstein cattle

We sequenced the whole-genome of a Danish Jutland bull to identify genetic variants (SNP/indel). Using UnifiedGenotyper from the Genome Analysis Toolkit (GATK), we identified 6,812,198 SNPs and 804,453 indels. There were 2,598,000 (38.1%) novel SNPs and 607,923(75.6%) novel indels while the remaining was annotated in dbSNP build 133. In-depth annotation of the variants revealed that 45,776 SNPs affected the coding sequences of 11,538 genes, 221 SNPs predicted to cause a premature stop codon, 17 to cause a gain in coding sequence and 20,828 predicted to be non-synonymous. We identified 1,122 indels in coding sequences, 832 predicted to cause frame shift, 89 predicted to be inframe insertion and 115 to be inframe deletion. We detected a higher level of genetic variation in the Jutland bull compared to similar data from Holstein cattle

CNVs are increasingly recognized as substantial source of genetic variation, fueling studies that assess their impact on complex traits. In particular rareCNVs have been suggested to potentially explain part of the missing heritability problem in genome wide association studies for complex traits. The objective of this study was to perform a high resolution genome scan for CNV, in a sample of 20 Brown Swiss dairy cattle bulls based on ~20x Illumina whole genome sesequencing data. Employing CNVnator for variant discovery, we present descriptive statistics for the CNVs detected and define consensus CNV regions at the population level. We identified 29,975 deletion-, 1,489 duplication- and 365 complex CNVRs, respectively, which cover 3.3% of the UMD3.1 autosome. We further compared NGS based CNV calls to CNV calls detected by PennCNV based on Illumina HD chip data for 17 bulls with high quality data for both platforms

Mastitis is a mammary disease that frequently affects dairy cattle. Despite considerable research on the development of effective prevention and treatment strategies, mastitis continues to be a significant issue in bovine veterinary medicine. To identify major genes that affect mastitis in dairy cattle, 6 chromosomal regions on Bos taurus autosome (BTA) 6, 13, 16, 19, and 20 were selected from a genome scan for 9 mastitis phenotypes using imputed high-density single nucleotide polymorphism arrays. Association analyses using sequence-level variants for the 6 targeted regions were carried out to map causal variants using whole-genome sequence data from 3 breeds. The quantitative trait loci (QTL) discovery population comprised 4,992 progeny-tested Holstein bulls, and QTL were confirmed in 4,442 Nordic Red and 1,126 Jersey cattle. The targeted regions were imputed to the sequence level. The highest association signal for clinical mastitis was observed on BTA 6 at 88.97 Mb in Holstein cattle and was confirmed in Nordic Red cattle. The peak association region on BTA 6 contained 2 genes: vitamin D-binding protein precursor (GC) and neuropeptide FF receptor 2 (NPFFR2), which, based on known biological functions, are good candidates for affecting mastitis. However, strong linkage disequilibrium in this region prevented conclusive determination of the causal gene. A different QTL on BTA 6 located at 88.32 Mb in Holstein cattle affected mastitis. In addition, QTL on BTA 13 and 19 were confirmed to segregate in Nordic Red cattle and QTL on BTA 16 and 20 were confirmed in Jersey cattle. Although several candidate genes were identified in these targeted regions, it was not possible to identify a gene or polymorphism as the causal factor for any of these regions.

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

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., PLoS ONE, 8 (3), ISSN:1932-6203

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.

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.

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.

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of nonmodel 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, micro-array development and for population genomic studies in the wild. ispartof: PLoS One vol:7 issue:8 pages:1-11 ispartof: location:United States status: published

The development in sequencing methods has made it possible to perform whole genome de novo sequencing of species without large commercial interests. Within the EU-financed QUANTOMICS project (KBBE-2A-222664), we have performed de novo sequencing of quail (Coturnix coturnix) and grey partridge (Perdix perdix) on a Genome Analyzer GAII (Illumina) using paired-end sequencing. The amount of generated sequences amounts to 8 to 9 Gb for each species. The analysis and assembly of the generated sequences is ongoing. Access to the whole genome sequence from these two species will enable enhanced comparative studies towards the chicken genome and will aid in identifying evolutionarily conserved sequences within the Galliformes. The obtained sequences from quail and partridge represent a beginning of generating the whole genome sequence for these species. The continuation of establishing the genome sequence of these species will require generation of more sequences from more diverse libraries. The QUANTOMICS group encourages other research groups to participate in a joint effort to establish the genome sequence of these two Galliformes species.

Transcriptome characterisation using next generation sequencing allows the global description of the genes expressed in a tissue or organ, the discovery of novel genes or alternative splicing events. In addition we can identify sequence variation (SNPs) and get information about transcript abundance. Our objective was to investigate animals previously not affected by lung disease and those that had been affected. To this end lung tissue samples were collected, separately pooled and tagged before sequencing using the Roche/454 FLX platform. We sequenced about one million reads that were clustered and mapped to the current pig genome reference sequence. Identified genes or clusters were annotated for functional classes and mined for singe nucleotide polymorphisms. In addition, we compared gene expression between sample groups in order to investigate possible changes in the lung transcriptome that may be related to previous infection

Actinobacillus pleuropneumoniae is a bacterial pathogen that causes highly contagious respiratory infection in pigs and has serious impact on production economy and animal welfare. As clear differences in virulence between serotypes have been observed, the genetic basis should be investigated at the genomic level. Here we present the draft genome sequences of the A. pleuropneumoniae serotypes 2 (strain 4226) and 6 (strain Femo) Udgivelsesdato: November Actinobacillus pleuropneumoniae is a bacterial pathogen that causes highly contagious respiratory infection in pigs and has serious impact on production economy and animal welfare. As clear differences in virulence between serotypes have been observed, the genetic basis should be investigated at the genomic level. Here we present the draft genome sequences of the A. pleuropneumoniae serotypes 2 (strain 4226) and 6 (strain Femo)

Actinobacillus pleuropneumoniae (Ap) is a gram-negative bacterium that causes porcine pleuropneumonia, which is a widespread, highly contagious and often fatal respiratory disease in swine. A total of 44 pigs were experimentally inoculated with Ap serotype 2 or 6 and samples of liver and tracheobronchial lung lymph nodes were collected 6, 12, 24 and 48 hours after experimental inoculation, as well as from six non-inoculated control pigs. Transcriptional profiles of the liver samples have been generated by preparation of 12-plexed mRNA-Seq libraries followed by sequencing on an Illumina GAIIx (51+7 cycles) obtaining more than 200 million tag sequences. The 12-plexed mRNA-Seq libraries of the lung lymph node samples have presently (April 2010) been prepared and are to be sequenced. The PCR amplicons of the liver libraries were quantified using both a fluorometer and a qPCR assay, including the use of a sequence-titrated, in-house control library. The libraries were diluted to 6 pM based on the qPCR assay, except for a single library set which was duplicated and diluted based on the fluorometer measurements as well. Analysis of the obtained sequences revealed that the qPCR based quantifications reduced the cluster density variability as compared to fluorometer based quantifications. Furthermore, it was found that the fluorometer based measurements tended to deviate for dilute as well as for more concentrated libraries. Following the sequencing of the lung lymph node samples analyses are to be conducted to study the time and serotype dependent transcriptional response to Ap infection.