Inclusive inheritance for residual feed intake in pigs and rabbits.

Non-genetic information (epigenetic, microbiota, behaviour) that results in different phenotypes in animals can be transmitted from one generation to the next and thus is potentially involved in the inheritance of traits. However, in livestock species, animals are selected based on genetic inheritance only. The objective of the present study was to determine whether non-genetic inherited effects play a role in the inheritance of residual feed intake (RFI) in two species: pigs and rabbits. If so, the path coefficients of the information transmitted from sire and dam to offspring would differ from the expected transmission factor of 0.5 that occurs if inherited information is of genetic origin only. Two pigs (pig1, pig2) and two rabbits (rabbit1, rabbit2) datasets were used in this study (1,603, 3,901, 5,213 and 4,584 records, respectively). The test of the path coefficients to 0.5 was performed for each dataset using likelihood ratio tests (null model: transmissibility model with both path coefficients equal to 0.5, full model: unconstrained transmissibility model). The path coefficients differed significantly from 0.5 for one of the pig datasets (pig2). Although not significant, we observed, as a general trend, that sire path coefficients of transmission were lower than dam path coefficients in three of the datasets (0.46 vs 0.53 for pig1, 0.39 vs 0.44 for pig2 and 0.38 vs 0.50 for rabbit1). These results suggest that phenomena other than genetic sources of inheritance explain the phenotypic resemblance between relatives for RFI, with a higher transmission from the dam's side than from the sire's side.

Genetic parameters of resistance to pasteurellosis using novel response traits in rabbits.

BACKGROUND: Pasteurellosis (Pasteurella infection) is one of the most common bacterial infections in rabbits on commercial farms and in laboratory facilities. Curative treatments using antibiotics are only partly efficient, with frequent relapses. Breeding rabbits for improved genetic resistance to pasteurellosis is a sustainable alternative approach. In this study, we infected 964 crossbred rabbits from six sire lines experimentally with Pasteurella multocida. After post-mortem examination and bacteriological analyses, abscess, bacteria, and resistance scores were derived for each rabbit based on the extent of lesions and bacterial dissemination in the body. This is the first study to use such an experimental design and response traits to measure resistance to pasteurellosis in a rabbit population. We investigated the genetic variation of these traits in order to identify potential selection criteria. We also estimated genetic correlations of resistance to pasteurellosis in the experimental population with traits that are under selection in the breeding populations (number of kits born alive and weaning weight). RESULTS: Heritability estimates for the novel response traits, abscess, bacteria, and resistance scores, ranged from 0.08 (± 0.05) to 0.16 (± 0.06). The resistance score showed very strong negative genetic correlation estimates with abscess (- 0.99 ± 0.05) and bacteria scores (- 0.98 ± 0.07). A very high positive genetic correlation of 0.99 ± 0.16 was estimated between abscess and bacteria scores. Estimates of genetic correlations of the resistance score with average daily gain traits for the first and second week after inoculation were 0.98 (± 0.06) and 0.70 (± 0.14), respectively. Estimates of genetic correlations of the disease-related traits with average daily gain pre-inoculation were favorable but with high standard errors. Estimates of genetic and phenotypic correlations of the disease-related traits with commercial selection traits were not significantly different from zero. CONCLUSIONS: Disease response traits are heritable and are highly correlated with each other, but do not show any significant genetic correlations with commercial selection traits. Thus, the prevalence of pasteurellosis could be decreased by selecting more resistant rabbits on any one of the disease response traits with a limited impact on the selection traits, which would allow implementation of a breeding program to improve resistance to pasteurellosis in rabbits.

Genome Sequences of 17 Pasteurella multocida Strains Involved in Cases of Rabbit Pasteurellosis.

This article reports draft genome sequences of 17 Pasteurella multocida strains isolated from naturally infected rabbits. The total lengths of the assembled contigs ranged between 2.21 and 2.48 Mb, and the total number of genes detected on the contigs ranged between 2,088 and 2,416.

Estimating direct genetic and maternal effects affecting rabbit growth and feed efficiency with a factorial design.

The aim of this experiment was to evaluate the significance of neonatal environment on feed efficiency. For that purpose, rabbits from a line selected for residual feed intake (RFI) during 10 generations (G10 kits) were cross-fostered with non-selected control does (i.e., G0 line), and reciprocally. In parallel, sibs were fostered by mothers from their original line. Nine hundred animals were raised in individual (N = 456) or collective (N = 320) cages. Traits analysed in this study were body weight at 32 days and at 63 days, average daily gain (ADG), feed intake between weaning and 63 days (FI), feed conversion ratio (FCR) and RFI. The maternal environment offered by does from the line selected for RFI deteriorated the FCR of the kits, independently of their line of origin, during fattening (+0.08 ± 0.02) compared to FCR of kits nursed by G0 does. The line, the type of housing and the batch were significant effects for all the measured traits: G10 kits were lighter than their G0 counterparts at 32 days (-82.9 ± 9 g, p < 0.0001) and at 63 days (-161 ± 16 g, p < 0.0001). They also had a lower ADG (-2.36 ± 0.36 g/day, p < 0.0001), RFI (-521 ± 24 g/day, p < 0.0001) and a lower FI (-855 ± 31 g, p < 0.0001), resulting in a more desirable feed efficiency (FCR: -0.35 ± 0.02). There was no significant difference in the contrast of G10 and G0 performances between collective and individual/digestive cages (p > 0.22): -2.35 g/day versus 2.94 g/day for ADG, -0.39 versus -0.40 for FCR, -577 g versus -565 g for RFI and -879 g versus -859 g for FI, respectively). Thus, no genotype-by-environment (housing) interaction is expected at the commercial level, that is, no re-ranking of the animals due to collective housing.

Genetic Parameters for Resistance to Non-specific Diseases and Production Traits Measured in Challenging and Selection Environments; Application to a Rabbit Case.

Breeding for disease resistance is a challenging but increasingly necessary objective to overcome the issues with the reduced use of antibiotics and growing concern for animal welfare while limiting economic losses. However, implementing such strategies is a complex process because animals face numerous diseases, and the environments on selection farms differ from those on commercial farms. We evaluated whether selection for resistance to non-specific diseases based on a single visual record in selection (S) and challenging (Ch) environments is possible. Records from 23,773 purebred rabbits born between 2012 and 2016 were used in this study. After weaning (at 32 days of age), 17,712 rabbits were raised in the S environment and 6,061 sibs were raised in the Ch environment. Clinical signs of disease were recorded for all animals at the end of the test, at a single time point, at 70 or 80 days of age. The causes of mortality occurring before the end of the test were also recorded. Three disease traits were analyzed: signs of respiratory disease, signs of digestive disease, and a composite trait (Resist) taking into account signs of digestive, respiratory and various infectious diseases. This latter composite trait is proposed to capture the global resistance to disease. All disease traits were binary, with 0 being the absence of symptoms. Two production traits were also recorded: the number of kits born alive (4,121 litters) and the weaning weight (13,090 rabbits). Disease traits were analyzed with animal threshold models, assuming that traits are different in the two environments. Bivariate analyses were carried out using linear animal models. The heritabilities of the disease traits ranged from 0.04 ± 0.01 to 0.11 ± 0.03. The genetic correlations between disease traits in both environments were below unity (≤ 0.84), indicating genotype by environment interactions. Most of the genetic correlations between disease and production traits were not significantly different from zero, except between the weaning weight and Resist_S, with a favorable correlation of -0.34 ± 0.12. Given these genetic parameters, for the same level of exposure of rabbits to pathogens, the expected response to selection is a reduction of disease incidence of 4-6% per generation.

Multiple-trait structured antedependence model to study the relationship between litter size and birth weight in pigs and rabbits.

BACKGROUND: Some genetic studies need to take into account correlations between traits that are repeatedly measured over time. Multiple-trait random regression models are commonly used to analyze repeated traits but suffer from several major drawbacks. In the present study, we developed a multiple-trait extension of the structured antedependence model (SAD) to overcome this issue and validated its usefulness by modeling the association between litter size (LS) and average birth weight (ABW) over parities in pigs and rabbits. METHODS: The single-trait SAD model assumes that a random effect at time [Formula: see text] can be explained by the previous values of the random effect (i.e. at previous times). The proposed multiple-trait extension of the SAD model consists in adding a cross-antedependence parameter to the single-trait SAD model. This model can be easily fitted using ASReml and the OWN Fortran program that we have developed. In comparison with the random regression model, we used our multiple-trait SAD model to analyze the LS and ABW of 4345 litters from 1817 Large White sows and 8706 litters from 2286 L-1777 does over a maximum of five successive parities. RESULTS: For both species, the multiple-trait SAD fitted the data better than the random regression model. The difference between AIC of the two models (AIC_random regression-AIC_SAD) were equal to 7 and 227 for pigs and rabbits, respectively. A similar pattern of heritability and correlation estimates was obtained for both species. Heritabilities were lower for LS (ranging from 0.09 to 0.29) than for ABW (ranging from 0.23 to 0.39). The general trend was a decrease of the genetic correlation for a given trait between more distant parities. Estimates of genetic correlations between LS and ABW were negative and ranged from -0.03 to -0.52 across parities. No correlation was observed between the permanent environmental effects, except between the permanent environmental effects of LS and ABW of the same parity, for which the estimate of the correlation was strongly negative (ranging from -0.57 to -0.67). CONCLUSIONS: We demonstrated that application of our multiple-trait SAD model is feasible for studying several traits with repeated measurements and showed that it provided a better fit to the data than the random regression model.

Dwarfism and Altered Craniofacial Development in Rabbits Is Caused by a 12.1 kb Deletion at the HMGA2 Locus.

The dwarf phenotype characterizes the smallest of rabbit breeds and is governed largely by the effects of a single dwarfing allele with an incompletely dominant effect on growth. Dwarf rabbits typically weigh under 1 kg and have altered craniofacial morphology. The dwarf allele is recessive lethal and dwarf homozygotes die within a few days of birth. The dwarf phenotype is expressed in heterozygous individuals and rabbits from dwarf breeds homozygous for the wild-type allele are normal, although smaller when compared to other breeds. Here, we show that the dwarf allele constitutes a ∼12.1 kb deletion overlapping the promoter region and first three exons of the HMGA2 gene leading to inactivation of this gene. HMGA2 has been frequently associated with variation in body size across species. Homozygotes for null alleles are viable in mice but not in rabbits and probably not in humans. RNA-sequencing analysis of rabbit embryos showed that very few genes (4-29 genes) were differentially expressed among the three HMGA2/dwarf genotypes, suggesting that dwarfism and inviability in rabbits are caused by modest changes in gene expression. Our results show that HMGA2 is critical for normal expression of IGF2BP2, which encodes an RNA-binding protein. Finally, we report a catalog of regions of elevated genetic differentiation between dwarf and normal-size rabbits, including LCORL-NCAPG, STC2, HOXD cluster, and IGF2BP2 Levels and patterns of genetic diversity at the LCORL-NCAPG locus further suggest that small size in dwarf breeds was enhanced by crosses with wild rabbits. Overall, our results imply that small size in dwarf rabbits results from a large effect, loss-of-function (LOF) mutation in HMGA2 combined with polygenic selection.

Levels and Patterns of Genetic Diversity and Population Structure in Domestic Rabbits.

Over thousands of years humans changed the genetic and phenotypic composition of several organisms and in the process transformed wild species into domesticated forms. From this close association, domestic animals emerged as important models in biomedical and fundamental research, in addition to their intrinsic economical and cultural value. The domestic rabbit is no exception but few studies have investigated the impact of domestication on its genetic variability. In order to study patterns of genetic structure in domestic rabbits and to quantify the genetic diversity lost with the domestication process, we genotyped 45 microsatellites for 471 individuals belonging to 16 breeds and 13 wild localities. We found that both the initial domestication and the subsequent process of breed formation, when averaged across breeds, culminated in losses of ~20% of genetic diversity present in the ancestral wild population and domestic rabbits as a whole, respectively. Despite the short time elapsed since breed diversification we uncovered a well-defined structure in domestic rabbits where the FST between breeds was 22%. However, we failed to detect deeper levels of structure, probably consequence of a recent and single geographic origin of domestication together with a non-bifurcating process of breed formation, which were often derived from crosses between two or more breeds. Finally, we found evidence for intrabreed stratification that is associated with demographic and selective causes such as formation of strains, colour morphs within the same breed, or country/breeder of origin. These additional layers of population structure within breeds should be taken into account in future mapping studies.

Rabbit genome analysis reveals a polygenic basis for phenotypic change during domestication.

The genetic changes underlying the initial steps of animal domestication are still poorly understood. We generated a high-quality reference genome for the rabbit and compared it to resequencing data from populations of wild and domestic rabbits. We identified more than 100 selective sweeps specific to domestic rabbits but only a relatively small number of fixed (or nearly fixed) single-nucleotide polymorphisms (SNPs) for derived alleles. SNPs with marked allele frequency differences between wild and domestic rabbits were enriched for conserved noncoding sites. Enrichment analyses suggest that genes affecting brain and neuronal development have often been targeted during domestication. We propose that because of a truly complex genetic background, tame behavior in rabbits and other domestic animals evolved by shifts in allele frequencies at many loci, rather than by critical changes at only a few domestication loci.

The genetic structure of domestic rabbits.

Understanding the genetic structure of domestic species provides a window into the process of domestication and motivates the design of studies aimed at making links between genotype and phenotype. Rabbits exhibit exceptional phenotypic diversity, are of great commercial value, and serve as important animal models in biomedical research. Here, we provide the first comprehensive survey of nucleotide polymorphism and linkage disequilibrium (LD) within and among rabbit breeds. We resequenced 16 genomic regions in population samples of both wild and domestic rabbits and additional 35 fragments in 150 rabbits representing six commonly used breeds. Patterns of genetic variation suggest a single origin of domestication in wild populations from France, supporting historical records that place rabbit domestication in French monasteries. Levels of nucleotide diversity both within and among breeds were ~0.2%, but only 60% of the diversity present in wild populations from France was captured by domestic rabbits. Despite the recent origin of most breeds, levels of population differentiation were high (F(ST) = 17.9%), but the majority of polymorphisms were shared and thus transferable among breeds. Coalescent simulations suggest that domestication began with a small founding population of less than 1,200 individuals. Taking into account the complex demographic history of domestication with two successive bottlenecks, two loci showed deviations that were consistent with artificial selection, including GPC4, which is known to be associated with growth rates in humans. Levels of diversity were not significantly different between autosomal and X-linked loci, providing no evidence for differential contributions of males and females to the domesticated gene pool. The structure of LD differed substantially within and among breeds. Within breeds, LD extends over large genomic distances. Markers separated by 400 kb typically showed r(2) higher than 0.2, and some LD extended up to 3,200 kb. Much less LD was found among breeds. This advantageous LD structure holds great promise for reducing the interval of association in future mapping studies.

Role of Gal11, a component of the RNA polymerase II mediator in stress-induced hyperphosphorylation of Msn2 in Saccharomyces cerevisiae.

In the yeast Saccharomyces cerevisiae, the Msn2 transcription factor is a key element in mediating the environmental stress response (ESR), leading to the induction of 100-200 genes through the cis-acting Stress Response Element (STRE) in response to various physico-chemical stresses and nutritional variations. This activation is accompanied by a stress-induced hyperphosphorylation of Msn2. By a systematic screening we identified two proteins essential in this process: (i) the cyclin-dependent Ssn3/Srb10 protein kinase, part of a module of the RNA polymerase II mediator, which has already been shown to be involved in hyperphosphorylation and degradation of Msn2 upon stress, and (ii) Gal11, a component of the mediator. In a gal11 mutant, stress-induced hyperphosphorylation of Msn2 is abolished, stress-induced transcription of Msn2-dependent genes is decreased and Msn2 degradation is impaired. Rgr1, another component of the mediator, is also critical for this hyperphosphorylation, indicating that the integrity of the mediator is required for this process. Moreover the transactivating region of Msn2 interacts in vitro with the N-terminal domain of Gal11. These results point out the role of the mediator, especially its Gal11 subunit, in the hyperphosphorylation and degradation of Msn2 during stress response.

The transcriptional activation region of Msn2p, in Saccharomyces cerevisiae, is regulated by stress but is insensitive to the cAMP signalling pathway.

Msn2p is a transcription factor that mediates a transient cellular response to multiple stresses and to changes in the nutritional environment. It was previously shown that the C-terminal half of Msn2p contains the DNA binding domain, a nuclear localization signal and nuclear export determinants which are activated by stress. In this report, we demonstrate that the N-terminal half of Msn2p contains the transcriptional activation domain(s). In addition, we present evidence that this region of Msn2p is able to mediate both the activation of transcription and export of the protein from the nucleus in response to stress. Interestingly, while the stress response integrated by the components of the C-terminal half that are involved in nucleocytoplasmic localization is reversed by elevated levels of cAMP, the effects of stress on the transcriptional activation domain and the localization determinants present in the N-terminal half of Msn2p are insensitive to variations in the intracellular cAMP concentrations.

Stress induces depletion of Cdc25p and decreases the cAMP producing capability in Saccharomyces cerevisiae.

In Saccharomyces cerevisiae the cAMP-dependent protein kinase A pathway antagonizes the cellular response to stress. It is shown here that the cellular content of Cdc25p, the upstream activator of Ras and adenylyl cyclase, decays upon various stresses such as heat shock and oxidative and ethanol shocks, whereas its phosphorylation level and its localization are unaffected. In parallel with the reduction of Cdc25p, the maximal capacity of the cell to accumulate cAMP decreases when its feedback regulation is abolished. A deletion of CDC25 prevents this decrease. Paradoxically, in wild-type cells, with normal feedback regulation, the level of cAMP, which is much lower, is not reduced but is rather increased upon stress. These observations are consistent with a role of Cdc25p in sensing and transducing stress to downstream targets, either through a cAMP-independent pathway or by large fluctuations in the cAMP content of the cell.

Hyperphosphorylation of Msn2p and Msn4p in response to heat shock and the diauxic shift is inhibited by cAMP in Saccharomyces cerevisiae.

In response to various stresses, as well as during the diauxic transition, the Msn2p and Msn4p transcription factors of Saccharomyces cerevisiae are activated and induce a large set of genes. This activation is inhibited by the Ras/cAMP/PKA (cAMP-dependent protein kinase) pathway. Here we show by immunoblotting experiments that Msn2p and Msn4p are phosphorylated in vivo during growth on glucose, and become hyperphosphorylated at the diauxic transition and upon heat shock. This hyperphosphorylation is correlated with activation of Msn2/4p-dependent transcription. An increased level of cAMP prevents and reverses these hyperphosphorylations, indicating that kinases other than PKA are involved. These results suggest that PKA and stress-activated kinases control Msn2/4p activity by antagonistic phosphorylation. It was also noted that Msn4p is transiently increased at the diauxic transition. Msn2p and Msn4p present different hyperphosphorylation patterns in response to different stresses.