Genomic characterization of the world's longest selection experiment in mouse reveals the complexity of polygenic traits.

BACKGROUND: Long-term selection experiments are a powerful tool to understand the genetic background of complex traits. The longest of such experiments has been conducted in the Research Institute for Farm Animal Biology (FBN), generating extreme mouse lines with increased fertility, body mass, protein mass and endurance. For >140 generations, these lines have been maintained alongside an unselected control line, representing a valuable resource for understanding the genetic basis of polygenic traits. However, their history and genomes have not been reported in a comprehensive manner yet. Therefore, the aim of this study is to provide a summary of the breeding history and phenotypic traits of these lines along with their genomic characteristics. We further attempt to decipher the effects of the observed line-specific patterns of genetic variation on each of the selected traits. RESULTS: Over the course of >140 generations, selection on the control line has given rise to two extremely fertile lines (>20 pups per litter each), two giant growth lines (one lean, one obese) and one long-distance running line. Whole genome sequencing analysis on 25 animals per line revealed line-specific patterns of genetic variation among lines, as well as high levels of homozygosity within lines. This high degree of distinctiveness results from the combined effects of long-term continuous selection, genetic drift, population bottleneck and isolation. Detection of line-specific patterns of genetic differentiation and structural variation revealed multiple candidate genes behind the improvement of the selected traits. CONCLUSIONS: The genomes of the Dummerstorf trait-selected mouse lines display distinct patterns of genomic variation harbouring multiple trait-relevant genes. Low levels of within-line genetic diversity indicate that many of the beneficial alleles have arrived to fixation alongside with neutral alleles. This study represents the first step in deciphering the influence of selection and neutral evolutionary forces on the genomes of these extreme mouse lines and depicts the genetic complexity underlying polygenic traits.

African and Asian leopards are highly differentiated at the genomic level.

Leopards are the only big cats still widely distributed across the continents of Africa and Asia. They occur in a wide range of habitats and are often found in close proximity to humans. But despite their ubiquity, leopard phylogeography and population history have not yet been studied with genomic tools. Here, we present population-genomic data from 26 modern and historical samples encompassing the vast geographical distribution of this species. We find that Asian leopards are broadly monophyletic with respect to African leopards across almost their entire nuclear genomes. This profound genetic pattern persists despite the animals' high potential mobility, and despite evidence of transfer of African alleles into Middle Eastern and Central Asian leopard populations within the last 100,000 years. Our results further suggest that Asian leopards originated from a single out-of-Africa dispersal event 500-600 thousand years ago and are characterized by higher population structuring, stronger isolation by distance, and lower heterozygosity than African leopards. Taxonomic categories do not take into account the variability in depth of divergence among subspecies. The deep divergence between the African subspecies and Asian populations contrasts with the much shallower divergence among putative Asian subspecies. Reconciling genomic variation and taxonomy is likely to be a growing challenge in the genomics era.

Ancestral mitogenome capture of the Southeast Asian banded linsang.

Utilising a reconstructed ancestral mitochondrial genome of a clade to design hybridisation capture baits can provide the opportunity for recovering mitochondrial sequences from all its descendent and even sister lineages. This approach is useful for taxa with no extant close relatives, as is often the case for rare or extinct species, and is a viable approach for the analysis of historical museum specimens. Asiatic linsangs (genus Prionodon) exemplify this situation, being rare Southeast Asian carnivores for which little molecular data is available. Using ancestral capture we recover partial mitochondrial genome sequences for seven banded linsangs (P. linsang) from historical specimens, representing the first intraspecific genetic dataset for this species. We additionally assemble a high quality mitogenome for the banded linsang using shotgun sequencing for time-calibrated phylogenetic analysis. This reveals a deep divergence between the two Asiatic linsang species (P. linsang, P. pardicolor), with an estimated divergence of ~12 million years (Ma). Although our sample size precludes any robust interpretation of the population structure of the banded linsang, we recover two distinct matrilines with an estimated tMRCA of ~1 Ma. Our results can be used as a basis for further investigation of the Asiatic linsangs, and further demonstrate the utility of ancestral capture for studying divergent taxa without close relatives.

Adaptive responses of animals to climate change are most likely insufficient.

Biological responses to climate change have been widely documented across taxa and regions, but it remains unclear whether species are maintaining a good match between phenotype and environment, i.e. whether observed trait changes are adaptive. Here we reviewed 10,090 abstracts and extracted data from 71 studies reported in 58 relevant publications, to assess quantitatively whether phenotypic trait changes associated with climate change are adaptive in animals. A meta-analysis focussing on birds, the taxon best represented in our dataset, suggests that global warming has not systematically affected morphological traits, but has advanced phenological traits. We demonstrate that these advances are adaptive for some species, but imperfect as evidenced by the observed consistent selection for earlier timing. Application of a theoretical model indicates that the evolutionary load imposed by incomplete adaptive responses to ongoing climate change may already be threatening the persistence of species.

Human-mediated introduction of introgressed deer across Wallace's line: Historical biogeography of Rusa unicolor and R. timorensis.

In this study we compared the phylogeographic patterns of two Rusa species, Rusa unicolor and Rusa timorensis, in order to understand what drove and maintained differentiation between these two geographically and genetically close species and investigated the route of introduction of individuals to the islands outside of the Sunda Shelf. We analyzed full mitogenomes from 56 archival samples from the distribution areas of the two species and 18 microsatellite loci in a subset of 16 individuals to generate the phylogeographic patterns of both species. Bayesian inference with fossil calibration was used to estimate the age of each species and major divergence events. Our results indicated that the split between the two species took place during the Pleistocene, ~1.8 Mya, possibly driven by adaptations of R. timorensis to the drier climate found on Java compared to the other islands of Sundaland. Although both markers identified two well-differentiated clades, there was a largely discrepant pattern between mitochondrial and nuclear markers. While nDNA separated the individuals into the two species, largely in agreement with their museum label, mtDNA revealed that all R. timorensis sampled to the east of the Sunda shelf carried haplotypes from R. unicolor and one Rusa unicolor from South Sumatra carried a R. timorensis haplotype. Our results show that hybridization occurred between these two sister species in Sundaland during the Late Pleistocene and resulted in human-mediated introduction of hybrid descendants in all islands outside Sundaland.

Fast and cost-effective single nucleotide polymorphism (SNP) detection in the absence of a reference genome using semideep next-generation Random Amplicon Sequencing (RAMseq).

Biodiversity has suffered a dramatic global decline during the past decades, and monitoring tools are urgently needed providing data for the development and evaluation of conservation efforts both on a species and on a genetic level. However, in wild species, the assessment of genetic diversity is often hampered by the lack of suitable genetic markers. In this article, we present Random Amplicon Sequencing (RAMseq), a novel approach for fast and cost-effective detection of single nucleotide polymorphisms (SNPs) in nonmodel species by semideep sequencing of random amplicons. By applying RAMseq to the Eurasian otter (Lutra lutra), we identified 238 putative SNPs after quality filtering of all candidate loci and were able to validate 32 of 77 loci tested. In a second step, we evaluated the genotyping performance of these SNP loci in noninvasive samples, one of the most challenging genotyping applications, by comparing it with genotyping results of the same faecal samples at microsatellite markers. We compared (i) polymerase chain reaction (PCR) success rate, (ii) genotyping errors and (iii) Mendelian inheritance (population parameters). SNPs produced a significantly higher PCR success rate (75.5% vs. 65.1%) and lower mean allelic error rate (8.8% vs. 13.3%) than microsatellites, but showed a higher allelic dropout rate (29.7% vs. 19.8%). Genotyping results showed no deviations from Mendelian inheritance in any of the SNP loci. Hence, RAMseq appears to be a valuable tool for the detection of genetic markers in nonmodel species, which is a common challenge in conservation genetic studies.

Homogenous Population Genetic Structure of the Non-Native Raccoon Dog (Nyctereutes procyonoides) in Europe as a Result of Rapid Population Expansion.

The extent of gene flow during the range expansion of non-native species influences the amount of genetic diversity retained in expanding populations. Here, we analyse the population genetic structure of the raccoon dog (Nyctereutes procyonoides) in north-eastern and central Europe. This invasive species is of management concern because it is highly susceptible to fox rabies and an important secondary host of the virus. We hypothesized that the large number of introduced animals and the species' dispersal capabilities led to high population connectivity and maintenance of genetic diversity throughout the invaded range. We genotyped 332 tissue samples from seven European countries using 16 microsatellite loci. Different algorithms identified three genetic clusters corresponding to Finland, Denmark and a large 'central' population that reached from introduction areas in western Russia to northern Germany. Cluster assignments provided evidence of long-distance dispersal. The results of an Approximate Bayesian Computation analysis supported a scenario of equal effective population sizes among different pre-defined populations in the large central cluster. Our results are in line with strong gene flow and secondary admixture between neighbouring demes leading to reduced genetic structuring, probably a result of its fairly rapid population expansion after introduction. The results presented here are remarkable in the sense that we identified a homogenous genetic cluster inhabiting an area stretching over more than 1500km. They are also relevant for disease management, as in the event of a significant rabies outbreak, there is a great risk of a rapid virus spread among raccoon dog populations.

Paternal heat exposure causes DNA methylation and gene expression changes of Stat3 in Wild guinea pig sons.

Epigenetic mechanisms convey environmental information through generations and can regulate gene expression. Epigenetic studies in wild mammals are rare, but enable understanding adaptation processes as they may occur in nature. In most wild mammal species, males are the dispersing sex and thus often have to cope with differing habitats and thermal changes more rapidly than the often philopatric females. As temperature is a major environmental selection factor, we investigated whether genetically heterogeneous Wild guinea pig (Cavia aperea) males adapt epigenetically to an increase in temperature, whether that response will be transmitted to the next generation(s), and whether it regulates mRNA expression. Five (F0) adult male guinea pigs were exposed to an increased ambient temperature for 2 months, corresponding to the duration of the species' spermatogenesis. To study the effect of heat, we focused on the main thermoregulatory organ, the liver. We analyzed CpG-methylation changes of male offspring (F1) sired before and after the fathers' heat treatment (as has recently been described in Weyrich et al. [Mol. Ecol., 2015]). Transcription analysis was performed for the three genes with the highest number of differentially methylated changes detected: the thermoregulation gene Signal Transducer and Activator of Transcription 3 (Stat3), the proteolytic peptidase gene Cathepsin Z (Ctsz), and Sirtuin 6 (Sirt6) with function in epigenetic regulation. Stat3 gene expression was significantly reduced (P < 0.05), which indicated a close link between CpG-methylation and expression levels for this gene. The two other genes did not show gene expression changes. Our results indicate the presence of a paternal transgenerational epigenetic effect. Quick adaptation to climatic changes may become increasingly relevant for the survival of wildlife species as global temperatures are rising.

Impact of enrichment conditions on cross-species capture of fresh and degraded DNA.

By combining high-throughput sequencing with target enrichment ('hybridization capture'), researchers are able to obtain molecular data from genomic regions of interest for projects that are otherwise constrained by sample quality (e.g. degraded and contamination-rich samples) or a lack of a priori sequence information (e.g. studies on nonmodel species). Despite the use of hybridization capture in various fields of research for many years, the impact of enrichment conditions on capture success is not yet thoroughly understood. We evaluated the impact of a key parameter--hybridization temperature--on the capture success of mitochondrial genomes across the carnivoran family Felidae. Capture was carried out for a range of sample types (fresh, archival, ancient) with varying levels of sequence divergence between bait and target (i.e. across a range of species) using pools of individually indexed libraries on Agilent SureSelect(™) arrays. Our results suggest that hybridization capture protocols require specific optimization for the sample type that is being investigated. Hybridization temperature affected the proportion of on-target sequences following capture: for degraded samples, we obtained the best results with a hybridization temperature of 65 °C, while a touchdown approach (65 °C down to 50 °C) yielded the best results for fresh samples. Evaluation of capture performance at a regional scale (sliding window approach) revealed no significant improvement in the recovery of DNA fragments with high sequence divergence from the bait at any of the tested hybridization temperatures, suggesting that hybridization temperature may not be the critical parameter for the enrichment of divergent fragments.

Whole genome sequencing and methylome analysis of the wild guinea pig.

BACKGROUND: DNA methylation is a heritable mechanism that acts in response to environmental changes, lifestyle and diseases by influencing gene expression in eukaryotes. Epigenetic studies of wild organisms are mandatory to understand their role in e.g. adaptational processes in the great variety of ecological niches. However, strategies to address those questions on a methylome scale are widely missing. In this study we present such a strategy and describe a whole genome sequence and methylome analysis of the wild guinea pig. RESULTS: We generated a full Wild guinea pig (Cavia aperea) genome sequence with enhanced coverage of methylated regions, benefiting from the available sequence of the domesticated relative Cavia porcellus. This new genome sequence was then used as reference to map the sequence reads of bisulfite treated Wild guinea pig sequencing libraries to investigate DNA-methylation patterns at nucleotide-specific level, by using our here described method, named 'DNA-enrichment-bisulfite-sequencing' (MEBS). The results achieved using MEBS matched those of standard methods in other mammalian model species. The technique is cost efficient, and incorporates both methylation enrichment results and a nucleotide-specific resolution even without a whole genome sequence available. Thus MEBS can be easily applied to extend methylation enrichment studies to a nucleotide-specific level. CONCLUSIONS: The approach is suited to study methylomes of not yet sequenced mammals at single nucleotide resolution. The strategy is transferable to other mammalian species by applying the nuclear genome sequence of a close relative. It is therefore of interest for studies on a variety of wild species trying to answer evolutionary, adaptational, ecological or medical questions by epigenetic mechanisms.

Rabbit haemorrhagic disease: virus persistence and adaptation in Australia.

In Australia, the rabbit haemorrhagic disease virus (RHDV) has been used since 1996 to reduce numbers of introduced European rabbits (Oryctolagus cuniculus) which have a devastating impact on the native Australian environment. RHDV causes regular, short disease outbreaks, but little is known about how the virus persists and survives between epidemics. We examined the initial spread of RHDV to show that even upon its initial spread, the virus circulated continuously on a regional scale rather than persisting at a local population level and that Australian rabbit populations are highly interconnected by virus-carrying flying vectors. Sequencing data obtained from a single rabbit population showed that the viruses that caused an epidemic each year seldom bore close genetic resemblance to those present in previous years. Together, these data suggest that RHDV survives in the Australian environment through its ability to spread amongst rabbit subpopulations. This is consistent with modelling results that indicated that in a large interconnected rabbit meta-population, RHDV should maintain high virulence, cause short, strong disease outbreaks but show low persistence in any given subpopulation. This new epidemiological framework is important for understanding virus-host co-evolution and future disease management options of pest species to secure Australia's remaining natural biodiversity.

Noninvasive molecular and morphological evidences for an undiscovered population of snow vole in Southern Spain.

Capturing wild animals can be time consuming and difficult or even impractical. Noninvasive sampling is potentially a cost-effective and efficient means to monitor wild animals, thereby avoiding the need of capture and disturb species in the wild. On the basis of the morphological and genetic analyses of owl pellet contents, a so far undetected European snow vole (Chionomys nivalis) population was discovered in the Sierra Segura mountain range (Southern Spain). The mtDNA sequence from the newly discovered haplotype was compared with sequences from vole populations of the Sierra Nevada and Sierra Peñalara mountain ranges (Spain) and from Churwalden (Switzerland). The nine recovered haplotypes clustered in four distinct lineages according to their geographical origin. The vole sequence from the Sierra Segura owl pellet belonged to a new haplotype, constituting a new lineage. The evolutionary divergence between sequences from the Sierra Segura and other Spanish populations was higher than that among other Spanish haplotypes. The new snow vole haplotype from this new locality duplicates the number of occurrence sites of this critically endangered species in Southern Spain, which is of great interest for further conservation and management plans of the European snow vole in the most southwestern area of its entire distribution range.

Applicability of major histocompatibility complex DRB1 alleles as markers to detect vertebrate hybridization: a case study from Iberian ibex × domestic goat in southern Spain.

BACKGROUND: Hybridization between closely related wild and domestic species is of great concern because it can alter the evolutionary integrity of the affected populations. The high allelic variability of Major Histocompatibility Complex (MHC) loci usually excludes them from being used in studies to detect hybridization events. However, if a) the parental species don't share alleles, and b) one of the parental species possesses an exceptionally low number of alleles (to facilitate analysis), then even MHC loci have the potential to detect hybrids. RESULTS: By genotyping the exon2 of the MHC class II DRB1 locus, we were able to detect hybridization between domestic goats (Capra hircus) and free-ranging Iberian ibex (Capra pyrenaica hispanica) by molecular means. CONCLUSIONS: This is the first documentation of a Capra pyrenaica × Capra hircus hybridization, which presented us the opportunity to test the applicability of MHC loci as new, simple, cost-effective, and time-saving approach to detect hybridization between wild species and their domesticated relatives, thus adding value to MHC genes role in animal conservation and management.

Epidemiology, control and management of an EBHS outbreak in captive hares.

Here we describe an outbreak of European brown hare syndrome (EBHS) in a captive hare population. The EBHS outbreak occurred in March 2009, at the beginning of the breeding season. Overall mortality was 53% out of an original population of 61 animals. Animals between five and eleven months showed a significantly higher mortality rate than other age classes. Pregnant females either aborted their foetuses and survived or died pregnant. All foetuses (n=10) of the pregnant hares were PCR positive for EBHSV. Only one offspring born during the outbreak survived. Shortly after the outbreak, the surviving hares developed a specific anti-EBHSV titre between 1:80 and 1:2560, which dropped to 1:10-1:160 nine months later. Hares between one and three years of age developed a significantly higher titre than hares younger than one year or older than four years. Offspring born after the outbreak showed a lower titre of 1:10, indicating passive antibody transfer via placenta and milk. After two months, the titre was not detectable any longer. In December 2009, the captive population was vaccinated against EBHS virus with inactivated virus prepared from the organs of infected hares. The titres after the first vaccination ranged from 1:10 to 1:640, and after the second vaccination from 1:10 to 1:320. To estimate the effect of EBHS on reproduction, we compared the breeding seasons 2008 and 2009. Several possible sources of infection of the colony are discussed, but the definite cause could not be determined.

A TaqMan real-time PCR-based assay for the identification of Fasciola spp.

Real time quantitative PCR (qPCR) is one of the key technologies of the post-genome era, with clear advantages compared to normal end-point PCR. In this paper, we report the first qPCR-based assay for the identification of Fasciola spp. Based on sequences of the second internal transcribed spacers (ITS-2) of the ribosomal rRNA gene, we used a set of genus-specific primers for Fasciola ITS-2 amplification, and we designed species-specific internal TaqMan probes to identify F. hepatica and F. gigantica, as well as the hybrid 'intermediate'Fasciola. These primers and probes were used for the highly specific, sensitive, and simple identification of Fasciola species collected from different animal host from China, Spain, Niger and Egypt. The novel qPCR-based technique for the identification of Fasciola spp. may provide a useful tool for the epidemiological investigation of Fasciola infection, including their intermediate snail hosts.

A fluorescence-based polymerase chain reaction-linked single-strand conformation polymorphism (F-PCR-SSCP) assay for the identification of Fasciola spp.

The present study aimed to establish a fluorescence-based polymerase chain reaction-linked single-strand conformation polymorphism (F-PCR-SSCP) assay for the identification of Fasciola spp. Based on the sequences of the second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA, we designed a set of genus-specific primers for the amplification of Fasciola ITS-2, with an estimated size of 140 bp. These primers were labelled by fluorescence dyes, and the PCR products were analyzed by capillary electrophoresis under non-denaturing conditions (F-PCR-SSCP). Capillary electrophoresis analysis of the fluorescence-labelled DNA fragments displayed three different peak profiles that allowed the accurate identification of Fasciola species: one single peak specific for either Fasciola hepatica or Fasciola gigantica and a doublet peak corresponding to the "intermediate" Fasciola. Validation of our novel method was performed using Fasciola specimens from different host animals from China, Spain, Nigeria, and Egypt. This F-PCR-SSCP assay provides a rapid, simple, and robust tool for the identification and differentiation between Fasciola spp.

Superconception in mammalian pregnancy can be detected and increases reproductive output per breeding season.

The concept of superfetation, a second conception during pregnancy, has been controversial for a long time. In this paper we use an experimental approach to demonstrate that female European brown hares (Lepus europaeus) frequently develop a second pregnancy while already pregnant and thereby increase their reproductive success. After a new, successful copulation, we confirmed additional ovulations before parturition in living, late-pregnant females by detecting a second set of fresh corpora lutea using high-resolution ultrasonography. The presence of early embryonic stages in the oviduct, demonstrated by oviduct flushing, next to fully developed fetuses in the uterus is best explained by passage of semen through the late-pregnant uterus; this was confirmed by paternity analysis using microsatellite profiling. Subsequent implantation occurred after parturition. This superfetation, categorized as superconception, significantly increased litter size and permitted females to produce up to 35.4% more offspring per breeding season. It is therefore most likely an evolutionary adaptation.

Investigations on reproductive physiology in the male Eurasian lynx (Lynx lynx).

This study characterized (in vivo) morphological and functional parameters of reproductive organs of adult male lynx (n = 3) prior to, during, and after the breeding season (n = 3). Size and morphology of the reproductive tract were monitored by transcutaneous (testes) and transrectal (accessory sex glands) ultrasonography. Semen was collected by electroejaculation. Ejaculate volume, sperm number, motility, and morphology of spermatozoa as well as testosterone concentrations in blood serum and feces were evaluated. The testes and prostate had seasonal changes in size and echotexture. The mean (+/- S.D.) maximum and minimum testicular volume were 2.8 +/- 0.8 cm3 and 1.5 +/- 0.3 cm3, respectively. Fecal testosterone concentrations were highest in February (1240 +/- 393 ng/g feces), with a second increase in May (971 +/- 202 ng/g feces), but concentrations were lowest in January (481 +/- 52.9 ng/g feces). Ejaculate volume, total sperm number and percentage of motile, and intact spermatozoa were maximal in March (the middle of the breeding season). In one of the eight litters, multiple paternity was proven; however, in the remaining seven litters, all 16 cubs were sired by the same male. This particular male had the most developed and active testes and best semen quality, which may be important for sperm competition.

Tannin-binding salivary proteins in three captive rhinoceros species.

Tannin-binding salivary proteins (TBSP) are considered to be counter-defences acquired in the course of evolution by animals whose natural forage contains such tannins. As tannins mostly occur in browse material but not in grasses, it is assumed that grazers do not have a need for TBSP. Whereas it has been shown in several non-ungulate species that TBSP can be induced by dietary tannins, their presence or absence in ungulates has, so far, been shown to be a species-specific characteristic independent of dietary manipulations. We investigated saliva from three rhinoceros species from zoological gardens fed comparable, conventional zoo diets. As expected, saliva from white rhinoceroses (Ceratotherum simum, grazer) had lower tannin-binding capacities than that from black rhinoceroses (Diceros bicornis, browser). Surprisingly, however, Indian rhinoceroses (Rhinoceros unicornis), commonly regarded as grazers as well, displayed the highest tannin-binding capacities of the three species investigated. It is speculated that this discrepancy might be a result of an evolutionarily recent switch to a grass-dominated diet in Indian rhinoceroses, and that the black rhinoceros, which is closer related to the white rhinoceros than the Indian species, has evolved an inducible mechanism of TBSP production. In separate trials during which the tannin content of the diets of black rhinoceroses was increased by the addition of either tannic acid or quebracho, the tannin-binding capacity of black rhinoceros saliva was increased to levels within the same range as that of Indian rhinoceroses on the conventional diets. While induction trials in white and Indian rhinoceroses remain to be performed for a full understanding of salivary anti-tannin defence in rhinoceroses, these results are the first report of an induced salivary response to increased dietary tannin levels in an ungulate species.