Construction of a cytogenetically anchored microsatellite map in rabbit

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.

A Genetic Evaluation System for New Zealand White Rabbit Germplasm Resources Based on SSR Markers

Simple Summary

The New Zealand white rabbit (Oryctolagus cuniculus) is one of the most important breeds of commercial and experimental rabbits in the world, and also one of the most raised rabbit breeds in China. Our goal was to develop a suite of microsatellite markers to aid future conservation genetics research for the Oryctolagus cuniculus breeds. Based on the genetic diversity of 130 New Zealand white rabbits, we obtained a set combination of 22 markers. Then, we performed a genetic analysis of 200 New Zealand white rabbits corresponding to two generations with this combination. It can be used to evaluate the breed conservation of rabbit germplasm resources.

Abstract

At present, there is an abundance of quality domestic rabbit breeds in China. However, due to the lack of technical standards for the genetic evaluation of rabbit germplasm resources, there have been a number of problems, such as poor breed conservation. By studying the genetic diversity of 130 New Zealand white rabbits (regardless of generation), we obtained the best simple sequence repeat (SSR) marker combination. We found that, when using microsatellite markers for the effective genetic evaluation of domestic rabbits, the number of records should be greater than 60 and the marker number more than 22. Through the comparative analysis of 30 combinations of 22 markers, the optimal combination of 22 markers was determined, and the 22 SSR polymorphic loci were distributed on different chromosomes. We performed a genetic analysis of 200 New Zealand white rabbits corresponding to two generations, using the best SSR polymorphic loci combination. There were no significant differences in the genetic diversity parameters between the two generations of rabbits (p > 0.05), indicating that the characteristics of this excellent rabbit germplasm have been effectively preserved. At the same time, we verified that the established method can be used to evaluate the breed conservation of rabbit germplasm resources.

Spatial Use and Survival of Sympatric Populations of New England and Eastern Cottontails in Connecticut

The New England cottontail (NEC) Sylvilagus transitionalis is strongly associated with shrubland and early successional habitat and is the only cottontail native to the U.S. Northeast. The distribution and abundance of young forest habitat and NEC populations have declined. The eastern cottontail (EC) Sylvilagus floridanus was introduced into the U.S. Northeast in the early 1900s and uses similar habitat as NEC, but is expanding in distribution and abundance. Little information exists on spatial use, survival, and competition in sympatric populations of NEC and EC. Understanding differences in population demographics may identify important factors or relationships influencing population trends and aid in developing effective management strategies. Our objectives were to quantify home range and core area sizes, annual survival rates, minimum population densities, and range overlap for sympatric populations of NEC and EC at four sites in Connecticut. We monitored spatial use and survival rates of 107 radio-collared rabbits over a 10-y period. Mean annual home ranges and core areas were 10.9 and 2.5 ha for NEC and 5.6 and 1.6 ha for EC. Overlap in home range and core areas was greater within species than between species (NEC-EC). For both species and sex, home range size expanded from winter to breeding seasons. Survival rates were greater for NEC than for EC at all four sites, with predation as the major cause of mortality for both species. Space-use patterns suggest that the potential for EC to interfere with NEC reproduction is limited and avoidance or resource partitioning between species in the same patch may be occurring.

Population structure and phylogenetic analysis of laboratory rabbits in Taiwan based on microsatellite markers

Laboratory rabbits used in Taiwan are primarily supplied by the Livestock Research Institute (LRI) and the Animal Drugs Inspection Branch (ADIB) of the Animal Health Research Institute. An analysis of the genetic characteristics and structure of these populations would thus be a fundamental step in building a long-term management programme for maintaining stable animal quality and preserving the genetic variation among the populations. In this study, DNA samples were isolated from founders of 5 populations: New Zealand White rabbits (NZW) and Japanese White rabbits (JPN) from the ADIB, NZW and Rex rabbits (REX) from the LRI, and NZW from a private rabbit breeding farm in Ban Ciao (BC). A set of microsatellite markers, 18 in total, was designed for genetic analysis. The average values for the allele number (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (H_E), and Wright’s fixation index (F_IS) were 5.50, 2.437, 0.442, 0.568 and 0.232, respectively. These results revealed that this set of microsatellite markers has high diversity and that the major local populations have a tendency toward inbreeding. At the same time, analysis of molecular variance results showed that the laboratory rabbits used in Taiwan have maintained a high level of within-population genetic differentiation (83%). The genetic differentiation among clusters was moderate (F_ST=0.18), and Bayesian cluster analysis showed that the most likely number of groups was 4 (K=4). Principal component analysis (PCA) also showed 4 divergent clusters. The LRI and BC NZW populations were not separated when K=4 was used in a Structure software analysis and were also hard to split until principal component 3 in PCA. The individual unrooted phylogenetic tree showed that the 5 populations were separated, except that some individuals from the LRI NZW population overlapped with the ADIB NZW and BC NZW populations. As such, in order to counteract the reduced F_IS (0.232) and maximise heterozygosity, the 3 NZW populations could be interbred or have new genes introduced into them. The set of microsatellite markers used herein was useful for studying the relationships and genetic diversities among these rabbit populations of Taiwan. Based on the resulting data, rabbit farms in Taiwan could select parental stocks for planned mating in the future as part of strategies to preserve and restore the rational breeding of laboratory rabbits.

Mitochondrial and nuclear DNA reveals reticulate evolution in hares (Lepus spp., Lagomorpha, Mammalia) from Ethiopia

For hares (Lepus spp., Leporidae, Lagomorpha, Mammalia) from Ethiopia no conclusive molecular phylogenetic data are available. To provide a first molecular phylogenetic model for the Abyssinian Hare (Lepus habessinicus), the Ethiopian Hare (L. fagani), and the Ethiopian Highland Hare (L. starcki) and their evolutionary relationships to hares from Africa, Eurasia, and North America, we phylogenetically analysed mitochondrial ATPase subunit 6 (ATP6; n = 153 / 416bp) and nuclear transferrin (TF; n = 155 / 434bp) sequences of phenotypically determined individuals. For the hares from Ethiopia, genotype composition at twelve microsatellite loci (n = 107) was used to explore both interspecific gene pool separation and levels of current hybridization, as has been observed in some other Lepus species. For phylogenetic analyses ATP6 and TF sequences of Lepus species from South and North Africa (L. capensis, L. saxatilis), the Anatolian peninsula and Europe (L. europaeus, L. timidus) were also produced and additional TF sequences of 18 Lepus species retrieved from GenBank were included as well. Median joining networks, neighbour joining, maximum likelihood analyses, as well as Bayesian inference resulted in similar models of evolution of the three species from Ethiopia for the ATP6 and TF sequences, respectively. The Ethiopian species are, however, not monophyletic, with signatures of contemporary uni- and bidirectional mitochondrial introgression and/ or shared ancestral polymorphism. Lepus habessinicus carries mtDNA distinct from South African L. capensis and North African L. capensis sensu lato; that finding is not in line with earlier suggestions of its conspecificity with L. capensis. Lepus starcki has mtDNA distinct from L. capensis and L. europaeus, which is not in line with earlier suggestions to include it either in L. capensis or L. europaeus. Lepus fagani shares mitochondrial haplotypes with the other two species from Ethiopia, despite its distinct phenotypic and microsatellite differences; moreover, it is not represented by a species-specific mitochondrial haplogroup, suggesting considerable mitochondrial capture by the other species from Ethiopia or species from other parts of Africa. Both mitochondrial and nuclear sequences indicate close phylogenetic relationships among all three Lepus species from Ethiopia, with L. fagani being surprisingly tightly connected to L. habessinicus. TF sequences suggest close evolutionary relationships between the three Ethiopian species and Cape hares from South and North Africa; they further suggest that hares from Ethiopia hold a position ancestral to many Eurasian and North American species.

REFERENCE AND BASELINE HEMATOCRIT MEASURES FOR THE THREATENED NEW ENGLAND COTTONTAIL (SYLVILAGUS TRANSITIONALIS) AND COMPARISON WITH SYMPATRIC EASTERN COTTONTAIL (SYLVILAGUS FLORIDANUS) RABBITS

From June 2014 to June 2015, capillary tube collections of blood were obtained concurrently with ear clips of trapped free-ranging, globally vulnerable New England cottontails (NEC; Sylvilagus transitionalis) and eastern cottontail rabbits (EC; Sylvilagus floridanus) in the Hudson Valley region of New York, United States. Species identification (NEC, EC) and sex (NEC) were determined genetically using a mitochondrial DNA assay and Y chromosome marker, respectively. Hematocrit values were obtained using a microhematocrit centrifuge. We provide the reference values 35.15-49.55 (2.5 and 97.5 percentiles) and 90% confidence intervals (CI) [lower: 33.00, 36.08; upper: 46.95, 51.00], for hematocrit of NEC. The mean hematocrit for NEC was 42.35% (SE = 0.58, n = 47) and a comparative contemporaneous mean in the same area for EC [39.96 (SE = 0.81, n = 26)], which was significantly different from NEC (P = 0.02). There was a significant sex difference for NEC [male: 43.99 (SE = 1.02, n = 28); female: 39.92 (SE = 0.78, n = 19), P < 0.0001], though not for EC.

A comprehensive linkage map and QTL map for carcass traits in a cross between Giant Grey and New Zealand White rabbits

BACKGROUND

Genomic resources for the rabbit are still limited compared to many other livestock species. The genomic sequence as well as linkage maps have gaps that hamper their use in rabbit genome research. Therefore, the aims of this study were the improvement of existing linkage maps and the mapping of quantitative trait loci (QTL) for carcass and meat quality traits. The study was performed in a F2 population of an initial cross between Giant Grey (GG) and New Zealand White (NZW) rabbits. The population consisted of 363 F2 animals derived from 9 F1 bucks and 33 F1 does. 186 microsatellite and three SNP markers were informative for mapping.

RESULTS

Out of 189 markers, which could be assigned to linkage groups, 110 markers were genetically mapped for the first time. The average marker distance was 7.8 cM. The map across all autosomes reached a total length of 1419 cM. The maternal linkage map was 1.4 times longer than the paternal. All linkage groups could be anchored to chromosomes. On the basis of the generated genetic map, we identified a highly significant QTL (genome-wide significance p < 0.01) for different carcass weights on chromosome 7 with a peak position at 91 cM (157 Mb), a significant QTL (p < 0.05) for bone mass on chromosome 9 at 61 cM (65 Mb), and another one for drip loss on chromosome 12 at 94 cM (128 Mb). Additional suggestive QTL were found on almost all chromosomes. Several genomic loci affecting the fore, intermediate and hind parts of the carcass were identified. The identified QTL explain between 2.5 to 14.6% of the phenotypic variance in the F2 population.

CONCLUSIONS

The results present the most comprehensive genetic map and the first genome-wide QTL mapping study for carcass and meat quality traits in rabbits. The identified QTL, in particular the major QTL on chromosome 7, provide starting points for fine mapping and candidate gene search. The data contribute to linking physical and genetic information in the rabbit genome.

Home-loving boreal hare mitochondria survived several invasions in Iberia: the relative roles of recurrent hybridisation and allele surfing

Genetic introgression from a resident species into an invading close relative can result from repeated hybridisation along the invasion front and/or allele surfing on the expansion wave. Cases where the phenomenon is massive and systematic, such as for hares (genus Lepus) in Iberia, would be best explained by recurrent hybridisation but this is difficult to prove because the donor populations are generally extinct. In the Pyrenean foothills, Lepus europaeus presumably replaced Lepus granatensis recently and the present species border is parallel to the direction of invasion, so that populations of L. granatensis in the contact zone represent proxies of existing variation before the invasion. Among three pairs of populations sampled across this border, we find less differentiation of mitochondrial DNA (mtDNA) across than along it, as predicted under recurrent hybridisation at the invasion front. Using autosomal microsatellite loci and X- and Y-linked diagnostic loci, we show that admixture across the border is quasi-absent, making it unlikely that lack of interspecific mtDNA differentiation results from ongoing gene flow. Furthermore, we find that the local species ranges are climatically contrasted, making it also unlikely that ongoing ecology-driven movement of the contact account for mtDNA introgression. The lack of mtDNA differentiation across the boundary is mostly due to sharing of mtDNA from a boreal species currently extinct in Iberia (Lepus timidus) whose mitochondria have thus remained in place since the last deglaciation despite successive invasions by two other species. Home-loving mitochondria thus witness past species distribution rather than ongoing exchanges across stabilised contact zones.

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.

Expressed sequence tags (ESTs) and single nucleotide polymorphisms (SNPs): What large-scale sequencing projects can tell us about ADME

To date over 800 complete genomes have been sequenced, with many more partially complete. Coupled with the large amount of mRNA transcript sequence data being produced from expression studies, there is now a daunting amount of information available to the research scientist. This review examines how this information may be best used, focusing on examples from sequences encoding absorption, distribution, metabolism and excretion (ADME)-related proteins in particular. Through the use of phylogenetic, splice variant and single nucleotide polymorphism (SNP) analysis, the review examines not only how insights into species-specific responses to drug exposure may be gained, but also how best to utilize this information to predict both individual human responses and the impact of population variance in response.

Localization and characterization of 170 BAC-derived clones and mapping of 94 microsatellites in the Hessian fly

Ninety-four microsatellites from enriched genomic libraries of Hessian fly (Hf, Mayetiola destructor [Say]) were localized to 170 cognate clones in an Hf bacterial artificial chromosome (BAC) library. These microsatellite-positive BAC clones were physically mapped to polytene chromosomes by fluorescent in situ hybridization. The mapped microsatellite loci can be used to study the genetic diversity and population structure of Hf.

Reduced introgression of the Y chromosome between subspecies of the European rabbit (Oryctolagus cuniculus) in the Iberian Peninsula

The role of the Y chromosome in speciation is unclear. Hybrid zones provide natural arenas for studying speciation, as differential introgression of markers may reveal selection acting against incompatibilities. Two subspecies of the European rabbit (Oryctolagus cuniculus) form a hybrid zone in the Iberian Peninsula. Previous work on mitochondrial DNA (mtDNA), Y- and X-linked loci revealed the existence of two divergent lineages in the rabbit genome and that these lineages are largely subspecies-specific for mtDNA and two X-linked loci. Here we investigated the geographic distribution of the two Y chromosome lineages by genotyping two diagnostic single nucleotide polymorphisms in a sample of 353 male rabbits representing both subspecies, and found that Y chromosome lineages are also largely subspecies-specific. We then sequenced three autosomal loci and discovered considerable variation in levels of differentiation at these loci. Finally, we compared estimates of population differentiation between rabbit subspecies at 26 markers and found a surprising bimodal distribution of F(ST)values. The vast majority of loci showed little or no differentiation between rabbit subspecies while a few loci, including the SRY gene, showed little or no introgression across the hybrid zone. Estimates of population differentiation for the Y chromosome were surprisingly high given that there is male-biased dispersal in rabbits. Taken together, these data indicate that there is a clear dichotomy in the rabbit genome and that some loci remain highly differentiated despite extensive gene flow following secondary contact.

A first-generation microsatellite-based integrated genetic and cytogenetic map for the European rabbit (Oryctolagus cuniculus) and localization of angora and albino

Although the European rabbit (Oryctolagus cuniculus) is used both in agronomics and in research, genomic resources for this species are still limited and no microsatellite-based genetic map has been reported. Our aim was to construct a rabbit genetic map with cytogenetically mapped microsatellites so as to build an integrated genetic and cytogenetic map. A reference population of 187 rabbits comprising eight three-generation families with 10-25 offspring per family was produced. One hundred and ninety-four of 305 previously identified microsatellites were included in this study. Of these, 158 were polymorphic with two to seven alleles. The map reported here comprises 111 markers, including 104 INRA microsatellites, five microsatellites from another source and two phenotypic markers (angora and albino). Ninety markers were integrated into 20 linkage groups. The remaining 21 microsatellites mapped to separate linkage groups, 19 with a precise cytogenetic position and two with only a chromosomal assignment. The genetic map spans 2766.6 cM and covers 20 rabbit chromosomes, excluding chromosomes 20, 21 and X. The density of this map is limited, but we used it to verify the location of angora and albino on chromosomes 15q and 1q, respectively, in agreement with previously published data. This first generation genetic/cytogenetic map will help gene identification and quantitative trait loci mapping projects in rabbit.