Search

A dataset was simulated and distributed to participants of the QTLMAS XII workshop who were invited to develop genomic selection models. Each contributing group was asked to describe the model development and validation as well as to submit genomic predictions for three generations of individuals, for which they only knew the genotypes. The organisers used these genomic predictions to perform the final validation by comparison to the true breeding values, which were known only to the organisers. Methods used by the 5 groups fell in 3 classes 1) fixed effects models 2) BLUP models, and 3) Bayesian MCMC based models. The Bayesian analyses gave the highest accuracies, followed by the BLUP models, while the fixed effects models generally had low accuracies and large error variance. The best BLUP models as well as the best Bayesian models gave unbiased predictions. The BLUP models are clearly sensitive to the assumed SNP variance, because they do not estimate SNP variance, but take the specified variance as the true variance. The current comparison suggests that Bayesian analyses on haplotypes or SNPs are the most promising approach for Genomic selection although the BLUP models may provide a computationally attractive alternative with little loss of efficiency. On the other hand fixed effect type models are unlikely to provide any gain over traditional pedigree indexes for selection. [ABSTRACT FROM AUTHOR]

Genome scans have identified quantitative trait loci (QTLs) affecting milk yield and composition in dairy cattle. For one QTL on bovine chromosome 6 (BTA6), previously fine-mapped to a 420-Kb region, mutations in two different genes (OPN and ABCG2) have been proposed as the underlying functional mutation. Comparing the arguments for each gene suggests that both mutations are equally probable. However, functional studies and/or additional populations are required to provide a definite answer. [Copyright &y& Elsevier]

A whole-genome scan was conducted to map quantitative trait loci (QTL) for BSE resistance or susceptibility. Cows from four half-sib families were included and 173 microsatellite markers were used to construct a 2835-cM (Kosambi) linkage map covering 29 autosomes and the pseudoautosomal region of the sex chromosome. Interval mapping by linear regression was applied and extended to a multiple-QTL analysis approach that used identified QTL on other chromosomes as cofactors to increase mapping power. In the multiple-QTL analysis, two genome-wide significant QTL (BTA17 and X/Yps) and four genome-wide suggestive QTL (BTA1, 6, 13, and 19) were revealed. The QTL identified here using linkage analysis do not overlap with regions previously identified using TDT analysis. One factor that may explain the disparity between the results is that a more extensive data set was used in the present study. Furthermore, methodological differences between TDT and linkage analyses may affect the power of these approaches. [ABSTRACT FROM AUTHOR]

Skeletal damage is a challenge for laying hens, as the adaptations for egg laying make them susceptible to osteoporosis. We know that activity, growth and puberty influence skeletal quality genetics. However, in this story, we follow the process of resolving a QTL for skeletal quality to a single gene. Providing both a nutritional and a genetic solution and illustrating another level of complexity. Using divergently selected populations in an F2 cross, a QTL for skeletal quality was described. Utilising a series of experiments in the commercial founder line the QTL was fine mapped on chromosome 1 and the differences in physiology of hens that segregated for markers at the locus defined e.g. tibia breaking strength 200.4 vs 218.1 Newton (p< 0.002). Transcriptome profiling of tibia from the high and low bone strength genotype revealed a highly differentially expressed gene co localising to the QTL, the enzyme cystathionine beta synthase (CBS). CBS is a component of one carbon metabolism. Plasma homocysteine, the substrate of CBS had a higher concentration in the genotype associated with poor bone quality. Homocysteine interferes with crosslinking of the collagen matrix, critical to bone formation. Homocysteine can be re-methylated into methionine, a limiting amino acid in laying hens, using betaine as a methyl donor. Testing the hypothesis that reducing plasma homocysteine by feeding betaine, a widely available feed additive, would improve bone quality in laying hens proved successful. This is just one component of a complex trait and its elucidation took us along unexpected routes. However, I am optimistic that with modern tools and an open mind we can progress more quickly. Alongside understanding the factors that determine skeletal quality genetics, the discovery of the mechanisms behind complex traits, although rarely easy, offer the potential for novel interventions to improve animal health and production. [ABSTRACT FROM AUTHOR]

The author comments on the report "Global plan of action for animal genetic resources (AnGR) and the interlaken declaration." AnGR are crucial for global food security and livelihoods. Livestock products have high densities of energy, protein, and other crucial nutrients, which are specifically of benefits for infants and mothers-to-be.

Background: Substantial variation in male fertility is regularly observed in farmed Arctic charr. However, detailed investigations of its fluctuation during a reproductive season and across years are lacking. Furthermore, information about the effect of underlying genetic factors influencing sperm quality is scarce. The current study focused on seasonal and age-related factors that may affect sperm quality characteristics in males reared in natural and delayed photoperiods. Animals were sampled three times for two consecutive years, and sperm quality parameters were recorded using a computer-assisted sperm analysis (CASA) system. Thereafter, high-throughput sequencing technologies were applied, aiming to identify genomic regions related to the variation of sperm quality throughout the reproductive season. Results: An across-season variation in the recorded sperm quality parameters was evident. Overall, 29% and 42% of males from the natural and delayed spawning groups had a highly variable total progressive motility. Males at four years of age showed significantly higher sperm motility and velocities during the early October and November recordings compared to the following year when the same animals were five years of age. On the other hand, the opposite was observed regarding sperm concentration during the last sampling. A genome-wide FST scan detected SNP differentiation among males with high and low variability in total progressive motility (PM) on eight chromosomes (FST > 0.17), Genome wide windows with the highest FST contained SNPs in proximity (within 250 kb up- and downstream distance) to 16 genes with sperm quality biological functions in mammalian species. Conclusion: Our findings provide a detailed view of seasonal, age-related, and genetic effects on sperm quality and can be used to guide decisions on broodstock selection and hatchery management. [ABSTRACT FROM AUTHOR]

Discusses the abstract of the study 'Molecular genetic mapping of imprinted QTL for gestation length in Bos taurus x Bos indicus cross cattle produced by embryo transfer,' presented at the American Society of Animal Science's Southern Section conference on January 31-February 4, 2003, in Mobile, Alabama.

Osteoporosis and bone fractures are a severe problem for the welfare of laying hens, with genetics and environment, such as housing system, each making substantial contributions to bone strength. In this work, we performed genetic analyses of bone strength, bone mineral density, and bone composition, as well as body weight, in 860 commercial crossbred laying hens from 2 different companies, kept in either furnished cages or floor pens. We compared bone traits between housing systems and crossbreds and performed a genome-wide association study of bone properties and body weight. As expected, the 2 housing systems produced a large difference in bone strength, with layers housed in floor pens having stronger bones. These differences were accompanied by differences in bone geometry, mineralization, and chemical composition. Genome scans either combining or independently analyzing the 2 housing systems revealed no genome-wide significant loci for bone breaking strength. We detected 3 loci for body weight that were shared between the housing systems on chromosomes 4, 6, and 27 (either genome-wide significant or suggestive) and these coincide with associations for bone length. In summary, we found substantial differences in bone strength, content, and composition between hens kept in floor pens and furnished cages that could be attributed to greater physical activity in pen housing. We found little evidence for large-effect loci for bone strength in commercial crossbred hens, consistent with a highly polygenic architecture for bone strength in the production environment. The lack of consistent genetic associations between housing systems in combination with the differences in bone phenotypes could be due to gene-by-environment interactions with housing system or a lack of power to detect shared associations for bone strength. [ABSTRACT FROM AUTHOR]

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is possible with models for inherited variability but may require a two‐step analysis. [ABSTRACT FROM AUTHOR]

Routine implementation of genomic information for guiding selection decisions is not yet common in the majority of aquaculture species. Reduced representation sequencing approaches offer a cost‐effective solution for obtaining genome‐wide information in species with a limited availability of genomic resources. In the current study, we implemented double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) on an Arctic charr strain with the longest known history of selection (approximately 40 years) aiming to improve selection decisions. In total, 1730 animals reared at four different farms in Sweden and spanning from year classes 2013–2017 were genotyped using ddRAD‐seq. Approximately 5000 single nucleotide polymorphisms (SNPs) were identified, genetic diversity‐related metrics were estimated, and genome‐wide association studies (GWAS) for body length at different time points and age of sexual maturation were conducted. Low genetic differentiation amongst animals from the different farms was observed based on both the results from pairwise Fst values and principal component analysis (PCA). The existence of associations was investigated between the mean genome‐wide heterozygosity of each full‐sib family (year class 2017) and the corresponding inbreeding coefficient or survival to the eyed stage. A moderate correlation (−0.33) was estimated between the mean observed heterozygosity of each full‐sib family and the corresponding inbreeding coefficient, while no linear association was obtained with the survival to the eyed stage. GWAS did not detect loci with major effect for any of the studied traits. However, genomic regions explaining more than 1% of the additive genetic variance for either studied traits were suggested across 14 different chromosomes. Overall, key insights valuable for future selection decisions of Arctic charr have been obtained, suggesting ddRAD as an attractive genotyping platform for obtaining genome‐wide information in a cost‐effective manner. [ABSTRACT FROM AUTHOR]

The mouse Collaborative Cross (CC) is a panel of eight-way recombinant inbred lines: eight diverse parental strains are intermated, followed by repeated sibling mating, many times in parallel, to create a new set of inbred lines whose genomes are random mosaics of the genomes of the original eight strains. Many generations are required to reach inbreeding, and so a number of investigators have sought to make use of phenotype and genotype data on mice from intermediate generations during the formation of the CC lines (so-called pre-CC mice). The development of a hidden Markov model for genotype reconstruction in such pre-CC mice, on the basis of incompletely informative genetic markers (such as single-nucleotide polymorphisms), formally requires the two-locus genotype probabilities at an arbitrary generation along the path to inbreeding. In this article, I describe my efforts to calculate such probabilities. While closed-form solutions for the two-locus genotype probabilities could not be derived, I provide a prescription for calculating such probabilities numerically. In addition, I present a number of useful quantities, including single-locus genotype probabilities, two-locus haplotype probabilities, and the fixation probability and map expansion at each generation along the course to inbreeding. [ABSTRACT FROM AUTHOR]

Tilapia hatcheries in Tanzania rely heavily on importing germplasm. Nevertheless, the genetic structure of the imported stocks is poorly understood. In the current study, the level of genetic diversity and differentiation of eight populations of Nile tilapia (Oreochromis niloticus) strains imported in Tanzania was investigated. Four of the studied strains originated from Thailand, three from Uganda, and one from the Netherlands. Double-digest restriction site–associated DNA sequencing (ddRAD-seq) was applied to identify and genotype single nucleotide polymorphisms (SNPs). In total, 2214 SNPs passed all the quality control steps and were utilized for downstream analysis. Mean heterozygosity estimates were higher for the Thailand strains (Ho, 0.23) compared with the strains from Uganda (Ho, 0.12). Low genetic distance was observed amongst populations from the same geographic origin (Fst, 0.01–0.04). However, genetic distance between populations from different geographic origins was substantial (Fst, 0.24–0.44). Bayesian model–based clustering (STRUCTURE) and discriminant analysis of principal components (DAPC) grouped the studied animals into three distinct clusters. A cross-validation approach (where 25% of animals from each population were considered of unknown origin) was conducted in order to test the efficiency of the SNP dataset for identifying the population of origin. The cross-validation procedure was repeated 10 times resulting in approximately 97% of the tested animals being allocated to the correct geographic population of origin. The breeding history and hatchery practices used to manage these stocks prior and after import appear to be the main factors for the genetic diversity observed in this study. Our study will help inform hatchery stock management and future breeding program designs in Tanzania. [ABSTRACT FROM AUTHOR]

Non-coagulating milk is a serious problem in the cheese industry, since it decreases cheese yield, resulting in decreased economic output. This study evaluated rennet-induced coagulation properties and composition of milk from individual Swedish Red Dairy Cattle. Milk samples from 679 individual cows were rheologically evaluated, of which 18.1% of the cows produced non-coagulating milk and 18.9% produced poor-coagulating milk. This resulted in 37% of the milk samples being non-optimal in cheese production, which is an alarmingly high figure. A comparison between non-coagulating and coagulating milk showed a significantly lower calcium content and less free Ca2+ in non-coagulating milk. The results provide more information about non- and poor-coagulating milk and will be further used to understand the genetic background of non-coagulating milk and breed against this undesired milk property. [ABSTRACT FROM AUTHOR]

Imputation is a cost-effective approach to augment marker data for genomic selection and genome-wide association studies. However, most imputation studies have focused on autosomes. Here, we assessed the imputation of markers on the X chromosome in Holstein cattle for nongenotyped animals and animals genotyped with low-density (Illumina BovineLD, Illumina Inc., San Diego, CA) chips, using animals genotyped with medium-density (Illumina BovineSNP50) chips. A total of 26,884 genotyped Holstein individuals genotyped with medium-density chips were used in this study. Imputation was carried out using FImpute V2.2. The following parameters were examined: treating the pseudoautosomal region as autosomal or as X specific, different sizes of reference groups, different male/female proportions in the reference group, and cumulated degree of relationship between the reference group and target group. The imputation accuracy of markers on the X chromosome was improved if the pseudoautosomal region was treated as autosomal. Increasing the proportion of females in the reference group improved the imputation accuracy for the X chromosome. Imputation for nongenotyped animals in general had lower accuracy compared with animals genotyped with the low-density single nucleotide polymorphism array. In addition, higher cumulative pedigree relationships between the reference group and the target animal led to higher imputation accuracy. In the future, better marker coverage of the X chromosome should be developed to facilitate genomic studies involving the X chromosome. [ABSTRACT FROM AUTHOR]

Using targeted genetical genomics, a quantitative trait locus (QTL) affecting the initial postmortem pH value of chicken breast muscle (Pectoralis major) on chromosome 1 (GGA1) recently was fine-mapped. Thirteen genes were present in the QTL region of approximately 1 Mb. In this study, 10 birds that were inferred to be homozygous for either the high (QQ) or low (qq) QTL allele were selected for resequencing. After enrichment for 1 Mb around the QTL region, >500 x coverage for the QTL region in each of the 10 birds was obtained. In total 5056 single-nucleotide polymorphisms (SNPs) were identified for which the genotypes were consistent with one of the QTL genotypes. We used custom tools to identify putative causal mutations in the mapped QTL region from these SNPs. Four nonsynonymous SNPs differentiating the two QTL genotype groups were identified within four local genes (PRDX4, EIF2S3, PCYT1B, and E1BTD2). Although these are likely candidate SNPs to explain the QTL effect, 54 additional consensus SNPs were detected within gene-related regions (untranslated regions, splicing sites CpG island, and promoter regions) for the QQ birds and 71 for the qq birds. These could also play a role explaining the observed QTL effect. The results provide an important step for prioritizing among a large amount of candidate mutations and significantly contribute to the understanding of the genetic mechanisms affecting the initial postmortem pH value of chicken muscle. [ABSTRACT FROM AUTHOR]

The JAX Diversity Outbred population is a new mouse resource derived from partially inbred Collaborative Cross strains and maintained by randomized outcrossing. As such, it segregates the same allelic variants as the Collaborative Cross but embeds these in a distinct population architecture in which each animal has a high degree of heterozygosity and carries a unique combination of alleles. Phenotypic diversity is striking and often divergent from phenotypes seen in the founder strains of the Collaborative Cross. Allele frequencies and recombination density in early generations of Diversity Outbred mice are consistent with expectations based on simulations of the mating design. We describe analytical methods for genetic mapping using this resource and demonstrate the power and high mapping resolution achieved with this population by mapping a serum cholesterol trait to a 2-Mb region on chromosome 3 containing only 11 genes. Analysis of the estimated allele effects in conjunction with complete genome sequence data of the founder strains reduced the pool of candidate polymorphisms to seven SNPs, five of which are located in an intergenic region upstream of the Foxo1 gene. [ABSTRACT FROM AUTHOR]

The classic diallel takes a set of parents and produces offspring from all possible mating pairs. Phenotype values among the offspring can then be related back to their respective parentage. When the parents are diploid, sexed, and inbred, the diallel can characterize aggregate effects of genetic background on a phenotype, revealing effects of strain dosage, heterosis, parent of origin, epistasis, and sex-specific versions thereof. However, its analysis is traditionally intricate, unforgiving of unplanned missing information, and highly sensitive to imbalance, making the diallel unapproachable to many geneticists. Nonetheless, imbalanced and incomplete diallels arise frequently, albeit unintentionally, as by-products of larger-scale experiments that collect F1 data, for example, pilot studies or multiparent breeding efforts such as the Collaborative Cross or the Arabidopsis MAGIC lines. We present a general Bayesian model for analyzing diallel data on dioecious diploid inbred strains that cleanly decomposes the observed patterns of variation into biologically intuitive components, simultaneously models and accommodates outliers, and provides shrinkage estimates of effects that automatically incorporate uncertainty due to imbalance, missing data, and small sample size. We further present a model selection procedure for weighing evidence for or against the inclusion of those components in a predictive model. We evaluate our method through simulation and apply it to incomplete diallel data on the founders and F1's of the Collaborative Cross, robustly characterizing the genetic architecture of 48 phenotypes. [ABSTRACT FROM AUTHOR]

We present full-genome genotype imputations for 100 classical laboratory mouse strains, using a novel method. Using genotypes at 549,683 SNP loci obtained with the Mouse Diversity Array, we partitioned the genome of 100 mouse strains into 40,647 intervals that exhibit no evidence of historical recombination. For each of these intervals we inferred a local phylogenetic tree. We combined these data with 12 million loci with sequence variations recently discovered by whole-genome sequencing in a common subset of 12 classical laboratory strains. For each phylogenetic tree we identified strains sharing a leaf node with one or more of the sequenced strains. We then imputed high- and medium-confidence genotypes for each of 88 nonsequenced genomes. Among inbred strains, we imputed 92% of SNPs genome-wide, with 71% in high-confidence regions. Our method produced 977 million new genotypes with an estimated per-SNP error rate of 0.083% in high-confidence regions and 0.37% genome-wide. Our analysis identified which of the 88 nonsequenced strains would be the most informative for improving full-genome imputation, as well as which additional strain sequences will reveal more new genetic variants. Imputed sequences and quality scores can be downloaded and visualized online. [ABSTRACT FROM AUTHOR]

There has been a great deal of interest in the development of methodologies to map quantitative trait loci (QTL) using experimental crosses in the last 2 decades. Experimental crosses in animal and plant sciences provide important data sources for mapping QTL through linkage analysis. The Collaborative Cross (CC) is a renewable mouse resource that is generated from eight genetically diverse founder strains to mimic the genetic diversity in humans. The recombinant inbred intercrosses (RIX) generated from CC recombinant inbred (RI) lines share similar genetic structures of F2 individuals but with up to eight alleles segregating at any one locus. In contrast to F2 mice, genotypes of RIX can be inferred from the genotypes of their RI parents and can be produced repeatedly. Also, RIX mice typically do not share the same degree of relatedness. This unbalanced genetic relatedness requires careful statistical modeling to avoid false-positive findings. Many quantitative traits are inherently complex with genetic effects varying with other covariates, such as age. For such complex traits, if phenotype data can be collected over a wide range of ages across study subjects, their dynamic genetic patterns can be investigated. Parametric functions, such as sigmoidal or logistic functions, have been used for such purpose. In this article, we propose a flexible nonparametric time-varying coefficient QTL mapping method for RIX data. Our method allows the QTL effects to evolve with time and naturally extends classical parametric QTL mapping methods. We model the varying genetic effects nonparametrically with the B-spline bases. Our model investigates gene-by-time interactions for RIX data in a very flexible nonparametric fashion. Simulation results indicate that the varying coefficient QTL mapping has higher power and mapping precision compared to parametric models when the assumption of constant genetic effects fails. We also apply a modified permutation procedure to control overall significance level. [ABSTRACT FROM AUTHOR]

Abstract: Fish pasteurellosis is an infectious disease that affects several fish species living in marine temperate waters. Its causative agent is the Gram-negative bacterium Photobacterium damselae subsp. piscicida (Phdp). Fish pasteurellosis represents a serious health problem for the majority of intensive sea bream hatcheries, with 90–100% mortality during disease outbreaks. Larvae and juveniles are the most susceptible stages. A potential strategy to prevent fish pasteurellosis is to select for animals that are genetically resistant to it. The aim of this work was to evaluate the genetic variation of disease resistance and growth in the gilthead sea bream (Sparus aurata L.). An experimental population, approximately 3500 animals originating from mass spawning of four broodstocks, was experimentally infected with a highly virulent strain of Phdp. Mortality was monitored daily for 19days. Upon completion of the challenge experiment, genetic profiles at nine microsatellite loci were obtained for 1753 animals and for all (256) broodstock fish, for parentage assignment. A subset of families defined by means of genetic analysis was used to calculate heritability for survival to infection. For body length (measured post-mortem), heritability was 0.38±0.07. The heritability for disease resistance varied depending on how resistance was defined. Heritability was 0.12±0.04 for days of survival post challenge, defined as a continuous trait, while it ranged from 0.45±0.04 to 0.18±0.08 for the binary trait dead/alive at a specific day. The genetic correlation between body length and survival was positive and significant (r =0.61±0.16). These results confirm the existence of genetic variation for growth and resistance against Phpd and highlight also the potential for selective breeding programs to improve these traits. [Copyright &y& Elsevier]

Abstract: Rapid development of genomics technology is providing new opportunities for genetic studies, including QTL mapping, in many aquaculture species. This paper investigates the strengths and limitations of QTL mapping designs for fish and shellfish under three different controlled breeding schemes. For each controlled breeding scheme, the potential and limitations are described for typical species and are illustrated by three different designs using interval mapping. The results show that, regardless of the species, the family structure is extremely important in experimental designs. The heritability of the QTL (controlled by its allele frequency and effect on the trait) also has an important impact on the power to detect QTL, while the overall polygenic heritability of the trait is less important. Marker density does not greatly affect the power when the distance between markers is less than 10 cM; but ideally spacing should not exceed 20 cM. For each of the systems studied, it is possible to design an experiment that would have an 80% power to detect a QTL of moderate effect (explaining between 1.5 and 5% of the trait variation) by genotyping 1000 or fewer individuals. [Copyright &y& Elsevier]

Combining global gene-expression profiling and genetic analysis of natural allelic variation (genetical genomics) has great potential in dissecting the genetic pathways underlying complex phenotypes. Efficient use of microarrays is paramount in experimental design as the cost of conducting this type of study is high. For those organisms where recombinant inbred lines are available for mapping, the "distant pair design" maximizes the number of informative contrasts over all marker loci. Here, we describe an extension of this design, named the "optimal pair design," for use with F2 crosses between outbred lines. The performance of this design is investigated by simulation and compared to several other two-color microarray designs. We show that, for a given number of microarrays, the optimal pair design outperforms all other designs considered for detection of expression quantitative trait loci (eQTL) with additive effects by linkage analysis. We also discuss the suitability of this design for outbred crosses in organisms with large genomes and for detection of dominance. [ABSTRACT FROM AUTHOR]

Field cress (Lepidium campestre L.), despite its potential as a sustainable alternative oilseed plant, has been underutilized, and no prior attempts to characterize the genome at the genetic or molecular cytogenetic level have been conducted. Genetic maps are the foundation for anchoring and orienting annotated genome assemblies and positional cloning of candidate genes. Our principal goal was to construct a genetic map using integrated approaches of genetic, comparative and cytogenetic map analyses. In total, 503 F2 interspecific hybrid individuals were genotyped using 7,624 single nucleotide polymorphism markers. Comparative analysis demonstrated that ~57% of the sequenced loci in L. campestre were congruent with Arabidopsis thaliana (L.) genome and suggested a novel karyotype, which predates the ancestral crucifer karyotype. Aceto-orcein chromosome staining and fluorescence in situ hybridization (FISH) analyses confirmed that L. campestre, L. heterophyllum Benth. and their hybrids had a chromosome number of 2n = 2x = 16. Flow cytometric analysis revealed that both species possess 2C roughly 0.4 picogram DNA. Integrating linkage and comparative maps with cytogenetic map analyses assigned two linkage groups to their particular chromosomes. Future work could incorporate FISH utilizing A. thaliana mapped BAC clones to allow the chromosomes of field cress to be identified reliably. [ABSTRACT FROM AUTHOR]

Most commercial laying hens suffer from sternum (keel) bone damage including deviations and fractures. X-raying hens, followed by segmenting and assessing the keel bone, is a key to automating the monitoring of keel bone condition. The aim of the current work is to train a deep learning model to segment the keel bone out of whole-body x-ray images. We obtained full-body x-ray images of laying hens (n = 1,051) and manually drew the outline of the keel bone on each image. Using the annotated images, a U-net model was then trained to segment the keel bone. The proposed model was evaluated using 5-fold cross validation. We obtained high segmentation accuracy (Dice coefficients of 0.88-0.90) repeatably over several validation folds. In conclusion, automatic segmentation of the keel bone from full-body x-ray images is possible with good accuracy. Segmentation is a requirement for automated measurements of keel geometry and density, which can subsequently be connected to susceptibility to keel deviations and fractures., Competing Interests: DISCLOSURES The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Keel bone damage, include deviations and fractures, is common in both white and brown laying hens, regardless of the housing system. Radiography for assessing birds' keel bones is was proposed by previous studies. However, radiographs show only 2 out of 3 dimensions of the dissected keel bones. The current study aimed to (1) investigate the association of radiographic optical density (keel and tibiotarsal) and geometry (keel) with dissected keel bone pathology. Previous studies suggested that keel bone fractures may result from internal pressure exerted by pelvic cavity contents. The current study also aimed to (2) investigate the potential associations between pelvic dimensions and measures of keel bone damage. A sample of 200 laying hens on a commercial farm were radiographed at 16, 29, 42, 55, and 68 weeks, and culled at the end of the laying period (week 74). The birds were examined post-mortem for pelvic dimensions and underwent whole-body radiography, followed by keel and tibiotarsal bone dissection and radiography, and keel bone scoring. The radiographs were used to estimate radiographic optical density (keel and tibiotarsal bone) and keel bone geometry (ratio of keel bone length to mid-depth). The method for on-farm radiography of laying hens, including live bird restraint, positioning for live keel imaging, and post-imaging measurements, was developed, tested, and found to be reproducible. The radiographs (1,116 images of 168 birds) and the respective measurements and post-mortem scores of keel bones are also provided for further development of radiographic metrics relevant to keel bone damage. Some longitudinal radiographic measurements of keel geometry (ratio of length to mid-depth) and optical density (keel and tibiotarsal) showed associations with the damage (deviations/fractures) observed on the dissected keel bones. The associations of keel damage were clearer with the radiographic keel geometry than with keel and tibiotarsal optical density, also clearer for the keel deviations than for keel fractures. The higher radiography ratio of keel length to mid-depth at weeks 42, 55 and 68 of age, the larger deviations size observed on the dissected keels at age of 74 weeks. The higher the tibiotarsal radiographic optical density at week 55 of age, the lower deviations size and fractures count observed on the dissected keels at age of 74 weeks. Pelvic dimensions showed a positive correlation with body weight, but a larger pelvic cavity was associated with increased keel bone damage. These findings lay the foundations for future use of on-farm radiography in identifying appropriate phenotypes for genetic selection for keel bone health., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Sallam, Göransson, Larsen, Alhamid, Johnsson, Wall, de Koning and Gunnarsson.)

African swine fever (ASF) is a global threat to animal health and food security. ASF is typically controlled by strict biosecurity, rapid diagnosis, and culling of affected herds. Much progress has been made in developing modified live virus vaccines against ASF. There is host variation in response to ASF infection in the field and under controlled conditions. To better understand the dynamics underlying this host differential morbidity, whole transcriptome profiling was carried out in twelve immunized and five sham immunized pigs. Seventeen MHC homozygous inbred Large white Babraham pigs were sampled at three time points before and after the challenge. The changes in the transcriptome profiles of infected animals were surveyed over time. In addition, the immunization effect on the host response was studied as well among the contrasts of all protection subgroups. The results showed two promising candidate genes to distinguish between recovered and non-recovered pigs after infection with a virulent African swine fever virus (ASFV) pre-infection: HTRA3 and GFPT2 (padj < 0.05). Variant calling on the transcriptome assemblies showed a two-base pair insertion into the ACOX3 gene closely located to HTRA3 that may regulate its expression as a putative genomic variant for ASF. Several significant DGEs, enriched gene ontology (GO) terms, and KEGG pathways at 1 day and 7 days post-infection, compared to the pre-infection, indicate a significant inflammation response immediately after ASF infection. The presence of the virus was confirmed by the mapping of RNA-Seq reads on two whole viral genome sequences. This was concordant with a higher virus load in the non-recovered animals 7 days post-infection. There was no transcriptome signature on the immunization at pre-infection and 1 day post-infection. More samples and data from additional clinical trials may support these findings., (© 2024. The Author(s).)

Impaired rennet coagulation properties in milk could lead to prolonged processing times and production losses. Heritability for milk coagulation has previously been estimated to be 0.28 to 0.45, indicating that genetic selection can be used to manipulate this trait. The CN proteins are expressed by the genes CSN1S1, CSN2, CSN1S2 , and CSN3 and are located on bovine chromosome 6. To better understand the effect of genetic variation in the CN genes on milk coagulation, blood and milk samples from 30 Swedish Red Dairy Cattle (RDC) with divergent coagulation properties were investigated. DNA from the 30 cows was sequenced for the CN genes to determine the theoretical AA sequence and to look for genetic variation in the untranslated regions. The aim is to confirm the protein genetic variants previously reported, while searching for additional genetic variation in the CN genes of 30 RDC. We observed genetic variation in 116 SNPs in the known CN genes where 10% of the SNPs are exon variants and the remaining 90% are intron variants. A total of 2.5% of the SNPs are found in the 5'- or 3'-untranslated region (UTR) regions of the exons; 2% are synonymous variants and 6% are missense variants that concurred with the known protein variants for CSN1S1, CSN2 , and CSN3 . Furthermore, 6% of the SNPs are splice polypyrimidine tract intron variants. The 2 genetic variants in the 5'- and 3'-UTR in CSN1S1 and CSN3 are found with protein variants CSN1S1C and CSN3B . Because both UTR variants are associated with gain and loss of micro RNA and transcription factors, this could explain differences in expression of the genetic protein variants. Preliminary chi-squared analysis and comparison with previous GWAS studies showed potential connections between the identified SNPs and coagulation properties of milk. By advancing the knowledge of the connection between the DNA sequence and the functional properties of the CN proteins, we hope to learn more about the cheese coagulation properties of milk from RDC., (© 2024.)

Rainbow trout is one of the most popular aquaculture species worldwide, with a long history of domestication. However, limited information exists about the genetic diversity of farmed rainbow trout populations globally, with most available reports relying on low-throughput genotyping technologies. Notably, no information exists about the genetic diversity status of farmed rainbow trout in Sweden. Double-digest restriction-site-associated DNA sequencing was performed on more than 500 broodfish from two leading producers in Sweden and from the country's national breeding program. Following the detection of single nucleotide polymorphisms (SNPs), genetic diversity was studied by using either individual SNPs (n = 8680; one SNP retained per 300 bp sequence reads) or through SNP haplotypes (n = 20 558; all SNPs retained in 300 bp sequence reads). Similar amounts of genetic diversity were found amongst the three populations when individual SNPs were used. Furthermore, principal component analysis and discriminant analysis of principal components suggested two genetic clusters with the two industry populations grouped together. Genetic differentiation based on the F ST fixation index was ~0.01 between the industry populations and ~0.05 when those were compared with the breeding program. Preliminary estimates of effective population size (N e ) and inbreeding (based on runs of homozygosity; F ROH ) were similar amongst the three populations (N e  ≈ 50-80; median F ROH  ≈ 0.11). Finally, the haplotype-based analysis suggested that animals from the breeding program had higher shared coancestry levels than those from the other two populations. Overall, our study provides novel insights into the genetic diversity and structure of Sweden's three main farmed rainbow trout populations, which could guide their future management., (© 2023 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Background: Bone damage has welfare and economic impacts on modern commercial poultry and is known as one of the major challenges in the poultry industry. Bone damage is particularly common in laying hens and is probably due to the physiological link between bone and the egg laying process. Previous studies identified and validated quantitative trait loci (QTL) for bone strength in White Leghorn laying hens based on several measurements, including bone composition measurements on the cortex and medulla of the tibia bone. In a previous pedigree-based analysis, bone composition measurements showed heritabilities ranging from 0.18 to 0.41 and moderate to strong genetic correlations with tibia strength and density. Bone composition was measured using infrared spectroscopy and thermogravimetry. The aim of this study was to combine these bone composition measurements with genotyping data via a genome-wide association study (GWAS) to investigate genetic markers that contribute to genetic variance in bone composition in Rhode Island Red laying hens. In addition, we investigated the genetic correlations between bone composition and bone strength., Results: We found novel genetic markers that are significantly associated with cortical lipid, cortical mineral scattering, medullary organic matter, and medullary mineralization. Composition of the bone organic matter showed more significant associations than bone mineral composition. We also found interesting overlaps between the GWAS results for tibia composition traits, particularly for cortical lipid and tibia strength. Bone composition measurements by infrared spectroscopy showed more significant associations than thermogravimetry measurements. Based on the results of infrared spectroscopy, cortical lipid showed the highest genetic correlations with tibia density, which was negative (- 0.20 ± 0.04), followed by cortical CO3/PO4 (0.18 ± 0.04). Based on the results of thermogravimetry, medullary organic matter% and mineral% showed the highest genetic correlations with tibia density (- 0.25 ± 0.04 and 0.25 ± 0.04, respectively)., Conclusions: This study detected novel genetic associations for bone composition traits, particularly those involving organic matter, that could be used as a basis for further molecular genetic investigations. Tibia cortical lipids displayed the strongest genetic associations of all the composition measurements, including a significantly high genetic correlation with tibia density and strength. Our results also highlighted that cortical lipid may be a key measurement for further avian bone studies., (© 2023. The Author(s).)

Background: Red dairy cattle breeds have an important role in the European dairy sector because of their functional characteristics and good health. Extensive pedigree information is available for these breeds and provides a unique opportunity to examine their population structure, such as effective population size, depth of the pedigree, and effective number of founders and ancestors, and inbreeding levels. Animals with the highest genetic contributions were identified. Pedigree data included 9,073,403 animals that were born between 1900 and 2019 from Denmark, Finland, Germany, Latvia, Lithuania, the Netherlands, Norway, Poland, and Sweden, and covered 32 breeds. The numerically largest breeds were Red Dairy Cattle and Meuse-Rhine-Yssel., Results: The deepest average complete generation equivalent (9.39) was found for Red Dairy Cattle in 2017. Mean pedigree completeness ranged from 0.6 for Finncattle to 7.51 for Red Dairy Cattle. An effective population size of 166 animals was estimated for the total pedigree and ranged from 35 (Rotes Höhenvieh) to 226 (Red Dairy Cattle). Average generation intervals were between 5 and 7 years. The mean inbreeding coefficient for animals born between 1960 and 2018 was 1.5%, with the highest inbreeding coefficients observed for Traditional Angler (4.2%) and Rotes Höhenvieh (4.1%). The most influential animal was a Dutch Meuse-Rhine-Yssel bull born in 1960. The mean inbreeding level for animals born between 2016 and 2018 was 2% and highest for the Meuse-Rhine-Yssel (4.64%) and Rotes Hohenvieh breeds (3.80%)., Conclusions: We provide the first detailed analysis of the genetic diversity and inbreeding levels of the European red dairy cattle breeds. Rotes Höhenvieh and Traditional Angler have high inbreeding levels and are either close to or below the minimal recommended effective population size, thus it is necessary to implement tools to monitor the selection process in order to control inbreeding in these breeds. Red Dairy Cattle, Vorderwälder, Swedish Polled and Hinterwälder hold more genetic diversity. Regarding the Meuse-Rhine-Yssel breed, given its decreased population size, increased inbreeding and low effective population size, we recommend implementation of a breeding program to prevent further loss in its genetic diversity., (© 2022. The Author(s).)

Background: A sufficient IgG content in the colostrum is essential for the newborn calf, as it provides passive immunity which substantially affects the probability of survival during rearing. Failure of passive transfer (FPT) occurs when a calf does not absorb enough antibodies from the colostrum and is defined by an IgG concentration in calf serum lower than 10 g/L. Apart from delayed access to colostrum, FPT can be due to a low production of IgG in the mother or poor IgG absorption by the calf. The aim of this study was to estimate the genetic background of antibody levels and indicator traits for antibodies in the colostrum and calf serum, and their correlation with milk production., Results: Colostrum data were available for 1340 dairy cows with at least one calving and calf serum data were available for 886 calves from these cows. Indicator traits for antibody concentrations were estimated using refractometry (a digital Brix refractometer for colostrum and an optical refractometer for serum), and enzyme-linked immunosorbent assays (ELISA) were used to determine the levels of total IgG and natural antibodies (NAb) of various antibody isotypes in the colostrum and calf serum. Colostrum traits had heritabilities ranging from 0.16 to 0.31 with repeatabilities ranging from 0.21 to 0.55. Brix percentages had positive genetic correlations with all colostrum antibody traits including total IgG (0.68). Calf serum antibody concentrations had heritabilities ranging from 0.25 to 0.59, with a significant maternal effect accounting for 17 to 27% of the variance. When later in life calves produced their first lactation, the lactation average somatic cell score was found to be negatively correlated with NAb levels in calf serum., Conclusions: Our results suggest that antibody levels in the colostrum and calf serum can be increased by means of selection., (© 2022. The Author(s).)

Nile tilapia is a key aquaculture species with one of the highest production volumes globally. Genetic improvement of feed efficiency via selective breeding is an important goal, and genomic selection may expedite this process. The aims of this study were to 1 ) dissect the genetic architecture of feed-efficiency traits in a Nile tilapia breeding population, 2 ) map the genomic regions associated with these traits and identify candidate genes, 3 ) evaluate the accuracy of breeding value prediction using genomic data, and 4 ) assess the impact of the genetic marker density on genomic prediction accuracies. Using an experimental video recording trial, feed conversion ratio (FCR), body weight gain (BWG), residual feed intake (RFI) and feed intake (FI) traits were recorded in 40 full-sibling families from the GIFT (Genetically Improved Farmed Tilapia ) Nile tilapia breeding population. Fish were genotyped with a ThermoFisher Axiom 65 K Nile tilapia SNP array. Significant heritabilities, ranging from 0.12 to 0.22, were estimated for all the assessed traits using the genomic relationship matrix. A negative but favourable genetic correlation was found between BWG and the feed-efficiency related traits; -0.60 and -0.63 for FCR and RFI, respectively. While the genome-wide association analyses suggested a polygenic genetic architecture for all the measured traits, there were significant QTL identified for BWG and FI on chromosomes seven and five respectively. Candidate genes previously found to be associated with feed-efficiency traits were located in these QTL regions, including ntrk3a , ghrh and eif4e3 . The accuracy of breeding value prediction using the genomic data was up to 34% higher than using pedigree records. A SNP density of approximately 5,000 SNPs was sufficient to achieve similar prediction accuracy as the full genotype data set. Our results highlight the potential of genomic selection to improve feed efficiency traits in Nile tilapia breeding programmes., Competing Interests: JB was employed by the company WorldFish. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Barría, Benzie, Houston, De Koning and de Verdal.)

Routine implementation of genomic information for guiding selection decisions is not yet common in the majority of aquaculture species. Reduced representation sequencing approaches offer a cost-effective solution for obtaining genome-wide information in species with a limited availability of genomic resources. In the current study, we implemented double-digest restriction site-associated DNA sequencing (ddRAD-seq) on an Arctic charr strain with the longest known history of selection (approximately 40 years) aiming to improve selection decisions. In total, 1730 animals reared at four different farms in Sweden and spanning from year classes 2013-2017 were genotyped using ddRAD-seq. Approximately 5000 single nucleotide polymorphisms (SNPs) were identified, genetic diversity-related metrics were estimated, and genome-wide association studies (GWAS) for body length at different time points and age of sexual maturation were conducted. Low genetic differentiation amongst animals from the different farms was observed based on both the results from pairwise F st values and principal component analysis (PCA). The existence of associations was investigated between the mean genome-wide heterozygosity of each full-sib family (year class 2017) and the corresponding inbreeding coefficient or survival to the eyed stage. A moderate correlation (-0.33) was estimated between the mean observed heterozygosity of each full-sib family and the corresponding inbreeding coefficient, while no linear association was obtained with the survival to the eyed stage. GWAS did not detect loci with major effect for any of the studied traits. However, genomic regions explaining more than 1% of the additive genetic variance for either studied traits were suggested across 14 different chromosomes. Overall, key insights valuable for future selection decisions of Arctic charr have been obtained, suggesting ddRAD as an attractive genotyping platform for obtaining genome-wide information in a cost-effective manner., Competing Interests: The authors declare no conflict of interest., (© 2021 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Background: Meiotic recombination results in the exchange of genetic material between homologous chromosomes. Recombination rate varies between different parts of the genome, between individuals, and is influenced by genetics. In this paper, we assessed the genetic variation in recombination rate along the genome and between individuals in the pig using multilocus iterative peeling on 150,000 individuals across nine genotyped pedigrees. We used these data to estimate the heritability of recombination and perform a genome-wide association study of recombination in the pig., Results: Our results confirmed known features of the recombination landscape of the pig genome, including differences in genetic length of chromosomes and marked sex differences. The recombination landscape was repeatable between lines, but at the same time, there were differences in average autosome-wide recombination rate between lines. The heritability of autosome-wide recombination rate was low but not zero (on average 0.07 for females and 0.05 for males). We found six genomic regions that are associated with recombination rate, among which five harbour known candidate genes involved in recombination: RNF212, SHOC1, SYCP2, MSH4 and HFM1., Conclusions: Our results on the variation in recombination rate in the pig genome agree with those reported for other vertebrates, with a low but nonzero heritability, and the identification of a major quantitative trait locus for recombination rate that is homologous to that detected in several other species. This work also highlights the utility of using large-scale livestock data to understand biological processes.

Five introduced strains of Nile tilapia (Oreochromis niloticus) were tested for growth performance both in fresh- and brackish-water (2 salinity units) environments for 56 days. The BIG NIN, GIFT, Chitralada, "Ruvu Farm" and Silver YY strains with initial mean average weight (± standard error) of 96.4 ± 6.90 g, 104.1 ± 7.19 g, 137.2 ± 7.21 g, 53.2 ± 6.98 g and 95.3 ± 7.11 g, respectively were used. Individuals were tagged and pooled in hapas (12 m × 8.5 m × 2 m each), aligned into different ponds (20 m × 20 m each). Stocking density of 5 fish/m 2 and 350 g/kg crude protein diet were used. Overall, the average weight gain for GIFT strain was 7.5%, 32%, 45% and 86.5% higher than BIG NIN, Chitralada, "Ruvu Farm" and Silver YY strains, respectively, across both environments. All strains performed significantly better (p < 0.05) when reared in brackish-water than their respective counterparts in freshwater, except for the BIG NIN strain. The morphometric correlations for all strains in both environments ranged from moderate (0.50) to strong positive (0.92). The GIFT strain demonstrated superior growth and genotype by environment interaction was weak and not important to be prioritized in breeding programs.

Background: The physiological adaptations that have evolved for egg laying make hens susceptible to bone fractures and keel bone damage. In modern laying hen breeds, longer periods of egg laying could result in a greater risk of poor bone quality, and selection for increased egg production has frequently been stated to be a cause. However, the existing literature does not support this hypothesis. To test the hypothesis that egg production is associated with quality, breaking strength and density of bone, genetic correlations between these traits were estimated in White Leghorn and Rhode Island Red breeds. Genetic correlations of cortical and medullary bone material chemical properties with bone quality were also estimated, in order to identify methods to improve bone quality with appropriately targeted measurement of key traits., Results: Estimates of heritability for bone quality traits were moderate (0.19-0.59) for both White Leghorn and Rhode Island Red breeds, except for the keel bone trait, which had a heritability estimate equal to zero. There was no evidence for genetic or phenotypic relationships between post-peak egg production and bone quality. In the White Leghorn breed, the estimate of the genetic correlation between pre-peak production/age at first egg and bone quality was significant and negative (- 0.7 to - 0.4). Estimates of heritability of thermogravimetric measurements of tibial medullary bone mineralisation were significant (0.18-0.41), as were estimates of their genetic correlations with tibia breaking strength and density (0.6-0.9)., Conclusions: The low genetic correlation of post-peak egg production with bone quality suggests that selection for increased persistency of egg production may not adversely affect bone quality. Onset of puberty and mineralisation of the medullary bone, which is a specialised adaptation for egg laying, were identified as important factors associated with the quality of the skeleton later during egg production. These are traits for which genetic, as well as environmental and management factors can positively impact the overall quality of the skeleton of laying hens.

Lepidium campestre (L.) or field cress is a multifaceted oilseed plant, which is not yet domesticated. Moreover, the molecular and genetic mechanisms underlying the domestication traits of field cress remain largely elusive. The overarching goal of this study is to identify quantitative trait loci (QTL) that are fundamental for domestication of field cress. Mapping and dissecting quantitative trait variation may provide important insights into genomic trajectories underlying field cress domestication. We used 7624 single nucleotide polymorphism (SNP) markers for QTL mapping in 428 F 2 interspecific hybrid individuals, while field phenotyping was conducted in F 2:3 segregating families. We applied multiple QTL mapping algorithms to detect and estimate the QTL effects for seven important domestication traits of field cress. Verification of pod shattering across sites revealed that the non-shattering lines declined drastically whereas the shattering lines increased sharply, possibly due to inbreeding followed by selection events. In total, 1461 of the 7624 SNP loci were mapped to eight linkage groups (LGs), spanning 571.9 cM map length. We identified 27 QTL across all LGs of field cress genome, which captured medium to high heritability, implying that genomics-assisted selection could deliver domesticated lines in field cress breeding. The use of high throughput genotyping can accelerate the process of domestication in novel crop species. This is the first QTL mapping analysis in the field cress genome that may lay a foundational framework for positional or functional QTL cloning, introgression as well as genomics-assisted breeding in field cress domestication.

Background: Skeletal damage is a challenge for laying hens because the physiological adaptations required for egg laying make them susceptible to osteoporosis. Previously, we showed that genetic factors explain 40% of the variation in end of lay bone quality and we detected a quantitative trait locus (QTL) of large effect on chicken chromosome 1. The aim of this study was to combine data from the commercial founder White Leghorn population and the F2 mapping population to fine-map this QTL and understand its function in terms of gene expression and physiology., Results: Several single nucleotide polymorphisms on chromosome 1 between 104 and 110 Mb (galGal6) had highly significant associations with tibial breaking strength. The alternative genotypes of markers of large effect that flanked the region had tibial breaking strengths of 200.4 vs. 218.1 Newton (P < 0.002) and, in a subsequent founder generation, the higher breaking strength genotype was again associated with higher breaking strength. In a subsequent generation, cortical bone density and volume were increased in individuals with the better bone genotype but with significantly reduced medullary bone quality. The effects on cortical bone density were confirmed in a further generation and was accompanied by increased mineral maturity of the cortical bone as measured by infrared spectrometry and there was evidence of better collagen cross-linking in the cortical bone. Comparing the transcriptome of the tibia from individuals with good or poor bone quality genotypes indicated four differentially-expressed genes at the locus, one gene, cystathionine beta synthase (CBS), having a nine-fold higher expression in the genotype for low bone quality. The mechanism was cis-acting and although there was an amino-acid difference in the CBS protein between the genotypes, there was no difference in the activity of the enzyme. Plasma homocysteine concentration, the substrate of CBS, was higher in the poor bone quality genotype., Conclusions: Validated markers that predict bone strength have been defined for selective breeding and a gene was identified that may suggest alternative ways to improve bone health in addition to genetic selection. The identification of how genetic variants affect different aspects of bone turnover shows potential for translational medicine.

Exploring dominance variance and loci contributing to dominance variation is important to understand the genetic architecture behind quantitative traits. The objectives of this study were i) to estimate dominance variances, ii) to detect quantitative trait loci (QTL) with dominant effects, and iii) to evaluate the power and the precision of identifying loci with dominance effect through post-hoc simulations, with applications for female fertility in Danish Holstein cattle. The female fertility records analyzed were number of inseminations (NINS), days from calving to first insemination (ICF), and days from the first to last insemination (IFL), covering both abilities to recycle and to get pregnant in the female reproductive cycle. There were 3,040 heifers and 4,483 cows with both female fertility records and Illumina BovineSNP50 BeadChip genotypes (35,391 single nucleotide polymorphisms (SNP) after quality control). Genomic best linear unbiased prediction (BLUP) models were used to estimate additive and dominance genetic variances. Linear mixed models were used for association analyses. A post-hoc simulation study was performed using genotyped heifers' data. In heifers, estimates of dominance genetic variances for female fertility traits were larger than additive genetic variances, but had large standard errors. The variance components for fertility traits in cows could not be estimated due to non-convergence of the statistical model. In total, five QTL located on chromosomes 9, 11 (2 QTL), 19, and 28 were identified and all of them showed both additive and dominance genetic effects. Among them, the SNP rs29018921 on chromosome 9 is close to a previously identified QTL in Nordic Holstein for interval between first and last insemination. This SNP is located in the 3' untranslated region of gene peptidylprolyl isomerase like 4 (PPIL4), which was shown to be associated with milk production traits in US Holstein cattle but not known for fertility-related functions. Simulations indicated that the current sample size had limited power to detect QTL with dominance effects for female fertility probably due to low QTL variance. More females need to be genotyped to achieve reliable mapping of QTL with dominance effects for female fertility.

Background: In this paper, we simulate deleterious load in an animal breeding program, and compare the efficiency of genome editing and selection for decreasing it. Deleterious variants can be identified by bioinformatics screening methods that use sequence conservation and biological prior information about protein function. However, once deleterious variants have been identified, how can they be used in breeding?, Results: We simulated a closed animal breeding population that is subject to both natural selection against deleterious load and artificial selection for a quantitative trait representing the breeding goal. Deleterious load was polygenic and was due to either codominant or recessive variants. We compared strategies for removal of deleterious alleles by genome editing (RAGE) to selection against carriers. When deleterious variants were codominant, the best strategy for prioritizing variants was to prioritize low-frequency variants. When deleterious variants were recessive, the best strategy was to prioritize variants with an intermediate frequency. Selection against carriers was inefficient when variants were codominant, but comparable to editing one variant per sire when variants were recessive., Conclusions: Genome editing of deleterious alleles reduces deleterious load, but requires the simultaneous editing of multiple deleterious variants in the same sire to be effective when deleterious variants are recessive. In the short term, selection against carriers is a possible alternative to genome editing when variants are recessive. Our results suggest that, in the future, there is the potential to use RAGE against deleterious load in animal breeding.

Bone fracture in egg laying hens is a growing welfare and economic concern in the industry. Although environmental conditions and management (especially nutrition) can exacerbate it, the primary cause of bone weakness and the resulting fractures is believed to have a genetic basis. To test this hypothesis, we performed a genome-wide association study to identify the loci associated with bone strength in laying hens. Genotype and phenotype data were obtained from 752 laying hens belonging to the same pure line population. These hens were genotyped for 580,961 SNPs, with 232,021 SNPs remaining after quality control. Each of the SNPs were tested for association with tibial breaking strength using the family-based score test for association. A total of 52 SNPs across chromosomes 1, 3, 8, and 16 were significantly associated with tibial breaking strength with the genome-wide significance threshold set as a corrected P value of 10e-5. Based on the local linkage disequilibrium around the significant SNPs, 5 distinct and novel QTLs were identified on chromosomes 1 (2 QTLs), 3 (1 QTL), 8 (1 QTL) and 16 (1 QTL). The strongest association was detected within the QTL region on chromosome 8, with the most significant SNP having a corrected P value of 4e-7. A number of candidate genes were identified within the QTL regions, including the BRD2 gene that is required for normal bone physiology. Bone-related pathways involving some of the genes were also identified including chloride channel activity, which regulates bone reabsorption, and intermediate filament organization, which plays a role in the regulation of bone mass. Our result supports previous studies that suggest that bone strength is highly regulated by genetics. It is therefore possible to reduce bone fractures in laying hens through genetic selection and ultimately improve hen welfare.

Large genomic deletions are potential candidate for loss-of-function, which could be lethal as homozygote. Analysing whole genome data of 175 cattle, we report 8,480 large deletions (199 bp-773 KB) with an overall false discovery rate of 8.8%; 82% of which are novel compared with deletions in the dbVar database. Breakpoint sequence analyses revealed that majority (24 of 29 tested) of the deletions contain microhomology/homology at breakpoint, and therefore, most likely generated by microhomology-mediated end joining. We observed higher differentiation among breeds for deletions in some genic-regions, such as ABCA12, TTC1, VWA3B, TSHR, DST/BPAG1, and CD1D. The genes overlapping deletions are on average evolutionarily less conserved compared with known mouse lethal genes (P-value = 2.3 × 10-6). We report 167 natural gene knockouts in cattle that are apparently nonessential as live homozygote individuals are observed. These genes are functionally enriched for immunoglobulin domains, olfactory receptors, and MHC classes (FDR = 2.06 × 10-22, 2.06 × 10-22, 7.01 × 10-6, respectively). We also demonstrate that deletions are enriched for health and fertility related quantitative trait loci (2-and 1.5-fold enrichment, Fisher's P-value = 8.91 × 10-10 and 7.4 × 10-11, respectively). Finally, we identified and confirmed the breakpoint of a ∼525 KB deletion on Chr23:12,291,761-12,817,087 (overlapping BTBD9, GLO1 and DNAH8), causing stillbirth in Nordic Red Cattle., (© The Author 2017. Published by Oxford University Press on behalf of Kazusa DNA Research Institute.)

The number of functional teats is an important selection criterion in pig breeding. Inherited defects of the udder, such as the inverted teat, do have a considerable negative impact on the nursing ability of the sow. To investigate the genetic background of this defect and the number of functional teats in Swedish maternal lines, samples from 230 Yorkshire pigs were selected for genotyping using the PorcineSNP60K BeadChip (Illumina Inc.), each pig with at least one inverted teat was matched with one non-affected pig (fullsib or pairs with matching herd and gender). A genome-wide association study on these 230 pigs was performed using the two-step approach implemented in GenABEL using 46,652 single nucleotide polymorphisms across all autosomes and the X chromosome. A number of significant regions were identified for the inverted teat defect on chromosomes 2, 10, and 18. Many of the regions associated with the number of functional teats were located in the same or close regions, except two associated markers on the X chromosome and one on chromosome 3. We identified some of the regions on chromosomes previously reported in one linkage and one gene expression study. We conclude, despite being able to suggest new candidate genes, that further studies are needed to better understand the biologic background of the teat development. Despite the in-depth comparison of identified regions for the inverted teat defect done here, more studies are required to allow a clear identification of genetic regions relevant for this defect across many pig populations.

Imputation is a cost-effective approach to augment marker data for genomic selection and genome-wide association studies. However, most imputation studies have focused on autosomes. Here, we assessed the imputation of markers on the X chromosome in Holstein cattle for nongenotyped animals and animals genotyped with low-density (Illumina BovineLD, Illumina Inc., San Diego, CA) chips, using animals genotyped with medium-density (Illumina BovineSNP50) chips. A total of 26,884 genotyped Holstein individuals genotyped with medium-density chips were used in this study. Imputation was carried out using FImpute V2.2. The following parameters were examined: treating the pseudoautosomal region as autosomal or as X specific, different sizes of reference groups, different male/female proportions in the reference group, and cumulated degree of relationship between the reference group and target group. The imputation accuracy of markers on the X chromosome was improved if the pseudoautosomal region was treated as autosomal. Increasing the proportion of females in the reference group improved the imputation accuracy for the X chromosome. Imputation for nongenotyped animals in general had lower accuracy compared with animals genotyped with the low-density single nucleotide polymorphism array. In addition, higher cumulative pedigree relationships between the reference group and the target animal led to higher imputation accuracy. In the future, better marker coverage of the X chromosome should be developed to facilitate genomic studies involving the X chromosome., (Copyright © 2016 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Non-coagulating (NC) milk, defined as milk not coagulating within 40 min after rennet-addition, can have a negative influence on cheese production. Its prevalence is estimated at 18% in the Swedish Red (SR) cow population. Our study aimed at identifying genomic regions and causal variants associated with NC milk in SR cows, by doing a GWAS using 777k SNP genotypes and using imputed sequences to fine map the most promising genomic region. Phenotypes were available from 382 SR cows belonging to 21 herds in the south of Sweden, from which individual morning milk was sampled. NC milk was treated as a binary trait, receiving a score of one in case of non-coagulation within 40 min. For all 382 SR cows, 777k SNP genotypes were available as well as the combined genotypes of the genetic variants of αs1-β-κ-caseins. In addition, whole-genome sequences from the 1000 Bull Genome Consortium (Run 3) were available for 429 animals of 15 different breeds. From these sequences, 33 sequences belonged to SR and Finish Ayrshire bulls with a large impact in the SR cow population. Single-marker analyses were run in ASReml using an animal model. After fitting the casein loci, 14 associations at -Log10(P-value) > 6 identified a promising region located on BTA18. We imputed sequences to the 382 genotyped SR cows using Beagle 4 for half of BTA18, and ran a region-wide association study with imputed sequences. In a seven mega base-pairs region on BTA18, our strongest association with NC milk explained almost 34% of the genetic variation in NC milk. Since it is possible that multiple QTL are in strong LD in this region, 59 haplotypes were built, genetically differentiated by means of a phylogenetic tree, and tested in phenotype-genotype association studies. Haplotype analyses support the existence of one QTL underlying NC milk in SR cows. A candidate gene of interest is the VPS35 gene, for which one of our strongest association is an intron SNP in this gene. The VPS35 gene belongs to the mammary gene sets of pre-parturient and of lactating cows.

Using targeted genetical genomics, a quantitative trait locus (QTL) affecting the initial postmortem pH value of chicken breast muscle (Pectoralis major) on chromosome 1 (GGA1) recently was fine-mapped. Thirteen genes were present in the QTL region of approximately 1 Mb. In this study, 10 birds that were inferred to be homozygous for either the high (QQ) or low (qq) QTL allele were selected for resequencing. After enrichment for 1 Mb around the QTL region, >500 × coverage for the QTL region in each of the 10 birds was obtained. In total 5056 single-nucleotide polymorphisms (SNPs) were identified for which the genotypes were consistent with one of the QTL genotypes. We used custom tools to identify putative causal mutations in the mapped QTL region from these SNPs. Four nonsynonymous SNPs differentiating the two QTL genotype groups were identified within four local genes (PRDX4, EIF2S3, PCYT1B, and E1BTD2). Although these are likely candidate SNPs to explain the QTL effect, 54 additional consensus SNPs were detected within gene-related regions (untranslated regions, splicing sites CpG island, and promoter regions) for the QQ birds and 71 for the qq birds. These could also play a role explaining the observed QTL effect. The results provide an important step for prioritizing among a large amount of candidate mutations and significantly contribute to the understanding of the genetic mechanisms affecting the initial postmortem pH value of chicken muscle., (Copyright © 2015 Li et al.)