Association of single nucleotide polymorphisms in Neuropeptide Y (NPY) and Phosphoglycerate Mutase 2 (PGAM2) genes with growth traits in rabbits.

Genetic improvement of local rabbit breeds using modern approaches such as marker-assisted selection requires accurate and precise information about marker‒trait associations in animals with different genetic backgrounds. Therefore, this study was designed to estimate the association between two mutations located in the Neuropeptide Y (NPY, g.1778G > C) and Phosphoglycerate Mutase 2 (PGAM2, c.195 C > T) genes in New Zealand White (NZW), Baladi (BR), and V-line rabbits. The first mutation was genotyped using high-resolution melting, and the second mutation was genotyped using the PCR-RFLP method. The results revealed significant associations between the NPY mutation and body weight at 10 (V-line) and 12 weeks of age (NZW, BR, and V-line), body weight gain (BWG) from 10 to 12 weeks of age (BR), BWG from 6 to 12 weeks of age (NZW, BR, and V-line), average daily gain (NZW, BR, and V-line, and BR), growth rate (GR) from 8 to10 weeks (V-line), 10 to 12 weeks (BR), and GR from 6 to 12 weeks of age (BR, and V-line). The PGAM2 mutation was associated with body weight at 10 (V-line) and 12 (NZW, and V-line) weeks of age, with significant positive additive effects at 12 weeks of age in all breeds, and was associated with BWG from 8 to 10 and 10 to 12 in BR, and BWG from 6 to 12 weeks of age (NZW, and BR), and average daily gain (NZW, and BR), and was associated with GR form 8 to 10 weeks (BR), from10 to 12 weeks (BR, and V-line) and from 6 to 12 weeks (BR). The results highlighted the importance of the two mutations in growth development, and the possibility of considering them as candidate genes for late growth in rabbits.

Selection against domestication alleles in introduced rabbit populations.

Humans have moved domestic animals around the globe for thousands of years. These have occasionally established feral populations in nature, often with devastating ecological consequences. To understand how natural selection shapes re-adaptation into the wild, we investigated one of the most successful colonizers in history, the European rabbit. By sequencing the genomes of 297 rabbits across three continents, we show that introduced populations exhibit a mixed wild-domestic ancestry. We show that alleles that increased in frequency during domestication were preferentially selected against in novel natural environments. Interestingly, causative mutations for common domestication traits sometimes segregate at considerable frequencies if associated with less drastic phenotypes (for example, coat colour dilution), whereas mutations that are probably strongly maladaptive in nature are absent. Whereas natural selection largely targeted different genomic regions in each introduced population, some of the strongest signals of parallelism overlap genes associated with neuronal or brain function. This limited parallelism is probably explained by extensive standing genetic variation resulting from domestication together with the complex mixed ancestry of introduced populations. Our findings shed light on the selective and molecular mechanisms that enable domestic animals to re-adapt to the wild and provide important insights for the mitigation and management of invasive populations.

Intragenic MFSD8 duplication and histopathological findings in a rabbit with neuronal ceroid lipofuscinosis.

Neuronal ceroid lipofuscinoses (NCL) are among the most prevalent neurodegenerative disorders of early life in humans. Disease-causing variants have been described for 13 different NCL genes. In this study, a refined pathological characterization of a female rabbit with progressive neurological signs reminiscent of NCL was performed. Cytoplasmic pigment present in neurons was weakly positive with Sudan black B and autofluorescent. Immunohistology revealed astrogliosis, microgliosis and axonal degeneration. During the subsequent genetic investigation, the genome of the affected rabbit was sequenced and examined for private variants in NCL candidate genes. The analysis revealed a homozygous ~10.7 kb genomic duplication on chromosome 15 comprising parts of the MFSD8 gene, NC_013683.1:g.103,727,963_103,738,667dup. The duplication harbors two internal protein coding exons and is predicted to introduce a premature stop codon into the transcript, truncating ~50% of the wild-type MFSD8 open reading frame encoding the major facilitator superfamily domain containing protein 8, XP_002717309.2:p.(Glu235Leufs*23). Biallelic loss-of-function variants in MFSD8 have been described to cause NCL7 in human patients, dogs and a single cat. The available clinical and pathological data, together with current knowledge about MFSD8 variants and their functional impact in other species, point to the MFSD8 duplication as a likely causative defect for the observed phenotype in the affected rabbit.

De novo transcriptome assembly database for 100 tissues from each of seven species of domestic herbivore.

Domesticated herbivores are an important agricultural resource that play a critical role in global food security, particularly as they can adapt to varied environments, including marginal lands. An understanding of the molecular basis of their biology would contribute to better management and sustainable production. Thus, we conducted transcriptome sequencing of 100 to 105 tissues from two females of each of seven species of herbivore (cattle, sheep, goats, sika deer, horses, donkeys, and rabbits) including two breeds of sheep. The quality of raw and trimmed reads was assessed in terms of base quality, GC content, duplication sequence rate, overrepresented k-mers, and quality score distribution with FastQC. The high-quality filtered RNA-seq raw reads were deposited in a public database which provides approximately 54 billion high-quality paired-end sequencing reads in total, with an average mapping rate of ~93.92%. Transcriptome databases represent valuable resources that can be used to study patterns of gene expression, and pathways that are related to key biological processes, including important economic traits in herbivores.

Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement.

Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.

Marker Density and Models to Improve the Accuracy of Genomic Selection for Growth and Slaughter Traits in Meat Rabbits.

The selection and breeding of good meat rabbit breeds are fundamental to their industrial development, and genomic selection (GS) can employ genomic information to make up for the shortcomings of traditional phenotype-based breeding methods. For the practical implementation of GS in meat rabbit breeding, it is necessary to assess different marker densities and GS models. Here, we obtained low-coverage whole-genome sequencing (lcWGS) data from 1515 meat rabbits (including parent herd and half-sibling offspring). The specific objectives were (1) to derive a baseline for heritability estimates and genomic predictions based on randomly selected marker densities and (2) to assess the accuracy of genomic predictions for single- and multiple-trait linear mixed models. We found that a marker density of 50 K can be used as a baseline for heritability estimation and genomic prediction. For GS, the multi-trait genomic best linear unbiased prediction (GBLUP) model results in more accurate predictions for virtually all traits compared to the single-trait model, with improvements greater than 15% for all of them, which may be attributed to the use of information on genetically related traits. In addition, we discovered a positive correlation between the performance of the multi-trait GBLUP and the genetic correlation between the traits. We anticipate that this approach will provide solutions for GS, as well as optimize breeding programs, in meat rabbits.

Cost-effectively dissecting the genetic architecture of complex wool traits in rabbits by low-coverage sequencing.

BACKGROUND: Rabbit wool traits are important in fiber production and for model organism research on hair growth, but their genetic architecture remains obscure. In this study, we focused on wool characteristics in Angora rabbits, a breed well-known for the quality of its wool. Considering the cost to generate population-scale sequence data and the biased detection of variants using chip data, developing an effective genotyping strategy using low-coverage whole-genome sequencing (LCS) data is necessary to conduct genetic analyses. RESULTS: Different genotype imputation strategies (BaseVar + STITCH, Bcftools + Beagle4, and GATK + Beagle5), sequencing coverages (0.1X, 0.5X, 1.0X, 1.5X, and 2.0X), and sample sizes (100, 200, 300, 400, 500, and 600) were compared. Our results showed that using BaseVar + STITCH at a sequencing depth of 1.0X with a sample size larger than 300 resulted in the highest genotyping accuracy, with a genotype concordance higher than 98.8% and genotype accuracy higher than 0.97. We performed multivariate genome-wide association studies (GWAS), followed by conditional GWAS and estimation of the confidence intervals of quantitative trait loci (QTL) to investigate the genetic architecture of wool traits. Six QTL were detected, which explained 0.4 to 7.5% of the phenotypic variation. Gene-level mapping identified the fibroblast growth factor 10 (FGF10) gene as associated with fiber growth and diameter, which agrees with previous results from functional data analyses on the FGF gene family in other species, and is relevant for wool rabbit breeding. CONCLUSIONS: We suggest that LCS followed by imputation can be a cost-effective alternative to array and high-depth sequencing for assessing common variants. GWAS combined with LCS can identify new QTL and candidate genes that are associated with quantitative traits. This study provides a cost-effective and powerful method for investigating the genetic architecture of complex traits, which will be useful for genomic breeding applications.

A single introduction of wild rabbits triggered the biological invasion of Australia.

Biological invasions are a major cause of environmental and economic disruption. While ecological factors are key determinants of their success, the role of genetics has been more challenging to demonstrate. The colonization of Australia by the European rabbit is one of the most iconic and devastating biological invasions in recorded history. Here, we show that despite numerous introductions over a 70-y period, this invasion was triggered by a single release of a few animals that spread thousands of kilometers across the continent. We found genetic support for historical accounts that these were English rabbits imported in 1859 by a settler named Thomas Austin and traced the origin of the invasive population back to his birthplace in England. We also find evidence of additional introductions that established local populations but have not spread geographically. Combining genomic and historical data we show that, contrary to the earlier introductions, which consisted mostly of domestic animals, the invasive rabbits had wild ancestry. In New Zealand and Tasmania, rabbits also became a pest several decades after being introduced. We argue that the common denominator of these invasions was the arrival of a new genotype that was better adapted to the natural environment. These findings demonstrate how the genetic composition of invasive individuals can determine the success of an introduction and provide a mechanism by which multiple introductions can be required for a biological invasion.

The complete mitochondrial genome of the 'Zacatuche' Volcano rabbit (Romerolagus diazi), an endemic and endangered species from the Volcanic Belt of Central Mexico.

BACKGROUND: The 'Zacatuche', 'Teporingo', or Volcano rabbit (Romerolagus diazi) belongs to the family Leporidae, is an endemic species restricted to the Central part of the Trans-Mexican Volcanic Belt, and is considered 'endangered' by the IUCN Red List of Threatened Species. METHODS AND RESULTS: This study reports, for the first time, the complete mitochondrial genome of R. diazi and examined the phylogenetic position of R. diazi among other closely related co-familiar species using mitochondrial protein-coding genes (PCGs). The mitogenome of R. diazi was assembled from short Illumina 150 bp pair-end reads with a coverage of 189x. The AT-rich mitochondrial genome of R. diazi is 17,400 bp in length and is comprised of 13 PCGs, two ribosomal RNA genes, and 22 transfer RNA genes. The gene order observed in the mitochondrial genome of R. diazi is identical to that reported for other leporids. Phylogenetic analyses based on PCGs support the basal position of Romerolagus within the Leporidae, at least when compared to the genera Oryctolagus and Lepus. Nonetheless, additional mitochondrial genomes from species belonging to the genera Bunolagus, Sylvilagus, and Pronolagus, among others, are needed before a more robust conclusion about the derived vs basal placement of Romerolagus within the family Leporidae can be reached based on mitochondrial PCGs. CONCLUSIONS: This is the first genomic resource developed for R. diazi and it represents a tool to improve our understanding about the ecology and evolutionary biology of this iconic and endangered species.

Selection for environmental variance of litter size in rabbits involves genes in pathways controlling animal resilience.

BACKGROUND: Environmental variance (VE) is partially under genetic control, which means that the VE of individuals that share the same environment can differ because they have different genotypes. Previously, a divergent selection experiment for VE of litter size (LS) during 13 generations in rabbit yielded a successful response and revealed differences in resilience between the divergent lines. The aim of the current study was to identify signatures of selection in these divergent lines to better understand the molecular mechanisms and pathways that control VE of LS and animal resilience. Three methods (FST, ROH and varLD) were used to identify signatures of selection in a set of 473 genotypes from these rabbit lines (377) and a base population (96). A whole-genome sequencing (WGS) analysis was performed on 54 animals to detect genes with functional mutations. RESULTS: By combining signatures of selection and WGS data, we detected 373 genes with functional mutations in their transcription units, among which 111 had functions related to the immune system, stress response, reproduction and embryo development, and/or carbohydrate and lipid metabolism. The genes TTC23L, FBXL20, GHDC, ENSOCUG00000031631, SLC18A1, CD300LG, MC2R, and ENSOCUG00000006264 were particularly relevant, since each one carried a functional mutation that was fixed in one of the rabbit lines and absent in the other line. In the 3'UTR region of the MC2R and ENSOCUG00000006264 genes, we detected a novel insertion/deletion (INDEL) variant. CONCLUSIONS: Our findings provide further evidence in favour of VE as a measure of animal resilience. Signatures of selection were identified for VE of LS in genes that have a functional mutation in their transcription units and are mostly implicated in the immune response and stress response pathways. However, the real implications of these genes for VE and animal resilience will need to be assessed through functional analyses.

A genome-wide association study for the number of teats in European rabbits (Oryctolagus cuniculus) identifies several candidate genes affecting this trait.

In the European rabbit (Oryctolagus cuniculus), a polytocous livestock species, the number of teats indirectly impacts the doe reproduction efficiency and, in turn, the sustainable production of rabbit meat. In this study, we carried out a genome-wide association study (GWAS) for the total number of teats in 247 Italian White does included in the Italian White rabbit breed selection program, by applying a selective genotyping approach. Does had either 8 (n = 121) or 10 teats (n = 126). All rabbits were genotyped with the Affymetrix Axiom OrcunSNP Array. Genomic data from the two extreme groups of rabbits were also analysed with the single-marker fixation index statistic and combined with the GWAS results. The GWAS identified 50 significant SNPs and the fixation index analysis identified a total of 20 SNPs that trespassed the 99.98th percentile threshold, 19 of which confirmed the GWAS results. The most significant SNP (P = 4.31 × 10-11 ) was located on OCU1, close to the NUDT2 gene, a breast carcinoma cells proliferation promoter. Another significant SNP identified as candidate gene NR6A1, which is well known to play an important role in affecting the correlated number of vertebrae in pigs. Other significant markers were close to candidate genes involved in determining body length in mice. Markers associated with increased number of teats could be included in selection programmes to speed up the improvement for this trait in rabbit lines that need to increase maternal performances.

Effect of host breeds on gut microbiome and serum metabolome in meat rabbits.

BACKGROUND: Gut microbial compositional and functional variation can affect health and production performance of farm animals. Analysing metabolites in biological samples provides information on the basic mechanisms that affect the well-being and production traits in farm animals. However, the extent to which host breeds affect the gut microbiome and serum metabolome in meat rabbits is still unknown. In this study, the differences in phylogenetic composition and functional capacities of gut microbiota in two commercial rabbit breeds Elco and Ira were determined by 16S rRNA gene and metagenomic sequencing. The alternations in serum metabolome in the two rabbit breeds were detected using ultra-performance liquid chromatography system coupled with quadrupole time of flight mass spectrometry (UPLC-QTOFMS). RESULTS: Sequencing results revealed that there were significant differences in the gut microbiota of the two breeds studied, suggesting that host breeds affect structure and diversity of gut microbiota. Numerous breed-associated microorganisms were identified at different taxonomic levels and most microbial taxa belonged to the families Lachnospiraceae and Ruminococcaceae. In particular, several short-chain fatty acids (SCFAs) producing species including Coprococcus comes, Ruminococcus faecis, Ruminococcus callidus, and Lachnospiraceae bacterium NK4A136 could be considered as biomarkers for improving the health and production performance in meat rabbits. Additionally, gut microbial functional capacities related to bacterial chemotaxis, ABC transporters, and metabolism of different carbohydrates, amino acids, and lipids varied greatly between rabbit breeds. Several fatty acids, amino acids, and organic acids in the serum were identified as breed-associated, where certain metabolites could be regarded as biomarkers correlated with the well-being and production traits of meat rabbits. Correlation analysis between breed-associated microbial species and serum metabolites revealed significant co-variations, indicating the existence of cross-talk among host-gut microbiome-serum metabolome. CONCLUSIONS: Our study provides insight into how gut microbiome and serum metabolome of meat rabbits are affected by host breeds and uncovers potential biomarkers important for breed improvement of meat rabbits.

Rabbits - their domestication and molecular genetics of hair coat development and quality.

The European rabbit (Oryctolagus cuniculus) is the only representative of its genus living in present-day Europe and North Africa, and all domestic rabbits are descendants of this one species, which is native to the Iberian Peninsula. There are over 300 breeds of rabbits that differ in size, coat color, length of ears and type of fur. Rabbits are bred for various reasons, such as for laboratory animals and a source of meat, wool and fur, as well as for pets and exhibition animals. The hair coat is a important economic trait of rabbits. Its development and quality are influenced by various factors, both environmental and genetic. The genetic mechanisms underlying its development have not been thoroughly researched. The aim of this review is to discuss the domestication of rabbits and the different aspects of rabbit genetics. A brief review of the properties of rabbit hair coat, hair coat development and hair cycle will be provided, followed by discussion of the factors regulating hair coat development, molecular control of hair coat development and the role of non-coding RNAs in the regulation of gene expression in the hair follicles of rabbits. Information about genetic regulation of pathways could provide useful tools for improving hair coat quality and be of practical use in rabbit breeding.

Brain Transcriptomics of Wild and Domestic Rabbits Suggests That Changes in Dopamine Signaling and Ciliary Function Contributed to Evolution of Tameness.

Domestication has resulted in immense phenotypic changes in animals despite their relatively short evolutionary history. The European rabbit is one of the most recently domesticated animals, but exhibits distinct morphological, physiological, and behavioral differences from their wild conspecifics. A previous study revealed that sequence variants with striking allele frequency differences between wild and domestic rabbits were enriched in conserved noncoding regions, in the vicinity of genes involved in nervous system development. This suggests that a large proportion of the genetic changes targeted by selection during domestication might affect gene regulation. Here, we generated RNA-sequencing data for four brain regions (amygdala, hypothalamus, hippocampus, and parietal/temporal cortex) sampled at birth and revealed hundreds of differentially expressed genes (DEGs) between wild and domestic rabbits. DEGs in amygdala were significantly enriched for genes associated with dopaminergic function and all 12 DEGs in this category showed higher expression in domestic rabbits. DEGs in hippocampus were enriched for genes associated with ciliary function, all 21 genes in this category showed lower expression in domestic rabbits. These results indicate an important role of dopamine signaling and ciliary function in the evolution of tameness during rabbit domestication. Our study shows that gene expression in specific pathways has been profoundly altered during domestication, but that the majority of genes showing differential expression in this study have not been the direct targets of selection.

Genome-wide association study for feed efficiency in collective cage-raised rabbits under full and restricted feeding.

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two-trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non-genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome-wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits.

Inbreeding depression for kit survival at birth in a rabbit population under long-term selection.

BACKGROUND: Accumulation of detrimental mutations in small populations leads to inbreeding depression of fitness traits and a higher frequency of genetic defects, thus increasing risk of extinction. Our objective was to quantify the magnitude of inbreeding depression for survival at birth, in a closed rabbit population under long-term selection. METHODS: We used an information theory-based approach and multi-model inference to estimate inbreeding depression and its purging with respect to the trait 'kit survival at birth' over a 25-year period in a closed population of Pannon White rabbits, by analysing 22,718 kindling records. Generalised linear mixed models based on the logit link function were applied, which take polygenic random effects into account. RESULTS: Our results indicated that inbreeding depression occurred during the period 1992-1997, based on significant estimates of the z-standardised classical inbreeding coefficient z.FL (CI95% - 0.12 to - 0.03) and of the new inbreeding coefficient of the litter z.FNEWL (CI95% - 0.13 to - 0.04) as well as a 59.2% reduction in contributing founders. Inbreeding depression disappeared during the periods 1997-2007 and 2007-2017. For the period 1992-1997, the best model resulted in a significantly negative standardised estimate of the new inbreeding coefficient of the litter and a significantly positive standardised estimate of Kalinowski's ancestral inbreeding coefficient of the litter (CI95% 0.01 to 0.17), which indicated purging of detrimental load. Kindling season and parity had effects on survival at birth that differed across the three periods of 1992-1997, 1997-2007 and 2007-2017. CONCLUSIONS: Our results support the existence of inbreeding depression and its purging with respect to kit survival at birth in this Pannon White rabbit population. However, we were unable to exclude possible confounding from the effects of parity and potentially other environmental factors during the study period, thus our results need to be extended and confirmed in other populations.

Identification of functional mutations associated with environmental variance of litter size in rabbits.

BACKGROUND: Environmental variance (VE) is partly under genetic control and has recently been proposed as a measure of resilience. Unravelling the genetic background of the VE of complex traits could help to improve resilience of livestock and stabilize their production across farming systems. The objective of this study was to identify genes and functional mutations associated with variation in VE of litter size (LS) in rabbits. To achieve this, we combined the results of a genome-wide association study (GWAS) and a whole-genome sequencing (WGS) analysis using data from two divergently selected rabbit lines for high and low VE of LS. These lines differ in terms of biomarkers of immune response and mortality. Moreover, rabbits with a lower VE of LS were found to be more resilient to infections than animals with a higher VE of LS. RESULTS: By using two GWAS approaches (single-marker regression and Bayesian multiple-marker regression), we identified four genomic regions associated with VE of LS, on chromosomes 3, 7, 10, and 14. We detected 38 genes in the associated genomic regions and, using WGS, we identified 129 variants in the splicing, UTR, and coding (missense and frameshift effects) regions of 16 of these 38 genes. These genes were related to the immune system, the development of sensory structures, and stress responses. All of these variants (except one) segregated in one of the rabbit lines and were absent (n = 91) or fixed in the other one (n = 37). The fixed variants were in the HDAC9, ITGB8, MIS18A, ENSOCUG00000021276 and URB1 genes. We also identified a 1-bp deletion in the 3'UTR region of the HUNK gene that was fixed in the low VE line and absent in the high VE line. CONCLUSIONS: This is the first study that combines GWAS and WGS analyses to study the genetic basis of VE. The new candidate genes and functional mutations identified in this study suggest that the VE of LS is under the control of functions related to the immune system, stress response, and the nervous system. These findings could also explain differences in resilience between rabbits with homogeneous and heterogeneous VE of litter size.

Population genetics of the European rabbit along a rural-to-urban gradient.

The European rabbit (Oryctolagus cuniculus) is declining in large parts of Europe but populations in some German cities remained so far unaffected by this decline. The question arises of how urbanization affects patterns of population genetic variation and differentiation in German rabbit populations, as urban habitat fragmentation may result in altered meta-population dynamics. To address this question, we used microsatellite markers to genotype rabbit populations occurring along a rural-to-urban gradient in and around the city of Frankfurt, Germany. We found no effect of urbanization on allelic richness. However, the observed heterozygosity was significantly higher in urban than rural populations and also the inbreeding coefficients were lower, most likely reflecting the small population sizes and possibly on-going loss of genetic diversity in structurally impoverished rural areas. Global FST and G'ST-values suggest moderate but significant differentiation between populations. Multiple matrix regression with randomization ascribed this differentiation to isolation-by-environment rather than isolation-by-distance. Analyses of migration rates revealed asymmetrical gene flow, which was higher from rural into urban populations than vice versa and may again reflect intensified agricultural land-use practices in rural areas. We discuss that populations inhabiting urban areas will likely play an important role in the future distribution of European rabbits.

Genetic variability of functional longevity in five rabbit lines.

The objectives of this study were to analyse the differences in the genetic determination of functional longevity in five Spanish lines of rabbits and to check how different systematic factors might affect this genetic determination. Four of the lines were maternal (lines A, V, H and LP), these lines were established selecting base generation animals according to different criteria, but in the subsequent generations all of them were selected for litter size at weaning. The other is the paternal line R, this line was constituted by selecting animals with an outstanding daily growth rate. The trait analysed, length of productive life, was the time in days between the date of the first positive pregnancy test and the date of culling or death of a doe. Four models extended from the Cox proportional hazard model were used to analyse data of each line separately and jointly. The complete model (Model 1) included the fixed effect of year-season (YS) combination, positive palpation order (OPP), that is, reproductive cycle, physiological status of the doe (PS) at service and number of kits born alive (NBA) in each kindling as time-dependent factors. The inbreeding coefficient was fitted as a continuous covariate and the animal's additive genetic effect was also fitted to the model (Model 1). The other models were identical to Model 1 but excluding OPP (Model 2) or PS (Model 3) or NBA (Model 4), which were explored to assess the consequence on additive variance estimates of not correcting for these animal-dependent factors. Estimated effective heritabilities of longevity were 0.07 ± 0.03, 0.03 ± 0.02, 0.14 ± 0.09, 0.05 ± 0.04, 0.02 ± 0.01 and 0.04 ± 0.01 for lines A, V, H, LP, R and for the merged data set, respectively. Removing the PS from the model led to an increase in the estimated additive genetic variance in all lines (0.17 ± 0.05, 0.05 ± 0.03, 0.29 ± 0.19, 0.29 ± 0.20, 0.07 ± 0.04 and 0.05 ± 0.02 for lines A, V, H, LP, R and the merged data set, respectively). The highest hazard of death and/or culling was observed during the first two parities and decreased as the order of parity progressed. Does non-pregnant-non-lactating had the highest risk of death or culling. The does that had zero kits born alive incurred the highest risk, and this risk decreased as the NBA increased. In conclusion, the consideration of longevity as selection criterion for the studied rabbit lines is not recommended.

Whole transcriptome sequencing reveals that non-coding RNAs are related to embryo morphogenesis and development in rabbits.

The roles of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in embryonic development remain unclear. We performed a comprehensive analysis of lncRNA and circRNA profiles in rabbit embryos at different stages by whole transcriptome sequencing. We identified 719 lncRNAs and 744 circRNAs that were differentially expressed between stages S1, S2 and S3. A total of 241 differentially expressed lncRNAs and 166 differentially expressed circRNAs were significantly involved in embryonic morphogenesis and development. An RNA network was established and of the embryonic development-associated RNAs, the lncRNAs TCONS_00009253 and TCONS_00010436 were persistently downregulated, while circRNA_07129, circRNA_15209, and circRNA_12526 were persistently upregulated, and their co-expressed mRNAs TBX1, WNT3 and FGFR2 were persistently downregulated during embryonic development. These candidate RNAs were mainly involved in the Wnt, PI3K-Akt, and calcium signaling pathways. This study reports candidate lncRNAs and circRNAs that may be indispensable for the morphogenesis and development of rabbit embryos.