Mirte Bosse

Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.

Genome sequencing reveals fine scale diversification and reticulation history during speciation in Sus

BackgroundElucidating the process of speciation requires an in-depth understanding of the evolutionary history of the species in question. Studies that rely upon a limited number of genetic loci do not always reveal actual evolutionary history, and often confuse inferences related to phylogeny and speciation. Whole-genome data, however, can overcome this issue by providing a nearly unbiased window into the patterns and processes of speciation. In order to reveal the complexity of the speciation process, we sequenced and analyzed the genomes of 10 wild pigs, representing morphologically or geographically well-defined species and subspecies of the genus Sus from insular and mainland Southeast Asia, and one African common warthog.ResultsOur data highlight the importance of past cyclical climatic fluctuations in facilitating the dispersal and isolation of populations, thus leading to the diversification of suids in one of the most species-rich regions of the world. Moreover, admixture analyses revealed extensive, intra- and inter-specific gene-flow that explains previous conflicting results obtained from a limited number of loci. We show that these multiple episodes of gene-flow resulted from both natural and human-mediated dispersal.ConclusionsOur results demonstrate the importance of past climatic fluctuations and human mediated translocations in driving and complicating the process of speciation in island Southeast Asia. This case study demonstrates that genomics is a powerful tool to decipher the evolutionary history of a genus, and reveals the complexity of the process of speciation.

Evolutionary dynamics of copy number variation in pig genomes in the context of adaptation and domestication.

BackgroundCopy number variable regions (CNVRs) can result in drastic phenotypic differences and may therefore be subject to selection during domestication. Studying copy number variation in relation to domestication is highly relevant in pigs because of their very rich natural and domestication history that resulted in many different phenotypes. To investigate the evolutionary dynamic of CNVRs, we applied read depth method on next generation sequence data from 16 individuals, comprising wild boars and domestic pigs from Europe and Asia.ResultsWe identified 3,118 CNVRs with an average size of 13 kilobases comprising a total of 39.2 megabases of the pig genome and 545 overlapping genes. Functional analyses revealed that CNVRs are enriched with genes related to sensory perception, neurological process and response to stimulus, suggesting their contribution to adaptation in the wild and behavioral changes during domestication. Variations of copy number (CN) of antimicrobial related genes suggest an ongoing process of evolution of these genes to combat food-borne pathogens. Likewise, some genes related to the omnivorous lifestyle of pigs, like genes involved in detoxification, were observed to be CN variable. A small portion of CNVRs was unique to domestic pigs and may have been selected during domestication. The majority of CNVRs, however, is shared between wild and domesticated individuals, indicating that domestication had minor effect on the overall diversity of CNVRs. Also, the excess of CNVRs in non-genic regions implies that a major part of these variations is likely to be (nearly) neutral. Comparison between different populations showed that larger populations have more CNVRs, highlighting that CNVRs are, like other genetic variation such as SNPs and microsatellites, reflecting demographic history rather than phenotypic diversity.ConclusionCNVRs in pigs are enriched for genes related to sensory perception, neurological process, and response to stimulus. The majority of CNVRs ascertained in domestic pigs are also variable in wild boars, suggesting that the domestication of the pig did not result in a change in CNVRs in domesticated pigs. The majority of variable regions were found to reflect demographic patterns rather than phenotypic.

Genomic analysis reveals selection for Asian genes in European pigs following human-mediated introgression

The independent domestication of local wild boar populations in Asia and Europe about 10,000 years ago led to distinct European and Asian pig breeds, each with very different phenotypic characteristics. During the Industrial Revolution, Chinese breeds were imported to Europe to improve commercial traits in European breeds. Here we demonstrate the presence of introgressed Asian haplotypes in European domestic pigs and selection signatures on some loci in these regions, using whole genome sequence data. The introgression signatures are widespread and the Asian haplotypes are rarely fixed. The Asian introgressed haplotypes are associated with regions harbouring genes involved in meat quality, development and fertility. We identify Asian-derived non-synonymous mutations in the AHR gene that associate with increased litter size in multiple European commercial lines. These findings demonstrate that increased fertility was an important breeding goal for early nineteenth century pig farmers, and that Asian variants of genes related to this trait were preferentially selected during the development of modern European pig breeds.

Transcriptomics reveals extensive inducible biotransformation in the soil-dwelling invertebrate Folsomia candida exposed to phenanthrene

BackgroundPolycyclic aromatic hydrocarbons are common pollutants in soil, have negative effects on soil ecosystems, and are potentially carcinogenic. The Springtail (Collembola) Folsomia candida is often used as an indicator species for soil toxicity. Here we report a toxicogenomic study that translates the ecological effects of the polycyclic aromatic hydrocarbon phenanthrene in soil to the early transcriptomic responses in Folsomia candida.ResultsMicroarrays were used to examine two different exposure concentrations of phenanthrene, namely the EC10 (24.95 mg kg-1 soil) and EC50 (45.80 mg kg-1 soil) on reproduction of this springtail, which evoked 405 and 251 differentially expressed transcripts, respectively. Fifty transcripts were differential in response to either concentration. Many transcripts encoding xenobiotic detoxification and biotransformation enzymes (phases I, II, and III) were upregulated in response to either concentration. Furthermore, indications of general and oxidative stress were found in response to phenanthrene. Chitin metabolism appeared to be disrupted particularly at the low concentration, and protein translation appeared suppressed at the high concentration of phenanthrene; most likely in order to reallocate energy budgets for the detoxification process. Finally, an immune response was evoked especially in response to the high effect concentration, which was also described in a previous transcriptomic study using the same effect concentration (EC50) of cadmium.ConclusionOur study provides new insights in the molecular mode of action of the important polluting class of polycyclic aromatic hydrocarbons in soil animals. Furthermore, we present a fast, sensitive, and specific soil toxicity test which enhances traditional tests and may help to improve current environmental risk assessments and monitoring of potentially polluted sites.

Using genome-wide measures of coancestry to maintain diversity and fitness in endangered and domestic pig populations

Conservation and breeding programs aim at maintaining the most diversity, thereby avoiding deleterious effects of inbreeding while maintaining enough variation from which traits of interest can be selected. Theoretically, the most diversity is maintained using optimal contributions based on many markers to calculate coancestries, but this can decrease fitness by maintaining linked deleterious variants. The heterogeneous patterns of coancestry displayed in pigs make them an excellent model to test these predictions. We propose methods to measure coancestry and fitness from resequencing data and use them in population management. We analyzed the resequencing data of Sus cebifrons, a highly endangered porcine species from the Philippines, and genotype data from the Pietrain domestic breed. By analyzing the demographic history of Sus cebifrons, we inferred two past bottlenecks that resulted in some inbreeding load. In Pietrain, we analyzed signatures of selection possibly associated with commercial traits. We also simulated the management of each population to assess the performance of different optimal contribution methods to maintain diversity, fitness, and selection signatures. Maximum genetic diversity was maintained using marker-by-marker coancestry, and least using genealogical coancestry. Using a measure of coancestry based on shared segments of the genome achieved the best results in terms of diversity and fitness. However, this segment-based management eliminated signatures of selection. We demonstrate that maintaining both diversity and fitness depends on the genomic distribution of deleterious variants, which is shaped by demographic and selection histories. Our findings show the importance of genomic and next-generation sequencing information in the optimal design of breeding or conservation programs.

Untangling the hybrid nature of modern pig genomes: a mosaic derived from biogeographically distinct and highly divergent Sus scrofa populations

The merging of populations after an extended period of isolation and divergence is a common phenomenon, in natural settings as well as due to human interference. Individuals with such hybrid origins contain genomes that essentially form a mosaic of different histories and demographies. Pigs are an excellent model species to study hybridization because European and Asian wild boars diverged ~1.2 Mya, and pigs were domesticated independently in Europe and Asia. During the Industrial Revolution in England, pigs were imported from China to improve the local pigs. This study utilizes the latest genomics tools to identify the origin of haplotypes in European domesticated pigs that are descendant from Asian and European populations. Our results reveal fine‐scale haplotype structure representing different ancient demographic events, as well as a mosaic composition of those distinct histories due to recently introgressed haplotypes in the pig genome. As a consequence, nucleotide diversity in the genome of European domesticated pigs is higher when at least one haplotype of Asian origin is present, and haplotype length correlates negatively with recombination frequency and nucleotide diversity. Another consequence is that the inference of past effective population size is influenced by the background of the haplotypes in an individual, but we demonstrate that by careful sorting based on the origin of haplotypes, both distinct demographic histories can be reconstructed. Future detailed mapping of the genomic distribution of variation will enable a targeted approach to increase genetic diversity of captive and wild populations, thus facilitating conservation efforts in the near future.

Copy number variation in the speciation of pigs: a possible prominent role for olfactory receptors

BackgroundUnraveling the genetic mechanisms associated with reduced gene flow between genetically differentiated populations is key to understand speciation. Different types of structural variations (SVs) have been found as a source of genetic diversity in a wide range of species. Previous studies provided detailed knowledge on the potential evolutionary role of SVs, especially copy number variations (CNVs), between well diverged species of e.g. primates. However, our understanding of their significance during ongoing speciation processes is limited due to the lack of CNV data from closely related species. The genus Sus (pig and its close relatives) which started to diverge ~4 Mya presents an excellent model for studying the role of CNVs during ongoing speciation.ResultsIn this study, we identified 1408 CNV regions (CNVRs) across the genus Sus. These CNVRs encompass 624 genes and were found to evolve ~2.5 times faster than single nucleotide polymorphisms (SNPs). The majority of these copy number variable genes are olfactory receptors (ORs) known to play a prominent role in food foraging and mate recognition in Sus. Phylogenetic analyses, including novel Bayesian analysis, based on CNVRs that overlap ORs retain the well-accepted topology of the genus Sus whereas CNVRs overlapping genes other than ORs show evidence for random drift and/or admixture.ConclusionWe hypothesize that inter-specific variation in copy number of ORs provided the means for rapid adaptation to different environments during the diversification of the genus Sus in the Pliocene. Furthermore, these regions might have acted as barriers preventing massive gene flow between these species during the multiple hybridization events that took place later in the Pleistocene suggesting a possible prominent role of ORs in the ongoing Sus speciation.

Whole-genome sequence analysis reveals differences in population management and selection of European low-input pig breeds

BackgroundA major concern in conservation genetics is to maintain the genetic diversity of populations. Genetic variation in livestock species is threatened by the progressive marginalisation of local breeds in benefit of high-output pigs worldwide. We used high-density SNP and re-sequencing data to assess genetic diversity of local pig breeds from Europe. In addition, we re-sequenced pigs from commercial breeds to identify potential candidate mutations responsible for phenotypic divergence among these groups of breeds.ResultsOur results point out some local breeds with low genetic diversity, whose genome shows a high proportion of regions of homozygosis (>50%) and that harbour a large number of potentially damaging mutations. We also observed a high correlation between genetic diversity estimates using high-density SNP data and Next Generation Sequencing data (r = 0.96 at individual level). The study of non-synonymous SNPs that were fixed in commercial breeds and also in any local breed, but with different allele, revealed 99 non-synonymous SNPs affecting 65 genes. Candidate mutations that may underlie differences in the adaptation to the environment were exemplified by the genes AZGP1 and TAS2R40. We also observed that highly productive breeds may have lost advantageous genotypes within genes involve in immune response – e.g. IL12RB2 and STAB1–, probably as a result of strong artificial in the intensive production systems in pig.ConclusionsThe high correlation between genetic diversity computed with the 60K SNP and whole genome re-sequence data indicates that the Porcine 60K SNP Beadchip provides reliable estimates of genomic diversity in European pig populations despite the expected bias. Moreover, this analysis gave insights for strategies to the genetic characterization of local breeds. The comparison between re-sequenced local pigs and re-sequenced commercial pigs made it possible to report candidate mutations to be responsible for phenotypic divergence among those groups of breeds. This study highlights the importance of low input breeds as a valuable genetic reservoir for the pig production industry. However, the high levels of ROHs, inbreeding and potentially damaging mutations emphasize the importance of the genetic characterization of local breeds to preserve their genomic variability.

Linking rhizosphere microbiome composition of wild and domesticated Phaseolus vulgaris to genotypic and root phenotypic traits

Plant domestication was a pivotal accomplishment in human history, but also led to a reduction in genetic diversity of crop species compared to their wild ancestors. How this reduced genetic diversity affected plant–microbe interactions belowground is largely unknown. Here, we investigated the genetic relatedness, root phenotypic traits and rhizobacterial community composition of modern and wild accessions of common bean (Phaseolus vulgaris) grown in agricultural soil from the highlands of Colombia, one of the centers of common bean diversification. Diversity Array Technology-based genotyping and phenotyping of local common bean accessions showed significant genetic and root architectural differences between wild and modern accessions, with a higher specific root length for the wild accessions. Canonical Correspondence Analysis indicated that the divergence in rhizobacterial community composition between wild and modern bean accessions is associated with differences in specific root length. Along the bean genotypic trajectory, going from wild to modern, we observed a gradual decrease in relative abundance of Bacteroidetes, mainly Chitinophagaceae and Cytophagaceae, and an increase in relative abundance of Actinobacteria and Proteobacteria, in particular Nocardioidaceae and Rhizobiaceae, respectively. Collectively, these results establish a link between common bean domestication, specific root morphological traits and rhizobacterial community assembly.