José Antonio Blanco-Aguiar

A Comparison of Brain Gene Expression Levels in Domesticated and Wild Animals

Domestication has led to similar changes in morphology and behavior in several animal species, raising the question whether similarities between different domestication events also exist at the molecular level. We used mRNA sequencing to analyze genome-wide gene expression patterns in brain frontal cortex in three pairs of domesticated and wild species (dogs and wolves, pigs and wild boars, and domesticated and wild rabbits). We compared the expression differences with those between domesticated guinea pigs and a distant wild relative (Cavia aperea) as well as between two lines of rats selected for tameness or aggression towards humans. There were few gene expression differences between domesticated and wild dogs, pigs, and rabbits (30–75 genes (less than 1%) of expressed genes were differentially expressed), while guinea pigs and C. aperea differed more strongly. Almost no overlap was found between the genes with differential expression in the different domestication events. In addition, joint analyses of all domesticated and wild samples provided only suggestive evidence for the existence of a small group of genes that changed their expression in a similar fashion in different domesticated species. The most extreme of these shared expression changes include up-regulation in domesticates of SOX6 and PROM1, two modulators of brain development. There was almost no overlap between gene expression in domesticated animals and the tame and aggressive rats. However, two of the genes with the strongest expression differences between the rats (DLL3 and DHDH) were located in a genomic region associated with tameness and aggression, suggesting a role in influencing tameness. In summary, the majority of brain gene expression changes in domesticated animals are specific to the given domestication event, suggesting that the causative variants of behavioral domestication traits may likewise be different.

Speciation in the European rabbit (Oryctolagus cuniculus): islands of differentiation on the X chromosome and autosomes

Studies of gene flow between closely related taxa can provide insight into the genetic basis of speciation. To evaluate the importance of the X chromosome in reproductive isolation between subspecies of the European rabbit and to study the genomic scale over which islands of differentiation extend, we resequenced a total of 34 loci distributed along the X chromosome and chromosome 14. Previous studies based on few markers suggested that loci in centromeric regions were highly differentiated between rabbit subspecies, whereas loci in telomeric regions were less differentiated. Here, we confirmed this finding but also discovered remarkable variation in levels of differentiation among loci, with FST values from nearly 0 to 1. Analyses using isolation‐with‐migration models suggest that this range appears to be largely explained by differential levels of gene flow among loci. The X chromosome was significantly more differentiated than the autosomes. On chromosome 14, differentiation decayed very rapidly at increasing distances from the centromere, but on the X chromosome distinct islands of differentiation encompassing several megabases were observed both at the centromeric region and along the chromosome arms. These findings support the idea that the X chromosome plays an important role in reproductive isolation between rabbit subspecies. These results also demonstrate the mosaic nature of the genome at species boundaries.

The genomic architecture of population divergence between subspecies of the european rabbit

The analysis of introgression of genomic regions between divergent populations provides an excellent opportunity to determine the genetic basis of reproductive isolation during the early stages of speciation. However, hybridization and subsequent gene flow must be relatively common in order to localize individual loci that resist introgression. In this study, we used next-generation sequencing to study genome-wide patterns of genetic differentiation between two hybridizing subspecies of rabbits (Oryctolagus cuniculus algirus and O. c. cuniculus) that are known to undergo high rates of gene exchange. Our primary objective was to identify specific genes or genomic regions that have resisted introgression and are likely to confer reproductive barriers in natural conditions. On the basis of 326,000 polymorphisms, we found low to moderate overall levels of differentiation between subspecies, and fewer than 200 genomic regions dispersed throughout the genome showing high differentiation consistent with a signature of reduced gene flow. Most differentiated regions were smaller than 200 Kb and contained very few genes. Remarkably, 30 regions were each found to contain a single gene, facilitating the identification of candidate genes underlying reproductive isolation. This gene-level resolution yielded several insights into the genetic basis and architecture of reproductive isolation in rabbits. Regions of high differentiation were enriched on the X-chromosome and near centromeres. Genes lying within differentiated regions were often associated with transcription and epigenetic activities, including chromatin organization, regulation of transcription, and DNA binding. Overall, our results from a naturally hybridizing system share important commonalities with hybrid incompatibility genes identified using laboratory crosses in mice and flies, highlighting general mechanisms underlying the maintenance of reproductive barriers.

Spatial variation in helminth community structure in the red-legged partridge ( Alectoris rufa L.): effects of definitive host density

Parasite community ecology has recently focused on understanding the forces structuring these communities. There are few surveys, however, designed to study the spatial repeatability and predictability of parasite communities at the local scale in one host. The purpose of our study was to address the relationship between infracommunity and component community richness, and to describe spatial variations on the local scale, of helminth parasite communities in an avian host, the red-legged partridge (Alectoris rufta). We sampled 235 wild partridges from 8 separate localities, with different partridge population densities, in the Ciudad Real and Toledo provinces of central Spain, and we determined their overall and intestinal helminth species. We found that habitat variables (mean temperature and land use) were not significantly associated with any component community. The partridge population abundance index was directly correlated with the prevalence and mean intensity of infection but not with component community species richness. There was a curvilinear relationship between infracommunity and component community species richness, as well as negative interspecific associations, for the helminth species assemblage parasitizing the intestine. A nestedness/anti-nestedness pattern, considered as part of a continuum, was associated with prevalence, mean intensity and partridge population abundance index, but not with component community richness. Increases in the partridge population abundance index and the prevalence and mean intensity of infection were associated with increases in helminth community nestedness. Although negative interactions between helminth species could not be ruled out as forces structuring helminth communities, our results suggest that parasite community structure in the red-legged partridge was primarily determined by the extrinsic influence of parasite habitat heterogeneity and its amplification of the differing probabilities of colonization of parasite species.

Evidence for widespread positive and purifying selection across the European rabbit (Oryctolagus cuniculus) genome

The nearly neutral theory of molecular evolution predicts that the efficacy of both positive and purifying selection is a function of the long-term effective population size (N(e)) of a species. Under this theory, the efficacy of natural selection should increase with N(e). Here, we tested this simple prediction by surveying ~1.5 to 1.8 Mb of protein coding sequence in the two subspecies of the European rabbit (Oryctolagus cuniculus algirus and O. c. cuniculus), a mammal species characterized by high levels of nucleotide diversity and N(e) estimates for each subspecies on the order of 1 × 10(6). When the segregation of slightly deleterious mutations and demographic effects were taken into account, we inferred that >60% of amino acid substitutions on the autosomes were driven to fixation by positive selection. Moreover, we inferred that a small fraction of new amino acid mutations (<4%) are effectively neutral (defined as 0 < N(e)s < 1) and that this fraction was negatively correlated with a genes expression level. Consistent with models of recurrent adaptive evolution, we detected a negative correlation between levels of synonymous site polymorphism and the rate of protein evolution, although the correlation was weak and nonsignificant. No systematic X chromosome-autosome difference was found in the efficacy of selection. For example, the proportion of adaptive substitutions was significantly higher on the X chromosome compared with the autosomes in O. c. algirus but not in O. c. cuniculus. Our findings support widespread positive and purifying selection in rabbits and add to a growing list of examples suggesting that differences in N(e) among taxa play a substantial role in determining rates and patterns of protein evolution.

Phylogeography and genetic structure of the red‐legged partridge (Alectoris rufa): more evidence for refugia within the Iberian glacial refugium

The Pleistocene climatic oscillations promoted the diversification in avian species during the last glacial period. The red‐legged partridge (Alectoris rufa, Family Phasianidae) has a large natural distribution extending from the Mediterranean to humid temperate zones. However, the genetic structure for this species is unknown. The present study investigates the phylogeography, genetic structure and demographic history of A. rufa across its distribution, employing both mitochondrial DNA control region sequences and nuclear microsatellite loci. Our results propose that this species was greatly affected by Pleistocene glaciations. The mismatch analyses suggest that the current populations resulted from post‐glacial expansion and subsequent differentiation resulting in five diagnosable genetic clusters: Southwestern, Central‐eastern, Northwestern, Balearic and French and Italian. Further, we found evidence of three glacial refugia within the currently recognized Iberian glacial refugium. The intraspecific structure revealed by both maternal and biparental phylogeographic analyses was not resolved in the phylogenetic analyses. Based on all considerations, we recommended that five management units be recognized.

Is the interaction between rabbit hemorrhagic disease and hyperpredation by raptors a major cause of the red-legged partridge decline in Spain?

Hyperpredation can be described as a restrictive case of apparent competition where an increased number of primary prey species indirectly induces the decrease of the secondary prey species through numerical response of predators to the primary prey dynamics. It has been proposed that rabbit hemorrhagic disease (RHD), which decimated populations of European wild rabbit (Oryctolagus cuniculus) in Spain, led to prey switching by raptors towards red-legged partridges (Alectoris rufa) causing declines in their populations as a peculiar case of hyperpredation.We have reviewed field studies that disagree with estimates (based on bag records) of recent increases of rabbit and partridge numbers in Spain. Because of an increase in releases of farm-reared animals, there is doubt about the use of bag records to accurately estimate population trends in recent decades. We also provide new data and discuss some analytical considerations related to the temporal and spatial scales that might affect the interpretation of data. Finally, we discuss why studies associated with raptor food habits, predator population dynamics and predator distribution suggest that diet data are not sufficient to link patterns and processes. Although we agree that the RHD outbreak has markedly affected rabbit populations, which has indirectly affected many other species in the Iberian Peninsula, we consider that hyperpredation mediated by raptors has not been clearly demonstrated. In contrast, endorsing the hypothesis of raptor-mediated hyperpredation without sufficient proof may have conflictive consequences if we consider the increased persecution of raptors in recent decades in Spain.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole‐genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein–protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome‐wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

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

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

Full genome sequences are key to disclose RHDV2 emergence in the Macaronesian islands

A recent publication by Carvalho et al. in “Virus Genes” (June 2017) reported the presence of the new variant of rabbit hemorrhagic disease virus (RHDV2) in the two larger islands of the archipelago of Madeira. Based on the capsid protein sequence, the authors suggested that the high sequence identity, along with the short time spanning between outbreaks, points to dissemination from Porto Santo to Madeira. By including information of the full RHDV2 genome of strains from Azores, Madeira, and the Canary Islands, we confirm the results obtained by Carvalho et al., but further show that several subtypes of RHDV2 circulate in these islands: non-recombinant RHDV2 in the Canary Islands, G1/RHDV2 in Azores, Porto Santo and Madeira, and NP/RHDV2 also in Madeira. Here we conclude that RHDV2 has been independently introduced in these archipelagos, and that in Madeira at least two independent introductions must have occurred. We provide additional information on the dynamics of RHDV2 in the Macaronesian archipelagos of Azores, Madeira, and the Canary Islands and highlight the importance of analyzing RHDV2 complete genome.