Genetic variants associated with two major bovine milk fatty acids offer opportunities to breed for altered milk fat composition.

BACKGROUND: Although bovine milk is regarded as healthy and nutritious, its high content of saturated fatty acids (FA) may be harmful to cardiovascular health. Palmitic acid (C16:0) is the predominant saturated FA in milk with adverse health effects that could be countered by substituting it with higher levels of unsaturated FA, such as oleic acid (C18:1cis-9). In this work, we performed genome-wide association analyses for milk fatty acids predicted from FTIR spectroscopy data using 1811 Norwegian Red cattle genotyped and imputed to a high-density 777k single nucleotide polymorphism (SNP)-array. In a follow-up analysis, we used imputed whole-genome sequence data to detect genetic variants that are involved in FTIR-predicted levels of C16:0 and C18:1cis-9 and explore the transcript profile and protein level of candidate genes. RESULTS: Genome-wise significant associations were detected for C16:0 on Bos taurus (BTA) autosomes 11, 16 and 27, and for C18:1cis-9 on BTA5, 13 and 19. Closer examination of a significant locus on BTA11 identified the PAEP gene, which encodes the milk protein ß-lactoglobulin, as a particularly attractive positional candidate gene. At this locus, we discovered a tightly linked cluster of genetic variants in coding and regulatory sequences that have opposing effects on the levels of C16:0 and C18:1cis-9. The favourable haplotype, linked to reduced levels of C16:0 and increased levels of C18:1cis-9 was also associated with a marked reduction in PAEP expression and ß-lactoglobulin protein levels. ß-lactoglobulin is the most abundant whey protein in milk and lower levels are associated with important dairy production parameters such as improved cheese yield. CONCLUSIONS: The genetic variants detected in this study may be used in breeding to produce milk with an improved FA health-profile and enhanced cheese-making properties.

SalmoBase: an integrated molecular data resource for Salmonid species.

BACKGROUND: Salmonids are ray-finned fishes which constitute 11 genera and at least 70 species including Atlantic salmon, whitefishes, graylings, rainbow trout, and char. The common ancestor of all Salmonidae experienced a whole genome duplication (WGD) ~80 million years ago, resulting in an autotetraploid genome. Genomic rediplodization is still going on in salmonid species, providing an unique system for studying evolutionary consequences of whole genome duplication. In recent years, high quality genome sequences of Atlantic salmon and Rainbow trout has been established, due to their scientific and commercial values. In this paper we introduce SalmoBase ( http://www.salmobase.org/ ), a tool for making molecular resources for salmonids public available in a framework of visualizations and analytic tools. RESULTS: SalmoBase has been developed as a part of the ELIXIR.NO project. Currently, SalmoBase contains molecular resources for Atlantic salmon and Rainbow trout. Data can be accessed through BLAST, Genome Browser (GBrowse), Genetic Variation Browser (GVBrowse) and Gene Expression Browser (GEBrowse). CONCLUSIONS: To the best of our knowledge, SalmoBase is the first database which integrates salmonids data and allow users to study salmonids in an integrated framework. The database and its tools (e.g., comparative genomics tools, synteny browsers) will be expanded as additional public resources describing other Salmonidae genomes become available.

Genome-wide association mapping for milk fat composition and fine mapping of a QTL for de novo synthesis of milk fatty acids on bovine chromosome 13.

BACKGROUND: Bovine milk is widely regarded as a nutritious food source for humans, although the effects of individual fatty acids on human health is a subject of debate. Based on the assumption that genomic selection offers potential to improve milk fat composition, there is strong interest to understand more about the genetic factors that influence the biosynthesis of bovine milk and the molecular mechanisms that regulate milk fat synthesis and secretion. For this reason, the work reported here aimed at identifying genetic variants that affect milk fatty acid composition in Norwegian Red cattle. Milk fatty acid composition was predicted from the nation-wide recording scheme using Fourier transform infrared spectroscopy data and applied to estimate heritabilities for 36 individual and combined fatty acid traits. The recordings were used to generate daughter yield deviations that were first applied in a genome-wide association (GWAS) study with 17,343 markers to identify quantitative trait loci (QTL) affecting fatty acid composition, and next on high-density and sequence-level datasets to fine-map the most significant QTL on BTA13 (BTA for Bos taurus chromosome). RESULTS: The initial GWAS revealed 200 significant associations, with the strongest signals on BTA1, 13 and 15. The BTA13 QTL highlighted a strong functional candidate gene for de novo synthesis of short- and medium-chained saturated fatty acids; acyl-CoA synthetase short-chain family member 2. However, subsequent fine-mapping using single nucleotide polymorphisms (SNPs) from a high-density chip and variants detected by resequencing showed that the effect was more likely caused by a second nearby gene; nuclear receptor coactivator 6 (NCOA6). These findings were confirmed with results from haplotype studies. NCOA6 is a nuclear receptor that interacts with transcription factors such as PPARγ, which is a major regulator of bovine milk fat synthesis. CONCLUSIONS: An initial GWAS revealed a highly significant QTL for de novo-synthesized fatty acids on BTA13 and was followed by fine-mapping of the QTL within NCOA6. The most significant SNPs were either synonymous or situated in introns; more research is needed to uncover the underlying causal DNA variation(s).

Fine mapping of a QTL on bovine chromosome 6 using imputed full sequence data suggests a key role for the group-specific component (GC) gene in clinical mastitis and milk production.

BACKGROUND: Clinical mastitis is an inflammation of the mammary gland and causes significant costs to dairy production. It is unfavourably genetically correlated to milk production, and, thus, knowledge of the mechanisms that underlie these traits would be valuable to improve both of them simultaneously through breeding. A quantitative trait locus (QTL) that affects both clinical mastitis and milk production has recently been fine-mapped to around 89 Mb on bovine chromosome 6 (BTA6), but identification of the gene that underlies this QTL was not possible due to the strong linkage disequilibrium between single nucleotide polymorphisms (SNPs) within this region. Our aim was to identify the gene and, if possible, the causal polymorphism(s) responsible for this QTL through association analysis of high-density SNPs and imputed full sequence data in combination with analyses of transcript and protein levels of the identified candidate gene. RESULTS: Associations between SNPs and the studied traits were strongest for SNPs that were located within and immediately upstream of the group-specific component (GC) gene. This gene encodes the vitamin D-binding protein (DBP) and has multiple roles in immune defense and milk production. A 12-kb duplication that was identified downstream of this gene covered its last exon and segregated with the QTL allele that is associated with increased mastitis susceptibility and milk production. However, analyses of GC mRNA levels on the available samples revealed no differences in expression between animals having or lacking this duplication. Moreover, we detected no differences in the concentrations of DBP and its ligand vitamin D between the animals with different GC genotypes that were available for this study. CONCLUSIONS: Our results suggest GC as the gene that underlies the QTL for clinical mastitis and milk production. However, since only healthy animals were sampled for transcription and expression analyses, we could not draw any final conclusion on the absence of quantitative differences between animals with different genotypes. Future studies should investigate GC RNA expression and protein levels in cows with different genotypes during an infection.

The gastric mucosa of Atlantic salmon (Salmo salar) is abundant in highly active chitinases.

Atlantic salmon (Salmo salar) possesses a genome containing 10 genes encoding chitinases, yet their functional roles remain poorly understood. In other fish species, chitinases have been primarily linked to digestion, but also to other functions, as chitinase-encoding genes are transcribed in a variety of non-digestive organs. In this study, we investigated the properties of two chitinases belonging to the family 18 glycoside hydrolase group, namely Chia.3 and Chia.4, both isolated from the stomach mucosa. Chia.3 and Chia.4, exhibiting 95% sequence identity, proved inseparable using conventional chromatographic methods, necessitating their purification as a chitinase pair. Biochemical analysis revealed sustained chitinolytic activity against ß-chitin for up to 24 h, spanning a pH range of 2 to 6. Moreover, subsequent in vitro investigations established that this chitinase pair efficiently degrades diverse chitin-containing substrates into chitobiose, highlighting the potential of Atlantic salmon to utilize novel chitin-containing feed sources. Analysis of the gastric matrix proteome demonstrates that the chitinases are secreted and rank among the most abundant proteins in the gastric matrix. This finding correlates well with the previously observed high transcription of the corresponding chitinase genes in Atlantic salmon stomach tissue. By shedding light on the secreted chitinases in the Atlantic salmon's stomach mucosa and elucidating their functional characteristics, this study enhances our understanding of chitinase biology in this species. Moreover, the observed capacity to effectively degrade chitin-containing materials implies the potential utilization of alternative feed sources rich in chitin, offering promising prospects for sustainable aquaculture practices.

Gene family expansion and functional diversification of chitinase and chitin synthase genes in Atlantic salmon (Salmo salar).

Chitin is one of the most abundant polysaccharides in nature, forming important structures in insects, crustaceans, and fungal cell walls. Vertebrates on the other hand are generally considered "nonchitinous" organisms, despite having highly conserved chitin metabolism-associated genes. Recent work has revealed that the largest group of vertebrates, the teleosts, have the potential to both synthesize and degrade endogenous chitin. Yet, little is known about the genes and proteins responsible for these dynamic processes. Here, we used comparative genomics, transcriptomics, and chromatin accessibility data to characterize the repertoire, evolution, and regulation of genes involved in chitin metabolism in teleosts, with a particular focus on Atlantic salmon. Reconstruction of gene family phylogenies provides evidence for an expansion of teleost and salmonid chitinase and chitin synthase genes after multiple whole-genome duplications. Analyses of multi-tissue gene expression data demonstrated a strong bias of gastrointestinal tract expression for chitin metabolism genes, but with different spatial and temporal tissue specificities. Finally, we integrated transcriptomes from a developmental time series of the gastrointestinal tract with chromatin accessibility data to identify putative transcription factors responsible for regulating chitin metabolism gene expression (CDX1 and CDX2) as well as tissue-specific divergence in the regulation of gene duplicates (FOXJ2). The findings presented here support the hypothesis that chitin metabolism genes in teleosts play a role in developing and maintaining a chitin-based barrier in the teleost gut and provide a basis for further investigations into the molecular basis of this barrier.

Characterization of wooden breast myopathy: a focus on syndecans and ECM remodeling.

Introduction: The skeletal muscle deformity of commercial chickens (Gallus gallus), known as the wooden breast (WB), is associated with fibrotic myopathy of unknown etiology. For future breeding strategies and genetic improvements, it is essential to identify the molecular mechanisms underlying the phenotype. The pathophysiological hallmarks of WB include severe skeletal muscle fibrosis, inflammation, myofiber necrosis, and multifocal degeneration of muscle tissue. The transmembrane proteoglycans syndecans have a wide spectrum of biological functions and are master regulators of tissue homeostasis. They are upregulated and shed (cleaved) as a regulatory mechanism during tissue repair and regeneration. During the last decades, it has become clear that the syndecan family also has critical functions in skeletal muscle growth, however, their potential involvement in WB pathogenesis is unknown. Methods: In this study, we have categorized four groups of WB myopathy in broiler chickens and performed a comprehensive characterization of the molecular and histological profiles of two of them, with a special focus on the role of the syndecans and remodeling of the extracellular matrix (ECM). Results and discussion: Our findings reveal differential expression and shedding of the four syndecan family members and increased matrix metalloproteinase activity. Additionally, we identified alterations in key signaling pathways such as MAPK, AKT, and Wnt. Our work provides novel insights into a deeper understanding of WB pathogenesis and suggests potential therapeutic targets for this condition.

Genetic Differentiation in Hatchery and Stocked Populations of Sea Trout in the Southern Baltic: Selection Evidence at SNP Loci.

The impacts and interactions between hatchery-bred fish and wild fish populations has been a topic of active investigation in recent decades. In some instances, the benefits of stocking can be overshadowed by negative effects such as genetic introgression with natural populations, loss of genetic diversity, and dilution of local adaptations. Methods that facilitate the identification of stocked fish enable us to estimate not only the effectiveness of stocking but also the level of natural reproduction and the degree of hybridization. The longest Baltic river, the Vistula, also has the second highest discharge. Historically, it hosted numerous populations of the anadromous form of brown trout (sea trout); however, dam construction has since interfered with and reduced spawning migration to a rate that is much lower than before. Reduced spawning has resulted in a population collapse and a negative flow-on effect on commercial catches. In response, Poland (along with many other Baltic countries) initiated an intensive stocking program which continues today and which sees the average annual release of 700,000 smolts. As a consequence, today's main-river and inshore catches come from stock-enhanced populations. High-throughput single-nucleotide polymorphism (SNP) genotyping was performed on samples of sea trout from southern Baltic populations; results suggest that a significant portion of the sea trout catches in the Vistula mouth region have direct hatchery origin and indicate the presence of Pomeranian specimens. SNP loci identified as outliers indicate a potential selection pressure that may be related with effects of hatchery breeding and mixing with natural populations. The brown trout SNP array applied in this study showed high effectiveness not only for population differentiation, but more importantly, it emerged as a sensitive tool to provide evidence of detection selection.

Northern European Salmo trutta (L.) populations are genetically divergent across geographical regions and environmental gradients.

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large-scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene-environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub-set of regions, suggesting that locus-specific variation was associated with local processes. A paired-population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection (pcadapt and bayescan) showed low overlap in loci identified as statistical outliers across sub-sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units.

Bos taurus genome assembly.

BACKGROUND: We present here the assembly of the bovine genome. The assembly method combines the BAC plus WGS local assembly used for the rat and sea urchin with the whole genome shotgun (WGS) only assembly used for many other animal genomes including the rhesus macaque. RESULTS: The assembly process consisted of multiple phases: First, BACs were assembled with BAC generated sequence, then subsequently in combination with the individual overlapping WGS reads. Different assembly parameters were tested to separately optimize the performance for each BAC assembly of the BAC and WGS reads. In parallel, a second assembly was produced using only the WGS sequences and a global whole genome assembly method. The two assemblies were combined to create a more complete genome representation that retained the high quality BAC-based local assembly information, but with gaps between BACs filled in with the WGS-only assembly. Finally, the entire assembly was placed on chromosomes using the available map information.Over 90% of the assembly is now placed on chromosomes. The estimated genome size is 2.87 Gb which represents a high degree of completeness, with 95% of the available EST sequences found in assembled contigs. The quality of the assembly was evaluated by comparison to 73 finished BACs, where the draft assembly covers between 92.5 and 100% (average 98.5%) of the finished BACs. The assembly contigs and scaffolds align linearly to the finished BACs, suggesting that misassemblies are rare. Genotyping and genetic mapping of 17,482 SNPs revealed that more than 99.2% were correctly positioned within the Btau_4.0 assembly, confirming the accuracy of the assembly. CONCLUSION: The biological analysis of this bovine genome assembly is being published, and the sequence data is available to support future bovine research.

Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod.

The genetic mechanisms determining sex in teleost fishes are highly variable and the master sex determining gene has only been identified in few species. Here we characterize a male-specific region of 9 kb on linkage group 11 in Atlantic cod (Gadus morhua) harboring a single gene named zkY for zinc knuckle on the Y chromosome. Diagnostic PCR test of phenotypically sexed males and females confirm the sex-specific nature of the Y-sequence. We identified twelve highly similar autosomal gene copies of zkY, of which eight code for proteins containing the zinc knuckle motif. 3D modeling suggests that the amino acid changes observed in six copies might influence the putative RNA-binding specificity. Cod zkY and the autosomal proteins zk1 and zk2 possess an identical zinc knuckle structure, but only the Y-specific gene zkY was expressed at high levels in the developing larvae before the onset of sex differentiation. Collectively these data suggest zkY as a candidate master masculinization gene in Atlantic cod. PCR amplification of Y-sequences in Arctic cod (Arctogadus glacialis) and Greenland cod (Gadus macrocephalus ogac) suggests that the male-specific region emerged in codfishes more than 7.5 million years ago.

Unravelling genetic variation underlying de novo-synthesis of bovine milk fatty acids.

The relative abundance of specific fatty acids in milk can be important for consumer health and manufacturing properties of dairy products. Understanding of genes controlling milk fat synthesis may contribute to the development of dairy products with high quality and nutritional value. This study aims to identify key genes and genetic variants affecting de novo synthesis of the short- and medium-chained fatty acids C4:0 to C14:0. A genome-wide association study using 609,361 SNP markers and 1,811 animals was performed to detect genomic regions affecting fatty acid levels. These regions were further refined using sequencing data to impute millions of additional genetic variants. Results suggest associations of PAEP with the content of C4:0, AACS with the content of fatty acids C4:0-C6:0, NCOA6 or ACSS2 with the longer chain fatty acids C6:0-C14:0, and FASN mainly associated with content of C14:0. None of the top-ranking markers caused amino acid shifts but were mostly situated in putatively regulating regions and suggested a regulatory role of the QTLs. Sequencing mRNA from bovine milk confirmed the expression of all candidate genes which, combined with knowledge of their roles in fat biosynthesis, supports their potential role in de novo synthesis of bovine milk fatty acids.

A genome scan for selection signatures comparing farmed Atlantic salmon with two wild populations: Testing colocalization among outlier markers, candidate genes, and quantitative trait loci for production traits.

Comparative genome scans can be used to identify chromosome regions, but not traits, that are putatively under selection. Identification of targeted traits may be more likely in recently domesticated populations under strong artificial selection for increased production. We used a North American Atlantic salmon 6K SNP dataset to locate genome regions of an aquaculture strain (Saint John River) that were highly diverged from that of its putative wild founder population (Tobique River). First, admixed individuals with partial European ancestry were detected using STRUCTURE and removed from the dataset. Outlier loci were then identified as those showing extreme differentiation between the aquaculture population and the founder population. All Arlequin methods identified an overlapping subset of 17 outlier loci, three of which were also identified by BayeScan. Many outlier loci were near candidate genes and some were near published quantitative trait loci (QTLs) for growth, appetite, maturity, or disease resistance. Parallel comparisons using a wild, nonfounder population (Stewiacke River) yielded only one overlapping outlier locus as well as a known maturity QTL. We conclude that genome scans comparing a recently domesticated strain with its wild founder population can facilitate identification of candidate genes for traits known to have been under strong artificial selection.

Population Structure, Genetic Variation, and Linkage Disequilibrium in Perennial Ryegrass Populations Divergently Selected for Freezing Tolerance.

Low temperature is one of the abiotic stresses seriously affecting the growth of perennial ryegrass (Lolium perenne L.), and freezing tolerance is a complex trait of major agronomical importance in northern and central Europe. Understanding the genetic control of freezing tolerance would aid in the development of cultivars of perennial ryegrass with improved adaptation to frost. The plant material investigated in this study was an experimental synthetic population derived from pair-crosses among five European perennial ryegrass genotypes, representing adaptations to a range of climatic conditions across Europe. A total number of 80 individuals (24 of High frost [HF]; 29 of Low frost [LF], and 27 of Unselected [US]) from the second generation of the two divergently selected populations and an unselected (US) control population were genotyped using 278 genome-wide SNPs derived from perennial ryegrass transcriptome sequences. Our studies investigated the genetic diversity among the three experimental populations by analysis of molecular variance and population structure, and determined that the HF and LF populations are very divergent after selection for freezing tolerance, whereas the HF and US populations are more similar. Linkage disequilibrium (LD) decay varied across the seven chromosomes and the conspicuous pattern of LD between the HF and LF population confirmed their divergence in freezing tolerance. Furthermore, two F st outlier methods; finite island model (fdist) by LOSITAN and hierarchical structure model using ARLEQUIN, both detected six loci under directional selection. These outlier loci are most probably linked to genes involved in freezing tolerance, cold adaptation, and abiotic stress. These six candidate loci under directional selection for freezing tolerance might be potential marker resources for breeding perennial ryegrass cultivars with improved freezing tolerance.