Genome-wide association studies for diarrhoea outcomes identified genomic regions affecting resistance to a severe enteropathy in suckling rabbits.

Selection and breeding strategies to improve resistance to enteropathies are essential to reaching the sustainability of the rabbit production systems. However, disease heterogeneity (having only as major visible symptom diarrhoea) and low disease heritability are two barriers for the implementation of these strategies. Diarrhoea condition can affect rabbits at different life stages, starting from the suckling period, with large negative economic impacts. In this study, from a commercial population of suckling rabbits (derived from 133 litters) that experienced an outbreak of enteropathy, we first selected a few animals that died with severe symptoms of diarrhoea and characterized their microbiota, using 16S rRNA gene sequencing data. Clostridium genus was consistently present in all affected specimens. In addition, with the aim to identify genetic markers in the rabbit genome that could be used as selection tools, we performed genome-wide association studies for symptoms of diarrhoea in the same commercial rabbit population. These studies were also complemented with FST analyses between the same groups of rabbits. A total of 332 suckling rabbits (151 with severe symptoms of diarrhoea, 42 with mild symptoms and 129 without any symptoms till the weaning period), derived from 45 different litters (a subset of the 133 litters) were genotyped with the Affymetrix Axiom OrcunSNP Array. In both genomic approaches, rabbits within litters were paired to constitute two groups (susceptible and resistant, including the mildly affected in one or the other group) and run case and control genome-wide association analyses. Genomic heritability estimated in the designed experimental structure integrated in a commercial breeding scheme was 0.19-0.21 (s.e. 0.09-0.10). A total of eight genomic regions on rabbit chromosome 2 (OCU2), OCU3, OCU7, OCU12, OCU13, OCU16 and in an unassembled scaffold had significant single nucleotide polymorphisms (SNPs) and/or markers that trespassed the FST percentile distribution. Among these regions, three main peaks of SNPs were identified on OCU12, OCU13 and OCU16. The QTL region on OCU13 encompasses several genes that encode members of a family of immunoglobulin Fc receptors (FCER1G, FCRLA, FCRLB and FCGR2A) involved in the immune innate system, which might be important candidate genes for this pathogenic condition. The results obtained in this study demonstrated that resistance to an enteropathy occurring in suckling rabbits is in part genetically determined and can be dissected at the genomic level, providing DNA markers that could be used in breeding programmes to increase resistance to enteropathies in meat rabbits.

Comparative analysis of genomic inbreeding parameters and runs of homozygosity islands in several fancy and meat rabbit breeds.

Runs of homozygosity (ROH) are defined as long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROH and their length are indicators of the level and origin of inbreeding. In this study, we analysed SNP chip datasets (obtained using the Axiom OrcunSNP Array) of a total of 702 rabbits from 12 fancy breeds and four meat breeds to identify ROH with different approaches and calculate several genomic inbreeding parameters. The highest average number of ROH per animal was detected in Belgian Hare (~150) and the lowest in Italian Silver (~106). The average length of ROH ranged from 4.001 ± 0.556 Mb in Italian White to 6.268 ± 1.355 Mb in Ermine. The same two breeds had the lowest (427.9 ± 86.4 Mb, Italian White) and the highest (921.3 ± 179.8 Mb, Ermine) average values of the sum of all ROH segments. More fancy breeds had a higher level of genomic inbreeding (as defined by ROH) than meat breeds. Several ROH islands contain genes involved in body size, body length, pigmentation processes, carcass traits, growth, and reproduction traits (e.g.: AOX1, GPX5, IFRD1, ITGB8, NELL1, NR3C1, OCA2, TRIB1, TRIB2). Genomic inbreeding parameters can be useful to overcome the lack of information in the management of rabbit genetic resources. ROH provided information to understand, to some extent, the genetic history of rabbit breeds and to identify signatures of selection in the rabbit genome.

Identification of Polymorphisms in the Rabbit Growth Hormone Receptor (GHR) Gene and Association with Finishing Weight in a Commercial Meat Rabbit Line.

A shortcut to identify DNA markers associated with economic traits is to use a candidate gene approach that is still useful in livestock species in which molecular tools and resources are not advanced or not well developed. Mutations in the growth hormone receptor (GHR) gene associated with production traits have been already described in several livestock species. For this reason GHR could be an interesting candidate gene in the rabbit. In this study we re-sequenced all exons and non-coding regions of the rabbit GHR gene in a panel of 10 different rabbits and identified 10 single nucleotide polymorphisms (SNPs). One of them (g.63453192C>G or c.106C>G), located in exon 3 was a missense mutation (p.L36V) substituting an amino acid in a highly conserved position across all mammals. This mutation was genotyped in 297 performance tested rabbits of a meat male line and association analysis showed that the investigated SNP was associated with weight at 70 days (P < 0.05). The most frequent genotype (GG) was in animals with higher weight at this age, suggesting that the high directional selection pressure toward this trait since the constitution of the genotyped line might have contributed to shape allele frequencies at this polymorphic site.

A missense mutation in the rabbit melanocortin 4 receptor (MC4R) gene is associated with finishing weight in a meat rabbit line.

In this study we resequenced 1729 bp of the rabbit melanocortin 4 receptor (MC4 R) gene in 31 rabbits from different breeds/lines and identified ten polymorphisms: one was an indel and 9 were single nucleotide polymorphisms (SNPs). The indel and 5 SNPs were in the 5'-flanking region, 3 were synonymous SNPs and one was a missense mutation (c.101G>A; p.G34D), located in a conserved position of the extracellular tail of the MC4 R protein. The missense mutation was analyzed in a panel of 74 rabbits of different breeds and in 516 performance tested rabbits of a commercial paternal line under selection for growth efficiency. Association analysis indicated that rabbits with the less frequent genotype in this population (DD) had a lighter weight at 70 postnatal days than animals with genotype GD (P < 0.10) and animals with genotype GG (P < 0.05). This is the third study on candidate genes, after those on GH1 and IGF2 that reported a marker associated with finishing weight. Therefore, it seems that a candidate gene approach in rabbit based on previous information accumulated in other livestock species could be useful to identify genes explaining a fraction of variability of performance traits with potential application on rabbit breeding and selection.