Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows.

BACKGROUND: Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS: In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION: Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Genetic and genomic analysis of reproduction traits in holstein cattle using SNP chip data and imputed sequence level genotypes.

BACKGROUND: Reproductive performance plays an important role in animal welfare, health and profitability in animal husbandry and breeding. It is well established that there is a negative correlation between performance and reproduction in dairy cattle. This relationship is being increasingly considered in breeding programs. By elucidating the genetic architecture of underlying reproduction traits, it will be possible to make a more detailed contribution to this. Our study followed two approaches to elucidate this area; in a first part, variance components were estimated for 14 different calving and fertility traits, and then genome-wide association studies were performed for 13 reproduction traits on imputed sequence-level genotypes with subsequent enrichment analyses. RESULTS: Variance components analyses showed a low to moderate heritability (h2) for the traits analysed, ranging from 0.014 for endometritis up to 0.271 for stillbirth, indicating variable degrees of variation within the reproduction traits. For genome-wide association studies, we were able to detect genome-wide significant association signals for nine out of 13 analysed traits after Bonferroni correction on chromosome 6, 18 and the X chromosome. In total, we detected over 2700 associated SNPs encircling more than 90 different genes using the imputed whole-genome sequence data. Functional associations were reviewed so far known and potential candidate regions in the proximity of reproduction events were hypothesised. CONCLUSION: Our results confirm previous findings of other authors in a comprehensive cohort including 13 different traits at the same time. Additionally, we identified new candidate genes involved in dairy cattle reproduction and made initial suggestions regarding their potential impact, with special regard to the X chromosome as a putative information source for further research. This work can make a contribution to reveal the genetic architecture of reproduction traits in context of trait specific interactions.

Genomic evidence for the suitability of Göttingen Minipigs with a rare seizure phenotype as a model for human epilepsy.

Epilepsy is a complex genetic disorder that affects about 2% of the global population. Although the frequency and severity of epileptic seizures can be reduced by a range of pharmacological interventions, there are no disease-modifying treatments for epilepsy. The development of new and more effective drugs is hindered by a lack of suitable animal models. Available rodent models may not recapitulate all key aspects of the disease. Spontaneous epileptic convulsions were observed in few Göttingen Minipigs (GMPs), which may provide a valuable alternative animal model for the characterisation of epilepsy-type diseases and for testing new treatments. We have characterised affected GMPs at the genome level and have taken advantage of primary fibroblast cultures to validate the functional impact of fixed genetic variants on the transcriptome level. We found numerous genes connected to calcium metabolism that have not been associated with epilepsy before, such as ADORA2B, CAMK1D, ITPKB, MCOLN2, MYLK, NFATC3, PDGFD, and PHKB. Our results have identified two transcription factor genes, EGR3 and HOXB6, as potential key regulators of CACNA1H, which was previously linked to epilepsy-type disorders in humans. Our findings provide the first set of conclusive results to support the use of affected subsets of GMPs as an alternative and more reliable model system to study human epilepsy. Further neurological and pharmacological validation of the suitability of GMPs as an epilepsy model is therefore warranted.

Resistance to chicken amyloid arthropathy is associated with a dysfunctional mutation in serum amyloid A.

Chicken amyloid arthropathy is a debilitating disease with a major impact on animal welfare. Since the disease is triggered by bacterial infection, preventative treatment also contributes to the widespread overuse of antibiotics. Bacterial infection initiates an acute phase response including increased serum amyloid A (SAA) production by the liver. SAA accumulates at sites of infection and in particular in large joints of affected birds. Interestingly, white egg-laying chickens (WL) are resistant to the disease whilst brown egg-laying chickens (BL) are most affected. Disease susceptibility has an immunological basis but the possible contribution of underlying genetic risk factors is not understood. Using a whole genome sequencing approach, we discovered a novel variant in the SAA gene in WL, which is predicted to result in an arginine to serine substitution at position 90 (SAA.R90S). Surprisingly, when overexpressed in chicken hepatocellular carcinoma cells, SAA.R90S was expressed at a higher rate and secreted to a greater degree than the wild-type SAA protein. Moreover, RNASeq analysis showed that the R90S mutant exerted a differential effect on the expression of core transcription factors linked to cell fate determination and cell differentiation. Comparative analysis of gene expression in murine CD4 T-cells stimulated with IL-6/SAA, suggests that SAA.R90S might block an induced cell fate change toward pro-inflammatory T helper 17 cells, which are required for immunological protection against pathogenic bacteria during an acute phase response. Our results provide first mechanistic insights into the genetic resistance of WL to amyloid arthropathy and could be applied to commercial layer breeding programs to improve animal welfare and reduce the negative effects of the overuse of antibiotics.

MG132 dramatically reduces SAA expression in chicken hepatocellular carcinoma cells at the transcript level independent of its endogenous promoter.

BACKGROUND: MG132, a proteasome inhibitor, is widely used to inhibit nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) activity by proteasome-mediated degradation of IκB. It has been marketed as a specific, reversible, cell-permeable and low-cost inhibitor. However, adverse effects of the compound have been reported in the literature. We recently discovered and characterised a point mutation in the acute phase protein serum amyloid A (SAA) in chickens, by overexpressing the protein in chicken hepatocellular carcinoma (LMH) cells. This serine to arginine exchange at amino acid position 90 (SAA.R90S) leads to intra- and extracellular accumulation of SAA, which is surprisingly counteracted by MG132 treatment, independent of SAA's intrinsic promoter. METHODS AND RESULTS: To test, whether low proteasomal degradation of SAA.R90S is responsible for the observed intra- and extracellular SAA accumulation, we intended to inhibit the proteasome in SAA wild type (SAA.WT) overexpressing cells with MG132. However, we observed an unexpected drastic decrease in SAA protein expression at the transcript level. NF-κB gene expression was unchanged by MG132 at the measured time point. CONCLUSIONS: The observed results demonstrate that MG132 inhibits SAA expression at the transcript level, independent of its endogenous promoter. Further, the data might indicate that NF-κB is not involved in the observed MG132-induced inhibition of SAA expression. We, consequently, question in this brief report whether MG132 should truly be categorised as a specific ubiquitin proteasome inhibitor and recommend the usage of alternative compounds.

Simulation of dual-purpose chicken breeding programs implementing gene editing.

BACKGROUND: In spite of being controversial and raising ethical concerns, the application of gene editing is more likely to be accepted when it contributes to improving animal welfare. One of the animal welfare and ethical issues in chicken breeding is chick culling, the killing of the male layer chicks after hatching due to the poor fattening performance. Although establishing dual-purpose chicken lines could solve this problem, unfavorable genetic correlations between egg and meat production traits hindered their competitiveness. Although it is also controversial in ethical terms, gene editing may accelerate genetic progress in dual-purpose chicken and alleviate the ethical concerns from chick culling. RESULTS: The simulation compared the utility improvement in dual-purpose use under two breeding schemes: one consisting in the improvement of the laying hens, and the second in the improvement of a synthetic line obtained from a layer broiler cross. In each breeding scheme, the breeding programs were simulated with and without gene editing. Polygenic breeding values and 500 simulated quantitative trait loci (QTL) with different levels of pleiotropy caused negative correlations between egg production, meat production, and overall health. The results of the simulation demonstrated that genetic gain could be accelerated by at most 81% for several generations if gene editing was used. The actual increase in genetic gain depended on the number of single nucleotide polymorphisms (SNPs) being edited per animal. The rate of genetic improvement became equal in scenarios with and without gene editing after 20 generations. This is because the remaining segregating QTL had small effects and their edition would have negative overall health effects from potential off-target edits. Although gene editing can improve genetic gain in quantitative traits, it can only be recommended as long as QTL with reasonable effect sizes are segregating and detectable. CONCLUSIONS: This simulation demonstrates the potential of gene editing to accelerate the simultaneous improvement of negatively correlated traits. When the risk of negative consequences from gene editing persists, the number of SNPs to be edited should be chosen carefully to obtain the optimal genetic gain.

Mendelian randomization analysis of 34,497 German Holstein cows to infer causal associations between milk production and health traits.

BACKGROUND: Claw diseases and mastitis represent the most important health issues in dairy cattle with a frequently mentioned connection to milk production. Although many studies have aimed at investigating this connection in more detail by estimating genetic correlations, they do not provide information about causality. An alternative is to carry out Mendelian randomization (MR) studies using genetic variants to investigate the effect of an exposure on an outcome trait mediated by genetic variants. No study has yet investigated the causal association of milk yield (MY) with health traits in dairy cattle. Hence, we performed a MR analysis of MY and seven health traits using imputed whole-genome sequence data from 34,497 German Holstein cows. We applied a method that uses summary statistics and removes horizontal pleiotropic variants (having an effect on both traits), which improves the power and unbiasedness of MR studies. In addition, genetic correlations between MY and each health trait were estimated to compare them with the estimates of causal effects that we expected. RESULTS: All genetic correlations between MY and each health trait were negative, ranging from - 0.303 (mastitis) to - 0.019 (digital dermatitis), which indicates a reduced health status as MY increases. The only non-significant correlation was between MY and digital dermatitis. In addition, each causal association was negative, ranging from - 0.131 (mastitis) to - 0.034 (laminitis), but the number of significant associations was reduced to five nominal and two experiment-wide significant results. The latter were between MY and mastitis and between MY and digital phlegmon. Horizontal pleiotropic variants were identified for mastitis, digital dermatitis and digital phlegmon. They were located within or nearby variants that were previously reported to have a horizontal pleiotropic effect, e.g., on milk production and somatic cell count. CONCLUSIONS: Our results confirm the known negative genetic connection between health traits and MY in dairy cattle. In addition, they provide new information about causality, which for example points to the negative energy balance mediating the connection between these traits. This knowledge helps to better understand whether the negative genetic correlation is based on pleiotropy, linkage between causal variants for both trait complexes, or indeed on a causal association.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes.

BACKGROUND: Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected. RESULTS: A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3. CONCLUSIONS: The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Structural variation and eQTL analysis in two experimental populations of chickens divergently selected for feather-pecking behavior.

Feather pecking (FP) is a damaging nonaggressive behavior in laying hens with a heritable component. Its occurrence has been linked to the immune system, the circadian clock, and foraging behavior. Furthermore, dysregulation of miRNA biogenesis, disturbance of the gamma-aminobutyric acid (GABAergic) system, as well as neurodevelopmental deficiencies are currently under debate as factors influencing the propensity for FP behavior. Past studies, which focused on the dissection of the genetic factors involved in FP, relied on single nucleotide polymorphisms (SNPs) and short insertions and deletions < 50 bp (InDels). These variant classes only represent a certain fraction of the genetic variation of an organism. Hence, we reanalyzed whole-genome sequencing data from two experimental populations, which have been divergently selected for FP behavior for over more than 15 generations, performed variant calling for structural variants (SVs) as well as tandem repeats (TRs), and jointly analyzed the data with SNPs and InDels. Genotype imputation and subsequent genome-wide association studies, in combination with expression quantitative trait loci analysis, led to the discovery of multiple variants influencing the GABAergic system. These include a significantly associated TR downstream of the GABA receptor subunit beta-3 (GABRB3) gene, two microRNAs targeting several GABA receptor genes, and dystrophin (DMD), a direct regulator of GABA receptor clustering. Furthermore, we found the transcription factor ETV1 to be associated with the differential expression of 23 genes, which points toward a role of ETV1, together with SMAD4 and KLF14, in the disturbed neurodevelopment of high-feather pecking chickens.

Defining valid breeding goals for animal breeds.

BACKGROUND: The objective of any valid breeding program is to increase the suitability of a breed for its future purposes. The approach most often followed in animal breeding for optimizing breeding goals assumes that the sole desire of the owners is profit maximization. As this assumption is often violated, a generalized approach is needed that does not rely on this assumption. RESULTS: The generalized approach is based on the niche concept. The niche of a breed is a set of environments in which a small population of the breed would have a positive population growth rate. Its growth rate depends on demand from prospective consumers and supply from producers. The approach involves defining the niche that is envisaged for the breed and identifying the trait optima that maximize the breed's adaptation to its envisaged niche within the set of permissible breeding goals. The set of permissible breeding goals is the set of all potential breeding goals that are compatible with animal welfare and could be reached within the planning horizon of the breeding program. In general, the breed's adaptation depends on the satisfaction of the producers with the animals and on the satisfaction of the consumers with the products produced by the animals. When consumers buy live animals, then the breed needs to adapt to both the environments provided by the producers, and the environments provided by the consumers. The profit function is replaced by a more general adaptedness function that measures the breed's adaptation to its envisaged niche. CONCLUSIONS: The proposed approach coincides with the traditional approach if the producers have the sole desire to maximize their income, and if consumer preferences are well reflected by the product prices. If these assumptions are not met, then the traditional approach to breeding goal optimization is unlikely to result in a valid breeding goal. Using the example of companion breeds, this paper shows that the proposed approach has the potential to fill the gap.

In vitro modelling of the influence of alternative feeds (Hermetia illucens, Arthrospira platensis) on the resistance of different rainbow trout populations (Oncorhynchus mykiss) against the viral haemorrhagic septicaemia virus and Yersinia ruckeri.

Replacing fishmeal, a finite resource with high market demand, in the diet of carnivorous rainbow trout with proteins from alternative sources may be a challenge for these fish. Therefore, this study investigated whether replacing fishmeal with protein derived from Hermetia illucens or Arthrospira platensis could promote disease susceptibility in local trout populations with different growth performance. This was assessed in vitro by measuring susceptibility to infection with the viral haemorrhagic septicaemia virus (VHSV) or the bacterium Yersinia ruckeri. Analysis of fin tissue explants and primary cell cultures from scales from the three trout populations infected in vitro with VHSV and gill explants infected with Y. ruckeri showed no significant differences in virus replication or bacterial counts. Evaluation of the virucidal or bactericidal effect of skin mucus showed a significant reduction in viral load and bacterial count for all samples with mucus addition, but no significant difference was observed between the experimental groups. This study documents no apparent impairment of innate immune mechanisms in the skin and gills of trout after feeding a diet replacing fishmeal with Arthrospira or Hermetia proteins. This underlines the potential of these alternative protein sources for the further development of sustainable trout aquaculture.

Genetic analysis of production traits and body size measurements and their relationships with metabolic diseases in German Holstein cattle.

This study sheds light on the genetic complexity and interplay of production, body size, and metabolic health in dairy cattle. Phenotypes for body size-related traits from conformation classification (130,166 animals) and production (101,562 animals) of primiparous German Holstein cows were available. Additionally, 21,992, 16,641, and 7,096 animals were from herds with recordings of the metabolic diseases ketosis, displaced abomasum, and milk fever in first, second, and third lactation. Moreover, all animals were genotyped. Heritabilities of traits and genetic correlations between all traits were estimated and GWAS were performed. Heritability was between 0.240 and 0.333 for production and between 0.149 and 0.368 for body size traits. Metabolic diseases were lowly heritable, with estimates ranging from 0.011 to 0.029 in primiparous cows, from 0.008 to 0.031 in second lactation, and from 0.037 to 0.052 in third lactation. Production was found to have negative genetic correlations with body condition score (BCS; -0.279 to -0.343) and udder depth (-0.348 to -0.419). Positive correlations were observed for production and body depth (0.138-0.228), dairy character (DCH) (0.334-0.422), and stature (STAT) (0.084-0.158). In first parity cows, metabolic disease traits were unfavorably correlated with production, with genetic correlations varying from 0.111 to 0.224, implying that higher yielding cows have more metabolic problems. Genetic correlations of disease traits in second and third lactation with production in primiparous cows were low to moderate and in most cases unfavorable. While BCS was negatively correlated with metabolic diseases (-0.255 to -0.470), positive correlations were found between disease traits and DCH (0.269-0.469) as well as STAT (0.172-0.242). Thus, the results indicate that larger and sharper animals with low BCS are more susceptible to metabolic disorders. Genome-wide association studies revealed several significantly associated SNPs for production and conformation traits, confirming previous findings from literature. Moreover, for production and conformation traits, shared significant signals on Bos taurus autosome (BTA) 5 (88.36 Mb) and BTA 6 (86.40 to 87.27 Mb) were found, implying pleiotropy. Additionally, significant SNPs were observed for metabolic diseases on BTA 3, 10, 14, 17, and 26 in first lactation and on BTA 2, 6, 8, 17, and 23 in third lactation. Overall, this study provides important insights into the genetic basis and interrelations of relevant traits in today's Holstein cattle breeding programs, and findings may help to improve selection decisions.

A genomic assessment of the correlation between milk production traits and claw and udder health traits in Holstein dairy cattle.

Claw diseases and mastitis represent the most important disease traits in dairy cattle with increasing incidences and a frequently mentioned connection to milk yield. Yet, many studies aimed to detect the genetic background of both trait complexes via fine-mapping of quantitative trait loci. However, little is known about genomic regions that simultaneously affect milk production and disease traits. For this purpose, several tools to detect local genetic correlations have been developed. In this study, we attempted a detailed analysis of milk production and disease traits as well as their interrelationship using a sample of 34,497 50K genotyped German Holstein cows with milk production and claw and udder disease traits records. We performed a pedigree-based quantitative genetic analysis to estimate heritabilities and genetic correlations. Additionally, we generated GWAS summary statistics, paying special attention to genomic inflation, and used these data to identify shared genomic regions, which affect various trait combinations. The heritability on the liability scale of the disease traits was low, between 0.02 for laminitis and 0.19 for interdigital hyperplasia. The heritabilities for milk production traits were higher (between 0.27 for milk energy yield and 0.48 for fat-protein ratio). Global genetic correlations indicate the shared genetic effect between milk production and disease traits on a whole genome level. Most of these estimates were not significantly different from zero, only mastitis showed a positive one to milk (0.18) and milk energy yield (0.13), as well as a negative one to fat-protein ratio (-0.07). The genomic analysis revealed significant SNPs for milk production traits that were enriched on Bos taurus autosome 5, 6, and 14. For digital dermatitis, we found significant hits, predominantly on Bos taurus autosome 5, 10, 22, and 23, whereas we did not find significantly trait-associated SNPs for the other disease traits. Our results confirm the known genetic background of disease and milk production traits. We further detected 13 regions that harbor strong concordant effects on a trait combination of milk production and disease traits. This detailed investigation of genetic correlations reveals additional knowledge about the localization of regions with shared genetic effects on these trait complexes, which in turn enables a better understanding of the underlying biological pathways and putatively the utilization for a more precise design of breeding schemes.

Genomic dominance variance analysis of health and milk production traits in German Holstein cattle.

Genomic analyses commonly explore the additive genetic variance of traits. The non-additive variance, however, is usually small but often significant in dairy cattle. This study aimed at dissecting the genetic variance of eight health traits that recently entered the total merit index in Germany and the somatic cell score (SCS), as well as four milk production traits by analysing additive and dominance variance components. The heritabilities were low for all health traits (between 0.033 for mastitis and 0.099 for SCS), and moderate for the milk production traits (between 0.261 for milk energy yield and 0.351 for milk yield). For all traits, the contribution of dominance variance to the phenotypic variance was low, varying between 0.018 for ovarian cysts and 0.078 for milk yield. Inbreeding depression, inferred from the SNP-based observed homozygosity, was significant only for the milk production traits. The contribution of dominance variance to the genetic variance was larger for the health traits, ranging from 0.233 for ovarian cysts to 0.551 for mastitis, encouraging further studies that aim at discovering QTLs based on their additive and dominance effects.

Single-step SNPBLUP evaluation in six German beef cattle breeds.

The implementation of genomic selection for six German beef cattle populations was evaluated. Although the multiple-step implementation of genomic selection is the status quo in most national dairy cattle evaluations, the breeding structure of German beef cattle, coupled with the shortcoming and complexity of the multiple-step method, makes single step a more attractive option to implement genomic selection in German beef cattle populations. Our objective was to develop a national beef cattle single-step genomic evaluation in five economically important traits in six German beef cattle populations and investigate its impact on the accuracy and bias of genomic evaluations relative to the current pedigree-based evaluation. Across the six breeds in our study, 461,929 phenotyped and 14,321 genotyped animals were evaluated with a multi-trait single-step model. To validate the single-step model, phenotype data in the last 2 years were removed in a forward validation study. For the conventional and single-step approaches, the genomic estimated breeding values of validation animals and other animals were compared between the truncated and the full evaluations. The correlation of the GEBVs between the full and truncated evaluations in the validation animals was slightly higher in the single-step evaluation. The regression of the full GEBVs on truncated GEBVs was close to the optimal value of 1 for both the pedigree-based and the single-step evaluations. The SNP effect estimates from the truncated evaluation were highly correlated with those from the full evaluation, with values ranging from 0.79 to 0.94. The correlation of the SNP effect was influenced by the number of genotyped animals shared between the full and truncated evaluations. The regression coefficients of the SNP effect of the full evaluation on the truncated evaluation were all close to the expected value of 1, indicating unbiased estimates of the SNP markers for the production traits. The Manhattan plot of the SNP effect estimates identified chromosomal regions harbouring major genes for muscling and body weight in breeds of French origin. Based on the regression intercept and slope of the GEBVs of validation animals, the single-step evaluation was neither inflated nor deflated across the six breeds. Overall, the single-step model resulted in a more accurate and stable evaluation. However, due to the small number of genotyped individuals, the single-step method only provided slightly better results when compared to the pedigree-based method.

Structural variants and tandem repeats in the founder individuals of four F2 pig crosses and implications to F2 GWAS results.

BACKGROUND: Structural variants and tandem repeats are relevant sources of genomic variation that are not routinely analyzed in genome wide association studies mainly due to challenging identification and genotyping. Here, we profiled these variants via state-of-the-art strategies in the founder animals of four F2 pig crosses using whole-genome sequence data (20x coverage). The variants were compared at a founder level with the commonly screened SNPs and small indels. At the F2 level, we carried out an association study using imputed structural variants and tandem repeats with four growth and carcass traits followed by a comparison with a previously conducted SNPs and small indels based association study. RESULTS: A total of 13,201 high confidence structural variants and 103,730 polymorphic tandem repeats (with a repeat length of 2-20 bp) were profiled in the founders. We observed a moderate to high (r from 0.48 to 0.57) level of co-localization between SNPs or small indels and structural variants or tandem repeats. In the association step 56.56% of the significant variants were not in high LD with significantly associated SNPs and small indels identified for the same traits in the earlier study and thus presumably not tagged in case of a standard association study. For the four growth and carcass traits investigated, many of the already proposed candidate genes in our previous studies were confirmed and additional ones were identified. Interestingly, a common pattern on how structural variants or tandem repeats regulate the phenotypic traits emerged. Many of the significant variants were embedded or nearby long non-coding RNAs drawing attention to their functional importance. Through which specific mechanisms the identified long non-coding RNAs and their associated structural variants or tandem repeats contribute to quantitative trait variation will need further investigation. CONCLUSIONS: The current study provides insights into the characteristics of structural variants and tandem repeats and their role in association studies. A systematic incorporation of these variants into genome wide association studies is advised. While not of immediate interest for genomic prediction purposes, this will be particularly beneficial for elucidating biological mechanisms driving the complex trait variation.

Immunological insights into the resistance of Nile tilapia strains to an infection with tilapia lake virus.

The emergence of viral diseases affecting fish and causing very high mortality can lead to the disruption of aquaculture production. Recently, this occurred in Nile tilapia aquaculture where a disease caused by a systemic infection with a novel virus named tilapia lake virus (TiLV) caused havoc in cultured populations. With mortality surpassing 90% in young tilapia, the disease caused by TiLV has become a serious challenge for global tilapia aquaculture. In order to partly mitigate the losses, we explored the natural resistance to TiLV-induced disease in three genetic strains of tilapia which were kept at the University of Göttingen, Germany. We used two strains originating from Nilotic regions (Lake Mansala (MAN) and Lake Turkana (ELM)) and one from an unknown location (DRE). We were able to show that the virus is capable of overcoming the natural resistance of tilapia when injected, providing inaccurate mortality results that might complicate finding the resistant strains. Using the cohabitation infection model, we found an ELM strain that did not develop any clinical signs of the infection, which resulted in nearly 100% survival rate. The other two strains (DRE and MAN) showed severe clinical signs and much lower survival rates of 29.3% in the DRE strain and 6.7% in the MAN strain. The disease resistance of tilapia from the ELM strain was correlated with lower viral loads both at the mucosa and internal tissues. Our results suggest that the lower viral load could be caused by a higher magnitude of a mx1-based antiviral response in the initial phase of infection. The lower pro-inflammatory responses also found in the resistant strain might additionally contribute to its protection from developing pathological changes related to the disease. In conclusion, our results suggest the possibility of using TiLV-resistant strains as an ad hoc, cost-effective solution to the TiLV challenge. However, as the fish from the disease-resistant strain still retained significant virus loads in liver and brain and thus could become persistent virus carriers, they should be used within an integrative approach also combining biosecurity, diagnostics and vaccination measures.\.

The light response in chickens divergently selected for feather pecking behavior reveals mechanistic insights towards psychiatric disorders.

BACKGROUND: Feather pecking is a serious behavioral disorder in chickens that has a considerable impact on animal welfare and poses an economic burden for poultry farming. To study the underlying genetics of feather pecking animals were divergently selected for feather pecking over 15 generations based on estimated breeding values for the behavior. METHODS AND RESULTS: By characterizing the transcriptomes of whole brains isolated from high and low feather pecking chickens in response to light stimulation we discovered a putative dysregulation of micro RNA processing caused by a lack of Dicer1. This results in a prominent downregulation of the GABRB2 gene and other GABA receptor transcripts, which might cause a constant high level of excitation in the brains of high feather pecking chickens. Moreover, our results point towards an increase in immune system-related transcripts that may be caused by higher interferon concentrations due to Dicer1 downregulation. CONCLUSION: Based on our results, we conclude that feather pecking in chickens and schizophrenia in humans have numerous common features. For instance, a Dicer1 dependent disruption of miRNA biogenesis and the lack of GABRB2 expression have been linked to schizophrenia pathogenesis. Furthermore, disturbed circadian rhythms and dysregulation of genes involved in the immune system are common features of both conditions.

Development and validation of a horse reference panel for genotype imputation.

BACKGROUND: Genotype imputation is a cost-effective method to generate sequence-level genotypes for a large number of animals. Its application can improve the power of genomic studies, provided that the accuracy of imputation is sufficiently high. The purpose of this study was to develop an optimal strategy for genotype imputation from genotyping array data to sequence level in German warmblood horses, and to investigate the effect of different factors on the accuracy of imputation. Publicly available whole-genome sequence data from 317 horses of 46 breeds was used to conduct the analyses. RESULTS: Depending on the size and composition of the reference panel, the accuracy of imputation from medium marker density (60K) to sequence level using the software Beagle 5.1 ranged from 0.64 to 0.70 for horse chromosome 3. Generally, imputation accuracy increased as the size of the reference panel increased, but if genetically distant individuals were included in the panel, the accuracy dropped. Imputation was most precise when using a reference panel of multiple but related breeds and the software Beagle 5.1, which outperformed the other two tested computer programs, Impute 5 and Minimac 4. Genome-wide imputation for this scenario resulted in a mean accuracy of 0.66. Stepwise imputation from 60K to 670K markers and subsequently to sequence level did not improve the accuracy of imputation. However, imputation from higher density (670K) was considerably more accurate (about 0.90) than from medium density. Likewise, imputation in genomic regions with a low marker coverage resulted in a reduced accuracy of imputation. CONCLUSIONS: The accuracy of imputation in horses was influenced by the size and composition of the reference panel, the marker density of the genotyping array, and the imputation software. Genotype imputation can be used to extend the limited amount of available sequence-level data from horses in order to boost the power of downstream analyses, such as genome-wide association studies, or the detection of embryonic lethal variants.

Multi-species transcriptome meta-analysis of the response to retinoic acid in vertebrates and comparative analysis of the effects of retinol and retinoic acid on gene expression in LMH cells.

BACKGROUND: Retinol (RO) and its active metabolite retinoic acid (RA) are major regulators of gene expression in vertebrates and influence various processes like organ development, cell differentiation, and immune response. To characterize a general transcriptomic response to RA-exposure in vertebrates, independent of species- and tissue-specific effects, four publicly available RNA-Seq datasets from Homo sapiens, Mus musculus, and Xenopus laevis were analyzed. To increase species and cell-type diversity we generated RNA-seq data with chicken hepatocellular carcinoma (LMH) cells. Additionally, we compared the response of LMH cells to RA and RO at different time points. RESULTS: By conducting a transcriptome meta-analysis, we identified three retinoic acid response core clusters (RARCCs) consisting of 27 interacting proteins, seven of which have not been associated with retinoids yet. Comparison of the transcriptional response of LMH cells to RO and RA exposure at different time points led to the identification of non-coding RNAs (ncRNAs) that are only differentially expressed (DE) during the early response. CONCLUSIONS: We propose that these RARCCs stand on top of a common regulatory RA hierarchy among vertebrates. Based on the protein sets included in these clusters we were able to identify an RA-response cluster, a control center type cluster, and a cluster that directs cell proliferation. Concerning the comparison of the cellular response to RA and RO we conclude that ncRNAs play an underestimated role in retinoid-mediated gene regulation.